Abstract:In capture–recapture experiments, the presence of overdispersion and heterogeneity necessitates the use of the negative binomial regression model for inferring population sizes. However, within this model, existing methods based on likelihood and ratio regression for estimating the dispersion parameter often face boundary and nonidentifiability issues. These problems can result in nonsensically large point estimates and unbounded upper limits of confidence intervals for the population size. We present a penalized empirical likelihood technique for solving these two problems by imposing a half-normal prior on the population size. Based on the proposed approach, a maximum penalized empirical likelihood estimator with asymptotic normality and a penalized empirical likelihood ratio statistic with asymptotic chi-square distribution are derived. To improve numerical performance, we present an effective expectation-maximization (EM) algorithm. In the M-step, optimization for the model parameters could be achieved by fitting a standard negative binomial regression model via the R basic function glm.nb(). This approach ensures the convergence and reliability of the numerical algorithm. Using simulations, we analyze several synthetic datasets to illustrate three advantages of our methods in finite-sample cases: complete mitigation of the boundary problem, more efficient maximum penalized empirical likelihood estimates, and more precise penalized empirical likelihood ratio interval estimates compared to the estimates obtained without penalty. These advantages are further demonstrated in a case study estimating the abundance of black bears (Ursus americanus) at the U.S. Army's Fort Drum Military Installation in northern New York.

Sample Size Estimation for Negative Binomial Regression Comparing Rates of Recurrent Events with Unequal Follow-Up Time

Sample size calculation for the Andersen-Gill model comparing rates of recurrent events

Sample size calculation for clinical trials with correlated count measurements based on the negative binomial distribution

Group sequential designs for negative binomial outcomes

Sample size calculation for studies with grouped survival data

Sample size estimation for recurrent event data using multifrailty and multilevel survival models

The design and analysis of randomized trials with recurrent events

Recurrent Events Analysis in the Presence of Terminal Event and Zero-recurrence Subjects

A Practical Simulation Method to Calculate Sample Size of Group Sequential Trials for Time-to-Event Data under Exponential and Weibull Distribution

A Penalized Empirical Likelihood Approach for Estimating Population Sizes under the Negative Binomial Regression Model

Sample size estimation for modified Poisson analysis of cluster randomized trials with a binary outcome

Estimating waiting time for compatible blood: A Negative Binomial approach

Recurrent event analysis in the presence of real-time high frequency data via random subsampling

An improved sample size calculation method for score tests in generalized linear models

Blinded continuous information monitoring of recurrent event endpoints with time trends in clinical trials

Simulation‐based sample size calculations of marginal proportional means models for recurrent events with competing risks

Closed-form REML estimators and sample size determination for mixed effects models for repeated measures under monotone missingness

Non-parametric cure rate estimation under insufficient follow-up using extremes

Left-censored recurrent event analysis in epidemiological studies: a proposal when the number of previous episodes is unknown

Growth Estimators and Confidence Intervals for the Mean of Negative Binomial Random Variables with Unknown Dispersion

Computation of the Distribution of the Sum of Independent Negative Binomial Random Variables