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

From |
"Paul Seed" <P.T.Seed@qmul.ac.uk> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
st: RE: statalist-digest V4 #970 |

Date |
Thu, 8 Aug 2002 13:42:51 +0100 |

Many thanks to Gary for this. Changing the saved result to a scalar did indeed take care of the problem. Although the term "control" is used, this is not a matched case-control study, rather a comparison of two groups of pregnant women, (with & without APS). An unpaied t-test (or non-parametric equivalent) is entirely appropriate. Date: Wed, 07 Aug 2002 10:47:41 -0700 From: Gary Longton <glongton@fhcrc.org> Subject: Re: st: Problems with -bs- command Paul Seed wrote: > I am interested in a bootstrapped confidence interval for the > difference between two means, and to save time > wish to work out the means using -summ, meanonly-. > > I have written a Stata program to do this, which runs perfectly well > on it's own. However, when trying to get -bs- to work out the confidence > intervals, I get an unexplained error message. > > Any thoughts on wher ethe problem lies, greatly appreciated. the program: > cap prog drop myprog > prog define myprog, rclass > summ centile if aps, mean > di "r(mean) = " r(mean) > local m1 = r(mean) > return local m1 = r(mean) > summ centile if control, mean > return local m2 = r(mean) > return local m2_m1 = `m1' - r(mean) > end myprog called with -bs- > set seed 35435141 > bs "myprog" "r(m1)", reps(50) or > bs "myprog" "r(m1) r(m2) r(m2_m1)", reps(50) and output: > . * Now use it for bootstrap" > . set seed 35435141 > > . bs "myprog" "r(m1) r(m2) r(m2_m1)", reps(50) > > command: myprog > statistics: r(m1) r(m2) r(m2_m1) > error in statistic: r(m1) > r(198); I believe the -bs- command is complaining about the use of macro results placed in the expression list, i.e. "r(m1)" It should work if these are instead scalar results. Eg. in the program - -myprog-, replace return local m1 = r(mean) with return scalar m1 = r(mean) ^^^^^^ and likewise for m2, and m2_m1 and it should work. However, am wondering whether this type of sampling is best for your design? The command > summ centile if control, mean hints that this might be a case-control design?, in which case you might want to be sampling separately from cases and controls? which can be done, but not with - -bs-. - - Gary * * 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/

- Prev by Date:
**Re:st:problems with outliers** - Next by Date:
**st: trivariate probit** - Previous by thread:
**Re:st:problems with outliers** - Next by thread:
**st: trivariate probit** - Index(es):

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