SAS PROC MIXED と JMP の出力が異なる<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />親愛なるみんな、
SAS PROC 混合モデルと JMP 混合モデルの結果を比較しようとしています。 結果は同じだと思っていましたが、そうではありませんでした。 JMP で SAS PROC MIXED と同じ結果が得られることを願っています。 誰でもそれを手伝ってもらえますか？
次のような私のSASコード：
PROC MIXED DATA = DATA;
CLASS サブ条件刺激。
モデル Y = 状態 / 解決策;
ランダムサブ;
ランダム刺激（条件）;
走る;
固定効果の結果は次のとおりです。
| 固定効果のソリューション |
|---|
| 効果 | 調子 | 見積もり | 標準
エラー | DF | t値 | Pr > |t| |
|---|
| インターセプト | | 4.4444 | 1.0943 | 4 | 4.06 | 0.0153 |
|---|
| 調子 | あ | 1.1111 | 1.5476 | 4 | 0.72 | 0.5125 |
|---|
| 調子 | B | 0 | . | . | . | . |
|---|
| 固定効果のタイプ 3 検定 |
|---|
| 効果 | DF数 | デンDF | F値 | Pr > F |
|---|
| 調子 | 1 | 4 | 0.52 | 0.5125 |
|---|
次に、JMP のインターフェイスを使用して次のコードを生成しました。
モデルの適合(
Y( :Y ),
効果( :状態 ),
Random Effects( :Sub, :Stimuli[:Condition] ),
制限なし( 0 ),
パーソナリティ（「混合モデル」）、
Run(反復効果共分散パラメータ推定値( 0 ) )
)
私は得た：
固定効果のタイプ 3 検定は JMP と SAS で同じように見えますが、私の固定効果パラメーターの推定値はまったく異なります。 SASと同じ結果になるようにJMPをコーディングする方法を知っている人がいるかどうか疑問に思っています。
皆様、どうもありがとうございました。
アリシア
Re: SAS PROC MIXED と JMP の出力が異なるねえ、ロン、アリシア: モデルの適合プラットフォームの「推定」サブメニュー (「実行」をクリックした後の分析の上部にある赤い三角形) で「指標パラメーター化推定」を選択する方が簡単な方法かもしれません。 これらの推定値をデフォルトで (JMP のデフォルトに加えて) 表示するように選択することもできます。これは、[設定] の [最小二乗法に適合] プラットフォームで行うことができます。
マーク
Re: SAS PROC MIXED と JMP の出力が異なるねえ、アリシア、
うまくいったことをうれしく思います。最後の方法はハックまたは配管工のトリックですが、基本に忠実です。以前は、カテゴリ変数を回帰モデルに入力する実際の方法でした。カテゴリを 0 と 1 の n-1 個の新しい列に再コーディングすることで、自分で対比を作成します。プログラムは、それらが推定目的で連続的であると考えるだけです。
一番、
ロン
Re: SAS PROC MIXED と JMP の出力が異なるやあロン、
どうもありがとう。 あなたが提案した3番目の方法を使用しましたが、うまくいきます!!! どうもありがとう。
一番、
アリシア
Re: SAS PROC MIXED と JMP の出力が異なる笑い
次の 3 つのことを試してください。
1) 他のケースでは、カテゴリ変数を名義型から順序型モデリング タイプに変更します。ただし、この場合、モデル プラットフォームはそれを受け入れないと思います。
2) 探しているコントラストを「明らかに」するために、これらのパラメーターのカスタム テストを実行できるかどうかを確認してください。
3) カテゴリが 2 つだけの場合は、カテゴリを 0 と 1 としてコーディングし、連続として実行して何が得られるかを確認してください。
ロン
Re: SAS PROC MIXED と JMP の出力が異なるやあロン、
どうもありがとう。 SAS と JMP が名義データの対比を異なる方法で処理する方法については、あなたが正しいと思います。 とにかく、JMP のデフォルトのコントラストを変更できることをご存知でしょうか? SAS の結果の方が妥当に見えると思います (標準誤差は JMP の固定効果係数ではすべて同じです)。そのため、JMP のデフォルトのコントラストを変更したいと思います。 両方を一致させたいのは、固定効果係数のベータ推定値を使用してモデルの R^2 を計算したいからです。 ベータ値が同じでない場合、R^2 は異なります。 本当にありがとうございました。
心から、
アリシア
Re: SAS PROC MIXED と JMP の出力が異なるこんにちは笑い、
p値は同じなので、分析は同じだと思います。したがって、この違いは、カテゴリの対比に対する異なるアプローチに基づくパラメータの JMP 推定によるものです。名義変数と順序変数の対比の定義の違いを確認してください。
このリンクに従ってくださいhttp://www.jmp.com/support/help/The_Factor_Models.shtml
名義因子を探します。この表は、パラメーターがどのように扱われるかを説明しています。
ロン
