=== 順序データの層別群間比較 ===                                   [DANS V7.4]

データ名:tck表4.4.6 薬剤と改善度の性別2x2分割表 (mtest_tck040406.tab 2024-10-01)

層別内容:男

群項目(縦)  :薬剤
順序項目(横):改善度判定

薬剤    非改善(  %  )  改善(  %  )  合計(  %  )
-----------------------------------------------
標準薬      25( 50.0)    25( 50.0)    50(100.0)
新薬         5( 10.0)    45( 90.0)    50(100.0)
-----------------------------------------------
合計        30( 30.0)    70( 70.0)   100(100.0)

・Mantel-Haenszel検定
  順位相関比:η=0.436436  寄与率:η2乗=0.190476
  順位相関比の検定:χ2乗=18.8571  自由度=1  p=1.40883e-05***
  リスク差(risk difference):RD=0.4
    95%信頼区間:下限=0.238377  上限=0.561623
  リスク比(risk ratio):RR=1.8  ln(RR)=0.587787
    95%信頼区間:下限=1.34533  上限=2.40832
  オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=3.16519  上限=14.8348
  共通性の検定(修正有):χ2乗=17.0186  自由度=1  p=3.7016e-05***

・ロジット検定
  オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=3.06406  上限=26.4356
  オッズ比の検定:χ2乗=15.9743  自由度=1  p=6.42074e-05***


=== 順序データの層別群間比較 ===                                   [DANS V7.4]

データ名:tck表4.4.6 薬剤と改善度の性別2x2分割表 (mtest_tck040406.tab 2024-10-01)

層別内容:女

群項目(縦)  :薬剤
順序項目(横):改善度判定

薬剤    非改善(  %  )  改善(  %  )  合計(  %  )
-----------------------------------------------
標準薬      45( 90.0)     5( 10.0)    50(100.0)
新薬        25( 50.0)    25( 50.0)    50(100.0)
-----------------------------------------------
合計        70( 70.0)    30( 30.0)   100(100.0)

・Mantel-Haenszel検定
  順位相関比:η=0.436436  寄与率:η2乗=0.190476
  順位相関比の検定:χ2乗=18.8571  自由度=1  p=1.40883e-05***
  リスク差(risk difference):RD=0.4
    95%信頼区間:下限=0.238377  上限=0.561623
  リスク比(risk ratio):RR=5  ln(RR)=1.60944
    95%信頼区間:下限=2.1357  上限=11.7058
  オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=2.1503  上限=15.8497
  共通性の検定(修正有):χ2乗=17.0186  自由度=1  p=3.7016e-05***

・ロジット検定
  オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=3.06406  上限=26.4356
  オッズ比の検定:χ2乗=15.9743  自由度=1  p=6.42074e-05***


=== 順序データの層別群間比較 ===                                   [DANS V7.4]

データ名:tck表4.4.6 薬剤と改善度の性別2x2分割表 (mtest_tck040406.tab 2024-10-01)

層別内容:全度数表合計

群項目(縦)  :薬剤-合計
順序項目(横):改善度判定-合計

薬剤    非改善(  %  )  改善(  %  )  合計(  %  )
-----------------------------------------------
標準薬      70( 70.0)    30( 30.0)   100(100.0)
新薬        30( 30.0)    70( 70.0)   100(100.0)
-----------------------------------------------
合計       100( 50.0)   100( 50.0)   200(100.0)

・Mantel-Haenszel検定
  順位相関比:η=0.4  寄与率:η2乗=0.16
  順位相関比の検定:χ2乗=31.84  自由度=1  p=1.67409e-08***
  リスク差(risk difference):RD=0.4
    95%信頼区間:下限=0.27298  上限=0.52702
  リスク比(risk ratio):RR=2.33333  ln(RR)=0.847298
    95%信頼区間:下限=1.68707  上限=3.22715
  オッズ比(odds ratio):OR=5.44444  ln(OR)=1.6946
    95%信頼区間:下限=2.64933  上限=8.23956
  共通性の検定(修正有):χ2乗=30.2679  自由度=1  p=3.76302e-08***

・ロジット検定
  オッズ比(odds ratio):OR=5.44444  ln(OR)=1.6946
    95%信頼区間:下限=2.97349  上限=9.96874
  オッズ比の検定:χ2乗=30.1524  自由度=1  p=3.99398e-08***

