Optimal sampling and estimation strategies under the linear model

In some cases model-based and model-assisted inferences canlead to very different estimators. These two paradigms are notso different if we search for an optimal strategy rather thanjust an optimal estimator, a strategy being a pair composedof a sampling design and an estimator. We show that, under alinear model, the optimal model-assisted strategy consists ofa balanced sampling design with inclusion probabilities thatare proportional to the standard deviations of the errors ofthe model and the Horvitz–Thompson estimator. If the heteroscedasticityof the model is 'fully explainable’ by the auxiliary variables,then this strategy is also optimal in a model-based sense. Moreover,under balanced sampling and with inclusion probabilities thatare proportional to the standard deviation of the model, thebest linear unbiased estimator and the Horvitz–Thompsonestimator are equal. Finally, it is possible to construct asingle estimator for both the design and model variance. Theinference can thus be valid under the sampling design and underthe model.     

Bayesian Estimation Using Ranked Set Sampling

Bayesian Estimation of the parameter of a distribution is considered using Ranked set sampling (RSS). It is shown that for at least one RSS plan, the Bayes estimator has smaller Bayes risk than the Bayes estimator using simple random sampling (SRS). Furthermore, for exponential family with conjugate prior, the Bayes estimator of the mean using balanced RSS dominates, in terms of its Bayes risk, the Bayes estimator of the mean using SRS. This procedure is used to estimate the average Milk yield of four hundreds and two sheep. The empirical efficiency supports the theoretical findings.     

Estimating Genetic Variability with Restriction Endonucleases

The estimation of the amount of sequence variation in samples of homologous DNA segments is considered. The data are assumed to have been obtained by restriction endonuclease digestion of the segments, from which the numbers and frequencies of the cleavage sites in the sample are determined. An estimator, p, of the proportion of sites that are polymorphic in the sample is derived without assuming any particular population genetic model for the evolution of the population. The estimator is very close to the EWENS, SPIELMAN and HARRIS (1981) estimator that was derived with the symmetric WRIGHT-FISHER neutral mode. ENGELS (1981) has also recently proposed an estimator of the same quantity, and he arrived at his estimator without assuming a particular population genetic model. The sampling variance of p and ENGELS' estimator are derived. It is found that the sampling variance of p is lower than the sampling variance of ENGELS' estimator. Also, the sampling variance of theta, an estimate of theta (=4Nu) is obtained for the symmetric WRIGHT-FISHER neutral model with free recombination and with no recombination.     

Combining multiple imputation and inverse-probability weighting

Two approaches commonly used to deal with missing data are multiple imputation (MI) and inverse-probability weighting (IPW). IPW is also used to adjust for unequal sampling fractions. MI is generally more efficient than IPW but more complex. Whereas IPW requires only a model for the probability that an individual has complete data (a univariate outcome), MI needs a model for the joint distribution of the missing data (a multivariate outcome) given the observed data. Inadequacies in either model may lead to important bias if large amounts of data are missing. A third approach combines MI and IPW to give a doubly robust estimator. A fourth approach (IPW/MI) combines MI and IPW but, unlike doubly robust methods, imputes only isolated missing values and uses weights to account for remaining larger blocks of unimputed missing data, such as would arise, e.g., in a cohort study subject to sample attrition, and/or unequal sampling fractions. In this article, we examine the performance, in terms of bias and efficiency, of IPW/MI relative to MI and IPW alone and investigate whether the Rubin's rules variance estimator is valid for IPW/MI. We prove that the Rubin's rules variance estimator is valid for IPW/MI for linear regression with an imputed outcome, we present simulations supporting the use of this variance estimator in more general settings, and we demonstrate that IPW/MI can have advantages over alternatives. IPW/MI is applied to data from the National Child Development Study.     

A Stochastic Regression Model for General Trend Analysis of Longitudinal Continuous Data

A predictive continuous time model is developed for continuous panel data to assess the effect of time‐varying covariates on the general direction of the movement of a continuous response that fluctuates over time. This is accomplished by reparameterizing the infinitesimal mean of an Ornstein–Uhlenbeck processes in terms of its equilibrium mean and a drift parameter, which assesses the rate that the process reverts to its equilibrium mean. The equilibrium mean is modeled as a linear predictor of covariates. This model can be viewed as a continuous time first‐order autoregressive regression model with time‐varying lag effects of covariates and the response, which is more appropriate for unequally spaced panel data than its discrete time analog. Both maximum likelihood and quasi‐likelihood approaches are considered for estimating the model parameters and their performances are compared through simulation studies. The simpler quasi‐likelihood approach is suggested because it yields an estimator that is of high efficiency relative to the maximum likelihood estimator and it yields a variance estimator that is robust to the diffusion assumption of the model. To illustrate the proposed model, an application to diastolic blood pressure data from a follow‐up study on cardiovascular diseases is presented. Missing observations are handled naturally with this model.     

