Dear all,

I run this model:
mod=lme(ratio_tige~Trait*Prlvt,random=~1|Plante,na.action=na.omit,data=dataAna)


I try to assess the assumptions on the within group error. I follow the
procedure described in Pinhero&Bates (2000).
To assess the normality and independence of errors, authors suggest to use
(page 175-176):

fm20rth.lm=distance~Sex*I(age-11),data=Orthodont,random=~1|Subject)
plot(fm20rth.lm,resid(.,type="p")~fitted(.)|*Sex*,id=0.05,adj=-0.03)
In their model:
distance=numeric
age=numeric
Subject=ordered factor
*Sex*=factor

When I run this procedure on my model, I get an error message:

plot(mod1,resid(.,type="p")~fitted(.)|*Trait*,id=0.05,adj=-0.03)
Erreur dans `[[<-.data.frame`(`*tmp*`, j, value = integer(0)) :
  replacement has 0 rows, data has 1900

In my model:
ratio_tige=numeric
*Trait*=factor
Prlvt=factor
Plante=ordered factor

Can someone please help ??

Thank you in advance,

 Cheers,

Amélie

-- 
Dr. Pinet Amélie
INRA, Plantes et Systèmes de culture Horticoles
Domaine St Paul, Site Agroparc, 84941 Avignon Cedex 9, France

	[[alternative HTML version deleted]]


