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.     

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.     

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.     

Randomly selected order statistics in ranked set sampling: A less expensive comparable alternative to simple random sampling

Rank-based sampling designs are powerful alternatives to simple random sampling (SRS) and often provide large improvements in the precision of estimators. In many environmental, ecological, agricultural, industrial and/or medical applications the interest lies in sampling designs that are cheaper than SRS and provide comparable estimates. In this paper, we propose a new variation of ranked set sampling (RSS) for estimating the population mean based on the random selection technique to measure a smaller number of observations than RSS design. We study the properties of the population mean estimator using the proposed design and provide conditions under which the mean estimator performs better than SRS and some existing rank-based sampling designs. Theoretical results are augmented with some numerical studies and a real-life example, where we also study the performance of our proposed design under perfect and imperfect ranking situations.     

Multivariate bootstrap confidence regions for abundance vector using

Abundance vector estimation is a well investigated problem in statistical ecology. The use of simple random sampling with replacement or replicated sampling ensures good asymptotic properties of the abundance vector estimators. However, real surveys are based on small sample sizes, and assuming any specific distribution of the abundance vector estimator may be hazardous.In this paper we focus our attention on situations where the population is not too large and the sample size is small. We propose bootstrap multivariate confidence regions based on data depth. Data depth is a geometrical concept of ordering data from the center outwardly in higher dimensions. The Simplicial depth, the Tukey's depth and the Mahalanobis depth are presented. In order to build confidence regions in the presence of a skewed distribution of the abundance vector estimator, the use of Tukey's depth is suggested. The proposed method has been applied to the benthic community of Lake Lesina. A comparison with Mahalanobis depth and standard existing methods is reported.     

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.     

Finite population properties of individual predictors based on spatial patterns

In this study a conceptual framework for assessing the statistical properties of a non-stochastic spatial interpolator is developed through the use of design-based finite population inference tools. By considering the observed locations as the result of a probabilistic sampling design, we propose a standardized weighted predictor for spatial data starting from a deterministic interpolator that usually does not provide uncertainty measures. The information regarding the coordinates of the spatial locations is known at the population level and is directly used in constructing the weighting system. Our procedure captures the spatial pattern by means of the Euclidean distances between locations, which are fixed and do not require any further assessment after the sample has been drawn. The predictor for any individual value turns in a ratio of design-based random quantities. We illustrate the predictor design-based statistical properties, i.e. asymptotically p-unbiasedness and p-consistency, for simple random sampling without replacement. An application to a couple of environmental datasets is presented, for assessing predictor performances in correspondence of different population characteristics. A comparison with the equivalent non-spatial predictor is presented.     

Ranked set sampling allocation models for multiple skewed variables: an application to agricultural data

The mean of a balanced ranked set sample is more efficient than the mean of a simple random sample of equal size and the precision of ranked set sampling may be increased by using an unbalanced allocation when the population distribution is highly skewed. The aim of this paper is to show the practical benefits of the unequal allocation in estimating simultaneously the means of more skewed variables through real data. In particular, the allocation rule suggested in the literature for a single skewed distribution may be easily applied when more than one skewed variable are of interest and an auxiliary variable correlated with them is available. This method can lead to substantial gains in precision for all the study variables with respect to the simple random sampling, and to the balanced ranked set sampling too.     

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.     

Spatially stratified sampling using auxiliary information for geostatistical mapping

This paper presents a method of spatial sampling based on stratification by Local Moran’s I _i calculated using auxiliary information. The sampling technique is compared to other design-based approaches including simple random sampling, systematic sampling on a regular grid, conditional Latin Hypercube sampling and stratified sampling based on auxiliary information, and is illustrated using two different spatial data sets. Each of the samples for the two data sets is interpolated using regression kriging to form a geostatistical map for their respective areas. The proposed technique is shown to be competitive in reproducing specific areas of interest with high accuracy.     

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.     

Habitat association studies in conjunction with adaptive cluster samples

Habitat association studies investigate the relationships between habitat characteristics and animal usage of study regions. These studies are often conducted in conjunction with surveys designed primarily to estimate population totals. This paper shows that habitat association studies may proceed from surveys using adaptive cluster sampling. The manner in which units appear in the sample turns out not be relevant to the habitat association study, which proceeds as though the units came from a simple random sample. However, it is also shown that the information about the habitat association parameters is greater than one would expect from a simple random sample of the same general size.     

Nonparametric Ranked-set Sampling Confidence Intervals for Quantiles of a Finite Population

