Ken,
> lme.fnl.REML <- lme(Max ~ S + Tr + Yr + Tr:Yr, random = ~1 |TID,
> method = "REML")
>
You have three variables in the model: S, Tr and Yr.
> tst <- contrast(lme.fnl.REML, a=list(Yr=levels(Yr), Tr="R"), b=list
> (Yr=levels(Yr, Tr="T97")))
> Error in gendata.default(fit = list(modelStruct = list(reStruct = list
> ( :
> not enough factors
>
I think that you need to specify a value for S in the lists for a and b.
Also, I think that there is a typo:
b=list(Yr=levels(Yr), Tr="T97")
instead of
b=list(Yr=levels(Yr, Tr="T97"))
> Trying to include the other factor to make R happy....
>
> > tst <- contrast(lme.fnl.REML, a=list(Yr=levels(Yr), Tr="R"), b=list
> (Yr=levels(Yr, Tr="T97")), c=list(Yr=levels(Yr, Tr="T98"))
> + )
> Error in contrastCalc(fit, ...) : argument 4 matches multiple formal
> arguments
Max
----------------------------------------------------------------------
LEGAL NOTICE\ Unless expressly stated otherwise, this messag...{{dropped}}
______________________________________________
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.
