公開日: 11/25/2021

名義ロジスティック回帰の例

インゴットの製造に関する実験として、加熱時間とソーキング時間を変え、仕様に適合した、圧延可能な状態に仕上がったかどうかを調べてみました。Cox and Snell(1989)を参照してください。データは、サンプルデータのフォルダにある「Ingots.jmp」データテーブルにまとめられています。この例では、「モデルのあてはめ」プラットフォームを使い、「状態」の確率に「加熱時間」「ソーキング時間」を説明変数としたロジスティック回帰モデルをあてはめます。

1. [ヘルプ]>[サンプルデータライブラリ]を選択し、「Ingots.jmp」を開きます。

カテゴリカル変数である「状態」は、「Ready」(適合)と「Not Ready」(不適合)の値を取ります。

2. [分析]>[モデルのあてはめ]を選択します。

3. 「状態」を選択し、[Y]をクリックします。

選択した列が名義尺度であるため、自動的に[名義ロジスティック]が手法として指定されます。

「状態」は、2水準の名義尺度の列であるため、[イベントを示す水準]オプションが表示されます。このオプションで、どの応答水準の確率をモデル化するかが指定できます。

4. 「イベントを示す水準」として[Ready]を選択します。

[イベントを示す水準]を「Ready」とすると、「Ready」が生じる確率のほうがロジットの分子に設定されます。

5. 「加熱時間」「ソーキング時間」を選択し、[追加]をクリックします。

6. 「度数」を選択し、[度数]をクリックします。

7. [実行]をクリックします。

反復計算が収束したら、名義ロジスティックの回帰レポートが表示されます。

図11.2 「名義ロジスティックのあてはめ」レポート 

Nominal Logistic Fit Report

「モデル全体の検定」レポートを見ると、カイ2乗統計量(11.64)のp値が0.0030と小さく、モデル全体が有意であることがわかります。しかし、「ソーキング時間」は、パラメータ推定値のp値が0.8639であり、統計的に有意ではないことが示されています。

8. 「名義ロジスティックのあてはめ 状態」の赤い三角ボタンをクリックし、[プロファイル]を選択します。

図11.3 予測プロファイル 

Prediction Profiler

「加熱時間」を45、「ソーキング時間」を2に設定したとき、「状態」が「Ready」である確率は0.85となります。

この例に対する追加の分析として、圧延可能な状態となる確率の値を指定した場合の、特定のソーキング時間に対応する加熱時間を、逆推定によって求めたくなるかもしれません。そのような逆推定の例については、逆推定の例を参照してください。

より詳細な情報が必要な場合や、質問があるときは、JMPユーザーコミュニティで答えを見つけましょう (community.jmp.com).