Salut à tous,
J'ai une question détaillée sur la différence entre R et JMP dans le modèle mixte.
1. Elle a créé des données de test avec R (attaché)
### Intention : Interception aléatoire SD 8
bibliothèque (lme4)
set.seed (3)
group.inter = rnorm (5, sd = 8)
groupe.slope = rnorm (5, sd = 5)
xx = 4 + rnorm (100, sd = 3)
yy = rep (group.inter, chaque = 20) + xx * rep (group.slope, chaque = 20) + rnorm (100, sd = 3)
group1 = rep (c ("A", "B", "C", "D", "E"), chacun = 20)
résumé (lmer (yy ~ xx + (xx | groupe1)))
2 Le résultat de l'analyse est comme ci-dessous dans R, le résultat est différent

3 Randow Effect (Intercept, Var) est trop différent
-R: 56
-JMP: 674
1) Je ne trouvais pas pourquoi cette différence se produisait
2) Quelle option doit être ajustée ou modifiée pour que le fabricant crée deux valeurs égales?
Merci
Ikju
Vous obtiendrez les mêmes estimations de variance et de covariance si vous désactivez le Polynômes de centre option (la valeur par défaut est).
Notez que pour ces estimations de variance et de covariance, JMP signale les erreurs-types tandis que R indique les écarts-types.
J'ai également comparé les résultats de JMP à ceux de SAS Proc Mixed. Ils correspondent.
voici les données avec mon script sauvegardé.