=== Coxの比例ハザードモデルによる生存時間解析 ===                  [DANS V7.4]

データ名:Coxの比例ハザードモデル用テストデータ (cox_00.dns 2005-5-23)

期間項目  :観察期間 (月)
転帰項目  :転帰 (0:死亡 1:生存)
共変数x  1:群 (0:A 1:B)
共変数x  2:重症度 (0:症状なし 1:軽症 2:重症)

・Kaplan-Meier法による累積生存率    死亡コード:0  +:打ち切り (+):脱落
番号  ケース番号  生存期間(転帰)  生存例数/観測例数  累積生存率    標準誤差
---------------------------------------------------------------------------
   1       1             1              69/70          0.985714   0.0141833
   2       2             2              68/69          0.971429   0.0199123
   3       3             2              67/68          0.957143   0.0242076
   4       4             3              66/67          0.942857   0.0277431
   5       5             3              65/66          0.928571   0.0307818
   6       6             3              64/65          0.914286   0.0334594
   7      37             3              63/64               0.9   0.0358569
   8       7             4              62/63          0.885714   0.0380272
   9       8             4              61/62          0.871429   0.0400073
  10       9             4              60/61          0.857143   0.0418243
  11      38             4              59/60          0.842857   0.0434986
  12      10             5              58/59          0.828571   0.0450461
  13      11             5              57/58          0.814286   0.0464796
  14      12             5              56/57               0.8   0.0478091
  15      13             5              55/56          0.785714   0.0490433
  16      39             5              54/55          0.771429   0.0501891
  17      40             5              53/54          0.757143   0.0512525
  18      14             6              52/53          0.742857   0.0522385
  19      41             7              51/52          0.728571   0.0531514
  20      15             8              50/51          0.714286   0.0539949
  21      16             8              49/50               0.7   0.0547723
  22      17             9              48/49          0.685714   0.0554862
  23      42             9              47/48          0.671429   0.0561391
  24      43            10              46/47          0.657143   0.0567332
  25      44            10              45/46          0.642857   0.0572703
  26      45            11              44/45          0.628571   0.0577519
  27      18            12              43/44          0.614286   0.0581794
  28      19            12              42/43               0.6    0.058554
  29      20            12              41/42          0.585714   0.0588768
  30      21            12              40/41          0.571429   0.0591485
  31      22            13              39/40          0.557143   0.0593699
  32      46            13              38/39          0.542857   0.0595415
  33      47            14              37/38          0.528571   0.0596638
  34      23            16              36/37          0.514286    0.059737
  35      48            18              35/36               0.5   0.0597614
  36      49            18              34/35          0.485714    0.059737
  37      50            19              33/34          0.471429   0.0596638
  38      51            19              32/33          0.457143   0.0595415
  39      52            21              31/32          0.442857   0.0593699
  40      53            23              30/31          0.428571   0.0591485
  41      54            25              29/30          0.414286   0.0588768
  42      55            26 +        (   29/29   )      0.414286   0.0588768
  43      24            27              27/28           0.39949   0.0586037
  44      56            27              26/27          0.384694    0.058271
  45      25            28              25/26          0.369898   0.0578778
  46      26            28              24/25          0.355102   0.0574228
  47      57            28              23/24          0.340306   0.0569045
  48      58            28 +        (   23/23   )      0.340306   0.0569045
  49      59            30              21/22          0.324838   0.0563811
  50      27            31              20/21          0.309369   0.0557779
  51      60            32              19/20          0.293901   0.0550921
  52      28            32 +        (   19/19   )      0.293901   0.0550921
  53      29            33              17/18          0.277573   0.0543972
  54      61            33 +        (   17/17   )      0.277573   0.0543972
  55      30            34              15/16          0.260225   0.0536925
  56      31            35 +        (   15/15   )      0.260225   0.0536925
  57      62            35 +        (   14/14   )      0.260225   0.0536925
  58      32            36 +        (   13/13   )      0.260225   0.0536925
  59      63            37              11/12          0.238539   0.0534181
  60      33            44 +        (   11/11   )      0.238539   0.0534181
  61      64            49               9/10          0.214685   0.0531361
  62      65            52 +        (    9/9    )      0.214685   0.0531361
  63      66            54               7/8            0.18785   0.0528378
  64      34            54 +        (    7/7    )       0.18785   0.0528378
  65      35            55               5/6           0.156541   0.0524939
  66      67            56               4/5           0.125233   0.0504753
  67      36            56 +        (    4/4    )      0.125233   0.0504753
  68      68            58 +        (    3/3    )      0.125233   0.0504753
  69      69            59 +        (    2/2    )      0.125233   0.0504753
  70      70            60 +        (    1/1    )      0.125233   0.0504753