Ranked-set sampling from a finite population is considered in this paper. Three sampling protocols are described, and procedures for constructing nonparametric confidence intervals for a population quantile are developed. Algorithms for computing coverage probabilities for these confidence intervals are presented, and the use of interpolated confidence intervals is recommended as a means to approximately achieve coverage probabilities that cannot be achieved exactly. A simulation study based on finite populations of sizes 20, 30, 40, and 50 shows that the three sampling protocols follow a strict ordering in terms of the average lengths of the confidence intervals they produce. This study also shows that all three ranked-set sampling protocols tend to produce confidence intervals shorter than those produced by simple random sampling, with the difference being substantial for two of the protocols. The interpolated confidence intervals are shown to achieve coverage probabilities quite close to their nominal levels. Rankings done according to a highly correlated concomitant variable are shown to reduce the level of the confidence intervals only minimally. An example to illustrate the construction of confidence intervals according to this methodology is provided.     

A robust weighted EM algorithm for use-availability data

Wildlife sampling for habitat selection often combines a random background sample with a random sample of used sites, because the background sample could contain too few used sites to be informative for rare species. This approach is referred to as use-availability sampling. Two variants are considered where there is: (1) a random background sample including used and unused sites augmented with a sample of used sites, and (2) a sample of used sites augmented with a contaminated background sample, i.e. use is not recorded. A weighted estimator first proposed by Manski and Lerman (Econometrica 45(8):1977?C1988, 1977) forms the basis for our suggested approach. The weighted estimator has been shown to perform better than the usual unweighted approach with uncontaminated data and mis-specified logit models (Xie and Manski in Sociol Methods Res 17(3):283?C302, 1989). A weighted EM algorithm is developed for use with contaminated background data. We show that the weighted estimator continues to perform well with contaminated data and maintains its robustness to model mis-specification. The weighted estimator has not been previously used for use-availability sampling due to reliance on the assumption that only the intercept is biased, which is valid for a correct logit model. We show that adjusting the intercept may not eliminate the bias with an incorrect logit model. In this case, the weighted estimator is a relatively simple and effective alternative.     

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.     

Will Observation Error and Biases Ruin the Use of Simple Extinction Models?

Abstract: Estimating the risk of extinction for populations of endangered species is an important component of conservation biology. These estimates must be made from data that contain both environmental noise in the year-to-year transitions in population size (so-called "process error"), random errors in sampling, and possible biases in sampling ( both forms of observation errors). To determine how much faith to place in estimated extinction rates, it is important to know how sensitive they are to observation error. We used three simple, commonly employed models of population dynamics to generate simulated population time series. We then combined random observation error or systematic biases with those data, fit models to the time series data, and observed how close the extinction dynamics of the fitted models compared with the dynamics of the underlying models. We found that systematic biases in sampling rarely affected estimates of extinction risk. We also found that even moderate levels of random observation error do not significantly affect extinction estimates except over a small range of process errors, corresponding to the region where extinction risk is most uncertain. With more substantial sampling error, estimates of extinction risk degraded rapidly. Field census techniques for a variety of taxa often involve observation errors within ±32% of actual population sizes. For typical time series used in conservation, therefore, we often may not need to be overly concerned about observation errors as an extra source of imperfection in our estimated extinction rates.     

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.     

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.     

Estimating a distance dependent contagion function using point sample data

Natural events and human activities cause changes in landscape structure. Landscape metrics are used as a useful tool to study landscape trends and ecological processes related to the landscape structure. These metrics are commonly calculated on wall-to-wall raster data from remote sensing. A recent trend is to use sample data to estimate landscape metrics. In this study, point sampling was used to estimate a vector-based and distance dependent contagion metric. The metric is an extension of the established contagion. The statistical properties, for both unconditional and conditional contagions, were assessed by a point (point pairs) sampling experiment in maps from the National Inventory of landscapes in Sweden. Random and systematic sampling designs were tested for nine point distances and five sample sizes and for two classification systems. The systematic design showed slightly smaller root mean square error (RMSE) and bias than the random design. Both true and estimated values were calculated using computer programs in FORTRAN, which was specifically written for the purpose of the study. For a given sample size, RMSE and bias increased with increasing point distance. The estimator of unconditional contagion had acceptable RMSE and bias for moderate sample sizes, but in the conditional case the bias (and thus the RMSE) was unacceptably large. The main reason for this is that small classes (by area) affect both the true value of the contagion and are often missing in the sample. The method proposed can be adopted in gradient-based model of landscape structure where no distinct border is assumed between polygons. The method can also be applied in field-based inventories.     

Distribution-free ranked-set sample procedures for umbrella alternatives in the m-sample setting

In settings where measurements are costly and/or difficult to obtain but ranking of the potential sample data is relatively easy and reliable, the use of statistical methods based on a ranked-set sampling approach can lead to substantial improvement over analogous methods associated with simple random samples. Previous nonparametric work in this area has been concentrated almost exclusively on the one- and two-sample location problems. In this paper we develop ranked-set sample procedures for the m-sample location setting where the treatment effect parameters follow a restricted umbrella pattern. Distribution-free testing procedures are developed for both the case where the peak of the umbrella is known and for the case where it is unknown. Small sample and asymptotic null distribution properties are provided for the peak-known test statistic.