Tengo problemas para encontrar una solución con respecto a cómo ejecutar una prueba post-hoc (Tukey HSD) después de un ANOVA de medidas repetidas de 2 factores (ambos dentro de los sujetos) en R. Para el ANOVA, he usado la función aov:
summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1))
Después de leer las respuestas a otras preguntas, deduje que primero tendría que volver a ejecutar el ANOVA usando alguna otra función (por ejemplo, lme). Esto es lo que se me ocurrió.
Lme.mod <- lme(dv ~ x1*x2, random=list(subject=pdBlocked(list(~1, pdIdent(~x1-1), pdIdent(~x2-1)))), data=df1)
anova(Lme.mod)
Ambos efectos principales fueron significativos, pero no hubo efectos de interacción. Luego, utilicé estas funciones para las comparaciones post-hoc:
summary(glht(Lme.mod, linfct=mcp(x1="Tukey")))
summary(glht(Lme.mod, linfct=mcp(x2="Tukey")))
Sin embargo, hubo algunos problemas:
En primer lugar, el archivo de Ayuda R dice que "La función mcp debe usarse con cuidado al definir parámetros de interés en los modelos ANOVA o ANCOVA de dos vías (...) multcomp versión 1.0-0 y superior genera comparaciones para los efectos principales solo, ignorando las covariables y las interacciones (las versiones anteriores promediaron automáticamente sobre los términos de interacción). Se da una advertencia ". Y efectivamente, recibí el siguiente mensaje de advertencia:
Warning message:
In mcp2matrix(model, linfct = linfct) :
covariate interactions found -- default contrast might be inappropriate
Otra cosa desconcertante fue que aunque ambos efectos principales fueron significativos, no hubo diferencias significativas en las comparaciones post-hoc para uno de los factores (x1). Nunca me he encontrado con esto antes. ¿Los guiones / análisis son correctos / apropiados, o hay algo que me falta? ¡Cualquier ayuda sería muy apreciada!