Additional Circular Systematic Sampling Methods

Several systematic sampling methods have been used to estimate the population mean when size of the population is a multiple of sample size. Among these, only few methods have been extended and used to estimate mean of the population when its size is not a multiple of sample size. In this paper, new methods called balanced circular systematic sampling and centered circular systematic sampling are introduced by extending balanced systematic sampling method and centered systematic sampling method respectively. These methods are compared with circular systematic sampling using average variance of corrected sample means for populations exhibiting approximate linear and parabolic trends. The suggested methods are found suitable to estimate the population mean.     

A Class of Ratio Cum Product-Type Estimator in Double Sampling

A class of ratio cum product-type estimator is proposed in case of double sampling in the present paper. Its bias and variance to the first order of approximation are obtained. For an appropriate weight ‘a’ and a good range of α-values, it is found that the proposed estimator is more efficient than the set of estimator viz., simple mean estimator, usual ratio and product estimators, SRIVASTAVA 's estimator (1967), CHAKARBARTY 's estimator and product-type estimator, which are in fact the particular cases of it. The proposed estimator is as efficient as linear regression estimator in double sampling at optimum value of α.     

Estimation of a Secondary Parameter in a Group Sequential Clinical Trial

In this article, we investigate estimation of a secondary parameter in group sequential tests. We study the model in which the secondary parameter is the mean of the normal distribution in a subgroup of the subjects. The bias of the naive secondary parameter estimator is studied. It is shown that the sampling proportions of the subgroup have a crucial effect on the bias: As the sampling proportion of the subgroup at or just before the stopping time increases, the bias of the naive subgroup parameter estimator increases as well. An unbiased estimator for the subgroup parameter and an unbiased estimator for its variance are derived. Using simulations, we compare the mean squared error of the unbiased estimator to that of the naive estimator, and we show that the differences are negligible. As an example, the methods of estimation are applied to an actual group sequential clinical trial, The Beta-Blocker Heart Attack Trial.     

A practical sampling design for acoustic surveys of bats

Acoustic surveys are widely used for describing bat occurrence and activity patterns and are increasingly important for addressing concerns for habitat management, wind energy, and disease on bat populations. Designing these surveys presents unique challenges, particularly when a probabilistic sample is required for drawing inference to unsampled areas. Sampling frame errors and other logistical constraints often require survey sites to be dropped from the sample and new sites added. Maintaining spatial balance and representativeness of the sample when these changes are made can be problematic. Spatially balanced sampling designs recently developed to support aquatic surveys along rivers provide solutions to a number of practical challenges faced by bat researchers and allow for sample site additions and deletions, support unequal-probability selection of sites, and provide an approximately unbiased local neighborhood-weighted variance estimator that is efficient for spatially structured populations such as is typical for bats. We implemented a spatially balanced design to survey canyon bat (Parastrellus hesperus) activity along a stream network. The spatially balanced design accommodated typical logistical challenges and yielded a 25% smaller estimated standard error for the mean activity level than the usual simple random sampling estimator. Spatially balanced designs have broad application to bat research and monitoring programs and will improve studies relying on model-based inference (e.g., occupancy models) by providing flexibility and protection against violations of the independence assumption, even if design-based estimators are not used. Our approach is scalable and can be used for pre- and post-construction surveys along wind turbine arrays and for regional monitoring programs. © 2011 The Wildlife Society.     

On Use of Distinct Respondents in Randomized Response Surveys

It is well known for direct response surveys (DR), where the responses are obtained from the respondents directly, that the sample mean, based on distinct units of a simple random sample selected with replacement (SRSWR) method, is more efficient than the sample mean based on all the units including repetition. In this paper, it is shown that a linear unbiased estimator based on distinct units is inadmissible for estimating a finite population mean when the sample is selected by an arbitrary with replacement (WR) sampling scheme and the responses are obtained independently by some RR technique. Efficiencies for a few linear unbiased estimators are compared under SRSWR sampling.     

Analysis of cytoplasmic and maternal effects I. A genetic model for diploid plant seeds and animals

