373441_zh-CN

liliconcom

JMP和R之间的VCA比较

你好，

我正在研究在R中执行的一项研究，我已经成功使用以下语法VisDil1 <-lmer（l2t〜（1 | Op / Day / MDCK），Dil1）复制了R中的VCA结果以下结果

随机效果：

组名称方差标准偏差

MDCK ：（天：操作）（拦截）0.02309 0.15195

日：运（拦截）0.00861 0.09279

运算（拦截）0.00000 0.00000

 

当我运行我认为与JMP中相同的模型时，我得到了

随机效应可变比率可变分量总数Pct
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中相同的模型来获得相同的结果。

 

有人有建议吗？？？

谢谢

 

分析工具回复：JMP和R之间的VCA比较当我发布它时，它已经从R上切除了剩余成分，这又是
随机效果：
组名称方差标准偏差
MDCK ：（天：操作）（拦截）0.02309 0.15195
日：运（拦截）0.00861 0.09279
运算（拦截）0.00000 0.00000
残留0.16670 0.40829

回复：JMP和R之间的VCA比较

请阐明设计的结构和随机化。

 

Op在层次结构的顶部。

 

两种操作的4天都不一样。

 

每天和Op的两个MDCK不同。

 

滴定度不是潜在的重要固定作用吗？

回复：JMP和R之间的VCA比较

嗨，马克，

 

多谢您的回覆，请随函附上嵌套设计的副本，l同意语法上的内容反映了实验的设计，l认为应该使用更合适的模型

红细胞
日[RBC]
Op [RBC，Day]

但是，这是另一个故事，我试图理解的是为什么当我运行我认为在R和JMP中都相同的模型时，我在Day中得到Day的方差成分，但是JMP似乎将此成分组合到了MDCK中多变的。

变量l2t是效价的对数log2，即独立变量。

 

谢谢