373441_zh-TW

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，即獨立變量。

 

謝謝