Empirical likelihood confidence intervals for adaptive cluster sampling

Adaptive cluster sampling (ACS) is an efficient sampling design for estimating parameters of rare and clustered populations. It is widely used in ecological research. The modified Hansen-Hurwitz (HH) and Horvitz-Thompson (HT) estimators based on small samples under ACS have often highly skewed distributions. In such situations, confidence intervals based on traditional normal approximation can lead to unsatisfactory results, with poor coverage properties. Christman and Pontius (Biometrics 56:503–510, 2000) showed that bootstrap percentile methods are appropriate for constructing confidence intervals from the HH estimator. But Perez and Pontius (J Stat Comput Simul 76:755–764, 2006) showed that bootstrap confidence intervals from the HT estimator are even worse than the normal approximation confidence intervals. In this article, we consider two pseudo empirical likelihood functions under the ACS design. One leads to the HH estimator and the other leads to a HT type estimator known as the Hájek estimator. Based on these two empirical likelihood functions, we derive confidence intervals for the population mean. Using a simulation study, we show that the confidence intervals obtained from the first EL function perform as good as the bootstrap confidence intervals from the HH estimator but the confidence intervals obtained from the second EL function perform much better than the bootstrap confidence intervals from the HT estimator, in terms of coverage rate.     

Statistical inference using stratified judgment post-stratified samples from finite populations

This paper develops statistical inference for population mean and total using stratified judgment post-stratified (SJPS) samples. The SJPS design selects a judgment post-stratified sample from each stratum. Hence, in addition to stratum structure, it induces additional ranking structure within stratum samples. SJPS is constructed from a finite population using either a with or without replacement sampling design. Inference is constructed under both randomization theory and a super population model. In both approaches, the paper shows that the estimators of population mean and total are unbiased. The paper also constructs unbiased estimators for the variance (mean square prediction error) of the sample mean (predictor of population mean), and develops confidence and prediction intervals for the population mean. The empirical evidence shows that the proposed estimators perform better than their competitors in the literature.     

Comparison of distribution function estimators for forestry applications following variable probability sampling

Practical considerations often motivate employing variable probability sampling designs when estimating characteristics of forest populations. Three distribution function estimators, the Horvitz-Thompson estimator, a difference estimator, and a ratio estimator, are compared following variable probability sampling in which the inclusion probabilities are proportional to an auxiliary variable, X. Relative performance of the estimators is affected by several factors, including the distribution of the inclusion probabilities, the correlation () between X and the response Y, and the position along the distribution function being estimated. Both the ratio and difference estimators are superior to the Horvitz-Thompson estimator. The difference estimator gains better precision than the ratio estimator toward the upper portion of the distribution function, but the ratio estimator is superior toward the lower end of the distribution function. The point along the distribution function at which the difference estimator becomes more precise than the ratio estimator depends on the sampling design, as well as the coefficient of variation of X and . A simple confidence interval procedure provides close to nominal coverage for intervals constructed from both the difference and ratio estimators, with the exception that coverage may be poor for the lower tail of the distribution function when using the ratio estimator.     

An estimating function approach to the inference of catch-effort models

A class of catch-effort models, which allows for heterogeneous removal probabilities, is proposed for closed populations. The model includes three types of removal probabilities: multiplicative, Poisson and logistic. The usual removal and generalized removal models then become special cases. The equivalence of the proposed model and a special type of capture-recapture model is discussed. A unified estimating function approach is used to estimate the initial population size. For the homogeneous model, the resulting population size estimator based on optimal estimating functions is asymptotically equivalent to the maximum likelihood estimator. One advantage for our approach is that it can be extended to handle the heterogeneous populations in which the maximum likelihood estimators do not exist. The bootstrap method is applied to construct variance estimators and confidence intervals. We illustrate the method by two real data examples. Results of a simulation study investigating the performance of the proposed estimation procedure are presented.     

Adaptive two-stage one-per-stratum sampling

