Adaptive cluster sampling with clusters selected without replacement and stopping rule

Adaptive cluster sampling (ACS) has received much attention in recent years since it yields more precise estimates than conventional sampling designs when applied to rare and clustered populations. These results, however, are impacted by the availability of some prior knowledge about the spatial distribution and the absolute abundance of the population under study. This prior information helps the researcher to select a suitable critical value that triggers the adaptive search, the neighborhood definition and the initial sample size. A bad setting of the ACS design would worsen the performance of the adaptive estimators. In particular, one of the greatest weaknesses in ACS is the inability to control the final sampling effort if, for example, the critical value is set too low. To overcome this drawback one can introduce ACS with clusters selected without replacement where one can fix in advance the number of distinct clusters to be selected or ACS with a stopping rule which stops the adaptive sampling when a predetermined sample size limit is reached or when a given stopping rule is verified. However, the stopping rule breaks down the theoretical basis for the unbiasedness of the ACS estimators introducing an unknown amount of bias in the estimates. The current study improves the performance of ACS when applied to patchy and clustered but not rare populations and/or less clustered populations. This is done by combining the stopping rule with ACS without replacement of clusters so as to further limit the sampling effort in form of traveling expenses by avoiding repeat observations and by reducing the final sample size. The performance of the proposed design is investigated using simulated and real data.     

A new design for sampling with adaptive sample plots

Adaptive cluster sampling (ACS) is a sampling technique for sampling rare and geographically clustered populations. Aiming to enhance the practicability of ACS while maintaining some of its major characteristics, an adaptive sample plot design is introduced in this study which facilitates field work compared to “standard” ACS. The plot design is based on a conditional plot expansion: a larger plot (by a pre-defined plot size factor) is installed at a sample point instead of the smaller initial plot if a pre-defined condition is fulfilled. This study provides insight to the statistical performance of the proposed adaptive plot design. A design-unbiased estimator is presented and used on six artificial and one real tree position maps to estimate density (number of objects per ha). The performance in terms of coefficient of variation is compared to the non-adaptive alternative without a conditional expansion of plot size. The adaptive plot design was superior in all cases but the improvement depends on (1) the structure of the sampled population, (2) the plot size factor and (3) the critical value (the minimum number of objects triggering an expansion). For some spatial arrangements the improvement is relatively small. The adaptive design may be particularly attractive for sampling in rare and compactly clustered populations with an appropriately chosen plot size factor.     

Designing an efficient adaptive cluster sample

The choice of neighborhood definition and critical value in adaptive cluster sampling is critical for designing an efficient survey. In designing an efficient adaptive cluster sample one should aim for a small difference between the initial and final sample size, and a small difference between the within-network and population variances. However, the two aims can be at odds with each other because small differences between initial and final sample size usually means small within-network variance. One way to help in designing an efficient survey is to think in terms of small network sizes since the network size is a function of both critical value and neighborhood definition. One should aim for networks that are small enough to ensure the final sample size is not excessively large compared with the initial sample size but large enough to ensure the within-network variance is a reasonable fraction of the population variance. In this study surveys that had networks that were two to four units in size were the most efficient.     

Two neighborhood-free plot designs for adaptive sampling of forests

Adaptive cluster sampling (ACS) has the potential of being superior for sampling rare and geographically clustered populations. However, setting up an efficient ACS design is challenging. In this study, two adaptive plot designs are proposed as alternatives: one for fixed-area plot sampling and the other for relascope sampling (also known as variable radius plot sampling). Neither includes a neighborhood search which makes them much easier to execute. They do, however, include a conditional plot expansion: at a sample point where a predefined condition is satisfied, sampling is extended to a predefined larger cluster-plot or a larger relascope plot. Design-unbiased estimators of population total and its variance are derived for each proposed design, and they are applied to ten artificial and one real tree position maps to estimate density (number of trees per ha) and basal area (the cross-sectional area of a tree stem at breast height) per hectare. The performances—in terms of relative standard error (SE%)—of the proposed designs and their non-adaptive alternatives are compared. The adaptive plot designs were superior for the clustered populations in all cases of equal sample sizes and in some cases of equal area of sample plots. However, the improvement depends on: (1) the plot size factor; (2) the critical value (the minimum number of trees triggering an expansion); (3) the subplot distance for the adapted cluster-plots, and (4) the spatial arrangement of the sampled population. For some spatial arrangements, the improvement is relatively small. The adaptive designs may be particularly attractive for sampling in rare and compactly clustered populations with critical value of 1, subplot distance equal to the diameter of initial circular plots, or plot size factor of 2.5 for an initial basal area factor of 2.     

