[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