○全体
・Mantel-Haenszel検定
  順位相関比:η=0.435338  寄与率:η2乗=0.189519
  順位相関比の検定:χ2乗=37.7143  自由度=1  p=8.19026e-10***
  調整リスク差(risk difference):RD=0.4
    95%信頼区間:下限=0.285715  上限=0.514285
  調整リスク比(risk ratio):RR=2.33333  ln(RR)=0.847298
    95%信頼区間:下限=1.72635  上限=3.15373
  調整オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=3.59166  上限=14.4083
  共通性の検定(修正有):χ2乗=35.8521  自由度=1  p=2.12874e-09***
  異質性の検定:χ2乗=0  自由度=1  p=1
  Breslow-Dayの検定(異質性):χ2乗=0  自由度=1  p=1

・ロジット検定
  調整オッズ比(odds ratio):OR=9  ln(OR)=2.19722
    95%信頼区間:下限=4.20102  上限=19.281
  調整オッズ比の検定:χ2乗=31.9486  自由度=1  p=1.58302e-08***
  異質性の検定:χ2乗=0  自由度=1  p=1


=== 共分散分析(analysis of covariance, ANCOVA) ===                 [DANS V7.4]

データ名:表4.4.61 性で層別した薬剤と改善度の分割表 (mtest_tck040406.dns 2024-10-01)

群項目    :薬剤 (0:標準薬 1:新薬)
集計項目y :改善 (0:非改善 1:改善)
共変数x  1:性別 (0:男 1:女)

・群    1:薬剤 (0:標準薬 1:新薬)=0
-----------------------------------------------------------------------------
y        :例数=100       平均値=0.3          SD=0.460566     SE=0.0460566    
x1       :例数=100       平均値=0.5          SD=0.502519     SE=0.0502519    
-----------------------------------------------------------------------------
群別回帰式:y=          0.5      -0.4x1 
共通回帰式:y=          0.5      -0.4x1 
群別回帰式の寄与率:r2乗=0.190476  r=-0.436436  有意確率:p=5.64554e-06***

・群    2:薬剤 (0:標準薬 1:新薬)=1
-----------------------------------------------------------------------------
y        :例数=100       平均値=0.7          SD=0.460566     SE=0.0460566    
x1       :例数=100       平均値=0.5          SD=0.502519     SE=0.0502519    
-----------------------------------------------------------------------------
群別回帰式:y=          0.9      -0.4x1 
共通回帰式:y=          0.9      -0.4x1 
群別回帰式の寄与率:r2乗=0.190476  r=-0.436436  有意確率:p=5.64554e-06***

・全体
-----------------------------------------------------------------------------
y        :例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
x1       :例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
-----------------------------------------------------------------------------
群別回帰式:y=          0.7      -0.4x1 
共通回帰式:y=          0.7      -0.4x1 
群別回帰式の寄与率:r2乗=0.16  r=-0.4  有意確率:p=4.40461e-09***

                          共分散分析表(ANCOVA table)
要因              平方和    自由度    平均平方和         F値   有意確率p値
------------------------------------------------------------------------------
群差                   8         1             8     46.1176   1.29709e-10***
共通回帰               8         1             8     46.1176   1.29709e-10***
------------------------------------------------------------------------------
修正群差               8         1             8     46.1176   1.29709e-10***
全体回帰               8         1             8     46.1176   1.29709e-10***
------------------------------------------------------------------------------
非平行性               0         1             0           0             1
残差                  34       196      0.173469
------------------------------------------------------------------------------
全体                  50       199

・修正群差のTukey型多重比較:共変数=全体平均値  群数=2  自由度=196
群 :  修正平均値 - 群 :  修正平均値    修正群差         q値   有意確率p値
-----------------------------------------------------------------------------
2  :         0.7 - 1  :         0.3         0.4     9.60392   1.29709e-10***
-----------------------------------------------------------------------------

・修正群差のTukey型95%同時信頼区間(simultaneous confidence interval)
群  - 群     修正群差    絶対精度        下限        上限
---------------------------------------------------------
2   - 1           0.4    0.116162    0.283838    0.516162
---------------------------------------------------------

                         共分散分析表(非平行性プール)
要因              平方和    自由度    平均平方和         F値   有意確率p値
------------------------------------------------------------------------------
群差                   8         1             8     46.3529   1.16419e-10***
共通回帰               8         1             8     46.3529   1.16419e-10***
------------------------------------------------------------------------------
修正群差               8         1             8     46.3529   1.16419e-10***
全体回帰               8         1             8     46.3529   1.16419e-10***
------------------------------------------------------------------------------
残差                  34       197      0.172589
------------------------------------------------------------------------------
全体                  50       199

・修正群差のTukey型多重比較:共変数=全体平均値  群数=2  自由度=197
群 :  修正平均値 - 群 :  修正平均値    修正群差         q値   有意確率p値
-----------------------------------------------------------------------------
2  :         0.7 - 1  :         0.3         0.4     9.62839   1.16419e-10***
-----------------------------------------------------------------------------