Estimator bias and efficiency for adaptive cluster sampling with order statistics and a stopping rule

Practical problems facing adaptive cluster sampling with order statistics (acsord) are explored using Monte Carlo simulation for three simulated fish populations and two known waterfowl populations. First, properties of an unbiased Hansen-Hurwitz (HH) estimator and a biased alternative Horvitz-Thompson (HT) estimator are evaluated. An increase in the level of population aggregation or the initial sample size increases the efficiencies of the two acsord estimators. For less aggregated fish populations, the efficiencies decrease as the order statistic parameter r (the number of units about which adaptive sampling is carried out) increases; for the highly aggregated fish and waterfowl populations, they increase with r. Acsord is almost always more efficient than simple random sampling for the highly aggregated populations. Positive bias is observed for the HT estimator, with the maximum bias usually occurring at small values of r. Secondly, a stopping rule at the Sth iteration of adaptive sampling beyond the initial sampling unit was applied to the acsord design to limit the otherwise open-ended sampling effort. The stopping rule induces relatively high positive bias to the HH estimator if the level of the population aggregation is high, the stopping level S is small, and r is large. The bias of HT is not very sensitive to the stopping rule and its bias is often reduced by the stopping rule at smaller values of r. For more aggregated populations, the stopping rule often reduces the efficiencies of the estimators compared to the non-stopping-rule scheme, but acsord still remains more efficient than simple random sampling. Despite its bias and lack of theoretical grounding, the HT estimator is usually more efficient than the HH estimator. In the stopping rule case, the HT estimator is preferable, because its bias is less sensitive to the stopping level.     

Adaptive cluster double sampling with post stratification with application to an epiphytic lichen community

The implementation of an adaptive cluster sampling design often becomes logistically challenging because variation in the final sampling effort introduces uncertainty in survey planning. To overcome this drawback, an inexpensive and easy to measure auxiliary variable could be used in a two-phase survey strategy, called adaptive cluster double sampling (Félix-Medina and Thompson in Biometrika 91:877–891, 2004). In this paper, a two-phase sampling strategy is proposed which combines the idea of adaptive cluster double sampling with the principle of post-stratification. In the first-phase an adaptive cluster sample is selected by means of an inexpensive auxiliary variable. Networks from the first phase sampling are then post-stratified according to their size. In the second-phase, the network structure is used to select a subsample of units by means of stratified random sampling. The proposed sampling strategy employs stratification without requiring an a priori delineation of the strata. Indeed, the strata sizes are estimated in the course of the two-phase sampling process. Therefore, it is suitable for situations where stratification is suspected to be efficient but strata cannot be easily delineated in advance. In this framework, a new type of estimator for the population mean which mimics the stratified sampling mean estimator and an estimator of the sampling variance are proposed. The results of a simulation study confirm, as expected, that the use of post-stratification leads to gain in precision for the estimator. The proposed sampling strategy is applied for targeting an epiphytic lichen community Lobarion pulmonariae in a forest area of the Northern Apennines (N-Italy), characterized by several species of conservation concern.     

Estimation of rare and clustered population mean using stratified adaptive cluster sampling

For many clustered populations, the prior information on an initial stratification exists but the exact pattern of the population concentration may not be predicted. Under this situation, the stratified adaptive cluster sampling (SACS) may provide more efficient estimates than the other conventional sampling designs for the estimation of rare and clustered population parameters. For practical interest, we propose a generalized ratio estimator with the single auxiliary variable under the SACS design. The expressions of approximate bias and mean squared error (MSE) for the proposed estimator are derived. Numerical studies are carried out to compare the performances of the proposed generalized estimator over the usual mean and combined ratio estimators under the conventional stratified random sampling (StRS) using a real population of redwood trees in California and generating an artificial population by the Poisson cluster process. Simulation results show that the proposed class of estimators may provide more efficient results than the other estimators considered in this article for the estimation of highly clumped population.

     

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.     

