je travaille actuellement sur une base de données correspondant à la mesure de la croissance en hauteur de plants (petits arbres), plantés dans différentes manières de travailler le sol. L'objectif est donc de mettre en lumière quelle est la meilleure manière de travailler le sol.
Pour ce faire, j'utilise un modèle linéaire mixte, que voici:
Code : Tout sélectionner
BL.lme=lme(haut~-1+mod, data=data, random=~1|parcelle)
-effet fixe: mod, correspondant à la modalité (manière de travailler le sol)
-effet aléatoire: parcelle (il y a eu différentes répétitions, sur différentes parcelles).
Ce que je veux savoir, c'est comment sortir une table des p-value de comparaison 2 à 2 de mes différentes modalités (mod) de travail du sol.
Ceci est possible avec la fonction "TurkeyHSD" (
Code : Tout sélectionner
TukeyHSD(aov(haut~mod, data=data))
Existe t-il une fonction similaire mais prenant en compte les modèles mixtes ?
Merci d'avance
Charles