こんにちは、
REMLを介して最小二乗平均を生成しようとしています。アイデアは、私の研究で複製された各個人に対して単一の値を取得することです(つまり、個人ごとに3つの担当者がいて、1つの数値で個人を表したい)。過去にJMP 14でこれを(同じモデルの同じデータで)実行したとき、lsmeansは算術平均に非常に近かった。現在、JMP 15では、すべての個人のlsmeansは互いに非常に近く、個人ごとの算術平均とは非常に異なります。これを説明するバージョン間で何か変更はありましたか?
助けや提案を事前にありがとう!
JMP 14とJMP 15の両方のFit Model起動ダイアログで同じ設定で分析された同じデータセットがありますが、LSMeansは同じではありませんか?この質問は明白なように思えますが、分析が両方のケースで同じではなかったことは過去に何度も発生しています。おそらく、プラットフォーム用のスクリプトを保存して、ここに投稿できます。その行為は再現性に向かって長い道のりを行くでしょう。
また、唯一の変数が個別の場合、変量効果モデルとREML推定を使用する必要はありません。ANOVAモデルは、有効なLSMeans推定も生成します。
早々に返事をくれてありがとう。恥ずかしいことに、分析は同じではなかったようです...
ただし、ANOVAを使用してlsmeansを生成することに関しては、ANOVAを使用するとREMLを使用した場合とは結果が大きく異なります。特性( "ADI")と効果( "UniqueID")はどちらの分析でも同じです。なぜそうなのかについての洞察はありますか?(これが別のディスカッショントピックでなければならない場合はお詫びします-初めてここに投稿します-また、それが私の側で別の愚かなエラーになった場合もお詫びします!)
REML:
Fit Model(
Y( :ADI ),
Effects,
Random Effects( :UniqueID ),
NoBounds( 1 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Emphasis( "Minimal Report" ),
Run(
:ADI << {Summary of Fit( 1 ), Analysis of Variance( 0 ),
Parameter Estimates( 1 ), Scaled Estimates( 0 ),
Plot Actual by Predicted( 0 ), Plot Regression( 0 ),
Plot Residual by Predicted( 0 ), Plot Studentized Residuals( 0 ),
Plot Effect Leverage( 0 ), Plot Residual by Normal Quantiles( 0 )}
),
SendToReport(
Dispatch( {"Response ADI"}, "Effect Details", OutlineBox, {Close( 0 )} ),
Dispatch(
{"Response ADI", "Effect Details", "UniqueID"},
"Least Squares Means Table",
OutlineBox,
{Close( 0 )}
)
)
);分散分析:
Oneway( Y( :ADI ), X( :UniqueID ), Means( 1 ), Mean Diamonds( 1 ) )
謝罪いたします。「モデルの起動」ダイアログの&UniqueID効果から&Random属性を削除すると、「最小二乗法」でLSMeansも得られるはずです。Onewayプラットフォームを使用する必要があるという意味ではありませんでした。
そして、私がそれについてもっと考えているように、個人が被験者の集団のサンプルを表す場合、個人による変動はランダムな効果として扱われるべきです。LSMeanとその標準エラーは、:UniqueIDを固定効果として扱うかランダム効果として扱うかによって異なります。