Complex national sampling design for long-term monitoring of protected dry grasslands in Switzerland

We describe a probabilistic sampling design of circular permanent plots for the long-term monitoring of protected dry grasslands in Switzerland. The population under study is defined by the perimeter of a national inventory. The monitoring focus is on the species composition of the protected grassland vegetation and derived conservation values. Efficient trend estimations are required for the whole country and for some predefined target groups (six biogeographical regions and eleven vegetation types). The target groups are equally important regardless of their size. Consequently, intensified sampling of the less frequent groups is essential for sample efficiency. The prior information needed to draw a targeted sample is obtained from the sampling frame and external databases. The logistics and generalized delineation of the target population may pose further problems. Thus, investments in fieldwork and travel time should be well balanced by selecting a cluster sample. Second, any access problems in the field and non-target units in the sample should be compensated for by selecting reserve plots as they otherwise may considerably reduce the effective sample size. Finally, the design has to be flexible as the sampling frame may change over time and sampling intensity might have to be adjusted to redefined budgets or requirements. Likewise, the variables and biological items of interest may change. To fulfil all these constraints and to optimally use the available prior information, we propose a multi-stage self-weighted unequal probability sampling design. The design uses modern techniques such as: balanced sampling, spreading, stratified balancing, calibration, unequal probability sampling and power allocation. This sampling design meets the numerous requirements of this study and provides a very efficient estimator.     

A framework for adapting survey design through time for wildlife population assessment

Sampling strategies for monitoring the status and trends in wildlife populations are often determined before the first survey is undertaken. However, there may be little information about the distribution of the population and so the sample design may be inefficient. Through time, as data are collected, more information about the distribution of animals in the survey region is obtained but it can be difficult to incorporate this information in the survey design. This paper introduces a framework for monitoring motile wildlife populations within which the design of future surveys can be adapted using data from past surveys whilst ensuring consistency in design-based estimates of status and trends through time. In each survey, part of the sample is selected from the previous survey sample using simple random sampling. The rest is selected with inclusion probability proportional to predicted abundance. Abundance is predicted using a model constructed from previous survey data and covariates for the whole survey region. Unbiased design-based estimators of status and trends and their variances are derived from two-phase sampling theory. Simulations over the short and long-term indicate that in general more precise estimates of status and trends are obtained using this mixed strategy than a strategy in which all of the sample is retained or all selected with probability proportional to predicted abundance. Furthermore the mixed strategy is robust to poor predictions of abundance. Estimates of status are more precise than those obtained from a rotating panel design.     

Ranked set sample design for environmental investigations

Ranked set sampling can provide an efficient basis for estimating parameters of environmental variables, particularly when sampling costs are intrinsically high. Various ranked set estimators are considered for the population mean and contrasted in terms of their efficiencies and useful- ness, with special concern for sample design considerations. Specifically, we consider the effects of the form of the underlying random variable, optimisation of efficiency and how to allocate sampling effort for best effect (e.g. one large sample or several smaller ones of the same total size). The various prospects are explored for two important positively skew random variables (lognormal and extreme value) and explicit results are given for these cases. Whilst it turns out that the best approach is to use the largest possible single sample and the optimal ranked set best linear estimator (ranked set BLUE), we find some interesting qualitatively different conclusions for the two skew distributions     

Regression-type estimators for adaptive two-stage sequential sampling

Adaptive two-stage sequential sampling (ATSSS) design was developed to observe more rare units and gain higher efficiency, in the sense of having a smaller variance estimator, than conventional sampling designs with equal effort for rare and spatially cluster populations. For certain rare populations, incorporating auxiliary variables into a sampling design can further improve the observation of rare units and increase efficiency. In this article, we develop regression-type estimators for ATSSS so that auxiliary variables can be incorporated into the ATSSS design when warranted. Simulation studies on two populations show that the regression-type estimators can significantly increase the efficiency of ATSSS and the detection of more rare units as compared to conventional sampling counterparts. Simulation of sampling of desert shrubs in Inner Mongolia (one of the two populations studied) showed that by incorporating a GIS auxiliary variable into ATSSS with the regression estimators resulted in a gain in efficiency over ATSSS without the auxiliary variable. Further, we found that the use of the GIS auxiliary variable in a conventional two-stage design with a regression estimator did not show a gain in efficiency.     

