R语言实战12.6.1章中队单个统计量使用自助法里,具体函数如下:
rsq <- function(formula, data, indices){
d <- data[indices,]
fit <- lm(formula, data = d)
return(summary(fit)$r.square)
}
library(boot)
set.seed(1234)
result <- boot(data = mtcars, statistic = rsq, R = 1000, formula = mpg ~ wt +disp)
print (result)
plot(result)
其中 d <- data[indices,]这一项书上写是boot函数中statistic中必须包括indices函数,不太明白有什么作用,而且在跑完数据后d也不能单独输出


雷达卡




京公网安备 11010802022788号







