Statalist


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

st: RE: issue with bootstrap


From   "Martin Weiss" <[email protected]>
To   <[email protected]>
Subject   st: RE: issue with bootstrap
Date   Tue, 25 Nov 2008 18:38:18 +0100

The issue here seems to be that -bootstrap- leaves behind "e(cmd)=bootstrap"
if the user specified additional statistics (like  e(rmse)) while it leaves
behind "e(cmd)=regress" if all you asked for was "_b". Along comes -estat
bootstrap- and asks " if "if "`e(cmd)'" != "regress" {" {" and exits with
error 301...

HTH
Martin


-----Original Message-----
From: [email protected]
[mailto:[email protected]] On Behalf Of Lachenbruch,
Peter
Sent: Tuesday, November 25, 2008 6:24 PM
To: [email protected]
Subject: st: issue with bootstrap

I was doing a bootstrap and wanted to get estat estimates and got "last
estimates not found"  Here are the commands and output

. bootstrap _b e(r2) e(rmse),bca reps(500) :reg tvas age anyapheresis
anyaza anyca anycyclosp anyhcq anyiv_cytoxan  anyivig anymtx
anypo_cytoxan anypred anyulc calcinosis ck252 cmasxt2 contracture
duractdz dysphagia dysphonia dzactvas extramuscleact finchaq funcclass
gender onsetseverity
(running regress on estimation sample)

Jackknife replications (84)
<snip]
Bootstrap replications (500)
<snip>
Linear regression                          Number of obs      =
84
                                           Replications       =
500
                                           R-squared          =
0.8502
                                           Adj R-squared      =
0.7856
                                           Root MSE           =
3.8019

      command:  regress tvas age anyapheresis anyaza anyca anycyclosp
anyhcq anyiv_cytoxan anyivig anymtx anypo_cytoxan anypred anyulc
calcinosis ck252 cmasxt2 contracture duractdz dysphagia dysphonia
dzactvas extramuscleact finchaq funcclass gender onsetseverity
  [_eq2]_bs_1:  e(r2)
  [_eq2]_bs_2:  e(rmse)

------------------------------------------------------------------------
----
             |   Observed   Bootstrap                      Normal-based
        tvas |      Coef.   Std. Err.      z    P>|z|   [95% Conf.
Interval]
-------------+----------------------------------------------------------
----
_eq1         |
         age |  -.0417141   .1664602    -0.25   0.802    -.36797
.2845418
<snip >
onsetsever~y |  -.4387078   .7729467    -0.57   0.570  -1.953655
1.07624
       _cons |   15.74078   7.812418     2.01   0.044   .4287244
31.05284
-------------+----------------------------------------------------------
---
_eq2         |
       _bs_1 |   .8501516    .053571    15.87   0.000   .7451544
.9551488
       _bs_2 |   3.801854   .5107574     7.44   0.000   2.800788
4.80292
------------------------------------------------------------------------
----

. estat boot , bca percentile
last estimates not found

This seems to follow the help file, but I'm mystified.  

Tony

Peter A. Lachenbruch
Department of Public Health
Oregon State University
Corvallis, OR 97330
Phone: 541-737-3832
FAX: 541-737-4001


 

*
*   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/



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