[R] bias and std. error from boot

Liaw, Andy andy_liaw at merck.com
Fri Jun 3 00:28:21 CEST 2005


They are calculated in the print() method for boot objects, but not
returned.  boot:::print.boot has:

[...]
            op <- cbind(t0, apply(t, 2, mean, na.rm = TRUE) - 
                t0, sqrt(apply(t, 2, function(t.st)
var(t.st[!is.na(t.st)]))))
            dimnames(op) <- list(rn, c("original", " bias  ", 
                " std. error"))
 
[...]

HTH,
Andy

> From: John Sorkin
> 
> I am running a bootstrap 
> bootResult<-boot(.......)
> and would like to obtain the bias and std. error produced by the
> bootstrap. 
> 
> I can get the values to print to the screen, I can get the original
> (non-bootstrap) estimates, OrgEst<-bootResult$t0
> 
> but I can not find any way to assign the bias and std. error values to
> variables, e.g. 
> bias<-bootResult$bias does not work. 
> 
> Any help would be appreciated.
> 
> Thank you, 
> John 
> 
> 
> John Sorkin M.D., Ph.D.
> Chief, Biostatistics and Informatics
> Baltimore VA Medical Center GRECC and
> University of Maryland School of Medicine Claude Pepper OAIC
> 
> University of Maryland School of Medicine
> Division of Gerontology
> Baltimore VA Medical Center
> 10 North Greene Street
> GRECC (BT/18/GR)
> Baltimore, MD 21201-1524
> 
> 410-605-7119 
> -- NOTE NEW EMAIL ADDRESS:
> jsorkin at grecc.umaryland.edu
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
> http://www.R-project.org/posting-guide.html
> 
> 
>




More information about the R-help mailing list