On two properties of an unequal probability sampling scheme

Summary For an inclusion probability proportional to size (IPPS) sampling scheme recently proposed by Saxena, Singh and Srivastava (1986), it is shown that under certain simple verifiable conditions (1) the Horvitz-Thompson (1952) estimator based on it has a smaller variance than the variance of the Hansen-Hurwitz (1943) estimator based on probability proportional to size (PPS) sampling with replacement (WR) both involving the same size-measures and the expected sample size in the former being equal to the number of draws in the latter and (2) the Yates-Grundy (1953) estimator for the variance of the Horvitz-Thompson estimator based on this IPPS scheme is uniformly non-negative.     

A note on Hájek’s theory of rejective sampling

Summary Supplementing results of Hájek (1960, 1964) sufficient conditions for Poisson convergence of the Horvitz-Thompson estimator in pps Rejective sampling are obtained.     

An invariance principle forU-statistics in simple random sampling without replacement

Summary Supplementing the well-known invariance principle forU-statistics based on i.i.d. observations (Miller and Sen 1972) we establish an invariance principle forU-statistics in case of simple random sampling from a sequence of finite populations. This generalizes the asymptotic normality-result of Nandi and Sen (1963), and permits e.g. to prove asymptotic normality ofU-statistics in the presence of random non-response.     

On the application of the Minkowski-Farkas theorem to sampling designs

In the paper, we consider the following problem: Let {π_k} be a sequence satisfying 0π_kΣ1 (k=1,…, N) and π=n.Tben, is there an unordered sampling design such that, for each k=1,…N, the inclusion probability of unit k is equal to π? It is shown that it can be solved by the straightforward application of the Minkowski-Farkas theorem.     

Variance estimation for unequal probability sampling

In this paper, we discuss the optimality of the variance estimator of the Horvitz-Thompson estimator proposed by Kott (1988) in the class of model-unbiased quadratic estimators. We also propose some improved estimators over Kott's estimator in the class of general quadratic estimators. Received: February 1999     

Model assisted survey sampling strategy in two phases

Postulating a super-population regression model connecting a size variable, a cheaply measurable variable and an expensively observable variable of interest, an asymptotically optimal double sampling strategy to estimate the survey population total of the third variable is specified. To render it practicable, unknown model-parameters in the optimal estimator are replaced by appropriate statistics. The resulting generalized regression estimator is then shown to have a model-cum-asymptotic design based expected square error equal to that of the asymptotically optimum estimator itself. An estimator for design variance of the estimator is also proposed.     

Multinomial distributions applied to random sampling of particulate materials

When sampling a batch consisting of particulate material, the distribution of a sample estimator can be characterized using knowledge about the sample drawing process. With Bernoulli sampling, the number of particles in the sample is binomially distributed. Because this is rarely realized in practice, we propose a sampling design in which the possible samples have a nearly equal mass. Expected values and variances of the sample estimator are calculated. It is shown that the sample estimator becomes identical to the Horvitz–Thompson estimator in the case of a large batch-to-sample mass ratio and a large sample mass. Simulations and experiments were performed to test the theory. Simulations confirm that the round-off error due to the discrete nature of particles is negligible for large sample sizes. Sampling experiments were carried out with a mixture of PolyPropylene (PP) and PolyTetraFluorEthylene (PTFE) spheres suspended in a viscous medium. The measured and theoretical variations are in good agreement.     

Estimation of finite population variance using random non-response in survey sampling

In this paper, an estimator of finite population variance proposed by Isaki (1983) is studied under the two different situations of random non-response suggested by Tracy and Osahan (1994). A distribution is proposed for the number of sampling units on which information could not be obtained due to random non-response. The estimators for the mean square errors of the proposed strategies are also suggested. This paper was written while both authors were members of the Dept. of Econometrics, Monash University, Clayton 3168, Australia. This paper was presented on SISC—1996, Sydney, Australia. The opinions and results discussed in this paper are of authors and not necessarily of their institutes.     

Two–step sequential sampling

Choosing the sample size in advance is a familiar problem: often, additional observations appear to be desirable. The final sample size then becomes a random variable, which has rather serious consequences.
Two such sample extension situations will be considered here. In the first situation, the observed sample variance determines whether or not to double the original sample size. In the second situation, the variances observed in two independent samples are compared; their ratio determines the number of additional observations.     

Mean square error estimation in multi-stage sampling

