A Bayesian-weighted approach to predicting the number of newly discovered rare species

In natural ecological communities, most species are rare and thus susceptible to extinction. Consequently, the prediction and identification of rare species are of enormous value for conservation purposes. How many newly found species will be rare in the next field survey? We took a Bayesian viewpoint and used observed species abundance information in an ecological sample to develop an accurate way to estimate the number of new rare species (e.g., singletons, doubletons, and tripletons) in an additional unknown sample. A similar method has been developed for incidence-based data sets. Five seminumerical tests (3 abundance cases and 2 incidence cases) showed that our proposed Bayesian-weight estimator accurately predicted the number of new rare species with low relative bias and low relative root mean squared error and, accordingly, high accuracy. Finally, we applied the proposed estimator to 6 conservation-directed empirical data sets (3 abundance cases and 3 incidence cases) and found the prediction of new rare species was quite accurate; the 95% CI covered the true observed value very well in most cases. Our estimator performed similarly to or better than an unweighted estimator derived from Chao et al. and performed consistently better than the naïve unweighted estimator. We recommend our Bayesian-weight estimator for conservation applications, although the unweighted estimator of Chao et al. may be better under some circumstances. We provide an R package RSE (r are s pecies e stimation) at https://github.com/ecomol/RSE for implementation of the estimators.     

A combination line transect and capture-recapture sampling model for multiple observers in aerial surveys

We present a robust sampling methodology to estimate population size using line transect and capture-recapture procedures for aerial surveys. Aerial surveys usually underestimate population density due to animals being missed. A combination of capture-recapture and line transect sampling methods with multiple observers allows violation of the assumption that all animals on the centreline are sighted from the air. We illustrate our method with an example of inanimate objects which shows evidence of failure of the assumption that all objects on the centreline have probability 1 of being detected. A simulation study is implemented to evaluate the performance of three variations of the Lincoln-Petersen estimator: the overall estimator, the stratified estimator, and the general stratified estimator based on the combined likelihood proposed in this paper. The stratified Lincoln-Petersen estimator based on the combined likelihood is found to be generally superior to the other estimators.     

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.     

An incidence-based richness estimator for quadrats sampled without replacement

Most richness estimators currently in use are derived from models that consider sampling with replacement or from the assumption of infinite populations. Neither of the assumptions is suitable for sampling sessile organisms such as plants where quadrats are often sampled without replacement and the area of study is always limited. In this paper, we propose an incidence-based parametric richness estimator that considers quadrat sampling without replacement in a fixed area. The estimator is derived from a zero-truncated binomial distribution for the number of quadrats containing a given species (e.g., species i) and a modified beta distribution for the probability of presence-absence of a species in a quadrat. The maximum likelihood estimate of richness is explicitly given and can be easily solved. The variance of the estimate is also obtained. The performance of the estimator is tested against nine other existing incidence-based estimators using two tree data sets where the true numbers of species are known. Results show that the new estimator is insensitive to sample size and outperforms the other methods as judged by the root mean squared errors. The superiority of the new method is particularly noticeable when large quadrat size is used, suggesting that a few large quadrats are preferred over many small ones when sampling diversity.     

Nonparametric estimation of Shannon’s index of diversity when there are unseen species in sample

A biological community usually has a large number of species with relatively small abundances. When a random sample of individuals is selected and each individual is classified according to species identity, some rare species may not be discovered. This paper is concerned with the estimation of Shannons index of diversity when the number of species and the species abundances are unknown. The traditional estimator that ignores the missing species underestimates when there is a non-negligible number of unseen species. We provide a different approach based on unequal probability sampling theory because species have different probabilities of being discovered in the sample. No parametric forms are assumed for the species abundances. The proposed estimation procedure combines the Horvitz–Thompson (1952) adjustment for missing species and the concept of sample coverage, which is used to properly estimate the relative abundances of species discovered in the sample. Simulation results show that the proposed estimator works well under various abundance models even when a relatively large fraction of the species is missing. Three real data sets, two from biology and the other one from numismatics, are given for illustration.     

On parametric estimation of population abundance for line transect sampling

Despite recent advances in nonparametric methods for estimating animal abundance, parametric methods are still used widely among biometricians due to their simplicity. In this paper, we propose an optimal shrinkage-type estimator and an empirical Bayes estimator for estimating animal density from line transect sampling data. The performances of the proposed estimators are compared with those of the maximum likelihood estimator and a bias-corrected maximum likelihood estimator both theoretically and numerically. Simulation results show that the optimal shrinkage-type estimator works the best if the detection function has a very thin tail (for example, the half normal detection function), while the maximum likelihood estimator is the best estimator if the detection function has relatively thick tail (for example, the polynomial detection function).     

The use of replicated plot,line and point sampling for estimating species abundance and ecological diversity

The paper deals with the problem of estimating diversity indexes for an ecological community. First the species abundances are unbiasedly and consistently estimated using designs based on n random and independent selections of plots, points or lines over the study area. The problem of sampling elusive populations is also considered. Finally, the diversity index estimates are obtained as functions of the abundance estimates. The resulting estimators turn out to be asymptotically (n large) unbiased, even if a considerable bias may occur for a small n. Accordingly, the method of jackknifing is made use of in order to reduce bias.     

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.     

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.     

Line transect sampling of Karner blue butterflies ( Lycaeides melissa samuelis)

Line transect sampling is an effective survey method for estimating butterfly densities because it provides unbiased estimates of site-density (provided key assumptions are met), and estimates are comparable among sites. For monitoring Karner blue butterflies in Wisconsin, USA, comparable estimates are required because each year a different selection of sites will be monitored. Annual state-wide indices of species abundance can be derived from the site-surveys and compared to previous year's indices to monitor trends. We advocate that line transect sampling is preferable to Pollard-Yates transects as a survey technique for monitoring Karner blue butter- flies. The Pollard-Yates surveys do not adjust for diferences in site detectability. As a consequence, estimates of among-site from Pollard-Yates surveys can be biased. © Rapid Science 1998