A genetic model for modified diallel crosses is proposed for estimating variance and covariance components of cytoplasmic, maternal additive and dominance effects, as well as direct additive and dominance effects. Monte Carlo simulations were conducted to compare the efficiencies of minimum norm quadratic unbiased estimation (MINQUE) methods. For both balanced and unbalanced mating designs, MINQUE (0/1), which has 0 for all the prior covariances and 1 for all the prior variances, has similar efficiency to MINQUE(), which has parameter values for the prior values. Unbiased estimates of variance and covariance components and their sampling variances could be obtained with MINQUE(0/1) and jackknifing. A t-test following jackknifing is applicable to test hypotheses for zero variance and covariance components. The genetic model is robust for estimating variance and covariance components under several situations of no specific effects. A MINQUE(0/1) procedure is suggested for unbiased estimation of covariance components between two traits with equal design matrices. Methods of unbiased prediction for random genetic effects are discussed. A linear unbiased prediction (LUP) method is shown to be efficient for the genetic model. An example is given for a demonstration of estimating variance and covariance components and predicting genetic effects.     

Statistics of natural populations. III. Sequential sampling plans for the estimation of gene frequencies

In previous work several models have been developed for genetic surveys of natural populations. Parents of unknown genotype are collected from a natural population, polymorphic at a single genetic locus. From each of these N cryptic parents a number of offspring are identified for their genotype. Our problem is to select an efficient offspring sampling plan for estimating the frequency of an allele in the cryptic adult population based on the N family profiles of juvenile genotypes. A criterion called the information per unit cost of observation is introduced to evaluate sequential sampling plans, in which the number of offspring per family examined is random. Some simple, practical schemes for stopping the sampling of offspring from a collected parent are introduced; one example is stopping when: (i) the offspring are definitive about the parental genotype(s) for the first time; (ii) a fixed number of one genotype only is seen; or (iii) a fixed maximum feasible number of offspring have been genotyped. This sampling scheme is recommended. For each sampling scheme, the best linear unbiased estimator and the sequential maximum likelihood estimator of the allele frequency are characterized. From the moments of these estimators, it is then possible to tabulate efficient sequential sampling plans, which are better (in the sense of information per unit cost), just as simple, and less costly than corresponding fixed sampling plans in use.     

A maximum-likelihood estimation of pairwise relatedness for autopolyploids

Relatedness between individuals is central to ecological genetics. Multiple methods are available to quantify relatedness from molecular data, including method-of-moment and maximum-likelihood estimators. We describe a maximum-likelihood estimator for autopolyploids, and quantify its statistical performance under a range of biologically relevant conditions. The statistical performances of five additional polyploid estimators of relatedness were also quantified under identical conditions. When comparing truncated estimators, the maximum-likelihood estimator exhibited lower root mean square error under some conditions and was more biased for non-relatives, especially when the number of alleles per loci was low. However, even under these conditions, this bias was reduced to be statistically insignificant with more robust genetic sampling. We also considered ambiguity in polyploid heterozygote genotyping and developed a weighting methodology for candidate genotypes. The statistical performances of three polyploid estimators under both ideal and actual conditions (including inbreeding and double reduction) were compared. The software package POLYRELATEDNESS is available to perform this estimation and supports a maximum ploidy of eight.     

Estimation of sample zooplankton abundance from Folsom splitter sub-samples

Procedures are described for estimating the abundance of eachzooplanklon species in a sample after counting at least twosub-samples from a Folsom splitter. A multinomiai model is assumedfor the splitting process. The method may be summarized as follows:two sub-samples, balanced with respect to the left — rightsplits, are counted. The counts for each species are testedfor homogeneity using a x² test. If the counts are homogeneous,an estimator is given which permits the estimation of the numberof animals in the sample, and the variance of this estimate.If the subsample counts are heterogeneous, it is assumed thatclumping has occurred. A procedure is described in which additionalsub-samples are counted to locate the clump, and estimatorsof abundance (and variance estimators) are derived. If the splitterbias is <5%, bias of the abundance estimator is negligiblecompared to the binomial sampling error. Comparisons of speciesabundance between plankton samples is made more rigorous usingthe estimates of mean abundance, and the variance of this mean,provided by the methods described here. ^*Present address: Dunoon Road, Dorroughby, NSW 2480, Australia.     

Estimation and comparison of changes in the presence of informative right censoring: conditional linear model

