[R] augPred and NAs

[Morten Gram Pedersen]

I have a problem using augPred from the nlme package on a data set with NAs.

Consider this example, modified from the augPred help page:

> library(nlme)
> Orthodont[100,]$age=NA     #### Insert NA somewhere in the data set
> fm1 <- lme(Orthodont, random = ~1,na.action=na.exclude)  #### Can still fit with lme
> augPred(fm1)    #### Error...!? Also with options "na.rm=T" or  "na.action=na.exclude"
Error in if (from == to) rep.int(from, length.out) else as.vector(c(from,  : 
  missing value where TRUE/FALSE needed

Any help on how to handle NAs in such cases is most appreciated. Neither Pinheiro & Bates nor Google helped...

Thanks,
Morten Pedersen