Я пытаюсь сохранить распределение значений R2 при загрузке модели, используя функции ols и bootcov в пакете rms.
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 приводит к той же ошибке. Как я могу решить эту проблему?