[R] augPred with lme(...,subset=...)

[Manuel Morales]

Hello,

Is there a way to get augPred to work with lme if a subset statement is
used? For example, if I modify the example from ?augPred.lme to include
a subset statement, I get the following error:


fm1 <- lme(Orthodont, random = ~1, subset=distance>19)
augPred(fm1, length.out = 2, level = c(0,1))
Error in tapply(object[[nm]], groups, FUN[["numeric"]], ...) :
        arguments must have same length

Thanks!
Manuel