・修正群差のTukey型95%同時信頼区間(simultaneous confidence interval)
群  - 群     修正群差    絶対精度        下限        上限
---------------------------------------------------------
2   - 1           0.4    0.115863    0.284137    0.515863
---------------------------------------------------------


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

データ名:表4.4.61 性で層別した薬剤と改善度の分割表 (mtest_tck040406.dns 2024-10-01)

目的変数y   :改善 (0:非改善 1:改善)
説明変数x  1:薬剤 (0:標準薬 1:新薬)
説明変数x  2:性別 (0:男 1:女)

・各変数の基礎統計量
------------------------------------------------------------------------
y   :例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
x  1:例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
x  2:例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
------------------------------------------------------------------------
反応有:コード=1 例数=100  反応無:コード=1以外 例数=100

・相関行列(correlation coefficient matrix)
        x  1    x  2    y  1
----------------------------
x  1    1      0.000   0.400
x  2   0.000    1     -0.400
y  1   0.400  -0.400    1   
----------------------------

・前進的変数増減法(stepwise forward selection method)による変数選択結果
  取り込み基準:χ2乗値＞2  追い出し基準:χ2乗値≦2
  反復計算回数=6
  ロジットモデル:p=1/{1+exp(-β0-Σβj・xj)}
  p:y=1(反応有)の確率  β0:切片  βj:変数xjの偏回帰係数
                                            標準                   有意確率
変数  偏回帰係数   標準誤差    オッズ比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片-1.65219e-14   0.263461                         3.93268e-27           1
x  1     2.19722    0.38873           9     1.10137     31.9486 1.58302e-08***
x  2    -2.19722    0.38873    0.111111    -1.10137     31.9486 1.58302e-08***
------------------------------------------------------------------------------

変数  偏回帰係数   95%CI下限        上限    オッズ比   95%CI下限        上限
----------------------------------------------------------------------------
切片-1.65219e-14   -0.516374    0.516374
x  1     2.19722     1.43533     2.95912           9     4.20102      19.281
x  2    -2.19722    -2.95912    -1.43533    0.111111   0.0518645    0.238037
----------------------------------------------------------------------------
対数尤度:L(モデル)=-101.823  L(切片)=-138.629  飽和Lf=-101.823
擬似寄与率R2乗=1  AIC(赤池の情報量基準)=209.646

                          回帰とズレの検定
要因         -(対数尤度比)    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰               36.8064         2       73.6128   1.04083e-16***
ズレ(LOF)       6.6791e-13         1             0             1
--------------------------------------------------------------------
全体               36.8064         3


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

データ名:表4.4.61 性で層別した薬剤と改善度の分割表 (mtest_tck040406.dns 2024-10-01)

目的変数y   :改善 (0:非改善 1:改善)
説明変数x  1:薬剤 (0:標準薬 1:新薬)

・各変数の基礎統計量
------------------------------------------------------------------------
y   :例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
x  1:例数=200       平均値=0.5          SD=0.501255     SE=0.0354441    
------------------------------------------------------------------------
反応有:コード=1 例数=100  反応無:コード=1以外 例数=100

・相関行列(correlation coefficient matrix)
        x  1    y  1
--------------------
x  1    1      0.400
y  1   0.400    1   
--------------------

・前進的変数増減法(stepwise forward selection method)による変数選択結果
  取り込み基準:χ2乗値＞2  追い出し基準:χ2乗値≦2
  反復計算回数=4
  ロジットモデル:p=1/{1+exp(-β0-Σβj・xj)}
  p:y=1(反応有)の確率  β0:切片  βj:変数xjの偏回帰係数
                                            標準                   有意確率
変数  偏回帰係数   標準誤差    オッズ比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片   -0.847298   0.218218                             15.0762 0.000103258***
x  1      1.6946   0.308607     5.44444    0.849424     30.1524 3.99398e-08***
------------------------------------------------------------------------------

変数  偏回帰係数   95%CI下限        上限    オッズ比   95%CI下限        上限
----------------------------------------------------------------------------
切片   -0.847298      -1.275   -0.419599
x  1      1.6946     1.08974     2.29945     5.44444     2.97349     9.96873
----------------------------------------------------------------------------
対数尤度:L(モデル)=-122.173  L(切片)=-138.629  飽和Lf=-122.173
擬似寄与率R2乗=1  AIC(赤池の情報量基準)=248.346

                          回帰とズレの検定
要因         -(対数尤度比)    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰               16.4566         1       32.9132   9.63691e-09***
ズレ(LOF)      5.68434e-14         0             0             1
--------------------------------------------------------------------
全体               16.4566         1


