(24)-chapter18-2描述与推断医学统计学

同学好，我们今天学习医学统计学中生存分析中的描述与推断方法。生存分析的主要内容描述生存过程比较生存过程影响生存事件的影响因素探索生存分析的基本方法非参数法半参数法生存率估计的非参数法寿命表法和Kaplan-Meier法原理：先求各时点生存概率，然后根据乘法原理计算生存率。

生存率估计的非参数法一、寿命表法

当样本量较大时，通常将资料先整理成频数表形式，用寿命表法计算样本资料的生存率及生存率的标准误。例1某人对食管癌患者手术后随访的资料如表1中的（1）~（4）栏，用寿命表法估计生存率。死亡概率：根据死亡概率的定义，相应时段内的死亡概率生存概率：生存率S(t+1)：表示各组段上限对应时点的生存率，即研究对象活满t+1年的概率生存率标准误SE[S(t+1)]：动态变化：前3年死亡危险性逐年增加，而后呈下降趋势，生存概率从反面说明了这一结果生存率的标准误都很小，说明生存率具有代表性；半数以上的病人术后活不到2年，提示此恶性肿瘤对生命威胁大；生存曲线（survivalcurve）

以观察（随访）时间为横轴，以生存率为纵轴，将各个时间点所对应的生存率连接在一起的曲线图。生存曲线是一条下降的曲线，分析时应注意曲线的高度和下降的坡度。平缓的生存曲线表示高生存率或较长生存期，陡峭的生存曲线表示低生存率或较短生存期。二、乘积极限法乘积极限法（product-limitmethod）是由Kaplan和Meier在1958年首先提出，故又称Kaplan-Meier法（K-M法）。主要适用于样本含量较小的资料。例：用某中药+化疗（中药组）和化疗（对照组）两种疗法治疗白血病后，随访记录患者死前存活月数，结果如下。试分别估计两组的生存率并绘制生存率曲线。中药组：102+12+13186+19+269+8+6+43+943124对照组：2+137+11+61113177用Kaplan-Meier法分别计算两组的生存率。（1）将含量为n的样本观察值（生存时间t）由小到大依次排列，秩次i=1,2…n。如遇非截尾值与截尾值相同时，将非截尾值排在前面。乘积极限法步骤如下（2）列出各时点（实为一短的时间单位）开始时的存活数，即期初观察单位数ni。（3）计算各时点死亡概率q及生存概率p(p=1-q）。（4）求活过各时点的生存率S(ti)，等于从开始观察时点到ti时点各生存概率的连乘积（5）按下式算生存率的标准误：（6）绘制生存率曲线。常绘制成阶梯形的曲线，方法是：将各非截尾值及其对应的生存率标在直角坐标纸上，然后将各点垂直向下再水平向右连成阶梯形。（7）必要时可按正态近似法估计总体生存率的可信区间。某时点ti的总体生存率的可信区间为：图2两种疗法治疗后白血病患者的生存率曲线(Kaplan-Meier法)10.900.7930.68670.56110.42130.21曲线特点：1.曲线阶梯形不能用直线或曲线连接相邻的两个生存率散点。随着死亡时点增多，曲线的阶梯形逐渐不明显。2.曲线左连续每一级台阶右端为断点，当前死亡(在该断点死亡)时点处的纵坐标值在下一个台阶。3.曲线尾部不稳定随着时间的增加，观察例数越来越少，误差越来越大，曲线尾部极不稳定。在多组比较时，常发生曲线尾部交叉现象，这很可能是因误差大而出现的一种假象。此时可适当提前终止日期，使得最后一个死亡时点仍有一定的观察例数。某医师收集20例脑瘤患者甲、乙两疗法治疗的生存时间（周）如下：甲疗法组13371015152330乙疗法组

57+1313233030+38424245+

三、中位生存期中位生存期（mediansurvivaltime）又称半数生存期，表示恰好有50％的个体尚存活的时间。中位生存期越长，表示疾病的预后越好；中位生存期越短，预后越差。估计中位生存期常用图解法或线性内插法。生存率的比较

（生存曲线比较的log-rank检验）生存率比较的假设检验方法有参数法、半参数法