373441_ja-JP

liliconcom

JMPとRのVCAの比較

こんにちは、

私はRで実行された研究を見ています、私は次の構文を使用してRでVCAの結果を複製することに成功しました、VisDil1 <-lmer（l2t〜（1 | Op / Day / MDCK）、Dil1）、次の結果

変量効果：

グループ名差異標準偏差

MDCK 日：操作）（切片）0.02309 0.15195

日：Op（切片）0.00861 0.09279

Op（切片）0.00000 0.00000

 

JMPで同じモデルだと思うものを実行すると、

変量効果変動率変動成分合計の割合
MDCK [Day、Op] 0.1904762 0.031746 16.000
日[操作] 0 0 0.000
Op 0 0 0.000
残差0.166666784.000
合計0.1984127100.000

 

l次のスクリプトを使用しました

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]の効果に結合しているので、私は感じます。同じ結果を得るために、Rと同じようにJMPでモデルを指定していません。

 

誰かアドバイスはありますか？

ありがとう

 

分析Re：JMPとRのVCAの比較私が投稿したとき、それはRから、残りの成分を切り落としました、ここに再びあります
変量効果：
グループ名差異標準偏差
MDCK 日：操作）（切片）0.02309 0.15195
日：Op（切片）0.00861 0.09279
Op（切片）0.00000 0.00000
残差0.166700.40829

Re：JMPとRのVCAの比較

デザインの構造とランダム化を明確にしてください。

 

Opは階層の最上位にあります。

 

4日間は両方のOpで異なります。

 

2つのMDCKは、毎日とOpで異なります。

 

力価は潜在的に重要な固定効果ではありませんか？

Re：JMPとRのVCAの比較

ハイマーク、

 

返信ありがとうございます。ネストされたデザインのコピーを添付してください。構文dosentが実験のデザインを反映していることに同意します。より適切なモデルは次のようになります。

RBC
日[RBC]
Op [RBC、Day]

しかし、それは別の話です。私が理解しようとしているのは、RとJMPの両方で同じモデルであると信じているものを実行すると、Day in Rの分散成分が得られるのですが、JMPはこの成分をMDCKに結合しているようです。変数。

変数l2tは、独立変数である力価のlog2です。

 

ありがとう