Pareto Sampling versus Sampford and Conditional Poisson Sampling

Abstract.  Pareto sampling was introduced by Rosén in the late 1990s. It is a simple method to get a fixed size π ps sample though with inclusion probabilities only approximately as desired. Sampford sampling, introduced by Sampford in 1967, gives the desired inclusion probabilities but it may take time to generate a sample. Using probability functions and Laplace approximations, we show that from a probabilistic point of view these two designs are very close to each other and asymptotically identical. A Sampford sample can rapidly be generated in all situations by letting a Pareto sample pass an acceptance–rejection filter. A new very efficient method to generate conditional Poisson ( CP ) samples appears as a byproduct. Further, it is shown how the inclusion probabilities of all orders for the Pareto design can be calculated from those of the CP design. A new explicit very accurate approximation of the second-order inclusion probabilities, valid for several designs, is presented and applied to get single sum type variance estimates of the Horvitz–Thompson estimator.     

Conditional and Restricted Pareto Sampling: Two New Methods for Unequal Probability Sampling

Abstract. Two new unequal probability sampling methods are introduced: conditional and restricted Pareto sampling. The advantage of conditional Pareto sampling compared with standard Pareto sampling, introduced by Rosén (J. Statist. Plann. Inference, 62, 1997, 135, 159), is that the factual inclusion probabilities better agree with the desired ones. Restricted Pareto sampling, preferably conditioned or adjusted, is able to handle cases where there are several restrictions on the sample and is an alternative to the recent cube method for balanced sampling introduced by Deville and Tillé (Biometrika, 91, 2004, 893). The new sampling designs have high entropy and the involved random numbers can be seen as permanent random numbers.     

Order Sampling Design with Prescribed Inclusion Probabilities

Order sampling with fixed distribution shape is a class of sampling schemes with inclusion probabilities approximately proportional to given size measures. In a recent article, methods were provided to compute the exact first and second order inclusion probabilities numerically when the distribution shape is of the Pareto type. In the same article, procedures were also provided for this case to adjust the parameters to get predetermined inclusion probabilities. In this paper we prove the existence and uniqueness of a solution for the latter problem, in general for any order sampling of fixed distribution shape.     

A New Sampling Design for Systematic Sampling

A sampling design called “Modified Systematic Sampling (MSS)” is proposed. In this design each unit has an equal probability of selection. Moreover, it works for both situations: N = nk or N ≠ nk. Consequently, the Linear Systematic Sampling (LSS) and Circular Systematic Sampling (CSS) become special cases of the proposed MSS design.     

Remainder Markov systematic sampling

Systematic sampling is the simplest and easiest of the most common sampling methods. However, when the population size N cannot be evenly divided by the sampling size n, systematic sampling cannot be performed. Not only is it difficult to determine the sampling interval k equivalent to the sampling probability of the sampling unit, but also the sample size will be inconstant and the sample mean will be a biased estimator of the population mean. To solve this problem, this paper introduces an improved method for systematic sampling: the remainder Markov systematic sampling method. This new method involves separately finding the first-order and second-order inclusion probabilities. This approach uses the Horvitz-Thompson estimator as an unbiased estimator of the population mean to find the variance of the estimator. This study examines the effectiveness of the proposed method for different super-populations.     

Robust Inference in Two‐phase Sampling Designs with Application to Unit Nonresponse

Influential units occur frequently in surveys, especially in business surveys that collect economic variables whose distributions are highly skewed. A unit is said to be influential when its inclusion or exclusion from the sample has an important impact on the sampling error of estimates. We extend the concept of conditional bias attached to a unit and propose a robust version of the double expansion estimator, which depends on a tuning constant. We determine the tuning constant that minimizes the maximum estimated conditional bias. Our results can be naturally extended to the case of unit nonresponse, the set of respondents often being viewed as a second‐phase sample. A robust version of calibration estimators, based on auxiliary information available at both phases, is also constructed.     

A Class of Sampling Two Units with Probability Proportional to Size

