生存分析Kaplan—meier法简介.pdfVIP

维普资讯 r froq ； 蝴 磁 … 6 技术政策 生存分析 Kaplan meier法简介 埔 桩 蜗 忠 ●0·--·-·—-一 …f… 摘要：本文扼要介绍了生存分析的三种方法：参数法、半参数法及非参数法 。 附实例 着重介绍了非参数法 中的Kaplan--meier法及时序检验法。最后，对生存分析时应注意 的概念作了解释 1、计算方法 (以l9例 肾移植病人生存率为例) 30～ 维普资讯 《卫生轭抖学》 1994年 第6期 1．1、将样本生存时间X 由小到大依次排列如表 1第 1栏。如遇非截尾值和截尾值相 同时，非截尾值排在前面，而截尾值排在后面 1．2、列出与X相应的死l亡人数 d，如表 1第 2栏 1．3、列出期初病例数 12，如表 1第 3栏。即生存期为某时间时尚存活的病例数 1．4、计算活过各时点的生存率 P (xt) 计算公式为： P ( =丌声= (n-d) 式 中 为连乘的符号。各时点生存率估计值 的连乘积即为活过该时点的概率。例 如 ： P (3)=声l= =0_9474 表 P (x9)=p ×声 = ( ) × ( ) 8947 l5)=p．× ×声，= ( ) × … _842l P (x3)表示病人活过 3天的概率，P (x9)表示病人活过 9天的概率，以此类 推 因截尾值 的 为 l，生存率必然同前一个非截尾值一样，故计算中只估计非截尾值的 P (xt)，上述生存率的Kaplan-meier计算法，又称乘积极限 (Product-limit)法，P (xt)有 PL估计量之称。 薄 2、生存率的标准误计算有两种方法 2．1、按式 2．1计算，如表第 6栏 sp ) =P r)√∑ ] (2．1) 2．2、按式 2．2计算，如 sp f) =P (x (2．2) 例数较多时，两法计算结果相同。但例数逐渐减少时，式 2．1计算的标准误偏小，而 式 2．2计算的标准误偏大。 3、生存率曲线，或称 Kaplan—meier曲线，它是以时间 t为横轴，生存率 P (xt) 为纵轴，表示时间与生存率关系的曲线 由图可直观地分析 比较各样本的生存曲线，也可 对某一病例任一时刻的生存率作出估计；反之可由随意生存率估计生存时间。本例存活天 数中的中位数，即P (xt)=0．5时的估计值在 46~64天间。用线性 内插法： (46_164) ： (46_t)= (05263--0．4737) ： (0．5263--0．50) 由上式即可解得： = 一 — 堕二{ _=舞 =s．9755(天)