A general linear regression model for the usual least squares estimated rate of change (slope) on censoring time is described as an approximation to account for informative right censoring in estimating and comparing changes of a continuous variable in two groups. Two noniterative estimators for the group slope means, the linear minimum variance unbiased (LMVUB) estimator and the linear minimum mean squared error (LMMSE) estimator, are proposed under this conditional model. In realistic situations, we illustrate that the LMVUB and LMMSE estimators, derived under a simple linear regression model, are quite competitive compared to the pseudo maximum likelihood estimator (PMLE) derived by modeling the censoring probabilities. Generalizations to polynomial response curves and general linear models are also described.     

On Double Sampling Ratio-Cum-Product-Type Estimator with Independent Samples

There are two cases in double sampling; case(i) when the second sample is a sub-sample from preliminary large sample, and case(ii) when the second sample is not a sub-sample from the preliminary large sample. Recently SISODIA and DWIVEDI (1981) proposed a ratio cum product-type estimator in double sampling in which they have studied the properties of this estimator under case (i). In this paper, we have made an attempt to study the properties of the same estimator under case (ii). It is found that the estimator is superior than double sampling linear regression estimator, usual ratio estimator, product estimator and among others. The estimator is also compared with simple mean per unit for a given cost of the survey.     

Global robust power-rate stability of delayed genetic regulatory networks with noise perturbations

In this paper, by using the Lyapunov method, Itô’s differential formula and linear matrix inequality (LMI) approach, the global robust power-rate stability in mean square is discussed for genetic regulatory networks with unbounded time-varying delay, noise perturbations and parameter uncertainties. Sufficient conditions are given to ensure the robust power-rate stability (in mean square) of the genetic regulatory networks. Meanwhile, the criteria ensuring global power-rate stability in mean square are a byproduct of the criteria guaranteeing global robust power-rate stability in mean square. The obtained conditions are derived in terms of linear matrix inequalities (LMIs) which are easy to be verified via the LMI toolbox. An illustrative example is given to show the effectiveness of the obtained result.     

Recommendations for estimating mark rate of cetaceans in photo-ID research: A critique of field sampling protocols and variance estimation

Mark rate, or the proportion of the population with unique, identifiable marks, must be determined in order to estimate population size from photographic identification data. In this study we address field sampling protocols and estimation methods for robust estimation of mark rate and its uncertainty in cetacean populations. We present two alternatives for estimating the variance of mark rate: (1) a variance estimator for clusters of unequal sizes (SRCS) and (2) a hierarchical Bayesian model (SRCS-Bayes), and compare them to the simple random sampling (SRS) variance estimator. We tested these variance estimators using a simulation to see how they perform at varying mark rates, number of groups sampled, photos per group, and mean group sizes. The hierarchical Bayesian model outperformed the frequentist variance estimators, with the true mark rate of the population held in its 95% HDI 91.9% of the time (compared with coverage of 79% for the SRS method and 76.3% for the SRCS-Cochran method). The simulation results suggest that, ideally, mark rate and its precision should be quantified using hierarchical Bayesian modeling, and researchers should attempt to sample as many unique groups as possible to improve accuracy and precision.     

A note on robust variance estimation for cluster-correlated data

There is a simple robust variance estimator for cluster-correlated data. While this estimator is well known, it is poorly documented, and its wide range of applicability is often not understood. The estimator is widely used in sample survey research, but the results in the sample survey literature are not easily applied because of complications due to unequal probability sampling. This brief note presents a general proof that the estimator is unbiased for cluster-correlated data regardless of the setting. The result is not new, but a simple and general reference is not readily available. The use of the method will benefit from a general explanation of its wide applicability.     

The probability of parallel genetic evolution from standing genetic variation

Parallel evolution is often assumed to result from repeated adaptation to novel, yet ecologically similar, environments. Here, we develop and analyse a mathematical model that predicts the probability of parallel genetic evolution from standing genetic variation as a function of the strength of phenotypic selection and constraints imposed by genetic architecture. Our results show that the probability of parallel genetic evolution increases with the strength of natural selection and effective population size and is particularly likely to occur for genes with large phenotypic effects. Building on these results, we develop a Bayesian framework for estimating the strength of parallel phenotypic selection from genetic data. Using extensive individual‐based simulations, we show that our estimator is robust across a wide range of genetic and evolutionary scenarios and provides a useful tool for rigorously testing the hypothesis that parallel genetic evolution is the result of adaptive evolution. An important result that emerges from our analyses is that existing studies of parallel genetic evolution frequently rely on data that is insufficient for distinguishing between adaptive evolution and neutral evolution driven by random genetic drift. Overcoming this challenge will require sampling more populations and the inclusion of larger numbers of loci.