2つの予測（重回帰モデル経由）が有意に異なるかどうかをテストする私は、JMP 16 で多重線形回帰モデルに取り組んでおり、以下の項目に基づいて脳卒中患者の臨床転帰 (「退院 mRS 」と呼ばれる数値スケールの形式) を予測しようとしています。
- 1つの連続予測因子（入院時のNIHSS 、これは入院時の脳卒中の重症度を測定する指標です）
- 1 つのカテゴリ予測子（ EVT タイミング、これは脳卒中患者が受ける可能性のある治療です）。次の 3 つの値のいずれかになります。
- 1 = EVT（治療）は病院到着後すぐに実施された
- 2 = EVTは遅れて実施された
- 0 = EVTは実施されなかった
- 相互作用項（入院時のNIHSS * EVTのタイミング）
現在、JMP でのモデルは次のようになります。
私の質問/問題は:
入院時の NIHSSの任意の値について、患者が直ちに EVT を受ける場合 (つまり、 EVT タイミング = 1 ) と EVT を受けない場合 (つまり、 EVT タイミング = 0 )の間で、予測される平均結果 (このモデルに基づく) が (p 値を使用して)有意に異なるかどうか、つまり、カテゴリ予測子の 2 つの異なる値の間で異なるかどうかを計算します。
JMP でこの分析を実行する方法はありますか?
(概念的には、私の目標は、このモデルの 2 つの予測因子を臨床的に有用な方法で利用できるようにすることです。たとえば、モデルが、治療による結果を予測して、どの患者 (この場合は脳卒中の重症度に基づいて) に治療を施すべきか、施すべきでないかを決定するのに役立つことを期待しています。)
ご協力いただきありがとうございます！
基本的なデータ分析とモデリングウィンドウズRe: 2 つの予測値 (重回帰モデル経由) が有意に異なるかどうかをテストするこんにちは@MRB3855 さん
本当にありがとうございます。うまくいったと思います。もちろん、自分でテストの出力を確認するための計算はしていません（やり方も知りません）が、カスタム テストの予測値がモデルの予測値と一致し、テストの出力が意味をなすようになりました。
本当に助かりました！
Re: 2 つの予測値 (重回帰モデル経由) が有意に異なるかどうかをテストするこんにちは@broccoli : 申し訳ありません、間違いを犯しました。時間 = 24 での予測については、以下を参照してください (私のデータの場合)。相互作用項には、(入院時の NIHSS - 14.7059) と入力する必要があります。したがって、入院時の NIHSS = 20 で予測する場合は、相互作用として 20-14.7059 = 5.2941 を入力する必要があります。
リストされていないバッチ（バッチが 4 つあります）については、次のようにパラメータ化する必要があります。
奇妙に見えますが、これは JMP がモデルをパラメータ化する方法によるものです。
https://www.jmp.com/support/help/en/18.0/index.shtml#page/jmp/nominal-factors.shtml#
また、テスト不可能なコントラストとは、すべてのコントラスト (列) を組み合わせたテストを提供する出力の下部にある平方和を指します。この場合、コントラストは直交していないため (非公式には、これは独立していると考えることができます)、意味がありません。ただし、3 つのコントラストのそれぞれを個別にテストすることは有効かつ正確です。
直交性についてはこちらをご覧ください。 https://stats.stackexchange.com/questions/171324/what-is-the-relationship-between-orthogonal-correlation-and-independence
Re: 2 つの予測値 (重回帰モデル経由) が有意に異なるかどうかをテストするこんにちは@broccoli切片を忘れました (最初の 2 列に 1 を入力してください)。
入学時の NIHSS が何であっても問題ないのでしょうか?
編集: うーん、共変量を含む相互作用項がある場合、期待どおりに機能するかどうかはわかりません。
乞うご期待。
Re: 2 つの予測値 (重回帰モデル経由) が有意に異なるかどうかをテストするこんにちは@MRB3855 、この目的のためにカスタム テストを設定しようとしましたが、うまくいかないようです。
たとえば、下の図では、入院時の NIHSS = 0で、 EVT タイミング = 0 (列 1) とEVT タイミング = 1 (列 2) を比較しようとしています。
テスト出力の最初の 2 列の「値」行にリストされている数値 ( 1.467と2.100 ) は、モデルがこれらの 2 つの条件で予測する値 (それぞれ0.0563と3.568 - 以下の予測プロファイラーを参照) と等しくありません。
同様に、 NIHSS on Admission = 20でこのカスタム テストを繰り返そうとすると (以下)、出力される値の最初の 2 つの列 ( 2.329と2.962 ) が、モデルの予測値 (それぞれ5.105と5.089 ) と一致しなくなります。
ここで何が欠けているのでしょうか? 明らかに、これは重要なことです。なぜなら、テストの「結果」は、私が正しく理解していれば、出力された 3 番目の列(値予測、標準誤差、 t 比、 p 値、およびSS の差) にリストされるはずであり、入学時に NIHSS をどのように設定したかに関係なく、まったく同じであるように見えるからです (これは絶対にそうではありません)。
ご迷惑をおかけして申し訳ございませんが、改めて感謝申し上げます。
追記ですが、「警告: テスト不可能なコントラスト」とはどういう意味ですか?
Re: 2 つの予測値 (重回帰モデル経由) が有意に異なるかどうかをテストするこんにちは@broccoli : はい、可能です。出力の上部にある赤い三角形をクリックし、見積もり/カスタム テストを選択します (下の図を参照)。例を示します。
カスタム テストの例 (jmp.com)
ただし、これは難しいかもしれません... 一度見て、助けが必要な場合は戻ってきてください。 2 つの列を追加する必要があります。次に、関心のある EVT タイミング (および相互作用) に 1 を入力し、最初の 2 つの列に、入院時の NIHSS の関心のある値を選択します。3 番目の列は、最初の 2 つの行の差です。
以下は、時間 = 3.14 でのバッチ 1_11 (最初の列 101.2) と時間 = 3.14 でのバッチ 3_13 (2 番目の列 101.4) の予測を比較した例です。3 番目の列は予測の差 (-0.202、p 値 0.6758) です。
