=== ロジスティック回帰分析(logistic regression analysis) ===       [DANS V7.4]

データ名:判別分析用テストデータ (disan_00.dns 2005-4-19)

目的変数y   :群 (0:正常 1:動脈硬化症)
説明変数x  1:TC (mg/dl)
説明変数x  2:TG (mg/dl)

・各変数の基礎統計量
------------------------------------------------------------------------
y   :例数=25        平均値=0.4          SD=0.5          SE=0.1          
x  1:例数=25        平均値=224.4        SD=28.2961      SE=5.65921      
x  2:例数=25        平均値=207.2        SD=60.1747      SE=12.0349      
------------------------------------------------------------------------
反応有:コード=1 例数=10  反応無:コード=1以外 例数=15

・相関行列(correlation matrix)
        x  1    x  2    y  1
----------------------------
x  1    1      0.499   0.783
x  2   0.499    1      0.025
y  1   0.783   0.025    1   
----------------------------

・前進的変数増減法(stepwise forward selection method)による変数選択結果
  取り込み基準:χ2乗値＞2  追い出し基準:χ2乗値≦2
  反復計算回数=18
  ロジットモデル:p=1/{1+exp(-β0-Σβj・xj)}
  p:y=1(反応有)の確率  β0:切片  βj:変数xjの偏回帰係数
  偏回帰係数初期値:β0=0 β1=0 β2=0
                                            標準                   有意確率
変数  偏回帰係数   標準誤差    オッズ比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片    -33.3481    13.9605                             5.70611   0.0169059*
x  1    0.144287  0.0608468     1.15522     4.08276     5.62315   0.0177248*
------------------------------------------------------------------------------

変数  偏回帰係数   95%CI下限        上限    オッズ比   95%CI下限        上限
----------------------------------------------------------------------------
切片    -33.3481    -60.7101    -5.98602
x  1    0.144287   0.0250295    0.263545     1.15522     1.02535     1.30154
----------------------------------------------------------------------------
対数尤度:L(モデル)=-6.0397  L(切片)=-16.8253  飽和Lf=-5.54518
擬似寄与率R2乗=0.95616  AIC(赤池の情報量基準)=16.0794

                          回帰とズレの検定
要因         -(対数尤度比)    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰               10.7856         1       21.5712   3.40935e-06***
ズレ(LOF)          0.49452         9      0.989039      0.999462
--------------------------------------------------------------------
全体               11.2801        10

・全変数を選択した結果  反復計算回数=14
  ロジットモデル:p=1/{1+exp(-β0-Σβj・xj)}
  p:y=1(反応有)の確率  β0:切片  βj:変数xjの偏回帰係数
  偏回帰係数初期値:β0=-61.569 β1=0.335192 β2=-0.074889
                                            標準                   有意確率
変数  偏回帰係数   標準誤差    オッズ比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片    -568.504    25515.6                         0.000496428    0.982224
x  1     2.86982    133.043     17.6338     81.2046 0.000465293     0.98279
x  2     -0.4266    28.2213    0.652725    -25.6705 0.000228501    0.987939
------------------------------------------------------------------------------

変数  偏回帰係数   95%CI下限        上限    オッズ比   95%CI下限        上限
----------------------------------------------------------------------------
切片    -568.504    -50578.1     49441.1
x  1     2.86982    -257.889     263.629     17.63381.00031e-1123.10855e+114
x  2     -0.4266    -55.7393     54.8861    0.652725 6.20513e-25 6.86609e+23
----------------------------------------------------------------------------
対数尤度:L(モデル)=-2.61293e-06  L(切片)=-16.8253  飽和Lf=0
擬似寄与率R2乗=1  AIC(赤池の情報量基準)=6.00001

                          回帰とズレの検定
要因         -(対数尤度比)    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰               16.8253         2       33.6506   4.93026e-08***
ズレ(LOF)      2.61293e-06        22   5.22586e-06             1
--------------------------------------------------------------------
全体               16.8253        24


