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.     

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.     

Restricted adaptive cluster sampling

Adaptive cluster sampling can be a useful design for sampling rare and patchy populations. With this design the initial sample size is fixed but the size of the final sample (and total sampling effort) cannot be predicted prior to sampling. For some populations the final sample size can be quite variable depending on the level of patchiness. Restricted adaptive cluster sampling is a proposed modification where a limit is placed on the sample size prior to sampling and quadrats are selected sequentially for the initial sample size. As a result there is less variation in the final sample size and the total sampling effort can be predicted with some certainty, which is impor- tant for many ecological studies. Estimates of density are biased with the restricted design but under some circumstances the bias can be estimated well by bootstrapping. © Rapid Science 1998     

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.     

A review of adaptive cluster sampling: 1990–2003

Adaptive cluster sampling (ACS) is an adaptive sampling scheme which operates under the rule that when the observed value of an initially selected sampling unit satisfies some condition of interest, C, other additional units in some pre-defined accompanying neighborhood are also added to the sample. In turn, if any of these additional units satisfy C, then their corresponding unit neighborhoods are added to the sample as well, and so on. This process stops when no additional units satisfying C are encountered. This paper will provide a review of the major developments and issues in ACS since its introduction by Thompson (1990) [Journal of the American Statistical Association, 85, 1050–1059].     

Using Niche-Based Models to Improve the Sampling of Rare Species

Abstract:  Because data on rare species usually are sparse, it is important to have efficient ways to sample additional data. Traditional sampling approaches are of limited value for rare species because a very large proportion of randomly chosen sampling sites are unlikely to shelter the species. For these species, spatial predictions from niche-based distribution models can be used to stratify the sampling and increase sampling efficiency. New data sampled are then used to improve the initial model. Applying this approach repeatedly is an adaptive process that may allow increasing the number of new occurrences found. We illustrate the approach with a case study of a rare and endangered plant species in Switzerland and a simulation experiment. Our field survey confirmed that the method helps in the discovery of new populations of the target species in remote areas where the predicted habitat suitability is high. In our simulations the model-based approach provided a significant improvement (by a factor of 1.8 to 4 times, depending on the measure) over simple random sampling. In terms of cost this approach may save up to 70% of the time spent in the field.     

The application of adaptive cluster sampling for rare subtidal macroalgae

Adaptive cluster sampling (ACS) is a targeting sampling method that provides unbiased abundance estimators for populations of rare species that may be inadequately sampled with simple random sampling (SRS). ACS has been used successfully to estimate abundances of rockfish and sardine larvae from shipboard surveys. In this study, we describe the application of ACS for subtidal macroalgae. Using SCUBA, we measured abundances of Codium mamillosum, C. pomoides, and Halimeda cuneata at three islands and two levels of wave exposure. The three species were relatively patchy and could be sampled with ACS at one site per dive. Their distributions differed among islands and with exposure to wave energy, with H. cuneata found at only one island. ACS is a useful tool for understanding the spatial distribution and abundance of populations of rare benthic species, but, as was the case in this study, may not be as efficient as sampling with SRS with comparable replication.     

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.     

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.

     

Application of adaptive cluster sampling to low-density populations of freshwater mussels

Freshwater mussels appear to be promising candidates for adaptive cluster sampling because they are benthic macroinvertebrates that cluster spatially and are frequently found at low densities. We applied adaptive cluster sampling to estimate density of freshwater mussels at 24 sites along the Cacapon River, WV, where a preliminary timed search indicated that mussels were present at low density. Adaptive cluster sampling increased yield of individual mussels and detection of uncommon species; however, it did not improve precision of density estimates. Because finding uncommon species, collecting individuals of those species, and estimating their densities are important conservation activities, additional research is warranted on application of adaptive cluster sampling to freshwater mussels. However, at this time we do not recommend routine application of adaptive cluster sampling to freshwater mussel populations. The ultimate, and currently unanswered, question is how to tell when adaptive cluster sampling should be used, i.e., when is a population sufficiently rare and clustered for adaptive cluster sampling to be efficient and practical? A cost-effective procedure needs to be developed to identify biological populations for which adaptive cluster sampling is appropriate.     

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.     

