Displaying 1 result from an estimated 1 matches for "metab2".
Did you mean:
metabb
2011 Nov 17
1
lme contrast Error in `$<-.data.frame`(`*tmp*`, "df", value = numeric(0)) :
...ceK tSac
tAceK cDMSO
Levels: cDMSO tAceK tSac
> time
[1] 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr 1hr
1hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr 4hr
15m 15m 15m 15m 15m 15m
[43] 15m 15m 15m 15m 15m 15m 15m 15m 15m 15m 15m
Levels: 15m 1hr 4hr
metab2<-data.frame(plate,lnY,gb,time)
fm1<-lme(lnY ~ time*gb,random=~1|plate,metab2,na.action=na.omit)
t1<-contrast(fm1,a = list(gb ="cDMSO" ,time="15m" ),b = list(gb =
"tAceK",time="15m"))
t2<-contrast(fm1,a = list(gb ="cDMSO" ,time="15m...