Bonjour à tous,
Je souhaiterais savoir quelle est la démarche à suivre pour obtenir la p-value associée au modèle linéaire mixte suivant (library lme4):
mod <- lmer(SES_alea~DistM_Prairies+DistM_Cultures+DistM_Bois + (1|data$ID1) + (1|data$ID2), REML=T)
Je ne cherche pas à obtenir les p-values associées à chaque effet fixe / à chaque effet aléatoire, mais à obtenir une p-value pour le modèle entier (comme on peut l'obtenir en faisant un summary d'un modèle établit avec la fonction lme () du package nlme)
J'ai essayé d'utiliser le package lmerTest, mais ce dernier ne me permet que d'obtenir la p-value associée à chaque variable du modèle (même resultat qu'en faisant une Anova(modele)... ce qui n'est pas ma requête)
Merci d'avance pour votre aide,
Léa Uroy