Re: [R] Extracting Data form Simpleboot Output

From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>
Date: Tue 14 Feb 2006 - 21:12:25 EST

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



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 Tue Feb 14 21:14:24 2006

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:42:29 EST