Evaluation of sampling methods for the estimation of structural indices in forest stands

The paper is about the accurate (i.e. unbiased and precise) and efficient estimation of structural indices in forest stands. We present SIAFOR, a computer programme for the calculation of four nearest-neighbour indices, which describe the spatial arrangement of tree positions, the distribution pattern of species, and the size differentiation between trees. The study uses SIAFOR as a sampling simulator in eight completely stem-mapped forest stands of varying area and structural complexity. We statistically evaluate two sample types (distance and plot sampling), comparing sampling error, bias and minimum sample size for index estimation. We introduce the concepts of measurement expansion factor (MEF) and design expansion factor (DEF) for the technical evaluation of sample type efficiency (optimal sample type). Results indicate that sampling error can reach high levels and that minimum sample sizes for index estimation often amply exceed the limit of 20% of tree density or 20 trees per species per hectare, that we set as the highest feasible sample size in normal situations. We found clear feasibility limits (in terms of minimal tree densities and reachable accuracy levels) for the estimation of all investigated indices. Generally, equal or higher sample sizes are needed for plot sampling than for distance sampling to reach equal accuracy levels. Nevertheless, plot sampling resulted more efficient for the estimation of tree size differentiation at low to medium accuracy levels. For all other investigated indices distance sampling resulted more efficient than plot sampling. Minimum sample size increases with accuracy and is negatively correlated with tree density. At a given accuracy level minimum sample size is highest for the estimation of relative mingling and lowest for tree size differentiation; furthermore it is generally lower in large stands than in small ones. Because of the consistency of our conclusions in all of the investigated stands, we think they apply in most stands of similar area (between 1 and 10 ha) and species diversity (not more than four species).     

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.     

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.     

Comparison of a two-stage sampling design and its composite sample alternative: An application to soil studies

The objective of a long-term soil survey is to determine the mean concentrations of several chemical parameters for the pre-defined soil layers and to compare them with the corresponding values in the past. A two-stage random sampling procedure is used to achieve this goal. In the first step, n subplots are selected from N subplots by simple random sampling without replacement; in the second step, m sampling sites are chosen within each of the n selected subplots. Thus n · m soil samples are collected for each soil layer. The idea of the composite sample design comes from the challenge of reducing very expensive laboratory analyses: m laboratory samples from one subplot and one soil layer are physically mixed to form a composite sample. From each of the n selected subplots, one composite sample per soil layer is analyzed in the laboratory, thus n per soil layer in total. In this paper we show that the cost is reduced by the factor m — 1 when instead of the two-stage sampling its composite sample alternative is used; however, the variance of the composite sample mean is increased. In the case of positive intraclass correlation the increase is less than 12.5%; in the case of negative intraclass correlation the increase depends on the properties of the variable as well. For the univariate case we derive the optimal number of subplots and sampling sites. A case study is discussed at the end.     

Comparison between Hansen-Hurwitz and Horvitz-Thompson estimators for adaptive cluster sampling

Thompson (1990) introduced the adaptive cluster sampling design. This sampling design has been shown to be a useful sampling method for parameter estimation of a clustered and scattered population (Roesch, 1993; Smith et al., 1995; Thompson and Seber, 1996). Two estimators, the modified Hansen-Hurwitz (HH) and Horvitz-Thompson (HT) estimators, are available to estimate the mean or total of a population. Empirical results from previous researches indicate that the modified HT estimator has smaller variance than the modified HH estimator. We analytically compare the properties of these two estimators. Some results are obtained in favor of the modified HT estimator so that practitioners are strongly recommended to use the HT estimator despite easiness of computations for the HH estimator.     

Variation among populations of the sailfin molly in the rate of concurrent multiple paternity and its implications for mating-system evolution

