《市政工程实战案例》第06讲生存分析I.ppt

两个指数分布模型的比较(似然比检验) H0：?1=?2(两组患者的生存率是相同的)； H1：?1??2(两组患者的生存率不相同)。 按自由度为1的?2界值作出决策。 例8.2 两组肿瘤患者的生存时间如下： 三期肿瘤患者的生存时间： 6 19 32 42 42 43+ 94 126+ 169+ 207 211+ 227+ 253 255+ 270+ 310+ 316+ 335+ 346+ 四期肿瘤患者的生存时间： 4 6 10 11 11 11 13 17 20 20 21 22 24 24 29 30 30 31 33 34 35 68 39 40 41+ 43+ 45 46 50 56 61+ 61+ 63 82 85 88 89 90 93 104 110 134 137 160+ 169 171 173 175 184 201 222 235+ 247+ 260+ 284+ 290+ 291+ 302+ 304+ 341+ 345+ Kaplan-Meier估计 n1=19，r1=8， T1=3303， n2=61，r2=46，T2=6415， ?1=0.002422, ?2=0.007171, ?C=0.005557 G = -2?[(0.002422+0.007171)?log(0.005557)- -0.002422?log(0.002422) -0.007171?log(0.007171) = 10.171342 P = 0.0014 可以认为三期病人比四期病人的生存时间长。 生存函数与风险函数： 三期： S(t)=e-0.002422 t h(t)=0.002422 四期： S(t)=e-0.007171 t h(t)=0.007171 四期与三期的风险比： h4/h3=2.96 肿瘤患者生存时间的Kaplan-Meier 估计与指数分布模型拟合 ? Weibull分布 如果生存资料服从Weibull分布，则 生存函数为： 风险函数为： m0为形状参数，决定了生存时间分布的形状。 m1 早夭现象 m1 老化现象 m=1 指数分布 ?=1时，不同m的Weibull分布生存率的变化规律 Weibull分布的均数和中位数 在已知?和m时，平均生存时间为： 中位生存时间 ?(?)为Gamma函数。 /wiki/%E9%9F%A6%E4%BC%AF%E5%88%86%E5%B8%83 Weibull分布的方差 其标准差为： 例8.2两组资料的Weibull分布拟合 例8.2资料的Weibull分布拟合 文献6导读(静脉血栓病人口服凝血酶抑制剂与安  慰剂组比较，两组静脉血栓再发生率) 文献6导读(静脉血栓病人口服凝血酶抑制剂与安 慰剂组比较，两组出血发生率) 正确应用 生存分析用于对随访研究资料的处理，描述生存过程、比较生存过程、对影响生存过程的因素进行分析。尤其是可以充分利用失访资料所提供的不完全信息。但一般控制失访率20%。 描述：Kaplan-Meier方法、指数分布法； 比较：logrank检验。 分析：回归模型。 * 生存过程的描述 非参数方法 参数方法 乘积-极限法 指数分布 (又称Kaplan-Meier法) Weibull分布 寿命表方法 Gamma分布 Turnbull估计 logistic分布 对数正态分布 生存率的估计: Kaplan-Meier法 又称：乘积极限法( product-limit，PL法) S(ti)=S(t1|0)S(t2|t1)S(t3|t2)…S(ti |ti-1 ) =p1? p2 ? p3 ???pi = S(ti-1 )S(ti |ti-1 ) 要在时刻ti生存，必须ti及ti 前的任何时刻均存活，而每一时刻存活的概率为pi，则在时刻ti的生存率为此前各时刻生存概率之乘积。