search for: sav1

...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, 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...