Resampling under Complex Sampling Designs: Comparison
Please note this is a comparison between Version 2 by Yvaine Wei and Version 1 by Pier Luigi Conti.

ResamIn pling methods have a long and honorable history. Srinciple, survey data are an ideal context to useapply resampling methods to approximate the (unknown) sampling distribution of statistics, due to both a generusually large sample size and data of typicontrolled quality. However, survey data cannot be generally assumed independent and identically good qualitydistributed (i.i.d.) so that any resampling methodologies to be used in sampling from finite populations must be adapted to account for the sample design effect. A principled appraisal is given and discussed here.

  • resampling
  • bootstrap
  • pseudo-population
  • design effect

1. Introduction

Generalities

Resampling methods have a long and honorable history, going back at least to the Efron's seminal paper by [1]. 

Iin exlatreme synthesis,e 70s [1]. vVirtually all resampling methodologies used in sampling from finite populations are based on the idea of accounting for the effect of the sampling design. T In fact, the main effect of the sampling design is that data cannot be generally assumed independent and identically distributed (i.i.d.).

The Amain approaches are large essentially two: the ad hoc approach and the plug in approach. The basic idea of the ad hoc apprortion of the ach consists in maintaining Efron’s bootstrap as a resampling procedure but in properly rescaling data in order to account for the dependence among units. This approach is used, among others, in [2][3], where the re-sampled data produced by the “usual” i.i.d. bootstrap are properly rescaled, as well as in [4][5]; cfr. also ther review in [6]. In [7] a “rescaled booture on strap process” based on asymptotic arguments is proposed. Among the ad hoc approaches, we also classify [8] (based on a rescaling of weights) and the “direct bootstrap” by [9]. Almost all ad hoc resampling ftechniques arom finite populations focuses one based on the same justification: in the case of linear statistics, the first two moments of the resampled statistic should match (at least approximately) the corresponding estimatiors; cfr., among the vothers, [9]. Cfr. also [8], where an analysiance of estimas in terms of the first three moments is performed for Poisson sampling.

Here the secorsnd approach based on pseudo-populations is considered. The main approachereasons beyond this choice are i) resampling based on pseudo-populations actually parallels Efron’s bootstrap for i.i.d. observations; ii) the basic ideas are essentially the ad hoc relatively simple to understand and to apply once the problem is approached and plug inin terms of an appropriate estimator of the finite population distribution function (f.p.d.f.); and iii) the main theoretical justification for resampling bapproachsed on pseudo-population is of asymptotic nature, similar in many respects, to the well known Bickel-Freedman results [10] for Efron’s bootstrap. 

The basic idea of the ad hoc approach consists in maintaining Efron’s bootstrap as a resampling procedure but in properly rescaling data in order to account for the dependence among units. This approach is used, among others, in [3,4], where the re-sampled data produced by the “usual” i.i.d. bootstrap are properly rescaled, as well as in [5,6]; cfr. also the review in [7]. In [8] a “rescaled bootstrap process” based on asymptotic arguments is proposed. Among the ad hoc approaches, we also classify [9] (based on a rescaling of weights) and the “direct bootstrap” by [10].
Almost all ad hoc resampling techniques are based on the same justification: in the case of linear statistics, the first two moments of the resampled statistic should match (at least approximately) the corresponding estimators; cfr., among the others, [10]. Cfr. also [9], where an analysis in terms of the first three moments is performed for Poisson sampling.

2. Accounting for the Sampling Design in Resampling: The Pseudo-Population Approach

Among several techniques that aim at accounting for the sampling design in resampling from finite populations, we consider here the approach based on pseudo-populations. The idea of pseudo-population goes back, at least, to [11] in the case of median estimation essentially under srs when the population size is a multiple of the sample size. Rather similar ideas are in [12] for srs, again under the condition that the ratio between population size and sample size is a ninteger, and in [13], for stratified random sampling. A major step forward is the paper by [14], where the construction of a pseudo-population is studied under a general πps sampling design, with general first order inclusion probabilities. In [19][15], a different approach to the construction of a pseudo-population, very interesting in many respects, is considered. The pseudo-population approach to resampling can be considered as a two-phase procedure. In the first phase, a pseudo-population (roughly speaking, a prediction of the population) is constructed. In the second phase, a (bootstrap) sample is drawn from the pseudo-population. Broadly speaking, this approach parallels the plug-in principle by Efron. The pseudo-population is plugged in the sampling process and is used as a “surrogate” of the actual finite population. In the second phase, a sample is drawn from the pseudo-population, according to a sampling design that mimics the original one. In this view, the pseudo-population mimics the real population, and the (re)sampling process from the pseudo-population mimics the (original) sampling process from the real population.

2.1. Resampling from Pseudo-Populations

