Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

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

From |
Stas Kolenikov <skolenik@gmail.com> |

To |
statalist@hsphsun2.harvard.edu |

Subject |
Re: st: Program syntax for svy bootstrap or bs4rw |

Date |
Tue, 26 Jun 2012 09:27:15 -0500 |

May be Jeff P can reply and clarify what's missing from this syntax. On Mon, Jun 25, 2012 at 5:45 PM, Mark McGovern <markemcgovern@gmail.com> wrote: > Thank you Steve and Stas for your suggestions. For now I can work with > bs4rw, although at some point I think I might need to use some of the > functionality of svy, so I may have to return to this in the future. > > As you pointed out Steve, from the Stata 12 programming manual (p.326): > > 1. Command must be eclass and allow iweights and accept the standard > estimation syntax. > 2. Command must save the model coefficients and ancillary parameters > in e(b) and the estimation sample size in e(N), and it must identify > the estimation subsample in e(sample); > 3. svy’s vce(bootstrap), vce(brr), and vce(sdr) require that command > have svyb as a property. > > I thought these were satisfied as the original code below does return > results in e() in the usual way from "regress". For example, > specifically overwriting e() does not help as in the below. > > Regards, > > Mark > > webuse nmihs_bs, clear > svyset _n [pweight=finwgt], bsrweight(bsrw1-bsrw5) vce(bootstrap) > singleunit(missing) > > capture program drop test1 > program define test1, eclass properties(svyb) > syntax [if] [in] [pw iw] > marksample touse > tempvar b V > qui xi:reg age race [`weight'`exp'] if `touse' > capture drop re1 > predict re1 > qui xi:reg birthwgt re1 [`weight'`exp'] if `touse' > local nobs = e(N) > matrix `b' = e(b) > matrix `V' = e(V) > ereturn post `b' `V', esample(`touse') > ereturn local cmd="test1" > ereturn scalar N = `nobs' > ereturn display > end > > svy : test1 > > On Sat, Jun 23, 2012 at 2:07 PM, Stas Kolenikov <skolenik@gmail.com> wrote: >> As a short cut, you can try to omit -eclass-, so that your program >> leaves the last estimates in the memory, treat the weights as [iw], >> and run -bs4rw _b- to produce the standard errors. >> >> On Sat, Jun 23, 2012 at 12:26 PM, Steve Samuels <sjsamuels@gmail.com> wrote: >>> Sorry for the double post. I had two drafts of the original message and forgotten that I had sent one. >>> >>> S. >>> >>> The Programming manual entry "Writing programs for use with svy" says that to use a command with a svy prefix it must return e(b), e(N), and e(sample). Even though you use none of these, perhaps their absence is causing the error. >>> >>> >>> Steve >>> sjsamuels@gmail.com >>> >>> On Jun 22, 2012, at 8:09 PM, Mark McGovern wrote: >>> >>> Thank you very much for the suggestion, I will look again at cmp. The >>> point of this is to use interviewer effects for selection. There are >>> some alternative approaches, however they all involve a similar two >>> step procedure which would require a similar bootstrap for the >>> standard errors. I have added weights to the program commands, however >>> I must still be omitting something else. Error: "last estimates not >>> found" Any further thoughts would be appreciated. >>> >>> webuse nmihs_bs >>> svyset _n [pweight=finwgt], bsrweight(bsrw1-bsrw5) vce(bootstrap) >>> singleunit(missing) >>> >>> program define test1, eclass properties(svyb) >>> syntax [if] [pw iw] [, NULLOPT] >>> xi:reg age race [`weight'`exp'] >>> capture drop re1 >>> predict re1 >>> xi:reg birthwgt re1 [`weight'`exp'] >>> ereturn local cmd="test1" >>> end >>> >>> svy : test1 >>> >>> Regards, >>> >>> Mark >>> >>> On Fri, Jun 22, 2012 at 3:35 PM, Stas Kolenikov <skolenik@gmail.com> wrote: >>>> Your program does not use the weights anywhere. It accepts them, as >>>> the -svy- requirements specify, but it just absorbs them without >>>> passing them on to the estimation commands. So it does run the same >>>> commands on the same data set, even though -svy bootstrap- does the >>>> best job it can in giving the right weights to you. >>>> >>>> In general, I would find it questionable that you make point >>>> predictions of the random effects and use them in any serious way >>>> later. The random effects must be integrated over, and you need to >>>> find a way to do that within your convoluted version of the Heckman >>>> model. May be you can cast this as a -cmp- model -- check David >>>> Roodman's work on that. >>>> >>>> On Fri, Jun 22, 2012 at 11:59 AM, Mark McGovern <markemcgovern@gmail.com> wrote: >>>>> Dear all, >>>>> >>>>> I am attempting to write a simple program for a bootstrap in which >>>>> predicted values (random effects) are estimated in a first stage and >>>>> then used in a heckman selection model in the second. >>>>> >>>>> This works fine when I use the standard bootstrap command, however I >>>>> am using survey data and would like to account for this with SVY as >>>>> discussed on this list previously. >>>>> >>>>> I have tried to meet the syntax requirements for SVY as outlined in >>>>> the programming manual, and have tried various additions but I am >>>>> having trouble with the program. I noted the discussion on >>>>> http://www.stata.com/statalist/archive/2011-03/msg01442.html and also >>>>> tried the bs4rw command. I would appreciate your advice. I am using >>>>> Stata 12 MP. >>>>> >>>>> Here is a simple illustration with the NHIMES data. I find that the >>>>> same sample is used in each bootstrap. >>>>> >>>>> I svyset the data and run the command (with just 5 replications for testing). >>>>> >>>>> webuse nmihs_bs >>>>> svyset _n [pweight=finwgt], bsrweight(bsrw1-bsrw5) vce(bootstrap) >>>>> singleunit(missing) >>>>> >>>>> program define test1, eclass properties(svyb) >>>>> syntax [if] [pw iw] [, NULLOPT] >>>>> xi:reg age race >>>>> capture drop re1 >>>>> predict re1 >>>>> xi:reg birthwgt re1 >>>>> ereturn local cmd="test1" >>>>> end >>>>> >>>>> svy : test1 >>>>> >>>>> bs4rw , rw(bsrw1-bsrw5) noisily: test1 [pweight=finwgt] >>>>> >>>>> Regards, >>>>> >>>>> Mark >>>>> * >>>>> * 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/ >>>> >>>> >>>> >>>> -- >>>> ---- Stas Kolenikov >>>> -- http://stas.kolenikov.name >>>> ---- Senior Survey Statistician, Abt SRBI >>>> -- Opinions stated in this email are mine only, and do not reflect the >>>> position of my employer >>>> >>>> * >>>> * 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/ >>> >>> * >>> * 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/ >>> >>> >>> * >>> * 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/ >> >> >> >> -- >> ---- Stas Kolenikov >> -- http://stas.kolenikov.name >> ---- Senior Survey Statistician, Abt SRBI >> -- Opinions stated in this email are mine only, and do not reflect the >> position of my employer >> >> * >> * 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/ > > * > * 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/ -- ---- Stas Kolenikov -- http://stas.kolenikov.name ---- Senior Survey Statistician, Abt SRBI -- Opinions stated in this email are mine only, and do not reflect the position of my employer * * 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/

**Follow-Ups**:**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Jeffrey Pitblado <jpitblado@stata.com>

**References**:**st: Program syntax for svy bootstrap or bs4rw***From:*Mark McGovern <markemcgovern@gmail.com>

**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Stas Kolenikov <skolenik@gmail.com>

**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Mark McGovern <markemcgovern@gmail.com>

**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Steve Samuels <sjsamuels@gmail.com>

**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Stas Kolenikov <skolenik@gmail.com>

**Re: st: Program syntax for svy bootstrap or bs4rw***From:*Mark McGovern <markemcgovern@GMAIL.COM>

- Prev by Date:
**Re: st: Re lincom not recognising macros** - Next by Date:
**Re: st: Regular expressions with locals** - Previous by thread:
**Re: st: Program syntax for svy bootstrap or bs4rw** - Next by thread:
**Re: st: Program syntax for svy bootstrap or bs4rw** - Index(es):