我正在尝试使用rms软件包中的ols和bootcov函数在引导模型时保存R2值的分布。
library(rms)
#Generate data
x = seq(1,100)
y = 0.2*x + rnorm(100)
#Model fitting and bootstrap attempt
mod = ols(y~x, myData, x=T, y=T)
mod.boot = bootcov(mod, B=1000, stat="R2")
但是这将返回错误:
Error in stats[b] <- f$stats[stat] : replacement has length zero
尝试保存mod$stats
的其他元素会产生相同的错误。我该如何解决?