Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: can I generate my own bootstrap weights for survey data?


From   "Stas Kolenikov" <skolenik@gmail.com>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: can I generate my own bootstrap weights for survey data?
Date   Mon, 17 Apr 2006 11:22:25 -0500

One thing you can try is to look into Nick Winter's -svr- code to see
how he created jackknife weights.

If I were to do this, I would take yat another look at Rao & Wu paper
on resampling with survey data (Journal of the American Statistical
Association, 83:231-- 241, 1988), figure out the correct number of
clusters to get resampled, check for other caveats, and then, in
Stata-like pseudocode,

g double bsweight = 0 if ~mi( actual_pweight )
foreach `h' in the list of strata {
  forvalues k=1/desired number of resampled clusters {
    local random_cluster = uniform()*(# clusters in the current stratum) + 1
    replace bsweight = bsweight + actual_pweight if stratum == `h'  &
cluster == `random_cluster'
  }
}

This should give you one bootstrap resample. Now, the question
remains, at least to me, what one is to do with the 0 weights (should
they indeed be zero or missing? I would keep them at zero), and what
happens if in a certain stratum, only a single cluster has been
randomly selected (I would probably discard such subsample). See if
that is at all addressed in the paper above, and try to find some
follow-ups -- and please let us know if you are successful.

Hope Nick Winter or Jeff Pitblado would comment on the above. It would
be nice if -svy- were to take in the bootstrap weights, too -- I don't
think it does currently. You can probably get around by feeding them
in as -jackknife- weights to the aforementioned -svr-.

On 4/14/06, Dongyi Du <ddu001@gmail.com> wrote:
> hi, all:
>
> I have information of weight, psu, and strata for a survey data, can
> anyone tell me if I can generate my own bootstrap weights? if I can,
> how?
>
> thanks.
>
> Dong
>
> *
> *   For searches and help try:
> *   http://www.stata.com/support/faqs/res/findit.html
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>


--
Stas Kolenikov
http://stas.kolenikov.name

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index