[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