search for: sav2

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