Ciao a tutti,
Ho alcune domande dettagliate su Difference of R e JMP in Mixed Model.
1. Ha fatto i dati del test usando R (allegato)
### Intenzione : Intercettazione casuale SD 8
biblioteca (lme4)
set.seed (3)
group.inter = rnorm (5, sd = 8)
group.slope = rnorm (5, sd = 5)
xx = 4 + rnorm (100, sd = 3)
yy = rep (group.inter, each = 20) + xx * rep (group.slope, each = 20) + rnorm (100, sd = 3)
group1 = rep (c ("A", "B", "C", "D", "E"), ciascuno = 20)
riepilogo (lmer (yy ~ xx + (xx | group1)))
2. Il risultato dell'analisi è come sotto in R, il risultato è diverso

3. L'effetto Randow (Intercept, Var) è troppo diverso
-R: 56
-JMP: 674
1) Non sono riuscito a trovare il motivo per cui questa differenza si verifica
2) Quale opzione deve essere regolata o modificata per creare due valori uguali?
Grazie
Ikju
Otterrai le stesse stime di varianza e covarianza se disattivi Polinomi centrali opzione (l'impostazione predefinita è attiva).
Si noti che per queste stime di varianza e covarianza JMP riporta gli errori standard mentre R mostra le deviazioni standard.
Ho anche confrontato i risultati di JMP con quelli di SAS Proc Mixed. Si abbinano.
ecco i dati con il mio script salvato.