Stratified Double Sampling with Continuous Outcomes: Design and Analysis

We provide a method for estimating the sample mean of a continuous outcome in a stratified population using a double sampling scheme. The stratified sample mean is a weighted average of stratum specific means. It is assumed that the fallible and true outcome data are related by a simple linear regression model in each stratum. The optimal stratified double sampling plan, i.e. , the double sampling plan that minimizes the cost of sampling for fixed variances, or alternatively, minimizes the variance for fixed costs, is found and compared to a standard sampling plan. The design parameters are the total sample size and the number of doubly sampled units in each stratum. We show that the optimal double sampling plan is a function of the between-strata and within-strata cost and variance ratios. The efficiency gains, relative to standard sampling plans, under broad set of conditions, are considerable.     

Paired double-ranked set sampling

Abstract

In environmental monitoring and assessment, the main focus is to achieve observational economy and to collect data with unbiased, efficient and cost-effective sampling methods. Ranked set sampling (RSS) is one traditional method that is mostly used for accomplishing observational economy. In this article, we propose an unbiased sampling scheme, named paired double RSS (PDRSS) for estimating the population mean. We study the performance of the mean estimators under PDRSS based on perfect and imperfect rankings. It is shown that, for perfect ranking, the variance of the mean estimator under PDRSS is always less than the variance of mean estimator based on simple random sampling, paired RSS and RSS. The mean estimators under RSS, median RSS, PDRSS, and double RSS are also compared with the regression estimator of population mean based on SRS. The procedure is also illustrated with a case study using a real data set.     

Efficient ascertainment schemes for maximum likelihood estimation

While well chosen sampling schemes may substantially increase efficiency of observational studies, some sampling schemes may instead decrease efficiency. Rules of thumb how to choose sampling schemes are only available for some special cases. In this paper we provide tools to compare efficiencies, and cost adjusted efficiencies, of different sampling schemes, in order to facilitate this choice. The method can be used for both categorical and continuous outcome variables. Some examples are presented, focusing on data from ascertainment sampling schemes. A Monte Carlo method is used to overcome computational issues wherever needed. The results are illustrated in graphs.     

A new sampling method for estimating the population mean

A double L ranked set sampling (DLRSS) method is suggested for estimating the population mean. The DLRSS is compared with the simple random sampling (SRS), ranked set sampling (RSS) and L ranked set sampling (LRSS) methods based on the same number of measured units. The conditions for which the suggested estimator performs better than the other estimators are derived. It is found that, the suggested DLRSS estimator is an unbiased of the population mean, and is more efficient than its counterparts using SRS, RSS, and LRSS methods. Real data sets are used for illustration.     

Estimation of mean based on modified robust extreme ranked set sampling

In this paper, double robust extreme ranked set sampling (DRERSS) and its properties for estimating the population mean are considered. It turns out that, when the underlying distribution is symmetric, DRERSS gives unbiased estimators of the population mean. Also, it is found that DRERSS is more efficient than the simple random sampling (SRS), ranked set sampling (RSS), and extreme ranked set sampling (ERSS) methods. For asymmetric distributions considered in this study, the DRERSS has a small bias and it is more efficient than SRS, RSS, and ERSS. A real data set is used to illustrate the DRERSS method.     

Characterization of MIL-STD-105D double sampling plans using normalized OC curves

A simple method of analyzing the MIL-STD-105D double sampling plans by ‘normalizing the OC curve’ is presented. Among a set of competing double sampling plans, procedures are given to select a particular plan when costs are unknown. Necessary tables are constructed and examples given.     

Regression analysis for a summed missing data problem under an outcome‐dependent sampling scheme

In this paper, we consider a regression analysis for a missing data problem in which the variables of primary interest are unobserved under a general biased sampling scheme, an outcome‐dependent sampling (ODS) design. We propose a semiparametric empirical likelihood method for accessing the association between a continuous outcome response and unobservable interesting factors. Simulation study results show that ODS design can produce more efficient estimators than the simple random design of the same sample size. We demonstrate the proposed approach with a data set from an environmental study for the genetic effects on human lung function in COPD smokers. The Canadian Journal of Statistics 40: 282–303; 2012 © 2012 Statistical Society of Canada     

