[R] augPred plot in nlme library

From: Afshartous, David <afshart_at_exchange.sba.miami.edu>
Date: Thu 07 Sep 2006 - 22:17:54 GMT


All,  

I'm trying to create an augPred plot in the nlme library, similar to the plot on
p.43 of Pinheiro & Bates (Mixed Effects Models in S and S-Plus) for their Pixel data.  

My data structure is the same as the example but I still get the error msg below.  

> comp.adj.UKV <- groupedData(adj.UKV ~ Time | Patient_no/Lisinopril,
 data = comp.adj.UKV.frm, order.groups = F)  

> fm1comp = lme(adj.UKV ~ Time + Time.sq, data = comp.adj.UKV, random =
list(Patient_no = ~ 1) )  

> plot(augPred(fm1comp, level= 1))

Error in model.frame(formula, rownames, variables, varnames, extras, extranames, :

        variable lengths differ

I've checked all the variale lengths, and have also made sure that factors are correctly defined
as factors. Is there anything special I need to be doing for augPred to work correctly? I checked
the help but didn't find much.  

cheers,
dave          

David Afshartous, PhD
University of Miami
School of Business
Rm KE-408
Coral Gables, FL 33124  

        [[alternative HTML version deleted]]



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 and provide commented, minimal, self-contained, reproducible code. Received on Fri Sep 08 08:25:33 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 07 Sep 2006 - 23:30:06 GMT.

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