ResampliIng based[16], onit pseudo-populations actually parallels Efron’s bootstrap for i.i.d. is thoroughly illustrated a class of resampling techniques for finite pobservpulations. The basic ideas are relatively simple, once the problem is approached in terms of an appropriate estimator under a general (complex) sampling design which is asymptotically ofcorrect thunde f.p.d.f.

2.2. Resampling Based on Pseudo-Populations: Basics Results for Descriptive Inference

Ther mild assumain theoreptions. Practical justification for resampling based on recommendations for finite sample sizes and on how to construct the pseudo-population is of asymptotic nature, similar, in many respects, to results in [17] for Efron’s bootstrapare also given.

3. Computational Issues

UseApplication of the pseudo-population approach, despite its many theoretical merits, is held backcan be limited in practice by its computational complexityburden. Real populations could contain millions of peopleunits, and thus the cactual construction of a pseudo-population could be computationally cumbersome. For this reason, it is of primary interest to develop shortcuts that, while possessing the fundamental theoretical properties described in the above sections, are computationally simple to implement because they avoid the physical construction of the pseudo-population. The above poiInts are[17] thoroughly discussed in [26], where the problem of resampling for finite populations is addressed as a problem of sampling with replacement directly from the sample data, the original sample, henceforth, with different drawing probabilities. An attempt to avoid complications related to integer-valued Ni bootstrap weights, i.e. the number of replications of each sampled is iunit into the pseudo-population, is offered in [27][18], where non-integer Ni bootstrap weights are allowed via  via the Horvitz–Thompson-based bootstrap (HTB) method. However, unless the sampling fraction n/N  tends to 0 as both population s N ize and n  sample size increase, HTB does not generally possess the good asymptotic properties outlined in the[16]. previous sections. An interesting computational shortcut is in [28][19], where the pseudo-population (again with possibly non-integer Nis) is only implicitly used, and a computational scheme based on drawings with replacements from the original sample is proposed. Unfortunately, although the main idea behind that paper is interesting, the proposed bootstrap method fails to possess good asymptotic properties. Computational shortcuts, based on ideas similar to those in [28][19], but based on correct approximations of first order inclusion probabilities, were developed in [29][20] for descriptive, design-based inference. In particular, in that paper, methodologies based on drawings with replacements from the original sample were proposed, and their merits, from both a theoretical and a computational point of view, were studied. As remarked by a referee, another drawback of the pseudo-population approach is the apparent necessity to generate and save a large number of bootstrap sample files. However, it is not necessary to save all the bootstrap sample files. Only the original sample file must be saved along with two additional variables for each bootstrap replicate: one variable that contains the number of times each sample unit is used to create the pseudo-population and another one containing the number of times each sample unit has been selected in the bootstrap sample. In other words, it can be implemented similar to methods that rescale the sampling weights.

Another practical drawback related to the pseudo-population approach is the seeming necessity to generate and store a large number of bootstrap sample files. However, it is not necessary to save all the bootstrap sample files. Only the original sample file should be saved along with two additional variables for each bootstrap replicate: one variable that contains the number of times each sample unit is used to create the pseudo-population and another one containing the number of times each sample unit has been selected in the bootstrap sample. In other words, it can be implemented similarly to methods that rescale the sampling weights.

4. Open Problems and Final Considerations

The pseudo-population approach, despite its merits, requires further development from both the theoretical and computational perspectives.

The pseudo-population app Froach, despite its merits, requires further development from both the theoretical and computational p a theoretical perspectives. From a theoretical point of view, the results obtained thus far only refer to non-informative single-stage designs. The consideration of multi-stage designs appears as a necessary dedevelopment as well as the consideration of non-respondent units. Again, from a t of theoretical perspective, a major issue is the development of theoretically sound resampling methodologies for informative sampling designs is a major issue calling for more research. The majorin drawback is that, apart from the exception of adaptive designs (cfr. [30][21]) and the references therein) first order inclusion probabilities can rarely be computed, as these might depend on unobserved quantities. This is what happens, for instance, with most of the network sampling designs that are actually used for hidden populations, where the inclusion probabilities are unknown and depend on unobserved/unknown network links (cfr. [30,31][21][22] and the references therein). F Again from athe computationtheoretical point of view, as indicated earlier, the computational shortcuts developed thus far only work in the case of descriptive inference. The the consideration of multistage designs appears as a further necessary development of theoretically well-founded computational schemes valid for analytic inference is an important issue that deserves further attentionas well as the consideration of non-respondent units.

From a computational point of view, the computational shortcuts developed thus far only apply to the case of descriptive inference. The development of theoretically well-founded  computational schemes valid for analytic inference is an important issue that deserves further attention.

Video Production Service