A class of sampling two units without replacement with inclusion probability proportional to size is proposed in this article. Many different well known probability proportional to size sampling designs are special cases from this class. The first and second inclusion probabilities of this class satisfy important properties and provide a non-negative variance estimator of the Horvitz and Thompson estimator for the population total. Suitable choice for the first and second inclusion probabilities from this class can be used to reduce the variance estimator of the Horvitz and Thompson estimator. Comparisons between different proportional to size sampling designs through real data and artificial examples are given. Examples show that the minimum variance of the Horvitz and Thompson estimator obtained from the proposed design is not attainable for the most cases at any of the well known designs.     

Optimal PPS Sampling with Vanishing Auxiliary Variables – with Applications in Microscopy

Recently, non‐uniform sampling has been suggested in microscopy to increase efficiency. More precisely, proportional to size (PPS) sampling has been introduced, where the probability of sampling a unit in the population is proportional to the value of an auxiliary variable. In the microscopy application, the sampling units are fields of view, and the auxiliary variables are easily observed approximations to the variables of interest. Unfortunately, often some auxiliary variables vanish, that is, are zero‐valued. Consequently, part of the population is inaccessible in PPS sampling. We propose a modification of the design based on a stratification idea, for which an optimal solution can be found, using a model‐assisted approach. The new optimal design also applies to the case where ‘vanish’ refers to missing auxiliary variables and has independent interest in sampling theory. We verify robustness of the new approach by numerical results, and we use real data to illustrate the applicability.     

An efficient non-rejective implementation of the πps sampling designs

Poisson sampling is a method for unequal probabilities sampling with random sample size. There exist several implementations of the Poisson sampling design, with fixed sample size, which almost all are rejective methods, that is, the sample is not always accepted. Thus, the existing methods can be time-consuming or even infeasible in some situations. In this paper, a fast and non-rejective method, which is efficient even for large populations, is proposed and studied. The method is a new design for selecting a sample of fixed size with unequal inclusion probabilities. For the population of large size, the proposed design is very close to the strict πps sampling which is similar to the conditional Poisson (CP) sampling design, but the implementation of the design is much more efficient than the CP sampling. And the inclusion probabilities can be calculated recursively.     

An improved cluster two-stage randomized response model

In this paper, Abdelfatah and Mazloum's (2015) two-stage randomized response model is extended to unequal probability sampling and stratified unequal probability sampling, both with and without replacement. The extended models result in more efficient estimators than Lee et al.'s (2014) estimators of the proportion of the population having a sensitive attribute.     

Randomized nomination sampling for finite populations

We propose a randomized minima–maxima nomination (RMMN) sampling design for use in finite populations. We derive the first- and second-order inclusion probabilities for both with and without replacement variations of the design. The inclusion probabilities for the without replacement variation are derived using a non-homogeneous Markov process. The design is simple to implement and results in simple and easy to calculate estimators and variances. It generalizes maxima nomination sampling for use in finite populations and includes some other sampling designs as special cases. We provide some optimality results and show that, in the context of finite population sampling, maxima nomination sampling is not generally the optimum design to follow. We also show, through numerical examples and a case study, that the proposed design can result in significant improvements in efficiency compared to simple random sampling without replacement designs for a wide choice of population types. Finally, we describe a bootstrap method for choosing values of the design parameters.     

On a generalization of Poisson sampling

In real-time sampling, the units of a population pass a sampler one by one. Alternatively the sampler may successively visit the units of the population. Each unit passes only once and at that time it is decided whether or not it should be included in the sample. The goal is to take a sample and efficiently estimate a population parameter. The list sequential sampling method presented here is called correlated Poisson sampling. The method is an alternative to Poisson sampling, where the units are sampled independently with given inclusion probabilities. Correlated Poisson sampling uses weights to create correlations between the inclusion indicators. In that way it is possible to reduce the variation of the sample size and to make the samples more evenly spread over the population. Simulation shows that correlated Poisson sampling improves the efficiency in many cases.     

Spatially correlated Poisson sampling

