Re: st: subpop() and over() with survey mean estimation using SDR weights

[Sam Schulhofer-Wohl <sschulh1.work@gmail.com>]

Jeff,

Thank you for checking this out and for the workaround.

Is -svy bootstrap- with -bsn(4)- equivalent to SDR weights in all cases?
For example, if I am using the IPUMS version of the Current Population
Survey with

svyset [iw=wtsupp], sdrweight(repwtp1-repwtp160) vce(sdr)

can I instead do

svyset [iw=wtsupp], bsrweight(repwtp1-repwtp160) bsn(4) vce(bootstrap)

(Comparing the formulas on pages 186 and 190 of [SVY] suggests that the two
settings are indeed always equivalent, but official confirmation would be
nice.)

Thanks,

Sam



On Mon, Jan 30, 2012 at 16:06 PM, Jeff Pitblado, StataCorp LP <
jpitblado@stata.com> wrote:

Sam Schulhofer-Wohl <sschulh1.work@gmail.com> is getting a syntax error from
> the -svy:- prefix when using SDR weights with a -subpop()- and -mean,
> over()-:
>
> > I get a syntax error message, r(198), when I try to estimate a survey
> > mean using both the subpop() and over() options and SDR weights.
> >
> > I can't find any indication in the manual that this combination of
> > options and weights is not allowed. Also, for other kinds of weights,
> > there is no error. Some examples are below.
> >
> > Has anyone run into this problem before? If so, any suggestions on how
> > to solve it? (Obviously, I can make a new variable that combines my
> > over() and subpop() categories, but in the particular analysis I'm
> > using, it would be more convenient to be able to use over() and
> > subpop() simultaneously.)
>
> We are able to reproduce the syntax error and believe that this is a bug in
> - -svy sdr-.
>
> We hope to have this fixed in a future update in Stata 12.
>
> In the mean time, Sam can trick the -svy bootstrap- prefix to produce the
> SDR
> standard errors by using option -bsn(4)- with the SDR weights -svyset- as
> - -bsrweight()- instead.
>
> Assuming that only the SDR weight variables are prefixed 'pwgtp' Sam's
> - -svyset- command would be
>
>        . svyset [pw=pwgtp], bsrweight(pwgtp?*) bsn(4) vce(bootstrap)
>
> - --Jeff
> jpitblado@stata.com
>


--
Sam Schulhofer-Wohl
Senior Economist
Research Department
Federal Reserve Bank of Minneapolis
90 Hennepin Ave.
Minneapolis MN 55480-0291
(612) 204-5484
wohls@minneapolisfed.org
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/