こんにちは、
私はRで行われた研究を見ています。私は、次の構文、VisDil1 <- lmer(l2t ~ (1|Op/Day/MDCK) , Dil1)を使用して、RでVCAの結果を再現することができました。次の結果が得られました。
ランダム効果:
グループ名 差異 標準偏差
MDCK:(日:Op) (切片) 0.02309 0.15195
日:Op (切片) 0.00861 0.09279
オン (インターセプト) 0.00000 0.00000
JMPで同じモデルを実行すると、
ランダム効果 分散率 分散成分 全体に占める割合
MDCK[日、操作] 0.1904762 0.031746 16,000
日[日] 0 0 0.000
0 0 0,000 時
残存 0.1666667 84,000
合計 0.1984127 100,000
私は次のスクリプトを使用しました
Fit Model(
Y( :l2t ),
Effects,
Random Effects( :MDCK[:Day, :Op], :Day[:Op], :Op ),
NoBounds( 0 ),
Estimate Only Variance Components( 1 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Run Model
);
これを 8 回の希釈で調べたところ、各希釈で同じ結果になりました。R には Day:Op の分散コンポーネントがありますが、JMP はこの効果を MDCK[Day,Op] の効果に組み合わせているため、同じ結果を得るために JMP で R と同じモデルを指定していないように感じます。
誰かアドバイスはありますか？
ありがとう
