=== 生存時間解析(survival time analysis) ===                       [DANS V7.4]

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

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

○Kaplan-Meier法による累積生存率    死亡コード:0  +:打ち切り  (+):脱落
・群1:重症度 (0:症状なし 1:軽症 2:重症)=0
例数=21  死亡数=13  観察期間合計=680
番号  ケース番号  生存期間(転帰)  生存例数/観測例数  累積生存率    標準誤差
---------------------------------------------------------------------------
   1      10             5              20/21          0.952381   0.0464714
   2      11             5              19/20          0.904762   0.0640564
   3      18            12              18/19          0.857143   0.0763604
   4      19            12              17/18          0.809524   0.0856891
   5      22            13              16/17          0.761905   0.0929429
   6      47            14              15/16          0.714286   0.0985808
   7      48            18              14/15          0.666667    0.102869
   8      56            27              13/14          0.619048    0.105971
   9      25            28              12/13          0.571429     0.10799
  10      26            28              11/12           0.52381    0.108985
  11      28            32 +        (   11/11   )       0.52381    0.108985
  12      30            34               9/10          0.471429    0.109956
  13      31            35 +        (    9/9    )      0.471429    0.109956
  14      62            35 +        (    8/8    )      0.471429    0.109956
  15      33            44 +        (    7/7    )      0.471429    0.109956
  16      66            54               5/6           0.392857    0.116364
  17      34            54 +        (    5/5    )      0.392857    0.116364
  18      67            56               3/4           0.294643    0.121865
  19      36            56 +        (    3/3    )      0.294643    0.121865
  20      68            58 +        (    2/2    )      0.294643    0.121865
  21      70            60 +        (    1/1    )      0.294643    0.121865
---------------------------------------------------------------------------
境界内平均生存時間:RMST(t=0-60)=36.9167  95%信頼区間=27.7546〜46.0788
観測データによる50%生存時間=32.9091  95%信頼区間=14〜60以上
指数分布モデル:最尤法
  ハザード:λ=0.0191176  95%信頼区間=0.00872537〜0.0295099
  RMST(t=0-60)=35.6964  95%信頼区間=0〜78.0425
  平均生存時間:RMST(t=0-∞)=52.3077  95%信頼区間=23.8734〜80.742
  50%生存時間=36.2569  95%信頼区間=16.5478〜55.9661
Weibull分布モデル:反復回数(微分法)=4
  対数尤度=-64.0598  AIC=132.12
  ln(a)=0.220355  a=1.24652  ln(λ)=-3.90249  λ=0.0201916
  hv0(0)=0.1  vhv0(0)=0  hv0(1)=-3.95714  vhv0(1)=0

・群2:重症度 (0:症状なし 1:軽症 2:重症)=1
例数=27  死亡数=22  観察期間合計=563
番号  ケース番号  生存期間(転帰)  生存例数/観測例数  累積生存率    標準誤差
---------------------------------------------------------------------------
   1       2             2              26/27          0.962963   0.0363447
   2       4             3              25/26          0.925926    0.050401
   3       7             4              24/25          0.888889   0.0604812
   4       8             4              23/24          0.851852   0.0683673
   5       9             4              22/23          0.814815   0.0747568
   6      12             5              21/22          0.777778   0.0800091
   7      39             5              20/21          0.740741    0.084337
   8      14             6              19/20          0.703704   0.0878772
   9      41             7              18/19          0.666667   0.0907218
  10      15             8              17/18           0.62963   0.0929349
  11      17             9              16/17          0.592593   0.0945607
  12      20            12              15/16          0.555556   0.0956292
  13      46            13              14/15          0.518519    0.096159
  14      23            16              13/14          0.481481    0.096159
  15      49            18              12/13          0.444444   0.0956292
  16      52            21              11/12          0.407407   0.0945607
  17      54            25              10/11           0.37037   0.0929349
  18      57            28               9/10          0.333333   0.0907218
  19      58            28 +        (    9/9    )      0.333333   0.0907218
  20      59            30               7/8           0.291667   0.0884338
  21      27            31               6/7               0.25   0.0850517
  22      61            33 +        (    6/6    )          0.25   0.0850517
  23      32            36 +        (    5/5    )          0.25   0.0850517
  24      64            49               3/4             0.1875   0.0836582
  25      65            52 +        (    3/3    )        0.1875   0.0836582
  26      35            55               1/2            0.09375    0.078385
  27      69            59 +        (    1/1    )       0.09375    0.078385
---------------------------------------------------------------------------
境界内平均生存時間:RMST(t=0-59)=23.3287  95%信頼区間=15.6406〜31.0168
観測データによる50%生存時間=14.5  95%信頼区間=6〜30
指数分布モデル:最尤法
  ハザード:λ=0.0390764  95%信頼区間=0.0227477〜0.0554051
  RMST(t=0-59)=23.0393  95%信頼区間=0〜59.5647
  平均生存時間:RMST(t=0-∞)=25.5909  95%信頼区間=14.8973〜36.2845
  50%生存時間=17.7383  95%信頼区間=10.3261〜25.1505
Weibull分布モデル:反復回数(微分法)=4
  対数尤度=-93.3292  AIC=190.658
  ln(a)=0.000888445  a=1.00089  ln(λ)=-3.24235  λ=0.0390718
  hv0(0)=0.1  vhv0(0)=0  hv0(1)=-3.24224  vhv0(1)=0

