Bonjour à tous,
Je suis en train d'effectuer des ANOVA suivies de tests post-hoc de comparaison de moyennes (tests de bonferroni).
Lors de mes comparaisons de moyennes, j'obtiens des pvalues très proches de 0.05, qui "m'empêchent" de regrouper certains traitements au sein de même groupes.
Ayant une puissance de 95%, j'aurais donc voulu le descendre à 90% afin de pouvoir tout de même exploiter ces différences entre groupes.
Cependant, je ne trouve pas la commande R me permettant de le faire... Quelqu'un aurait-il une piste pour moi ?
Voilà mon script actuel :
lm4 <- lm(Structure.maintenance~Treatment, data=dataSQI)
shapiro.test(resid(lm4))
leveneTest(resid(lm4)~Treatment,data=dataSQI)
Anova(lm4, type=3)
summary(glht(lm4,linfct=mcp(Treatment="Tukey"),adjust='bonferroni'))
grp1 = cld(glht(lm4,linfct=mcp(Treatment="Tukey"),adjust='bonferroni'));grp1
D'avance merci,
Clara