Displaying 1 result from an estimated 1 matches for "mass_f".
Did you mean:
mass_
2011 Mar 18
0
predict.nlme
Hi folks,
I am having trouble to plot a mixed model analysis of covariance (ANCOVA).
To do so I use the function predict in nlme but the line that is being drawn
is totally out of control!!!
here is my script (where MASS_S is dry mass and MASS_F is fresh mass):
MEN<-read.table("Mentha_lme2.txt", h=T)
attach(MEN)
lme1<-lme(log(MASS_S)~log(MASS_F)*TREAT, random=~1|INDIV)
summary(lme1)
anova.lme(lme1)
plot(log(MASS_F), log(MASS_S), type="n")
points((log(MASS_F)[TREAT=="A"]), (log(MASS_S)[TREAT=="A&quo...