[R] Extracting Data form Simpleboot Output
Uwe Ligges
ligges at statistik.uni-dortmund.de
Tue Feb 14 11:12:25 CET 2006
Terry W. Schulz wrote:
> Hello,
> I am new to R, so forgive my ignorance on what is probably simple.
> I find package simpleboot quite useful for LOESS bootstrapping and
> generation of plots.
> I want to calculate the standard error for x=60 of the 100 grid points
> and 50 bootstraps.
> The code below gives the first fitted value.
> How can I grab the other 49 values easily?
> I could do it one at a time for 50 bootstraps, but this becomes tedious
> for say 2000 bootstraps.
> Thanks for any help.
>
> library(simpleboot)
> library(bootstrap)
> attach(cholost)
> set.seed(13579)
> lo <- loess(y ~ z, data=cholost, span=.5, degree=2, family="gaussian",
> method="loess")
> lo.b <- loess.boot(lo, R=50, rows = TRUE, new.xpts = NULL, ngrid = 100,
> weights = NULL)
> lo.b$boot.list$"1"$fitted[60]
>
sapply(lo.b$boot.list, function(x) x$fitted)
gives you the corresponding 100x50 matrix ...
Uwe Ligges
More information about the R-help
mailing list