We examined patterns of concurrent multiple mating in a live-bearing poeciliid fish, the sailfin molly (Poecilia latipinna). We tested whether the probability of multiple paternity was related to female body size or fertility and whether the rate of multiple paternity varied among four populations that differed in their distributions of female body size and fertility. We analyzed data on mother and offspring genotypes for three polymorphic allozymes by three techniques, including a maximum-likelihood estimator that accounts for sampling error in both parental and offspring allele frequencies. The estimated rate of multiple paternity varied between 0.09 and 0.85, and the rate in one population varied seasonally between 0.33 (spring) and 0.85 (autumn). The variation in these rates was not associated with variation in body-size distributions among populations but was closely associated with variation in size-specific fertility: populations with greater variation in female fertility had higher multiple-paternity rates. Within two populations, logistic regression revealed that individual females of larger body size and greater size-specific fertility were more likely to carry multiply sired broods. This result is consistent with observations made in one of the populations 5 years earlier. In general, the results strongly suggest that the mating system varies markedly among conspecific populations of sailfin mollies and that larger, more fertile females are the objects of intermale competition. Received: 6 May 1996 / Accepted in revised form: 5 December 1996     

A hierarchical model for spatial capture-recapture data

Estimating density is a fundamental objective of many animal population studies. Application of methods for estimating population size from ostensibly closed populations is widespread, but ineffective for estimating absolute density because most populations are subject to short-term movements or so-called temporary emigration. This phenomenon invalidates the resulting estimates because the effective sample area is unknown. A number of methods involving the adjustment of estimates based on heuristic considerations are in widespread use. In this paper, a hierarchical model of spatially indexed capture-recapture data is proposed for sampling based on area searches of spatial sample units subject to uniform sampling intensity. The hierarchical model contains explicit models for the distribution of individuals and their movements, in addition to an observation model that is conditional on the location of individuals during sampling. Bayesian analysis of the hierarchical model is achieved by the use of data augmentation, which allows for a straightforward implementation in the freely available software WinBUGS. We present results of a simulation study that was carried out to evaluate the operating characteristics of the Bayesian estimator under variable densities and movement patterns of individuals. An application of the model is presented for survey data on the flat-tailed horned lizard (Phrynosoma mcallii) in Arizona, USA.     

Sample based estimation of landscape metrics; accuracy of line intersect sampling for estimating edge density and Shannon’s diversity index

A recent trend is to estimate landscape metrics using sample data and cost-efficiency is one important reason for this development. In this study, line intersect sampling (LIS) was used as an alternative to wall-to-wall mapping for estimating Shannon’s diversity index and edge length and density. Monte Carlo simulation was applied to study the statistical performance of the estimators. All combinations of two sampling designs (random and systematic distribution of transects), four sample sizes, five transect configurations (straight line, L, Y, triangle, and quadrat), two transect orientations (fixed and random), and three configuration lengths were tested, each with a large number of simulations. Reference was 50 photos of size 1 km², already manually delineated in vector format by photo interpreters using GIS environment. The performance was compared by root mean square error (RMSE) and bias. The best combination for all three metrics was found to be the systematic design and as response design the straight line configuration with random orientation of transects, with little difference between the fixed and random orientation of transects. The rate of decrease of RMSE for increasing sample size and line length was studied with a mixed linear model. It was found that the RMSE decreased to a larger degree with the systematic design than the random one, especially with increasing sample size. Due to the nonlinearity in the definition of Shannon diversity estimator its estimator has a small and negative bias, decreasing with sample size and line length. Finally, a time study was conducted, measuring the time for registration of line intersections and their lengths on non-delineated aerial photos. The time study showed that long sampling lines were more cost-efficient than short ones for photo-interpretation.     

Using the bootstrap with two-phase adaptive stratified samples from multiple populations at multiple locations

Recently the two-phase adaptive stratified sampling design proposed by Francis (1984) has been extended by Manly et al. (2002) for situations where several biological populations are sampled simultaneously, and where this is done at several different geographical locations in order to estimate population totals or means. The method uses the results from a first phase sample to decide how best to allocate a second phase sample to locations and strata, in order to maximise a criterion (based on estimated coefficients of variation) that measures the accuracy of estimation for population totals, for all variables at all locations. One potential problem with this method is bias in the estimators of the population totals and means. In this paper bootstrapping is considered as a means of overcoming these biases. It is shown using model populations of Pacific walrus and shellfish, based on real data, that bootstrapping is a useful tool for removing about half of the bias. This is also confirmed from some simulations using artificial data.