R help - Nov 2004 - augPred with lme(...,subset=...)

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

Manuel Morales wrote:
> 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
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
Well you could debug the augPred method and send the package maintainer 
a patch :-)