你好，
我正在研究一项在 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 中不同，以获得相同的结果。
有人有什么建议吗？？？
谢谢
