Hi everyone, I'm running the following analysis for two meadows and two
periods. I need to get the post hoc comparisons for the combination of the two
factors (meadowA*PeriodA, meadowA*PeriodB, meadowB*PeriodB, meadowB*PeriodA). I
use the function glht glht(model, linfct=mcp(Meadow="Tukey") in order to get
results just for the meadow factor. However I need to get the results for both
Meadow*Period. Does anyone know how to do this? Thank you in advance, Vasillis
lme(Yield~1+Meadow*Period, random=~1|Site/Area/Quadrate, control =
lmeControl( maxIter = 10000, # Of course msMaxIter = 5000, # it
helps niterEM = 500, # to be msMaxEval = 500, # patient!
msVerbose = TRUE, # Loquacity ... opt = "optim" # Optional, and #
often useful ),data=Data, weights=varPower(),
correlation=corRatio(form=~1|Site/Area/Quadrate,nugget=T),
method="ML")->model
Vasillis Papathanasiou,
Marine SciencesPalaion Fokon 35-37Thessaloniki, 54454Greecetel:+306945874855
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org 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.
