Имам следния модел
x <- rep(seq(0, 100, by=1), 10)
y <- 15 + 2*rnorm(1010, 10, 4)*x + rnorm(1010, 20, 100)
id <- NULL
for(i in 1:10){ id <- c(id, rep(i,101)) }
dtfr <- data.frame(x=x,y=y, id=id)
library(nlme)
with(dtfr, summary( lme(y~x, random=~1+x|id, na.action=na.omit)))
model.mx <- with(dtfr, (lme(y~x, random=~1+x|id, na.action=na.omit)))
pd <- predict( model.mx, newdata=data.frame(x=0:100), level=0)
with(dtfr, plot(x, y))
lines(0:100, predict(model.mx, newdata=data.frame(x=0:100), level=0), col="darkred", lwd=7)
с predict
и level=0
мога да начертая средния отговор на населението. Как мога да извлека и начертая 95% доверителни интервали / ленти за прогнозиране от обекта nlme за цялата популация?
curve(predict(model.lm, data.frame(x=x),interval ='confidence'),add=T)
, където model.lm например е lm(y~x) - person agstudy   schedule 16.01.2013intervals .lme
, но тя не дава на групата увереност само една точка. - person agstudy   schedule 16.01.2013intervals
получава CI на оценките/коефициентите на съвпаденията. Какви са ми необходими CI на y за всяко дадено x. - person ECII   schedule 16.01.2013