Spatial modelling of individual-level parasite counts using the negative binomial distribution

We present a spatial model for the mean and correlation of highly dispersed count data, and apply it to individual-level counts of the nematode Wuchereria bancrofti, a parasite of humans which causes the disease lymphatic filariasis. Our model uses the negative binomial distribution, whose shape parameter is a convenient index of over-dispersion. Spatial association is quantified in terms of a characteristic length, which has an intuitive interpretation as the distance over which correlation decreases by half. Demographic surveillance and mapping enable us to include individual-level covariates such as age and sex. We discuss the distinctive features of our model and interpret the results in terms of the epidemiology of lymphatic filariasis and possible implications for control programmes.     

Estimating negative binomial parameters from occurrence data with detection times

The negative binomial distribution is a common model for the analysis of count data in biology and ecology. In many applications, we may not observe the complete frequency count in a quadrat but only that a species occurred in the quadrat. If only occurrence data are available then the two parameters of the negative binomial distribution, the aggregation index and the mean, are not identifiable. This can be overcome by data augmentation or through modeling the dependence between quadrat occupancies. Here, we propose to record the (first) detection time while collecting occurrence data in a quadrat. We show that under what we call proportionate sampling, where the time to survey a region is proportional to the area of the region, that both negative binomial parameters are estimable. When the mean parameter is larger than two, our proposed approach is more efficient than the data augmentation method developed by Solow and Smith ( 2010 , Am. Nat. 176 , 96–98), and in general is cheaper to conduct. We also investigate the effect of misidentification when collecting negative binomially distributed data, and conclude that, in general, the effect can be simply adjusted for provided that the mean and variance of misidentification probabilities are known. The results are demonstrated in a simulation study and illustrated in several real examples.     

Bias-corrected maximum likelihood estimator of the negative binomial dispersion parameter

We derive a first-order bias-corrected maximum likelihood estimator for the negative binomial dispersion parameter. This estimator is compared, in terms of bias and efficiency, with the maximum likelihood estimator investigated by Piegorsch (1990, Biometrics46, 863-867), the moment and the maximum extended quasi-likelihood estimators investigated by Clark and Perry (1989, Biometrics45, 309-316), and a double-extended quasi-likelihood estimator. The bias-corrected maximum likelihood estimator has superior bias and efficiency properties in most instances. For ease of comparison we give results for the two-parameter negative binomial model. However, an example involving negative binomial regression is given.     

A score test for testing a zero-inflated Poisson regression model against zero-inflated negative binomial alternatives

Count data often show a higher incidence of zero counts than would be expected if the data were Poisson distributed. Zero-inflated Poisson regression models are a useful class of models for such data, but parameter estimates may be seriously biased if the nonzero counts are overdispersed in relation to the Poisson distribution. We therefore provide a score test for testing zero-inflated Poisson regression models against zero-inflated negative binomial alternatives.     

Optimizing sampling design to infer the number of marine turtles nesting on low and high density sea turtle rookeries using convolution of negative binomial distribution

Reliable monitoring of wildlife populations represents a non-negligible cost, and in a limited-resource world, resources allocated to monitoring are not devoted to actions to solve identified problems.I explore resource efficient survey designs based on a negative binomial distribution including variable survey intervals for marine turtles using track counts as an index of female activity. In the modified procedure, all new tracks between two monitoring patrols are recorded. These data are analyzed by statistical models that take advantage of the statistical properties of the sum of counts.The outputs of models with different lagged monitoring dates (3–10 days) are compared with the outputs of daily surveys using extrapolations from high and low density populations. Results show that the quality of the estimates is similar when total time series analysis is compared with situations in which only a fourth, a seventh, or a tenth of monitoring daily during the season are used.This solution permits the reallocation of funds from monitoring to other conservation activities. Furthermore, the efficient sampling design and the statistical methods allow getting similar information with less effort.     

Small-sample estimation of negative binomial dispersion, with applications to SAGE data

We derive a quantile-adjusted conditional maximum likelihood estimator for the dispersion parameter of the negative binomial distribution and compare its performance, in terms of bias, to various other methods. Our estimation scheme outperforms all other methods in very small samples, typical of those from serial analysis of gene expression studies, the motivating data for this study. The impact of dispersion estimation on hypothesis testing is studied. We derive an "exact" test that outperforms the standard approximate asymptotic tests.