A new method for sampling from a finite population that is spread in one, two or more dimensions is presented. Weights are used to create strong negative correlations between the inclusion indicators of nearby units. The method can be used to produce unequal probability samples that are well spread over the population in every dimension, without any spatial stratification. Since the method is very general there are numerous possible applications, especially in sampling of natural resources where spatially balanced sampling has proven to be efficient. Two examples show that the method gives better estimates than other commonly used designs.     

A List Sequential Sampling Method Suitable for Real-Time Sampling

Abstract.  A flexible list sequential π ps sampling method is introduced and studied. It can reproduce any given sampling design without replacement, of fixed or random sample size. The method is a splitting method and uses successive updating of inclusion probabilities. The main advantage of the method is in real-time sampling situations where it can be used as a powerful alternative to Bernoulli and Poisson sampling and can give any desired second-order inclusion probabilities and thus considerably reduce the variability of the sample size.     

On use of symmetrical balanced incomplete block design in construction of sampling design realising pre-assigned sets of inclusion probabilities of first two orders.

Gupta, Nigam and Kumar (1982) proposed a sampling scheme using certain combinatorial properties of balanced incomplete block design (BIBD) which realises first order inclusion probabilities proportional to measure of size(IPPS). Here their results have been extended by presenting samplng schemes realising pre-assigned sets of inclusion probabilities of first two orders.     

Optimal sample coordination using controlled selection

Sample coordination maximizes or minimizes the overlap of two or more samples selected from overlapping populations. It can be applied to designs with simultaneous or sequential selection of samples. We propose a method for sample coordination in the former case. We consider the case where units are to be selected with maximum overlap using two designs with given unit inclusion probabilities. The degree of coordination is measured by the expected sample overlap, which is bounded above by a theoretical bound, called the absolute upper bound, and which depends on the unit inclusion probabilities. If the expected overlap equals the absolute upper bound, the sample coordination is maximal. Most of the methods given in the literature consider fixed marginal sampling designs, but in many cases, the absolute upper bound is not achieved. We propose to construct optimal sampling designs for given unit inclusion probabilities in order to realize maximal coordination. Our method is based on some theoretical conditions on joint selection probability of two samples and on the controlled selection method with linear programming implementation. The method can also be applied to minimize the sample overlap.     

Simple random sample equivalent survey designs reducing undesirable units from a finite population

In this paper, we consider fixed size sampling plans for which the first order inclusion probabilities are identical for all units and the second order inclusion probabilities are constant for every pair-wise unit. The statistical conditions are identified under which these plans are equivalent to the usual simple random sampling plan. These sampling plans are constructed to reduce undesirable units.     

Examination of Influential Observations to Improve Discriminant Analysis

Three simple transformations are proposed in the context of ratio and product methods of estimation, based on any probability sampling design, and the usual unbiased estimation under varying probability sampling. These transformations may be effected

after the data are collected in a survey. The objective is to obtain improved estimators of the population total     

HORVITZ-THOMPSON AND DES RAJ ESTIMATORS REVISITED

In this paper it is shown that the generalized πPS sampling strategy consisting of the design with π_i, the probability of inclusion of the ith unit in the sample, proportional to the modified size together with the corresponding Horvitz-Thompson estimator (Rao, 1971), is superior to the symmetrized Des Raj strategy under a general super-population set-up for all values of the super-population parameter g, when the samples are of size two.     

On inclusion probabilities and relative estimator bias for Pareto πps sampling

A means for utilizing auxiliary information in surveys is to sample with inclusion probabilities proportional to given size values, to use a πps design, preferably with fixed sample size. A novel candidate in that context is Pareto πps. This sampling scheme was derived by limit considerations and it works with a degree of approximation for finite samples. Desired and factual inclusion probabilities do not agree exactly, which in turn leads to some estimator bias. The central topic in this paper is to derive conditions for the bias to be negligible.Practically useful information on small sample behavior of Pareto πps can, to the best of our understanding, be gained only by numerical studies. Earlier investigations to that end have been too limited to allow general conclusions, while this paper reports on findings from an extensive numerical study. The chief conclusion is that the estimator bias is negligible in almost all situations met in survey practice.