We briefly describe adaptive cluster sampling designs in which the initial sample is taken according to a Markov chain one-per-stratum design (Breidt, 1995) and one or more secondary samples are taken within strata if units in the initial sample satisfy a given condition C. An empirical study of the behavior of the estimation procedure is conducted for three small artificial populations for which adaptive sampling is appropriate. The specific sampling strategy used in the empirical study was a single random-start systematic sample with predefined systematic samples within strata when the initially sampled unit in that stratum satisfies C. The bias of the Horvitz-Thompson estimator for this design is usually very small when adaptive sampling is conducted in a population for which it is suited. In addition, we compare the behavior of several alternative estimators of the standard error of the Horvitz-Thompson estimator of the population total. The best estimator of the standard error is population-dependent but it is not unreasonable to use the Horvitz-Thompson estimator of the variance. Unfortunately, the distribution of the estimator is highly skewed hence the usual approach of constructing confidence intervals assuming normality cannot be used here.     

Multiple data sources improve DNA-based mark-recapture population estimates of grizzly bears.

A fundamental challenge to estimating population size with mark-recapture methods is heterogeneous capture probabilities and subsequent bias of population estimates. Confronting this problem usually requires substantial sampling effort that can be difficult to achieve for some species, such as carnivores. We developed a methodology that uses two data sources to deal with heterogeneity and applied this to DNA mark-recapture data from grizzly bears (Ursus arctos). We improved population estimates by incorporating additional DNA "captures" of grizzly bears obtained by collecting hair from unbaited bear rub trees concurrently with baited, grid-based, hair snag sampling. We consider a Lincoln-Petersen estimator with hair snag captures as the initial session and rub tree captures as the recapture session and develop an estimator in program MARK that treats hair snag and rub tree samples as successive sessions. Using empirical data from a large-scale project in the greater Glacier National Park, Montana, USA, area and simulation modeling we evaluate these methods and compare the results to hair-snag-only estimates. Empirical results indicate that, compared with hair-snag-only data, the joint hair-snag-rub-tree methods produce similar but more precise estimates if capture and recapture rates are reasonably high for both methods. Simulation results suggest that estimators are potentially affected by correlation of capture probabilities between sample types in the presence of heterogeneity. Overall, closed population Huggins-Pledger estimators showed the highest precision and were most robust to sparse data, heterogeneity, and capture probability correlation among sampling types. Results also indicate that these estimators can be used when a segment of the population has zero capture probability for one of the methods. We propose that this general methodology may be useful for other species in which mark-recapture data are available from multiple sources.     

Bootstrap confidence intervals for adaptive cluster sampling design based on Horvitz–Thompson type estimators

Perez and Pontius (J Stat Comput Simul 76:755–764, 2006) introduced several bootstrap methods under adaptive cluster sampling using a Horvitz–Thompson type estimator. Using a simulation study, they showed that their proposed methods provide confidence intervals with highly understated coverage rates. In this article, we first show that their bootstrap methods provide biased bootstrap estimates. We then define two bootstrap methods, based on the method of Gross (Proceeding of the survey research methods section. American Statistical Association, Alexandria, VA, pp 181–184, 1980) and Bootstrap With Replacement, that provide unbiased bootstrap estimates of the population mean with bootstrap variances matching the corresponding unbiased variance estimator. Using a simulation study, we show that the bootstrap confidence intervals based on our proposed methods have better performance than those based on available bootstrap methods, in the sense of having coverage proportion closer to the nominal coverage level. We also compare the proposed intervals to empirical likelihood based intervals in small samples.     

Direct calculation of likelihood-based benchmark dose levels for quantitative responses

A benchmark dose (BMD) for quantitative responses is a lower confidence limit (LCL) on the effective dose corresponding to a specified risk level r. A commonly adopted method for calculating the BMD is to obtain a pointwise upper confidence curve U(d) on the risk function and then invert this relationship by solving the equation U(d)=r. The solution d is taken to be the BMD. Sciullo et al. (2000) have shown that the coverage achieved by this inversion method is at least as great as the coverage achieved by U (·) but that there is otherwise no general relationship between the two coverage probabilities. The present paper develops a method for direct calculation of the BMD based on the asymptotic distribution of the likelihood ratio statistic. It is further shown that the direct method and the inversion method are equivalent when U (·) is also based on the likelihood ratio. Since the direct method is known to be asymptotically correct, it follows that the LR-based inversion method is also asymptotically correct. However, the direct method is computationally faster and easier to program. Finally, some simulation studies are conducted to assess the small sample coverage probabilities of the direct method when responses follow either a normal or a gamma distribution.     

