Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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

From   Sam Schulhofer-Wohl <>
Subject   Re: st: subpop() and over() with survey mean estimation using SDR weights
Date   Tue, 31 Jan 2012 08:13:36 -0600


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



On Mon, Jan 30, 2012 at 16:06 PM, Jeff Pitblado, StataCorp LP <> wrote:

Sam Schulhofer-Wohl <> 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
> 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

Sam Schulhofer-Wohl
Senior Economist
Research Department
Federal Reserve Bank of Minneapolis
90 Hennepin Ave.
Minneapolis MN 55480-0291
(612) 204-5484
*   For searches and help try:

© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index