[R] Bootstrap and Jackknife Bias using Survey Package

Thomas Lumley tlumley at u.washington.edu
Tue Apr 11 17:14:24 CEST 2006

On Tue, 11 Apr 2006, Carlos Creva Singano (M2004078) wrote:

> Dear R users,
> I´m student of Master in Statistic and Data analysis, in New University 
> of Lisbon. And now i´m writting my dissertation in variance 
> estimation.So i´m using Survey Package to compute the principal 
> estimators and theirs variances.
> My data is from Incoming and Expendire Survey. This is stratified 
> Multi-stage Survey care out by National Statistic Institute of 
> Mozambique. My domain of analysis is Maputo City, the Capital of 
> Mozambique. I just compute the sampling errors using Survey Package, but 
> i have some droubles:

> 1. How to compute Bootstrap and Jackknife Bias of estimates, like mean?

I don't know how to do this for survey estimates, but if you do, you can 
compute it from the replicates.

> 2. How to see each replicate estimate of the parameter?

All the survey functions for replicate designs will return the replicates 
with the return.replicates=TRUE option

R> a <- svymean(~api00,rclus1,return.replicates=TRUE)
R> summary(a$replicates)
    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
   636.1   640.4   642.9   644.3   645.7   667.7
R> str(a)
List of 2
  $ mean      : atomic [1:1] 644
   ..- attr(*, "var")= num 693
   ..- attr(*, "statistic")= chr "mean"
  $ replicates: num [1:15] 643 648 646 643 645 ...
  - attr(*, "class")= chr "svrepstat"

> 3. Is it possible to use the Sitter algoritm in Bootstrap for stratified 
> and multi-stage sampling in Survey Package? Or is possible to use the 
> pseudo-population (obtained by replication the sample) in estimation 
> using Bootstrap?

as.svrepdesign( ,type="bootstrap")
does create a pseudo-population in designs with a finite population 
correction. This is as described in the reference by Davison and Canty 
given on the help page.

I don't know what you mean by the "Sitter algorithm".  Prof Sitter has 
written several papers on bootstrapping survey data -- can you give a more 
precise reference?


Thomas Lumley			Assoc. Professor, Biostatistics
tlumley at u.washington.edu	University of Washington, Seattle

More information about the R-help mailing list