Testing the effectiveness of capture mark recapture population estimation techniques using a computer simulation with known population size

Estimation of small mammal population sizes is important for monitoring ecosystem condition and for conservation. Here, we test the accuracy of standard methods of population size estimation using Capture-Mark-Recapture (CMR) on a simulated population of agents. The use of a computer simulation allows complete control of population sizes and behaviors, thereby avoiding assumptions that may be violated in real populations. We find that the recommended protocol for CMR sampling, using uniformly distributed traps, consistently overestimates population sizes by as much as 100% when studies are conducted over only two trapping periods. More than 20 trapping periods are required before this method, or that of placing traps randomly, gives an accurate estimation of population size (i.e., within a 95% confidence limit of the actual value). Non-random sampling, by placing traps on runways used by small mammals, produces the most accurate, and least variable, estimates of population. However, we show that around 10 trapping periods are still required to produce an accurate population estimate using this method. Given that most real populations do not comply with the ‘ideal’ assumptions made by CMR, we suggest that population estimates based on CMR may be fundamentally flawed, and recommend that protocols for CMR population estimation methods may need revising.     

Spatial methods for plot-based sampling of wildlife populations

Classical sampling methods can be used to estimate the mean of a finite or infinite population. Block kriging also estimates the mean, but of an infinite population in a continuous spatial domain. In this paper, I consider a finite population version of block kriging (FPBK) for plot-based sampling. The data are assumed to come from a spatial stochastic process. Minimizing mean-squared-prediction errors yields best linear unbiased predictions that are a finite population version of block kriging. FPBK has versions comparable to simple random sampling and stratified sampling, and includes the general linear model. This method has been tested for several years for moose surveys in Alaska, and an example is given where results are compared to stratified random sampling. In general, assuming a spatial model gives three main advantages over classical sampling: (1) FPBK is usually more precise than simple or stratified random sampling, (2) FPBK allows small area estimation, and (3) FPBK allows nonrandom sampling designs.     

Evaluation of Local Ecological Knowledge as a Method for Collecting Extensive Data on Animal Abundance

Abstract:  The use of local ecological knowledge (LEK) has been advocated for biodiversity monitoring and management. To date, however, it has been underused in studying wild populations of animals and, particularly, in obtaining quantitative abundance estimates. We evaluated LEK as a tool for collecting extensive data on local animal abundance and population trends. We interviewed shepherds in southeastern Spain, asking them to estimate the local abundance of the terrestrial tortoise Testudo graeca . We quantified reliability of abundance estimates derived from interviews by comparing them with those obtained from standard field-sampling protocols (distance sampling). We also explored the complementarity of these 2 approaches. LEK provided high-quality and low-cost information about both distribution and abundance of T. graeca . Interviews with shepherds yielded abundance estimates in a much wider range than linear transects, which only detected the species in the upper two-thirds of its abundance range. Abundance estimates from both methodologies showed a close relationship. Analysis of confidence intervals indicated local knowledge could be used to estimate mean local abundances and to detect mean population trends. A cost analysis determined that the information derived from LEK was 100 times cheaper than that obtained through linear-transect surveys. Our results should further the use of LEK as a standard tool for sampling the quantitative abundance of a great variety of taxa, particularly when population densities are low and traditional sampling methods are expensive or difficult to implement.     

Design based estimation for ranked set sampling in finite populations

In this paper, we consider design-based estimation using ranked set sampling (RSS) in finite populations. We first derive the first and second-order inclusion probabilities for an RSS design and present two Horvitz–Thompson type estimators using these inclusion probabilities. We also develop an alternate Hansen–Hurwitz type estimator and investigate its properties. In particular, we show that this alternate estimator always outperforms the usual Hansen–Hurwitz type estimator in the simple random sampling with replacement design with comparable sample size. We also develop formulae for ratio estimator for all three developed estimators. The theoretical results are augmented by numerical and simulation studies as well as a case study using a well known data set. These show that RSS design can yield a substantial improvement in efficiency over the usual simple random sampling design in finite populations.     