Detecting Rare Species with Random or Subjective Sampling: a Case Study of Red-Listed Saproxylic Beetles in Boreal Sweden

Abstract:  Efficient sampling design in field studies is important for economical and statistical reasons. We compared two ways to distribute sampling effort over an area, either randomly or subjectively. We searched for red-listed saproxylic (wood-living) beetles in 30 spruce stands in boreal Sweden by sifting wood from dead trees. We randomly selected positions within each stand with a geographic positioning system and sampled the nearest dead tree (random sample). In the same stand we also sampled dead trees that, based on literature, were likely to host such species (subjective sampling). The subjective sampling (two to five samples per stand, depending on stand size) was compared with the higher, random sampling effort (fixed level of 12 samples/stand). Subjective sampling was significantly more efficient. Red-listed species were found in 36% of the subjective samples and in 16% of the random samples. Nevertheless, the larger random effort resulted in a comparable number of red-listed species per stand and in 13 detected species in total (vs. 12 species with subjective sampling). Random sampling was less efficient, but provided an unbiased alternative more suitable for statistical purposes, as needed in, for example, monitoring programs. Moreover, new species-specific knowledge can be gained through random searches.     

Untested assumptions about within-species sample size and missing data in interspecific studies

Phylogenetic comparative studies rely on species-specific data that often contain missing values and/or differ in sample size among species. These phenomena may violate statistical assumptions about the non-random variance component in sampling effort. A major reason why this assumption is often not fulfilled is because the probability of being sampled (i.e., being captured or observed) may depend on species-specific characteristics. Here, we test this assumption by using information on within-species sample sizes and missing data from five independent comparative datasets of European birds. First, we show that the two estimates of data availability (missing values and within-species sample size) are positively correlated and are associated with research effort in general (the number of papers published). Second, we demonstrate biologically meaningful relationships between data availability and phenotypic traits. For example, population size, risk-taking, and habitat specialization independently predicted within-species sample size. The key determinants of missing data were population size and distribution range. However, data availability was not structured by phylogenetic relationships. These results indicate that the accuracy of sampling is repeatable and distributed non-randomly among species, as several species-specific attributes determined the probability of observation. Therefore, data availability seems to be a species-specific trait that can be shaped by ecology, life history, and behavior. Such relationships raise issues about non-random sampling, which requires attention in comparative studies.     

Adjusted two-stage adaptive cluster sampling

An adjusted two-stage sampling procedure is discussed for adaptive cluster sampling where some networks are large and others are small. A two-stage sample is drawn from the large networks and a single-stage sample is drawn from the rest. The simple random sampling (SRS) procedure without replacement is used at the initial stage. An estimator for the population mean along with its properties is discussed.     

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.     

Asymptotics in adaptive cluster sampling

In this article we consider asymptotic properties of the Horvitz-Thompson and Hansen-Hurwitz types of estimators under the adaptive cluster sampling variants obtained by selecting the initial sample by simple random sampling without replacement and by unequal probability sampling with replacement. We develop an asymptotic framework, which basically assumes that the number of units in the initial sample, as well as the number of units and networks in the population tend to infinity, but that the network sizes are bounded. Using this framework we prove that under each of the two variants of adaptive sampling above mentioned, both the Horvitz-Thompson and Hansen-Hurwitz types of estimators are design-consistent and asymptotically normally distributed. In addition we show that the ordinary estimators of their variances are also design-consistent estimators.     

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.     

Contributions to composite sampling

The initial use of composite sampling involved the analysis of many negative samples with relatively high laboratory cost (Dorfman sampling). We propose a method of double compositing and compare its efficiency with Dorfman sampling. The variability of composite measurement samples has environmental interest (hot spots). The precision of these estimates depends on the kurtosis of the distribution; leptokurtic distributions (₂ > 0) have increased precision as the number of field samples is increased. The opposite effect is obtained for platykurtic distributions. In the lognormal case, coverage probabilities are reasonable for < 0.5. The Poisson distribution can be associated with temporal compositing, of particular interest where radioactive measurements are taken. Sample size considerations indicate that the total sampling effort is directly proportional to the length of time sampled. If there is background radiation then increasing levels of this radiation require larger sample sizes to detect the same difference in radiation.     

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.