Exemple pour la fonction Ajuster le modèle à effets mixtes

Un chercheur teste les rendements de six variétés de luzerne sur quatre terrains choisis aléatoirement. Le rendement de chaque variété a été enregistré pour chaque terrain.

Le chercheur souhaite savoir si la variété de la luzerne influe sur le rendement moyen. Il dispose de 4 terrains où collecter les données. Cependant, il souhaite modéliser la manière dont pousse la luzerne sur des terrains qui ne sont pas pris en compte dans l'expérience. C'est pourquoi il décide que le terrain où pousse la luzerne est un facteur aléatoire. Il utilise un modèle à effets mixtes pour évaluer à la fois les effets fixes et les effets aléatoires.

  1. Ouvrez le fichier de données échantillons, Alfalfa.MTW.
  2. Sélectionnez Stat > ANOVA > Modèle à effets mixtes > Ajuster le modèle à effets mixtes.
  3. Dans la zone Réponses , saisissez Rendement.
  4. Dans la zone Facteurs aléatoires (obligatoires) , saisissez Terrain.
  5. Dans la zone Facteurs fixes , saisissez Variété.
  6. Cliquez sur Graphiques.
  7. Dans Valeurs résiduelles des graphiques , sélectionnez Conditionnel, normalisé.
  8. Dans Graphiques des valeurs résiduelles, sélectionnez Quatre en un.
  9. Cliquez sur OK dans chaque boîte de dialogue.

Interprétation des résultats

Dans le tableau des composantes de la variance, la valeur de p pour Terrain est égale à 0,124. Le test d'hypothèse ne permet pas de prouver que la composante de la variance est différente de 0. La valeur de p de la composante de variance pour l'erreur est 0,003. Sachant que la valeur de p est inférieure au seuil de signification de 0,05, le chercheur peut conclure que la composante de variance pour l'erreur n'est pas égale à 0.

La valeur de p d'environ 0 pour le facteur fixe, Variété, indique qu'au moins un type de luzerne a un effet sur le rendement significativement différent des cinq autres types.

Les coefficients des effets principaux représentent la différence entre chaque moyenne de niveau et la moyenne générale. Par exemple, le facteur Variété 1 est associé à un rendement de la luzerne d'environ 0,385 unité de plus que la moyenne générale. La valeur de p d'environ 0 pour ce coefficient indique que l'effet de Variété 1 sur le rendement est significativement différent de l'effet d'un autre niveau du terme Variété. Pour déterminer les effets de niveaux qui sont statistiquement égaux et ceux qui sont statistiquement différents, le chercheur envisage d'effectuer une analyse de comparaisons multiples pour le terme.

La valeur R2 indique que le modèle explique environ 92 % de la variation du rendement. La valeur R carré (ajust) aussi est élevée, avec environ 90,2 %. Le chercheur utilise cette valeur pour comparer des modèles ayant un nombre de prédicteurs différent.

Les observations 1 et 5 sont aberrantes car elles contiennent des valeurs résiduelles normalisées supérieures à 2. Le chercheur examine les données pour vérifier que les valeurs de réponse pour ces observations sont correctes.

Les valeurs résiduelles sur la droite de Henry se rapprochent d'une ligne droite et les points semblent être répartis aléatoirement autour de 0 sur le diagramme des valeurs résiduelles en fonction des valeurs ajustées.

Méthode

Estimation de la varianceMaximum de vraisemblance restreint
DL pour effets fixesKenward-Roger

Informations sur les facteurs

FacteurTypeNiveauxValeurs
TerrainAléatoire41; 2; 3; 4
VariétéFixe61; 2; 3; 4; 5; 6

Composantes de la variance

SourceVar% du totalVar. ErTValeur de ZValeur de p
Terrain0,07791972,93%0,0675801,1529960,124
Erreur0,02892427,07%0,0105622,7386130,003
Total0,106843       
Log de vraisemblance -2 = 7,736012

Tests pour effets fixes

TermeDL NumDL DénValeur FValeur de p
Variété5,0015,0026,290,000

Récapitulatif du modèle

SR carréR carré
(ajust)
AICcBIC
0,17007192,33%90,20%12,5413,52

Coefficients

TermeCoeffCoef ErTDLValeur de TValeur de p
Constante3,0945830,1438223,0021,5166920,000
Variété         
  10,3854170,07762615,004,9650160,000
  20,1454170,07762615,001,8732870,081
  30,1079170,07762615,001,3902050,185
  4-0,3195830,07762615,00-4,1169380,001
  50,3954170,07762615,005,0938380,000

Ajustements et diagnostics marginaux pour les observations aberrantes

ObservationRendementValeur
ajustée
RésiduelleVal. résid.
norm.
14,1000003,4800000,6200002,190221R
54,2200003,4900000,7300002,578808R
R : Valeur résiduelle élevée

Ajustements et diagnostics conditionnels pour les observations aberrantes

ObservationRendementValeur
ajustée
RésiduelleVal. résid.
norm.
54,2200003,8953390,3246612,400733R
R : Valeur résiduelle élevée