Displaying 1 result from an estimated 1 matches for "fm1comp".
Did you mean:
fcomp
2006 Sep 07
5
augPred plot in nlme library
...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 als...