Network inclusion probabilities and Horvitz-Thompson estimation for adaptive simple Latin square sampling

Consider a survey of a plant or animal species in which abundance or presence/absence will be recorded. Further assume that the presence of the plant or animal is rare and tends to cluster. A sampling design will be implemented to determine which units to sample within the study region. Adaptive cluster sampling designs Thompson (1990) are sampling designs that are implemented by first selecting a sample of units according to some conventional probability sampling design. Then, whenever a specified criterion is satisfied upon measuring the variable of interest, additional units are adaptively sampled in neighborhoods of those units satisfying the criterion. The success of these adaptive designs depends on the probabilities of finding the rare clustered events, called networks. This research uses combinatorial generating functions to calculate network inclusion probabilities associated with a simple Latin square sample. It will be shown that, in general, adaptive simple Latin square sampling when compared to adaptive simple random sampling will (i) yield higher network inclusion probabilities and (ii) provide Horvitz-Thompson estimators with smaller variability.     

Variances and variance estimators of the improved ratio estimators under adaptive cluster sampling

Although not design-unbiased, the ratio estimator is recognized as more efficient when a certain degree of correlation exists between the variable of primary interest and the auxiliary variable. Meanwhile, the Rao–Blackwell method is another commonly used procedure to improve estimation efficiency. Various improved ratio estimators under adaptive cluster sampling (ACS) that make use of the auxiliary information together with the Rao–Blackwellized univariate estimators have been proposed in past research studies. In this article, the variances and the associated variance estimators of these improved ratio estimators are proposed for a thorough framework of statistical inference under ACS. Performance of the proposed variance estimators is evaluated in terms of the absolute relative percentage bias and the empirical mean-squared error. As expected, results show that both the absolute relative percentage bias and the empirical mean-squared error decrease as the initial sample size increases for all the variance estimators. To evaluate the confidence intervals based on these variance estimators and the finite-population Central Limit Theorem, the coverage rate and the interval width are used. These confidence intervals suffer a disadvantage similar to that of the conventional ratio estimator. Hence, alternative confidence intervals based on a certain type of adjusted variance estimators are constructed and assessed in this article.     

Optimization of capture–recapture monitoring of elusive species illustrated with a threatened grasshopper

Information on population sizes and trends of threatened species is essential for their conservation, but obtaining reliable estimates can be challenging. We devised a method to improve the precision of estimates of population size obtained from capture–recapture studies for species with low capture and recapture probabilities and short seasonal activity, illustrated with population data of an elusive grasshopper (Prionotropis rhodanica). We used data from 5 capture–recapture studies to identify methodological and environmental factors affecting capture and recapture probabilities and estimates of population size. In a simulation, we used the population size and capture and recapture probability estimates obtained from the field studies to identify the minimum number of sampling occasions needed to obtain unbiased and robust estimates of population size. Based on these results we optimized the capture–recapture design, implemented it in 2 additional studies, and compared their precision with those of the nonoptimized studies. Additionally, we simulated scenarios based on thresholds of population size in criteria C and D of the International Union for Conservation of Nature (IUCN) Red List to investigate whether estimates of population size for elusive species can reliably inform red-list assessments. Identifying parameters that affect capture and recapture probabilities (for the grasshopper time since emergence of first adults) and optimizing field protocols based on this information reduced study effort (−6% to −27% sampling occasions) and provided more precise estimates of population size (reduced coefficient of variation) compared with nonoptimized studies. Estimates of population size from the scenarios based on the IUCN thresholds were mostly unbiased and robust (only the combination of very small populations and little study effort produced unreliable estimates), suggesting capture–recapture can be considered reliable for informing red-list assessments. Although capture–recapture remains difficult and costly for elusive species, our optimization procedure can help determine efficient protocols to increase data quality and minimize monitoring effort.     

Estimation methods for nonlinear state-space models in ecology

