[R] Parsing results from boot
Scott Norton
nortonsm at verizon.net
Tue Jun 15 06:09:03 CEST 2004
This probably has a super easy answer...but I claim newbie status! (I did
search help lists but this question is hard to isolate keyword-wise)
Basically, I'm trying to figure out how to parse the results from executing
boot(). I'm mainly interested in assigning the standard error estimate to a
scalar variable.
For example:
--------------------------------------------------------------
> b<-c(100,100,120,130,1000,1200,1100,1150,125)
> b.boot <- boot(b, function(b,i) median(b[i]), R=1000)
> b.boot
ORDINARY NONPARAMETRIC BOOTSTRAP
Call:
boot(data = b, statistic = function(b, i) median(b[i]), R = 1000)
Bootstrap Statistics :
original bias std. error
t1* 130 354.75 450.9763
>
------------------------------------------------------------------
I'm interested in the value for std.error (i.e. 450.9763).
Now executing the folling:
------------------------------------------------------------------
> summary(b.boot)
Length Class Mode
t0 1 -none- numeric
t 1000 -none- numeric
R 1 -none- numeric
data 9 -none- numeric
seed 626 -none- numeric
statistic 1 -none- function
sim 1 -none- character
call 4 -none- call
stype 1 -none- character
strata 9 -none- numeric
weights 9 -none- numeric
>
-------------------------------------------------------------------
seems to imply that it's not actually assigned to an output variable...but
it is shown when I just type "b.boot"
Where is the std.error being assigned? and how to I access it?
Thanks in advance...
-Scott
More information about the R-help
mailing list