你好，
我正在研究在 R 中進行的一項研究，我已成功使用以下語法在 R 中複製 VCA 的結果，VisDil1 <- lmer(l2t ~ (1|Op/Day/MDCK) , Dil1)，對於以下結果
隨機效果：
組名稱變異數標準差
MDCK：（日：Op）（截距）0.02309 0.15195
日：Op（截距）0.00861 0.09279
開啟（截距） 0.00000 0.00000
當我在 JMP 中運行我認為是相同的模型時，我得到
隨機效應 Var 比率 Var 成分 佔總數的百分比
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 中相同的模型以獲得相同的結果。
有人有什麼建議嗎？
謝謝
