Displaying 1 result from an estimated 1 matches for "someweek".
Did you mean:
someweeks
2006 Apr 06
1
polynomial predict with lme
...<- 1:nweeks
mcurve <- 25 + 0.819*week - 0.0588*week^2 + 0.000686*week^3
cow.eff <- rnorm(ncows)
week <- rep(week, ncows)
cow <- gl(ncows,nweeks)
yield <- mcurve + cow.eff[cow] + rnorm(ncows*nweeks)
lmefit <- lme(yield ~ poly(week,3), random = ~1|cow)
summary(lmefit) # seems OK
someweeks <- seq(5,45,by=5)
new <- data.frame(week=someweeks)
predicts <- predict(lmefit, new, level=0)
print(predicts) # not even close
#plot(week, yield, las=1)
#lines(someweeks, predicts)
###############
--
************************************************
* I.White...