[R] Multiple comparisons when interacction

[Jorge Lampurlanes Castel]

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: jlampur at eagrof.udl.es