・群3:重症度 (0:症状なし 1:軽症 2:重症)=2
例数=22  死亡数=21  観察期間合計=302
番号  ケース番号  生存期間(転帰)  生存例数/観測例数  累積生存率    標準誤差
---------------------------------------------------------------------------
   1       1             1              21/22          0.954545   0.0444095
   2       3             2              20/21          0.909091   0.0612909
   3       5             3              19/20          0.863636    0.073165
   4       6             3              18/19          0.818182   0.0822304
   5      37             3              17/18          0.772727   0.0893461
   6      38             4              16/17          0.727273   0.0949514
   7      13             5              15/16          0.681818   0.0993026
   8      40             5              14/15          0.636364    0.102559
   9      16             8              13/14          0.590909    0.104824
  10      42             9              12/13          0.545455    0.106159
  11      43            10              11/12               0.5      0.1066
  12      44            10              10/11          0.454545    0.106159
  13      45            11               9/10          0.409091    0.104824
  14      21            12               8/9           0.363636    0.102559
  15      50            19               7/8           0.318182   0.0993026
  16      51            19               6/7           0.272727   0.0949514
  17      53            23               5/6           0.227273   0.0893461
  18      55            26 +        (    5/5    )      0.227273   0.0893461
  19      24            27               3/4           0.170455   0.0831355
  20      60            32               2/3           0.113636   0.0722771
  21      29            33               1/2          0.0568182   0.0540384
  22      63            37               0/1                  0           0
---------------------------------------------------------------------------
境界内平均生存時間:RMST(t=0-37)=14.0114  95%信頼区間=9.19757〜18.8252
観測データによる50%生存時間=10  95%信頼区間=4〜19
指数分布モデル:最尤法
  ハザード:λ=0.0695364  95%信頼区間=0.0397958〜0.0992771
  RMST(t=0-37)=13.2834  95%信頼区間=0〜35.143
  平均生存時間:RMST(t=0-∞)=14.381  95%信頼区間=8.23023〜20.5317
  50%生存時間=9.96812  95%信頼区間=5.70476〜14.2315
Weibull分布モデル:反復回数(微分法)=4
  対数尤度=-76.6413  AIC=157.283
  ln(a)=0.149905  a=1.16172  ln(λ)=-2.70975  λ=0.0665532
  hv0(0)=0.1  vhv0(0)=0  hv0(1)=-2.6659  vhv0(1)=0

・ログランク検定:χ2乗=12.9511  自由度=2  有意確率:p=0.00154065**

・一般化拡張Mantel検定  :χ2乗=11.5617  自由度=2  有意確率:p=0.00308609**
  交互作用(異質性)の検定:χ2乗=65.1031  自由度=54  有意確率:p=0.143203

○重症度 (0:症状なし 1:軽症 2:重症)のBonferroni型多重比較
・群1:0 - 群2:1
  パラメトリックHR検定:χ2乗=4.17634  自由度=1  有意確率:p=0.122976
  β=0.714906  標準誤差=0.349825  ハザード比(群2/群1)=2.044
  HRの95%同時信頼区間=0.884646〜4.7227  β=-0.122568〜1.55238
・群1:0 - 群3:2
  パラメトリックHR検定:χ2乗=13.3874  自由度=1  有意確率:p=0.000759953***
  β=1.29124  標準誤差=0.352905  ハザード比(群3/群1)=3.63729
  HRの95%同時信頼区間=1.56266〜8.46624  β=0.446391〜2.13609
・群2:1 - 群3:2
  パラメトリックHR検定:χ2乗=3.56878  自由度=1  有意確率:p=0.176627
  β=0.576333  標準誤差=0.30508  ハザード比(群3/群2)=1.7795
  HRの95%同時信頼区間=0.857253〜3.69392  β=-0.154022〜1.30669

○重症度 (0:症状なし 1:軽症 2:重症)のBonferroni型多重比較
・群1:0 - 群2:1
  Cox-Mantelの検定(修正有):χ2乗=3.81784  自由度=1  有意確率:p=0.152128
  交互作用(異質性)の検定  :χ2乗=16.8019  自由度=23  有意確率:p=1
  Coxのβ=0.742665  標準誤差=0.348932  ハザード比(群2/群1)=2.10153
  HRの95%同時信頼区間=0.911494〜4.84526  β=-0.0926706〜1.578
・群1:0 - 群3:2
  Cox-Mantelの検定(修正有):χ2乗=13.7236  自由度=1  有意確率:p=0.000635324***
  交互作用(異質性)の検定  :χ2乗=24.0976  自由度=20  有意確率:p=0.714464
  Coxのβ=1.52242  標準誤差=0.39039  ハザード比(群3/群1)=4.5833
  HRの95%同時信頼区間=1.80008〜11.6698  β=0.587832〜2.45701
・群2:1 - 群3:2
  Cox-Mantelの検定(修正有):χ2乗=2.53718  自由度=1  有意確率:p=0.33358
  交互作用(異質性)の検定  :χ2乗=26.7754  自由度=25  有意確率:p=1
  Coxのβ=0.579613  標準誤差=0.329751  ハザード比(群3/群2)=1.78535
  HRの95%同時信頼区間=0.810742〜3.93154  β=-0.209805〜1.36903