The use of nonlinear state-space models for analyzing ecological systems is increasing. A wide range of estimation methods for such models are available to ecologists, however it is not always clear, which is the appropriate method to choose. To this end, three approaches to estimation in the theta logistic model for population dynamics were benchmarked by Wang (2007). Similarly, we examine and compare the estimation performance of three alternative methods using simulated data. The first approach is to partition the state-space into a finite number of states and formulate the problem as a hidden Markov model (HMM). The second method uses the mixed effects modeling and fast numerical integration framework of the AD Model Builder (ADMB) open-source software. The third alternative is to use the popular Bayesian framework of BUGS. The study showed that state and parameter estimation performance for all three methods was largely identical, however with BUGS providing overall wider credible intervals for parameters than HMM and ADMB confidence intervals.     

Triangulation based inclusion probabilities: a design-unbiased sampling approach

A probabilistic sampling approach for design-unbiased estimation of area-related quantitative characteristics of spatially dispersed population units is proposed. The developed field protocol includes a fixed number of 3 units per sampling location and is based on partial triangulations over their natural neighbors to derive the individual inclusion probabilities. The performance of the proposed design is tested in comparison to fixed area sample plots in a simulation with two forest stands. Evaluation is based on a general approach for areal sampling in which all characteristics of the resulting population of possible samples is derived analytically by means of a complete tessellation of the areal sampling frame. The example simulation shows promising results. Expected errors under this design are comparable to sample plots including a much greater number of trees per plot.     

Spatially-balanced adaptive web sampling

The paper deals with sampling from a finite population that is distributed over space and has a highly uneven spatial distribution. It suggests a sampling design that allocates a portion of the sample units that are well spread over the population and sequentially selects the remaining units in sub-areas that appear to be of more interest according to the study variable values observed during the survey. In order to estimate the population mean while using this sampling design, a computationally intense estimator, obtained via the Rao–Blackwell approach, is proposed and a resampling method is used that makes the inference computationally feasible. The whole sampling strategy is evaluated through several Monte Carlo experiments.     

Animal abundance estimation in independent observer line transect surveys

The theory of conventional line transect surveys is based on an essential assumption that 100% detection of animals right on the transect lines can be achieved. When this assumption fails, independent observer line transect surveys are used. This paper proposes a general approach, based on a conditional likelihood, which can be carried out either parametrically or nonparametrically, to estimate the abundance of non-clustered biological populations using data collected from independent observer line transect surveys. A nonparametric estimator is specifically proposed which combines the conditional likelihood and the kernel smoothing method. It has the advantage that it allows the data themselves to dictate the form of the detection function, free of any subjective choice. The bias and the variance of the nonparametric estimator are given. Its asymptotic normality is established which enables construction of confidence intervals. A simulation study shows that the proposed estimator has good empirical performance, and the confidence intervals have good coverage accuracy.     

Assessing and testing prediction uncertainty for single tree-based models: A case study applied to northern hardwood stands in southern Québec, Canada

Estimating prediction uncertainty for a single tree-based model is hindered by the complex structure of these models. In this paper, we addressed this issue with a case study applied to northern hardwood stands in Québec, Canada. SaMARE is a stochastic single tree-based model that was designed for these types of stands. Using a Monte Carlo approach, the model can provide a mean predicted value and its confidence limits for some plot-level attributes.The mean predicted values were compared to observed values in terms of bias and accuracy. In addition to these common statistics, we compared nominal coverage of Monte Carlo-simulated confidence intervals with real (observed) coverage to verify the adequacy of the simulated uncertainty. A comparison was made using several plot-level attributes, which exhibited an increasing discriminative complexity. This complexity ranges from coarse attributes, such as all-species basal area, up to more complex ones, such as basal area for stems of a particular species and with sawlog potential.The results showed that in terms of absolute value, biases were small, but could be relatively high with respect to the average observed value when the discriminative complexity of the attribute increased. The comparison between nominal and real coverage of confidence intervals gave satisfactory results for all-species plot-level attributes. However, for some species-specific attributes, the Monte Carlo-simulated confidence intervals overestimated the real coverage.