混合モデルで列を保存オプションみなさん、こんにちは。
Fit Modelプラットフォームを使用して混合モデルを実行すると、残差と予測値を保存できますが、条件付き残差と条件付き予測値(固定効果とランダム効果の両方を考慮に入れる)を保存することもできます。ただし、スチューデント化残差列とStdErr Pred Formula列には、同様のオプションはありません。
これら2つの列の「条件付きの理由」を取得または再計算するための最善の方法は何ですか?
私はJMP Proを使用していないことに注意してください。
ご協力ありがとうございます。
混合モデルで列を保存オプションはい、その通りです。
混合モデルで列を保存オプションこの機能はJMPコミュニティのウィッシュリストの部分で提案できます。
混合モデルで列を保存オプションGMさん、こんにちは。
そうです、私たちはそのような変量効果の標準誤差を抽出することができますが、私が探しているのは
残差と予測値の標準誤差 私はそれがここで同じ情報だとは思わない。
混合モデルで列を保存オプションこの詳細な回答をありがとうキャメロン、私はそれをテストします!
しかし確かにそれは私たちが欲しいものを手に入れるための簡単な方法ではありません。列保存メニューから直接利用できないのはなぜでしょうか。
混合モデルで列を保存オプションPersonality = Standard Least SquaresおよびMethod = REMLでFit Modelプラットフォームを使用して変量効果を推定する場合、変量効果予測テーブルからデータテーブルを作成することによって変量効果の標準誤差を抽出できます。 混合モデルパーソナリティは、レポートにランダム効果係数を含む標準誤差を出力しません。
混合モデルで列を保存オプションあなたが正しい。 ええと。 JMP Proにアクセスできないので、チェックアウトすることはできませんが、標準化された条件付き残差を計算するために必要なすべてのコンポーネントをレポートから直接取得することはできないと思います。 PredSE列はVec Quadratic()の2番目の引数としてX行列を示し、レポートはMSEを示します。 あなたはGとZを構築しなければならないでしょう。
Zは次のようなものから得られます。
Design Nom( :RandomEffectCol << Get Values, {list of levels}) 複数のランダム効果がある場合は、それらを連結するだけで構いません。
Gを得るには、それはあなたが選んだ相関構造に依存します。 非対角上にすべて0がある単純な相関構造がある場合は、分散成分のベクトルを単位行列で乗算することによってGを構築できます。
X、G、Z、およびMSEを取得したら、次のように誤差行列の分散を計算できます。
n = N Row(dt);
V = Z*G*Z` + MSE*Identity(n);
K = Identity(n) - Z*G*Z`*Inv(V);
Q = X*Inv(X`*Inv(V)*X)*X`;
Var_e = K*(V-Q)*K`;
そこから、Var_eの対角要素のsqrtを取り、条件付き残差の条件付き標準誤差を取得します。 JMPには行列の対角要素だけを返す関数はないと思いますが、ここでは簡単なことをする簡単なものを紹介します。
getDiag = function({X},
n = ncol(X);
diags = {};
for(i=1, i<=N Col(X),i++,
Insert Into(diags,X[i,i])
)
);これは大変な作業です。すみません、もっと簡単な方法を知っていたらいいのにと思います。
混合モデルで列を保存オプションご回答ありがとうございますキャメロン!
あなたが残差の標準誤差を保存した場合しかし私はそれがランダム効果を含まないものと条件付き残差ではないものに関連していると思う。残差の条件付き標準誤差を保存することはできないようです。Hatsの値に基づいて計算をやり直すことができると思いますが、問題は同じです。それらはおそらく「基本的な残差」に関連しています。
混合モデルで列を保存オプション標準誤差で除算することで条件付き残差をスチューデント化できます。 両方の列をテーブルに保存してから、列式を使用してスチューデント化された値を計算します。 Std Error Pred Formulaでどのようにあなたを助けることができるかすぐにわかりません。