Displaying 1 result from an estimated 1 matches for "r76192".
Did you mean:
r76198
2019 Mar 05
0
Weird problem / bug with augPred() from nlme
...z.ch/pipermail/r-sig-mixed-models/2019q1/027620.html) but this might rather need to go to r-devel anyway, so reposting here:
I came across a weird problem / bug with augPred() from nlme (nlme_3.1-137). Here is a reproducible example to illustrate the issue (tested on R 3.5.2 and R-devel 2019-03-03 r76192):
library(nlme)
dat <- data.frame(id = c(1,1,1,2,3,3,4,4,4,4),
xi = c(1,2,3,1,1,2,1,2,3,4),
yi = c(2,1,4,2,3,2,5,4,6,8),
zi = c(rep("a",9), NA))
res1 <- lme(yi ~ xi, random = ~ 1 | id, data=dat)
sav1 <- augPred(res1, primar...