In the model: lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera)
I found TIL*YEAR interaction significant. Then I am trying to compare means of the different levels of TIL inside every YEAR using: mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey")) summary(mc.2, test = univariate()) but it does not work. There is any way of doing this, like the SLICE option in PROC GLM (SAS)? Thanks a lot, Jorge -- ************************************************** Jorge Lampurlanés Castel Departament d'Enginyeria Agroforestal Escola Tècnica Superior d'Enginyeria Agrària Universitat de Lleida Avinguda Rovira Roure, 191 25198-LLEIDA SPAIN Tl.: +34 973 70 25 37 Fax.:+34 073 70 26 73 e-mail: [EMAIL PROTECTED] ______________________________________________ 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.
