Re: [R] Extracting components from a 'boot' class output in R
Sarath pointed out (privately) that the standard error is not actually stored as part of the bootstrap object. Instead, it is calculated on the fly when you print the bootstrap object; this is done by boot:::print.boot Similarly for bias. (The 'boot:::' is necessary because "'print.boot' is not an exported object from 'namespace:boot'"). Tim Hesterberg >Do > names(bootObj) >to find out what the components are, and use $ or [[ to extract >components. >Do > help(boot) >for a description of components of the object (look in the Value section). > >That is general advice in R, applying to all kinds of objects - >boot, and many other functions such as lm(), return lists with >a class added, and you can operate on the object as a list using >names(), $, etc. > >Tim Hesterberg > >>Dear R user, >> >>I used the following to do a bootstrap. >> >> >>>bootObj<-boot(data=DAT, statistic=Lp.est, >>R=1000,x0=3) >> >>I have the following output from the above bootstrap. How >>can I extract components of the output. >>For example, how can I extract the std.error? >> >> >>> bootObj >> >>ORDINARY NONPARAMETRIC BOOTSTRAP >> >>Call: >>boot(data = DAT, statistic = Lp.est, R = 1000, x0 = 3) >> >>Bootstrap Statistics : >> original bias std. error >>t1* 794.9745 -0.341 4.042099 >> >>Any help is greatly appreciated. >> >>Thank you >> >> >>Sarath Banneheka __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Extracting components from a 'boot' class output in R
Do names(bootObj) to find out what the components are, and use $ or [[ to extract components. Do help(boot) for a description of components of the object (look in the Value section). That is general advice in R, applying to all kinds of objects - boot, and many other functions such as lm(), return lists with a class added, and you can operate on the object as a list using names(), $, etc. Tim Hesterberg >Dear R user, > >I used the following to do a bootstrap. > > >>bootObj<-boot(data=DAT, statistic=Lp.est, >R=1000,x0=3) > >I have the following output from the above bootstrap. How >can I extract components of the output. >For example, how can I extract the std.error? > > >> bootObj > >ORDINARY NONPARAMETRIC BOOTSTRAP > >Call: >boot(data = DAT, statistic = Lp.est, R = 1000, x0 = 3) > >Bootstrap Statistics : > original bias std. error >t1* 794.9745 -0.341 4.042099 > >Any help is greatly appreciated. > >Thank you > > >Sarath Banneheka __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Extracting components from a 'boot' class output in R
Sarath, Maybe someone can show a direct access to the statistic. One way to get around is to access 't' component of the boot object which contains individual estimates. so you can extract standard error by sqrt(var(bootObj$t)) Weidong Gu On Fri, Jul 22, 2011 at 11:36 AM, Sarath Gamini Banneheka wrote: > Dear R user, > > I used the following to do a bootstrap. > > >>bootObj<-boot(data=DAT, statistic=Lp.est, > R=1000,x0=3) > > I have the following output from the above bootstrap. How > can I extract components of the output. > For example, how can I extract the std.error? > > >> bootObj > > ORDINARY NONPARAMETRIC BOOTSTRAP > > Call: > boot(data = DAT, statistic = Lp.est, R = 1000, x0 = 3) > > Bootstrap Statistics : > original bias std. error > t1* 794.9745 -0.341 4.042099 > > Any help is greatly appreciated. > > Thank you > > > Sarath Banneheka > > [[alternative HTML version deleted]] > > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Extracting components from a 'boot' class output in R
Dear R user, I used the following to do a bootstrap. >bootObj<-boot(data=DAT, statistic=Lp.est, R=1000,x0=3) I have the following output from the above bootstrap. How can I extract components of the output. For example, how can I extract the std.error? > bootObj ORDINARY NONPARAMETRIC BOOTSTRAP Call: boot(data = DAT, statistic = Lp.est, R = 1000, x0 = 3) Bootstrap Statistics : original bias std. error t1* 794.9745 -0.341 4.042099 Any help is greatly appreciated. Thank you Sarath Banneheka [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.