Hi R-experts, I want to use the jackknife function from the bootstrap package onto a linear model. I can't figure out how to do that. The manual says the following:
# To jackknife functions of more complex data structures, # write theta so that its argument x # is the set of observation numbers # and simply pass as data to jackknife the vector 1,2,..n. # For example, to jackknife # the correlation coefficient from a set of 15 data pairs: xdata <- matrix(rnorm(30),ncol=2) n <- 15 theta <- function(x,xdata){ cor(xdata[x,1],xdata[x,2]) } results <- jackknife(1:n,theta,xdata) So I tried: DF <- as.data.frame(matrix(rnorm(250), ncol=5)) attach(DF) model.lm <- formula(V1 ~ V2 + V3 + V4) n <- 50 theta <- function(x, xdata){ lm(model.lm, data=xdata) } results <- jackknife(1:n,theta,xdata=DF) I was not able to find a similar example on the web, so I would really appreciate some help. TIA, Mark ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.