Summary: Suppose for a homogeneous linear unbiased function of the sampled first stage unit (fsu)-values taken as an estimator of a survey population total, the sampling variance is expressed as a homogeneous quadratic function of the fsu-values. When the fsu-values are not ascertainable but unbiased estimators for them are separately available through sampling in later stages and substituted into the estimator, Raj (1968) gave a simple variance estimator formula for this multi-stage estimator of the population total. He requires that the variances of the estimated fsu-values in sampling at later stages and their unbiased estimators are available in certain `simple forms'. For the same set-up Rao (1975) derived an alternative variance estimator when the later stage sampling variances have more ‘complex forms’. Here we pursue with Raj's (1968) simple forms to derive a few alternative variance and mean square error estimators when the condition of homogeneity or unbiasedness in the original estimator of the total is relaxed and the variance of the original estimator is not expressed as a quadratic form.  We illustrate a particular three-stage sampling strategy and present a simulation-based numerical exercise showing the relative efficacies of two alternative variance estimators. Received: 19 February 1999     

Simultaneous estimation of hazard rates of several exponential populations

Suppose independent random samples are drawn from k (2) populations with a common location parameter and unequal scale parameters. We consider the problem of estimating simultaneously the hazard rates of these populations. The analogues of the maximum likelihood (ML), uniformly minimum variance unbiased (UMVU) and the best scale equivariant (BSE) estimators for the one population case are improved using Rao‐Blackwellization. The improved version of the BSE estimator is shown to be the best among these estimators. Finally, a class of estimators that dominates this improved estimator is obtained using the differential inequality approach.     

On Invariant Post‐randomization for Statistical Disclosure Control

In this paper, we investigate certain operational and inferential aspects of invariant Post‐randomization Method (PRAM) as a tool for disclosure limitation of categorical data. Invariant PRAM preserves unbiasedness of certain estimators, but inflates their variances and distorts other attributes. We introduce the concept of strongly invariant PRAM, which does not affect data utility or the properties of any statistical method. However, the procedure seems feasible in limited situations. We review methods for constructing invariant PRAM matrices and prove that a conditional approach, which can preserve the original data on any subset of variables, yields invariant PRAM. For multinomial sampling, we derive expressions for variance inflation inflicted by invariant PRAM and variances of certain estimators of the cell probabilities and also their tight upper bounds. We discuss estimation of these quantities and thereby assessing statistical efficiency loss from applying invariant PRAM. We find a connection between invariant PRAM and creating partially synthetic data using a non‐parametric approach, and compare estimation variance under the two approaches. Finally, we discuss some aspects of invariant PRAM in a general survey context.     

A note on the distribution of the product of zero‐mean correlated normal random variables

The problem of finding an explicit formula for the probability density function of two zero‐mean correlated normal random variables dates back to 1936. Perhaps, surprisingly, this problem was not resolved until 2016. This is all the more surprising given that a very simple proof is available, which is the subject of this note; we identify the product of two zero‐mean correlated normal random variables as a variance‐gamma random variable, from which an explicit formula for the probability density function is immediate.     

Using Complex Surveys to Estimate the L1‐Median of a Functional Variable: Application to Electricity Load Curves

Mean profiles are widely used as indicators of the electricity consumption habits of customers. Currently, in Électricité De France (EDF), class load profiles are estimated using point‐wise mean profiles. Unfortunately, it is well known that the mean is highly sensitive to the presence of outliers, such as one or more consumers with unusually high‐levels of consumption. In this paper, we propose an alternative to the mean profile: the L ₁ ‐ median profile which is more robust. When dealing with large data sets of functional data (load curves for example), survey sampling approaches are useful for estimating the median profile avoiding storing the whole data. We propose here several sampling strategies and estimators to estimate the median trajectory. A comparison between them is illustrated by means of a test population. We develop a stratification based on the linearized variable which substantially improves the accuracy of the estimator compared to simple random sampling without replacement. We suggest also an improved estimator that takes into account auxiliary information. Some potential areas for future research are also highlighted.     

Estimation from two-stage unequal probability sampling with missing units

A formula is presented for an unbiased estimator for the variance of an unbiased estimator of a survey population total as well as for an unbiased estimator of its variance based on sampling in two-stages following Rao et al. J Roy Stat Soc B 24: 482–491 (1962) scheme in both stages when the originally selected units in both stages cannot be fully covered in the survey but are to be randomly sub-sampled. The development is helpful to tackle non-responses if assumed to have occurred at random in either or both the stages     

A note on estimation based on record data

In data-processing standpoint, an efficient algorithm for identifying the minimum value among a set of measurements are record statistics. From a sequence of n independent identically distributed continuous random variables only about log(n) records are expected, so we expect to have little data, hence any prior information is welcome (Houchens, Record value theory and inference, Ph.D. thesis, University of California, Riverside, 1984). In this paper, non-Bayesian and Bayesian estimates are derived for the two parameters of the Exponential distribution based on record statistics with respect to the squared error and Linear-Exponential loss functions and then compared with together. The admissibility of some estimators is discussed.     

Randomized response: Estimating mean square errors of linear estimators and finding optimal unbiased strategies

Summary General procedures are described to generate quantitative randomized response (RR) required to estimate the finite population total of a sensitive variable. Permitting sample selection with arbitrary probabilities a formula for the mean square error (MSE) of a linear estimator of total based on RR is noted indicating the simple modification over one that might be based on direct response (DR) if the latter were available. A general formula for an unbiased estimator of the MSE is presented. A simple approximation is proposed in case the RR ratio estimator is employed based on a simple random sample (SRS) taken without replacement (WOR). Among sampling strategies employing unbiased but not necessarily linear estimators based on RR, certain optimal ones are identified under two alternative models analogously to well-known counterparts based on DR, if available. Unlike Warner’s (1965) treatment of categorical RR we consider quantitative RR here.     

Optimum allocation in multivariate stratified random sampling: stochastic matrix mathematical programming

The allocation problem for multivariate stratified random sampling as a problem of stochastic matrix integer mathematical programming is considered, minimizing the estimated covariance matrix of estimated means subject to fixed cost or fixed total sample size. With these aims the asymptotic normality of sample covariance matrices for each strata is established. Some alternative approaches are suggested for its solution. An example is solved by applying the proposed techniques.     

Developments in Survey Research over the Past 60 Years: A Personal Perspective

Many developments have occurred in the practice of survey sampling and survey methodology in the past 60 years or so. These developments have been partly driven by the emergence of computers and the continuous growth in computer power over the years and partly by the increasingly sophisticated demands from the users of survey data. The paper reviews these developments with a main emphasis on survey sampling issues for the design and analysis of social surveys. Design‐based inference based on probability samples was the predominant approach in the early years, but over time, that predominance has been eroded by the need to employ model‐dependent methods to deal with missing data and to satisfy analysts' demands for survey estimates that cannot be met with design‐based methods. With the continuous decline in response rates that has occurred in recent years, much current research has focused on the use of non‐probability samples and data collected from administrative records and web surveys.