Displaying 1 result from an estimated 1 matches for "sav2".
Did you mean:
sav
2019 Mar 05
0
Weird problem / bug with augPred() from nlme
...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, primary = ~ xi) # WORKS
dat$zi <- "a"
res2 <- lme(yi ~ xi, random = ~ 1 | id, data=dat)
sav2 <- augPred(res2, primary = ~ xi) # WORKS
dat$zi[10] <- NA
res3 <- lme(yi ~ xi, random = ~ 1 | id, data=dat)
sav3 <- augPred(res3, primary = ~ xi) # ERROR
The error message:
Error in `[[<-.data.frame`(`*tmp*`, nm, value = c(1L, 1L, 1L, 2L, 3L, :
replacement has 10 rows, data h...