こんにちは、みんな、
混合モデルでのRとJMPの違いについて詳細な質問がありました。
1。彼女はR(添付)を使用してテストデータを作成しました
###意図 :ランダムインターセプトSD 8
ライブラリ(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")、each = 20)
summary(lmer(yy〜xx +(xx | group1)))
2。分析結果はRで次のようになり、結果は異なります

3。ランダム効果(切片、変数)が大きすぎる
-R:56
-JMP:674
1)なぜこの違いが起こるのかわからなかった
2)2つの値が等しいメーカーに調整または変更するオプションはどれですか?
ありがとう
一樹
をオフにすると、同じ分散と共分散の推定値が得られます中心多項式オプション(デフォルトはオン)。
これらの分散と共分散の推定値について、JMPは標準誤差を報告し、Rは標準偏差を示していることに注意してください。
また、JMPの結果をSAS Proc Mixedの結果と比較しました。彼らは一致します。
これは、保存されたスクリプトを使用したデータです。