Abstract
When detection or occupancy probability is small or when the number of sites and number of visits per site is small, maximum likelihood estimators (MLE) of site occupancy parameters have large biases, are numerically unstable, and the corresponding confidence intervals have smaller than nominal coverage. We propose an alternative method of estimation, based on penalized likelihood. This method is numerically stable, the estimators have smaller mean square error than the MLE, and associated confidence intervals have close to nominal coverage.