---------------------------------------------------------------------------
  50%生存時間推定値=18    95%信頼区間=12〜27


=== Coxの比例ハザードモデルによる生存時間解析 ===                  [DANS V7.4]

データ名:Coxの比例ハザードモデル用テストデータ (cox_00.dns 2005-5-23)

期間項目  :観察期間 (月)
転帰項目  :転帰 (0:死亡 1:生存)
共変数x  1:群 (0:A 1:B)
共変数x  2:重症度 (0:症状なし 1:軽症 2:重症)

・共変数の基礎統計量:y=ln(観察期間 (月))
------------------------------------------------------------------------
y   :例数=70        平均値=2.67845      SD=1.02104      SE=0.122038     
x  1:例数=70        平均値=0.485714     SD=0.503405     SE=0.0601684    
x  2:例数=70        平均値=1.01429      SD=0.789292     SE=0.0943384    
------------------------------------------------------------------------
死亡コード=0  死亡数=56  打ち切り数=14  平均瞬間死亡率=0.036246

・相関行列:y1=ln(観察期間 (月))
        x  1    x  2    y  1
----------------------------
x  1    1      0.165   0.307
x  2   0.165    1     -0.404
y  1   0.307  -0.404    1   
----------------------------

・前進的変数増減法による変数選択結果  反復回数=5
  取り込み基準:χ2乗＞2  追い出し基準:χ2乗≦2
  Coxの比例ハザードモデル(有限修正有):S(t)=S0(t)^HR=S0(t)^exp(β0+Σβjxj)
  S(t):補正生存関数  HR:ハザード比  β0:切片  βj:共変数xjの偏回帰係数
  S0(t):基準生存関数(全共変数の値が平均値の時の生存関数)
  初期値:β1=0 β2=0 
                                            標準                   有意確率
共変数 偏回帰係数  標準誤差  ハザード比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片   -0.420242
x  1   -0.669109   0.279439    0.512164   -0.336833     5.73352   0.0166441*
x  2    0.734742   0.185005     2.08494    0.579926     15.7725 7.14326e-05***
------------------------------------------------------------------------------

共変数  偏回帰係数   95%CI下限        上限  ハザード比   95%CI下限        上限
------------------------------------------------------------------------------
x  1     -0.669109     -1.2168    -0.12142    0.512164    0.296177    0.885662
x  2      0.734742    0.372138     1.09735     2.08494     1.45083      2.9962
------------------------------------------------------------------------------

対数尤度:L(モデル)=-191.79  L(切片)=-201.434  AIC=387.58
回帰の尤度比検定:χ2乗=19.2873  自由度=2  有意確率p=6.48375e-05***
C-index=1505/2218=0.678539  1-C-index=0.321461

・全変数を選択した結果  反復回数=4
  Coxの比例ハザードモデル(有限修正有):S(t)=S0(t)^HR=S0(t)^exp(β0+Σβjxj)
  S(t):補正生存関数  HR:ハザード比  β0:切片  βj:共変数xjの偏回帰係数
  S0(t):基準生存関数(全共変数の値が平均値の時の生存関数)
  初期値:β1=-0.77911 β2=0.60483 
                                            標準                   有意確率
共変数 偏回帰係数  標準誤差  ハザード比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片   -0.420249
x  1   -0.669106   0.279439    0.512166   -0.336831     5.73346   0.0166446*
x  2    0.734747   0.185005     2.08495     0.57993     15.7727 7.14246e-05***
------------------------------------------------------------------------------

共変数  偏回帰係数   95%CI下限        上限  ハザード比   95%CI下限        上限
------------------------------------------------------------------------------
x  1     -0.669106     -1.2168   -0.121416    0.512166    0.296178    0.885665
x  2      0.734747    0.372143     1.09735     2.08495     1.45084     2.99622
------------------------------------------------------------------------------

対数尤度:L(モデル)=-191.79  L(切片)=-201.434  AIC=387.58
回帰の尤度比検定:χ2乗=19.2873  自由度=2  有意確率p=6.48375e-05***
C-index=1505/2218=0.678539  1-C-index=0.321461


