[R] nlme plot residuals per group

From: Osman Al-Radi <osman.al.radi_at_gmail.com>
Date: Mon 08 May 2006 - 00:51:32 EST


dear list:

I used the nlme library according to the great Pinheiro/Bates book, on R2.3, WinXp

Lac.lme is an lme object with unbalanced data, group is a factor variable with three levels, when I tried to plot the residuals by group I got this error msg:

>plot(Lac.lme,resid(.,type='p')~fitted(.)|group)
Error in limits.and.aspect(prepanel.default.xyplot, prepanel = prepanel, :

        need at least one panel

Also When I try to use the auPred() function I get the follwoing error msg:
>plot(augPred(Lac.lme))

Error in tapply(as.character(object[[nm]]), groups, FUN[[dClass]]) :

        arguments must have same length

Any suggestions?

Thanks

--
Osman O. Al-Radi, MD, MSc, FRCSC
Fellow, Cardiovascular Surgery
The Hospital for Sick Children
University of Toronto, Canada

______________________________________________
R-help@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
Received on Mon May 08 00:56:36 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 11 May 2006 - 12:10:06 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.