RE: [R] bias and std. error from boot

From: Liaw, Andy <andy_liaw_at_merck.com>
Date: Fri 03 Jun 2005 - 08:28:21 EST


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@grecc.umaryland.edu
>
> ______________________________________________
> R-help@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
>
>
>



R-help@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 Received on Fri Jun 03 08:32:37 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:22 EST