Извиняюсь за, наверное, действительно простой вопрос. Я должен повторить следующий код N = 2000 раз, составить гистограмму оценок N = 2000 b1 и вычислить долю рассчитанных доверительных интервалов N = 2000, которые фактически покрывают истинное значение для X = 14. У меня нет проблем с репликацией данных 2000 раз, но я сталкиваюсь с проблемами, пытаясь повторить что-то еще.
X <- c(4,8,10,12,16,20)
Yhat <- c(36,52,60,68,84,100)
error <- (rnorm(6,0,5))
Y <- Yhat+error
mydata <- data.frame(X, Yhat, error, Y)
mydata
attach(mydata)
mydata.lm <- lm(Y~X)
mydata.lm
newdata <- data.frame(X=14)
predict(mydata.lm, newdata, level=0.95, interval="confidence")
detach(mydata)
Заранее спасибо за помощь. Я полный новичок в области R и статистики.