作者gsuper (数理统计-九阳真经)
看板Statistics
标题[问题] 两种hazard的算法 , 哪个对?
时间Fri Sep 2 23:39:59 2011
10个月的 example
以月为"瞬时"的时间单位
瞬时存活率为条件机率 = 月底存活/月初存活 = 当月存活率
并假定无 censor
-----------------------------------------------------------------------
月份 月初存活 月底存活 当月死亡 St 瞬时存活率 |
hazard_A hazard_B
1 1000 975 25 0.975 0.975 | 0.025 0.025
2 975 950 25 0.950 0.974 | 0.026 0.026
3 950 900 50 0.900 0.947 | 0.053 0.054
4 900 800 100 0.800 0.889 | 0.111 0.118
5 800 600 200 0.600 0.750 |
0.250 0.288
6 600 480 120 0.480 0.800 | 0.200 0.223
7 480 420 60 0.420 0.875 | 0.125 0.134
8 420 380 40 0.380 0.905 | 0.095 0.100
9 380 350 30 0.350 0.921 | 0.079 0.082
10 350 320 30 0.320 0.914 | 0.086 0.090
-------------------------------------------------------------------------
第一种算法 : hazard_A
1 - 瞬时存活率
= 当月死亡人数 / 当月月初存活人数
= 瞬时死亡率
= hazard
----------------------------------------------------
第二种算法 : hazard_B
-log(瞬时存活率)
= hazard
----------------------------------------------------
依公式来看应该是第二种算法
但我看到一篇 paper 举的计算实例是第一种算法
想请问哪个是正确的?
因为这个例子在 5月 的时候 (标记绿色)
两者的差异直达到顶峰
也就是说 , 虽然看起来很像 , 但还是不一样
请问这是怎麽回事?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.239.247
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:40)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:40)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:43)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:44)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:46)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:48)
※ 编辑: gsuper 来自: 140.113.239.247 (09/02 23:49)
1F：→ bmka:The are different estimators for the same parameter. 09/03 12:17
2F：→ bmka:一个是从product-limit estimator 来的, 一个则是利用 09/03 12:31
3F：→ bmka:survival function 跟hazard function的关系. 09/03 12:31
4F：→ bmka:但是严格讲後者的关系式要continuous random variable才对 09/03 12:33
5F：→ gsuper:所以每个月电访一次的资料应该属於不连续资料 , 选方法一吧 09/03 19:58