Maximum likelihood estimation for outcome‐dependent samples

In outcome‐dependent sampling, the continuous or binary outcome variable in a regression model is available in advance to guide selection of a sample on which explanatory variables are then measured. Selection probabilities may either be a smooth function of the outcome variable or be based on a stratification of the outcome. In many cases, only data from the final sample is accessible to the analyst. A maximum likelihood approach for this data configuration is developed here for the first time. The likelihood for fully general outcome‐dependent designs is stated, then the special case of Poisson sampling is examined in more detail. The maximum likelihood estimator differs from the well‐known maximum sample likelihood estimator, and an information bound result shows that the former is asymptotically more efficient. A simulation study suggests that the efficiency difference is generally small. Maximum sample likelihood estimation is therefore recommended in practice when only sample data is available. Some new smooth sample designs show considerable promise.     

On triple sampling scremes for estimating from binomial data with misclassification errors ∗

Previous work has been carried out on the use of double sampling schemes for inference from binomial data which are subject to misclassification. The double sampling scheme utilizes a sample of n units which are classified by both a fallible and a true device and another sample of n₂ units which are classified only by a fallible device. A triple sampljng scheme incorporates an additional sample of n_l units which are classified only by the true device. In this paper we apply this triple sampling to estimation from binomialdata. First estimation of a binomial proportion is discussed under different misclassification structures. Then, the problem of optimal allocation of sample sizes is discussed.     

Improved quality control charts for monitoring the process mean,using double-ranked set sampling methods

Statistical control charts are widely used in the manufacturing industry. The Shewhart-type control charts are developed to improve the monitoring process mean by using the double quartile-ranked set sampling, quartile double-ranked set sampling, and double extreme-ranked set sampling methods. In terms of the average run length, the performance of the proposed control charts are compared with the existing control charts based on simple random sampling, ranked set sampling and extreme-ranked set sampling methods. An application of real data is also considered to investigate the performance of the suggested process mean control charts. The findings of the study revealed that the newly suggested control charts are superior to the existing counterparts.     

Recent progresses in outcome-dependent sampling with failure time data

An outcome-dependent sampling (ODS) design is a retrospective sampling scheme where one observes the primary exposure variables with a probability that depends on the observed value of the outcome variable. When the outcome of interest is failure time, the observed data are often censored. By allowing the selection of the supplemental samples depends on whether the event of interest happens or not and oversampling subjects from the most informative regions, ODS design for the time-to-event data can reduce the cost of the study and improve the efficiency. We review recent progresses and advances in research on ODS designs with failure time data. This includes researches on ODS related designs like case–cohort design, generalized case–cohort design, stratified case–cohort design, general failure-time ODS design, length-biased sampling design and interval sampling design.     

Double and group acceptance sampling plan for truncated life test based on inverse log-logistic distribution

This paper introduces a double and group acceptance sampling plans based on time truncated lifetimes when the lifetime of an item follows the inverse log-logistic (ILL) distribution with known shape parameter. The operating characteristic function and average sample number (ASN) values of the double acceptance sampling inspection plan are provided. The values of the minimum number of groups and operating characteristic function for various quality levels are obtained for a group acceptance sampling inspection plan. A comparative study between single acceptance sampling inspection plan and double acceptance sampling inspection plan is carried out in terms of sample size. One simulated example and four real-life examples are discussed to show the applicability of the proposed double and group acceptance sampling inspection plans for ILL distributed quality parameters.     

Inference in generalized additive mixed modelsby using smoothing splines

Generalized additive mixed models are proposed for overdispersed and correlated data, which arise frequently in studies involving clustered, hierarchical and spatial designs. This class of models allows flexible functional dependence of an outcome variable on covariates by using nonparametric regression, while accounting for correlation between observations by using random effects. We estimate nonparametric functions by using smoothing splines and jointly estimate smoothing parameters and variance components by using marginal quasi-likelihood. Because numerical integration is often required by maximizing the objective functions, double penalized quasi-likelihood is proposed to make approximate inference. Frequentist and Bayesian inferences are compared. A key feature of the method proposed is that it allows us to make systematic inference on all model components within a unified parametric mixed model framework and can be easily implemented by fitting a working generalized linear mixed model by using existing statistical software. A bias correction procedure is also proposed to improve the performance of double penalized quasi-likelihood for sparse data. We illustrate the method with an application to infectious disease data and we evaluate its performance through simulation.     

