Re: st: svy bootstrap

[Stas Kolenikov <skolenik@gmail.com>]

Your -twopart- program expected weights, but you did not supply any
when you called it for the first time without any arguments. A more
protective code would have something like

if "`weight'" == "" {
  display as error "twopart was expecting [pweight=variable], but none
was specified"
  exit 198
}

to make sure all the required inputs were provided by the user (or by
the -bootstrap-).

On Wed, Mar 9, 2011 at 9:13 PM, Michael Palmer
<Michael.Palmer@anu.edu.au> wrote:
> Thanks so much for your reply. I went ahead with the suggested code in the single program.
> However, I was met with the error message: '=unknown weight type'.
>
> The revised do file is:
>
> capture program drop twopart
> program twopart, rclass
> syntax [pweight /]
> logit inpat_public $xlist1 if age > 5 [`weight'=`exp']
> predict `pr_y', pr
> reg ipexp_public_log $xlist1 if age > 5 & inpat_public==1 [`weight'=`exp']
> predict `yhat', xb
> predict `ehat', residuals

etc.

> return scalar bpwd = `m_1'*`m_3' - `m_2'*`m_4'
>
> end
>
> set more off
> twopart
> return list
> bsweights bsw, reps(100) n(-1)
> bs4rw (bpwd=r(pwd)), rw(bsw*): twopart [pw=psindwt]
>


-- 
Stas Kolenikov, also found at http://stas.kolenikov.name
Small print: I use this email account for mailing lists only.
*
*   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/