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

データ名:tck表5.3.6 マウスの生後日数分類と尾長分類 (logreg_tck050306.dns 2024-10-09)

目的変数y   :尾長分類 (0:8cm未満 1:8cm以上)
説明変数x  1:生後日数分類

・各変数の基礎統計量
------------------------------------------------------------------------
y   :例数=11        平均値=0.454545     SD=0.522233     SE=0.157459     
x  1:例数=11        平均値=9            SD=3.07409      SE=0.926872     
------------------------------------------------------------------------
反応有:コード=1 例数=5  反応無:コード=0 例数=6

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

・前進的変数増減法(stepwise forward selection method)による変数選択結果
  取り込み基準:χ2乗値＞2  追い出し基準:χ2乗値≦2
  反復計算回数=5
  ロジットモデル:p=1/{1+exp(-β0-Σβj・xj)}
  p:y=1(反応有)の確率  β0:切片  βj:変数xjの偏回帰係数
                                            標準                   有意確率
変数  偏回帰係数   標準誤差    オッズ比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片    -4.25831    2.68974                             2.50642    0.113383
x  1    0.445178   0.280664     1.56077     1.36852     2.51591    0.112703
------------------------------------------------------------------------------

変数  偏回帰係数   95%CI下限        上限    オッズ比   95%CI下限        上限
----------------------------------------------------------------------------
切片    -4.25831    -9.53011     1.01348
x  1    0.445178   -0.104913     0.99527     1.56077    0.900403     2.70545
----------------------------------------------------------------------------
対数尤度:L(モデル)=-5.8793  L(切片)=-7.5791  飽和Lf=-5.6144
擬似寄与率R2乗=0.865169  AIC(赤池の情報量基準)=15.7586

                          回帰とズレの検定
要因         -(対数尤度比)    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰                1.6998         1        3.3996     0.0652122+
ズレ(LOF)         0.264902         1      0.529805      0.466689
--------------------------------------------------------------------
全体                1.9647         2


=== 判別分析(discriminant analysis) ===                            [DANS V7.4]

データ名:tck表5.3.6 マウスの生後日数分類と尾長分類 (logreg_tck050306.dns 2024-10-09)

判別項目:尾長分類 (0:8cm未満 1:8cm以上)
変数x  1:生後日数分類

○群1:尾長分類 (0:8cm未満 1:8cm以上)=1
・各変数の基礎統計量
------------------------------------------------------------------------
x  1:例数=5         平均値=10.7         SD=2.46475      SE=1.10227      
------------------------------------------------------------------------

○群2:尾長分類 (0:8cm未満 1:8cm以上)=0
・各変数の基礎統計量
------------------------------------------------------------------------
x  1:例数=6         平均値=7.58333      SD=2.95663      SE=1.20704      
------------------------------------------------------------------------

○全体(SD・SE:群内)
・各変数の基礎統計量
------------------------------------------------------------------------
x  1:例数=11        平均値=9            SD=2.74891      SE=0.828826     
------------------------------------------------------------------------

・群内相関行列
        x  1
------------
x  1    1   
------------

・全変数を選択した結果
                                             標準                  有意確率
変数    判別係数   標準誤差 exp(判別係数)  判別係数       偏F値       p値
------------------------------------------------------------------------------
切片    -3.77048
x  1    0.412449    0.22028     1.51051     1.13378     3.50582   0.0939409+
------------------------------------------------------------------------------

変数    判別係数   95%CI下限        上限 exp(判別係数) 95%CI下限        上限
----------------------------------------------------------------------------
x  1    0.412449  -0.0858595    0.910758     1.51051    0.917723     2.48621
----------------------------------------------------------------------------
・判別効率(Mahalanobisの平方距離):D2乗=1.28547  相関比:η=0.529467
  D/2(正規分布z値)=0.566892  判別確率:p=0.714606(誤判別確率:1-p=0.285394)
・HotellingのT2乗検定:T2乗=D2乗×5×6/11=3.50582
  F=3.50582  第1自由度=1  第2自由度=9  有意確率:p=0.0939409+


=== 出現率の傾向分析 ===                                           [DANS V7.4]

データ名:tck表5.3.6 マウスの生後日数分類と尾長分類 (logreg_tck050306.tab 2024-10-09)

説明変数項目(縦):生後日数分類
目的変数項目(横):尾長分類 (0:8cm未満 1:8cm以上)

縦＼横     0(  %  )     1(  %  )  合計(  %  )
---------------------------------------------
4.5        2(100.0)     0(  0.0)     2(100.0)
8          3( 60.0)     2( 40.0)     5(100.0)
12.5       1( 25.0)     3( 75.0)     4(100.0)
---------------------------------------------
合計       6( 54.5)     5( 45.5)    11(100.0)

・横のコード=1の出現率直線:p=-0.354978+0.0899471x
  回帰係数の95%信頼区間=0.0899471±0.100392(-0.0104451-0.190339)
  寄与率r2乗=0.983632  r=0.991782

                     Cochran-Armitageの傾向分析
要因                平方和    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
直線性             3.08369         1       3.08369     0.0790811+
異質性(ズレ)     0.0513139         1     0.0513139      0.820793
--------------------------------------------------------------------
全体                 3.135         2         3.135      0.208566

・横のコード=1の出現率曲線:p=1/{1+exp(3.07941-0.334242x)}
  p0=1/{1+exp(3.07941)}=0.0439646  (繰り返し補完回数=10)
  回帰係数の95%信頼区間=0.334242±0.487992(-0.15375-0.822234)
  オッズ比:OR=1.39688  95%信頼区間=0.857487-2.27558
  寄与率r2乗=1  r=1

                            ロジット分析
要因                平方和    自由度        χ2乗値   有意確率p値
--------------------------------------------------------------------
回帰               1.80215         1       1.80215      0.179453
異質性(ズレ)   4.29021e-11         1   4.29021e-11      0.999995
--------------------------------------------------------------------
全体               1.80215         2


