=== パラメトリック生存時間解析 ===                                 [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


=== パラメトリック生存時間解析 ===                                 [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   
----------------------------

・前進的変数増減法による変数選択結果  反復回数(微分法)=9
  取り込み基準:χ2乗＞2  追い出し基準:χ2乗≦2
  指数分布モデル:S(t)=exp(-λt)=exp{-exp(β0+Σβjxj)t}
  S(t):生存関数  λ:ハザード  β0:切片  βj:共変数xjの偏回帰係数
  初期値:β0=0 β1=0 β2=0 
                                            標準                   有意確率
共変数 偏回帰係数  標準誤差  ハザード比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片    -3.63879   0.254483                             204.454           0***
x  1   -0.666585   0.271464    0.513459   -0.335562     6.02957   0.0140682*
x  2    0.708172   0.172818     2.03028    0.558954     16.7919 4.17115e-05***
------------------------------------------------------------------------------

共変数  偏回帰係数   95%CI下限        上限  ハザード比   95%CI下限        上限
------------------------------------------------------------------------------
切片      -3.63879    -4.13757    -3.14001
x  1     -0.666585    -1.19864   -0.134525    0.513459    0.301603    0.874131
x  2      0.708172    0.369455     1.04689     2.03028     1.44695     2.84878
------------------------------------------------------------------------------

対数尤度:L(モデル)=-231.78  L(切片)=-241.776  飽和Lf=-231.183
擬似寄与率:R2乗=0.943587  AIC=469.561
C-index=1505/2218=0.678539  1-C-index=0.321461

                          回帰とズレの検定
要因         -(対数尤度比)    自由度          χ2乗   有意確率p値
--------------------------------------------------------------------
回帰               9.99547         2       19.9909   4.56062e-05***
ズレ(LOF)         0.597586         3       1.19517      0.754163
--------------------------------------------------------------------
全体               10.5931         5

・全変数を選択した結果  反復回数(微分法)=4
  指数分布モデル:S(t)=exp(-λt)=exp{-exp(β0+Σβjxj)t}
  S(t):生存関数  λ:ハザード  β0:切片  βj:共変数xjの偏回帰係数
  初期値:β0=-3.55247 β1=-0.77911 β2=0.60483 
                                            標準                   有意確率
共変数 偏回帰係数  標準誤差  ハザード比  偏回帰係数  Waldのχ2乗       p値
------------------------------------------------------------------------------
切片    -3.63879   0.254483                             204.454           0***
x  1   -0.666585   0.271464    0.513459   -0.335562     6.02957   0.0140682*
x  2    0.708172   0.172818     2.03028    0.558954     16.7919 4.17115e-05***
------------------------------------------------------------------------------

共変数  偏回帰係数   95%CI下限        上限  ハザード比   95%CI下限        上限
------------------------------------------------------------------------------
切片      -3.63879    -4.13757    -3.14001
x  1     -0.666585    -1.19864   -0.134525    0.513459    0.301603    0.874131
x  2      0.708172    0.369455     1.04689     2.03028     1.44695     2.84878
------------------------------------------------------------------------------

対数尤度:L(モデル)=-231.78  L(切片)=-241.776  飽和Lf=-231.183
擬似寄与率:R2乗=0.943587  AIC=469.561
C-index=1505/2218=0.678539  1-C-index=0.321461

                          回帰とズレの検定
要因         -(対数尤度比)    自由度          χ2乗   有意確率p値
--------------------------------------------------------------------
回帰               9.99547         2       19.9909   4.56062e-05***
ズレ(LOF)         0.597586         3       1.19517      0.754163
--------------------------------------------------------------------
全体               10.5931         5