A case-base sampling method for estimating recurrent event intensities

Case-base sampling provides an alternative to risk set sampling based methods to estimate hazard regression models, in particular when absolute hazards are also of interest in addition to hazard ratios. The case-base sampling approach results in a likelihood expression of the logistic regression form, but instead of categorized time, such an expression is obtained through sampling of a discrete set of person-time coordinates from all follow-up data. In this paper, in the context of a time-dependent exposure such as vaccination, and a potentially recurrent adverse event outcome, we show that the resulting partial likelihood for the outcome event intensity has the asymptotic properties of a likelihood. We contrast this approach to self-matched case-base sampling, which involves only within-individual comparisons. The efficiency of the case-base methods is compared to that of standard methods through simulations, suggesting that the information loss due to sampling is minimal.     

A note on calibration weightings for stratified double sampling with equal probability

This study proposes a more efficient calibration estimator for estimating population mean in stratified double sampling using new calibration weights. The variance of the proposed calibration estimator has been derived under large sample approximation. Calibration asymptotic optimum estimator and its approximate variance estimator are derived for the proposed calibration estimator and existing calibration estimators in stratified double sampling. Analytical results showed that the proposed calibration estimator is more efficient than existing members of its class in stratified double sampling. Analysis and evaluation are presented.     

A general procedure for estimating the mean using double sampling for stratification and multi-auxiliary information

This paper defines a general procedure for estimating the population mean of the study variate based on double sampling for stratification in presence of multi-auxiliary information. Classes of combined and separate estimators have been suggested and their properties are studied under large sample approximation. A class of unstratified double sampling estimators is also proposed with its properties. Asymptotic optimum estimators in the classes are identified with their approximate variance formulae. Further the proposed classes of estimators are compared with the corresponding class of estimators based on un-stratified double sampling. All findings are encouraging and support the soundness of the proposed procedure for mean estimation.     

New acceptance sampling plans based on life tests for Birnbaum–Saunders distributions

In this paper, we consider the Birnbaum–Saunders distribution as a life model to develop various acceptance sampling schemes based on the truncated life tests. We develop the double sampling plan and determine the design parameters satisfying both the producer's and consumer's risks simultaneously for the specified reliability levels in terms of the mean ratio to the specified life. We also propose a group sampling plan and determine the parameters by the above-mentioned two-point method. Tables are constructed for the proposed sampling plans and results are explained with examples.     

Estimation of mean using double sampling for stratification and multivariate auxiliary information

Several estimators for estimating the mean of a principal variable are proposed based on double sampling for stratification (DSS) and multivariate auxiliary information. The general properties of the proposed estimators are studied, search for optimum estimators is made and the proposed estimators are compared with the corresponding estimators based on unstratified double sampling (USDS).     

An Extension of Sampford's Method for Unequal Probability Sampling

Abstract. Sampford's unequal probability sampling method is extended to the case that the inclusion probabilities do not sum to an integer. In this case, the sampling outcome is left open for exactly one randomly chosen unit and that unit gets a new inclusion probability. Three applications are presented. Two of them challenge traditional sampling routines. The simple Pareto sampling design, which was introduced by Rosén in 1997, is also extended. The extended Pareto design is shown to be close to the extended Sampford design.     

A unified empirical likelihood approach for testing MCAR and subsequent estimation

For an estimation with missing data, a crucial step is to determine if the data are missing completely at random (MCAR), in which case a complete‐case analysis would suffice. Most existing tests for MCAR do not provide a method for a subsequent estimation once the MCAR is rejected. In the setting of estimating means, we propose a unified approach for testing MCAR and the subsequent estimation. Upon rejecting MCAR, the same set of weights used for testing can then be used for estimation. The resulting estimators are consistent if the missingness of each response variable depends only on a set of fully observed auxiliary variables and the true outcome regression model is among the user‐specified functions for deriving the weights. The proposed method is based on the calibration idea from survey sampling literature and the empirical likelihood theory.