华西统计第十五章.ppt

1、华西公共卫生学院卫生统计教研室,1,随访资料的生存分析,朱 彩 蓉,华西公共卫生学院卫生统计教研室,2,120例食管癌手术后生存情况研究结果如下，若观察期间无人失访，请计算1年生存率和2年生存率。,实例,华西公共卫生学院卫生统计教研室,3,根据传统率的概念计算n年生存率：,华西公共卫生学院卫生统计教研室,4,120例食管癌手术后生存情况研究结果如下，其中术后1年到2年期间有20例患者失访，请计算1年生存率和2年生存率。,实例,华西公共卫生学院卫生统计教研室,5,用该方法计算的n年生存率损失了部分信息，不能对生存率作出正确的估计。,华西公共卫生学院卫生统计教研室,6,华西公共卫生学院卫生统计教研

2、室,7,某医生将22例肺癌患者随机分为两组，分别采用化疗和放化疗联合治疗，从缓解出院日开始随访， 22名患者生存时间(月)如下，试比较化疗和放化疗联合治疗肺癌的疗效有无差别。,实例,华西公共卫生学院卫生统计教研室,8,化疗组： 1，2，3，5，6，9，11，13， 16，26，37 放化疗联合组： 10，11，14，18，22，22，26，32， 38，40，42,华西公共卫生学院卫生统计教研室,9,某医生将22例肺癌患者随机分为两组，分别采用化疗和放化疗联合治疗，从缓解出院日开始随访，22名患者生存时间(月)如下，试比较化疗和放化疗联合治疗肺癌的疗效有无差别。,实例,华西公共卫生学院卫生统计

3、教研室,10,化疗组： 1，2，3，5，6，9+，11，13， 16，26，37+ 放化疗联合组： 10，11+，14，18，22，22，26，32， 38，40+，42+,注：带“+”号的数据表示患者至少存活了多少个月,华西公共卫生学院卫生统计教研室,11,隐含的研究目的,两组生存率是否有差别 两组生存时间是否有差别,华西公共卫生学院卫生统计教研室,12,随访资料具有其特殊性，观察对象既有随访结果又有随访时间 ，不能套用前面所学的统计方法； 当比较不同治疗方案的治疗效果时，仅比较某个时点的生存率，有时可能出现不正确的结论。,华西公共卫生学院卫生统计教研室,13,生存分析(survival a

4、nalysis),将研究对象的随访结局和随访时间两个因素同时结合起来考虑的一种统计分析方法； 能充分利用所获得的信息，达到较准确、全面地评价和比较随访资料的目的。,华西公共卫生学院卫生统计教研室,14,生存资料的特点,蕴涵有结局和时间两方面信息； 结局为两分类互斥事件； 通过随访收集资料； 常因失访等原因造成生存时间数据不完整，分布类型复杂。,华西公共卫生学院卫生统计教研室,15,例15.1 某医师采用手术疗法治疗12例宫颈癌患者，随访时间(月)记录如下：1，2，4，5，7，8+，11，15，18，33+，36，38+。试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时

5、间，并绘制生存曲线。,华西公共卫生学院卫生统计教研室,16,生存分析中的几个基本概念,“死亡”事件 失败事件(failure event) 生存时间(survival time) 死亡概率(mortality probability) 生存概率(survival probability) 生存率(survival rate) 中位生存时间(median survival time),华西公共卫生学院卫生统计教研室,17,死亡事件(death event),又称失效事件(failure event) 死亡事件是一个广义概念，不单是指通常意义下的生物体死亡，而是泛指标志某种处理措施失败或失效的特征

6、事件。,华西公共卫生学院卫生统计教研室,18,失效事件实例,例15.1 宫颈癌患者死亡; 接受健康教育后青少年戒烟一段时间后又吸烟； 投保某种健康险后中途退保。,华西公共卫生学院卫生统计教研室,19,生存时间(survival time),指观察到的存活时间，可用天、周、月、年等时间单位记录，常用符号t表示； 生存时间是一个广义概念，不单是指通常意义下生物体的存活时间，而是泛指研究者所关心的某现象的持续时间。,华西公共卫生学院卫生统计教研室,20,如：,例15.1 宫颈癌患者生存的时间; 白血病患者化疗后的缓解时间； 肾移植患者的存活时间； 接受健康教育后青少年的戒烟时间； 投保某种健康险至中

7、途退保的时间等。,华西公共卫生学院卫生统计教研室,21,华西公共卫生学院卫生统计教研室,22,生存时间的两种类型,完全数据(complete data)： 指从观察起点到发生死亡事件所经历的时间，如表15.1中2和6号患者对应的生存天数89天和85天。 截尾数据(censored data)： 简称截尾值(censored value)，又称删失值或终检值。习惯上在生存时间右上标注“”表示。 生存时间观察过程的截止不是由于死亡事件或失效事件，而是由于其他原因引起的，称截尾(censored)。,华西公共卫生学院卫生统计教研室,23,月 份（如1998年）,生存时间原始记录示意图(+为仍存活,

8、下同),华西公共卫生学院卫生统计教研室,24,天 数,图2 生存时间排序整理数据示意图,华西公共卫生学院卫生统计教研室,25,截尾的主要原因,失访(withdrawal)：指失去联系，如信访无回音、电话采访不应答、上门采访找不到人、搬迁没留地址等； 退出：指死于非研究因素或非处理因素而退出研究，如死于车祸等意外事件、死于其它疾病等； 终止：指设计时规定的研究时限已到而终止观察，但研究对象仍然存活。,华西公共卫生学院卫生统计教研室,26,死亡概率与生存概率,死亡概率(mortality probability)：记为 ，是指在某单位时段开始时存活的个体在该时段内死亡的可能性大小。 生存概率(su

9、rvival probability)：记为 ，与死亡概率相对立，表示在某单位时段开始时存活的个体到该时段结束时仍存活的可能性大小。,华西公共卫生学院卫生统计教研室,27,生存率(survival rate),指某个观察对象活过t 时刻的概率，常用p(xt)表示； 根据不同随访资料的失败事件 ，生存率可以是缓解率、某种状态的维持率等。,华西公共卫生学院卫生统计教研室,28,生存曲线(survival curve),将各个时点的生存率在坐标轴上连接在一起的曲线图； 用以描述生存过程。,华西公共卫生学院卫生统计教研室,29,中位生存时间(median survival time),又称半数生存期，

10、指生存率为0.5时对应的生存时间，表示有50%的观察对象可以活这么长时间； 一般采用内插法进行估计。,华西公共卫生学院卫生统计教研室,30,图15.1乘积极限法生存曲线(Kaplan-Meier曲线),华西公共卫生学院卫生统计教研室,31,生存资料的2种类型,未分组资料 分组资料,华西公共卫生学院卫生统计教研室,32,例15.1 某医师采用手术疗法治疗12例宫颈癌患者，随访时间(月)记录如下：1，2，4，5，7，8+，11，15，18，33+，36，38+。试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。,华西公共卫生学院卫生统计教研室,33,华西

11、公共卫生学院卫生统计教研室,34,例15.2某医师对110例原发性肺癌患者确诊后进行随访，得资料见表15.3第(2)栏第(4)栏，试估计各时点生存率及其标准误、各时点总体生存率的95%置信区间、中位生存时间，并绘制生存曲线。,华西公共卫生学院卫生统计教研室,35,华西公共卫生学院卫生统计教研室,36,未分组资料的生存分析,未分组资料生存率的估计：主要是用乘积极限法(product-limit estimates)，又称Kaplan-Meier法； 主要适用于小样本，是一种非参数方法。,华西公共卫生学院卫生统计教研室,37,例15.1 某医师采用手术疗法治疗12例宫颈癌患者，随访时间(月)记录如

12、下：1，2，4，5，7，8+，11，15，18，33+，36，38+。试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。,华西公共卫生学院卫生统计教研室,38,华西公共卫生学院卫生统计教研室,39,生存率及其标准误的计算,本例生存时间以“月”为单位，并将t月当作一个时点看待； 具体计算过程见表15.2；,华西公共卫生学院卫生统计教研室,40,第(1)栏为序号：本例k = 1，2，312； 第(2)栏是将生存时间t由小到大顺次排列； 如某时间点既有完全数据又有截尾数据，将截尾数据排在后面。,15.2中各栏的含义,华西公共卫生学院卫生统计教研室,41,第

13、(3)栏为t月的死亡人数dk，如生存时间t为5个月时有1例死亡，相应的d4=1。 截尾患者即便是已死亡，也非死于处理措施(手术疗法)，如生存时间t为8月时有1例截尾，相应的d6=0。,华西公共卫生学院卫生统计教研室,42,第(4)栏为期初病例数nk，即恰好在t时点以前尚存活的病例数 如t为5月时对应的n4=9，表示恰好在5月时点前有9人存活,华西公共卫生学院卫生统计教研室,43,第(5)栏计算各时点死亡概率qk，即恰好在t时点以前尚存活的患者在t时点上(第t个月)死亡的概率，其计算公式为： 如“q4”表示恰好在5月时点前尚存活的9例患者在5月时点上(实指治疗后第五个月的第一天到第三十天)的死亡

14、概率为,华西公共卫生学院卫生统计教研室,44,第(6)栏计算各时点生存概率pk，即恰好在t时点以前尚存活的患者在t时点上(第t个月)继续存活的概率，其计算公式为 如t为5月时对应：,华西公共卫生学院卫生统计教研室,45,第(7)栏计算各时点生存率S(tk)，即恰好在t时点以前尚存活的患者活过t时点的概率，按式(15.5)计算： 如生存时间t为4月的生存率为：,华西公共卫生学院卫生统计教研室,46,第(8)栏为各时点生存率的标准误，其计算公式为：,华西公共卫生学院卫生统计教研室,47,如S(t3)的标准误SES(t3)为：,华西公共卫生学院卫生统计教研室,48,各时点总体生存率的可信区间,采用正

15、态近似原理估计生存率的可信区间,华西公共卫生学院卫生统计教研室,49,本例4月总体生存率的95%可信区间 下限 上限 宫颈癌患者术后4月生存率的95%可信区间为50.5099.50%。,华西公共卫生学院卫生统计教研室,50,采用内插法 由表15.2可见，中位生存时间估计在711月之间，计算如下：,中位生存时间的计算,华西公共卫生学院卫生统计教研室,51,华西公共卫生学院卫生统计教研室,52,图15.1乘积极限法生存曲线(Kaplan-Meier曲线),华西公共卫生学院卫生统计教研室,53,(711): (7t) = (0.58330.4861) : (0.58330.5),即宫颈癌患者术后中位

16、生存时间为10.4个月。,华西公共卫生学院卫生统计教研室,54,生存曲线,它是以生存时间t为横轴，生存率 为纵轴，绘制而成的连续型的阶梯形曲线，说明生存时间与生存率之间的关系； 未分组资料的生存曲线也称Kaplan-Meier曲线。,华西公共卫生学院卫生统计教研室,55,图15.1乘积极限法生存曲线(Kaplan-Meier曲线),华西公共卫生学院卫生统计教研室,56,例15.2某医师对110例原发性肺癌患者确诊后进行随访，得资料见表15.3第(2)栏至第(4)栏，试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。,华西公共卫生学院卫生统计教研室,5

17、7,寿命表法 (life-table method),当样本例数足够多时，乘积极限法可按时间分组计算，这就是寿命表法。 寿命表法的提出早于乘积极限法，但实质上是乘积极限法的一种近似(频数表法)。,华西公共卫生学院卫生统计教研室,58,生存率及其标准误的计算,本例生存时间以“月”为单位，具体计算过程见表15.3；,华西公共卫生学院卫生统计教研室,59,第(1)栏为序号：本例k = 1，2，310； 第(2)栏是为确诊后月数tk。“0”表示从确诊日起不满1月，“1”表示确诊后满1月但未满2月，其他依此类推。,15.3中各栏的含义,华西公共卫生学院卫生统计教研室,60,第(3)栏为期内截尾人数Ck，

18、表示确诊后满t月但未满t+1月期间截尾的人数。 如C3=1，表示确诊后满2月但未满3月期间有1例截尾。,华西公共卫生学院卫生统计教研室,61,第(4)栏为期内死亡人数 ，表示确诊后满t月但未满t+1月期间发生死亡事件的人数。 如 =25，表示确诊后未满1月有25例死于原发性肺癌； =18表示确诊后满3月但未满4月有18例死于原发性肺癌。,华西公共卫生学院卫生统计教研室,62,第(5)栏为期初观察人数 ，指t月初尚存活的病例数。 此栏自下而上累计求得，见式 (15.8) 如本例 =4，则 余类推,华西公共卫生学院卫生统计教研室,63,第(6)栏为校正人数Nk，即实际观察人月数; 凡在t月内截尾的

19、病例,按平均观察半月计算; 因此校正人数Nk等于月初观察人数减去当月截尾人数的一半，其计算公式为:,华西公共卫生学院卫生统计教研室,64,第(7)栏为死亡概率 ，表示确诊后活满t月的病人在今后一个月内死亡的概率，其计算公式为：,华西公共卫生学院卫生统计教研室,65,第(8)栏为生存概率 ，表示确诊后活满t月的病人在今后一个月内存活的概率，其计算公式为,华西公共卫生学院卫生统计教研室,66,第(9)栏为t月生存率，表示确诊后活满t月的概率，采用式(15.5)计算得到。 如1月生存率 2月生存率 3月生存率 余类推,华西公共卫生学院卫生统计教研室,67,第(10)栏为生存率的标准误，其计算公式为：

20、,华西公共卫生学院卫生统计教研室,68,分组资料与未分组资料的生存率计算比较,分组资料与未分组资料生存率的计算相似 。 未分组资料生存率的计算是以观察人数为计算单位，生存时间未分组； 而分组资料的生存率的计算是以观察人时数为单位，生存时间等距离分组。,华西公共卫生学院卫生统计教研室,69,统计描述主要包括各组生存率及其标准误的计算; 统计推断包括总体生存率的可信区间估计、两样本或多样本生存率比较等。,生存率分析的内容,华西公共卫生学院卫生统计教研室,70,1. 描述生存过程,研究生存时间的分布特点，估计生存率、生存率曲线等。 例如 根据白血病化疗后的缓解年数资料，可以估计不同年数的缓解率p(X

21、t)，如p(X3)， p(X5)等，也可以获得这些病人的缓解率曲线。,华西公共卫生学院卫生统计教研室,71,2. 比较生存过程,在获得生存率及其标准误的估计值后，进行两组或多组生存率比较。 例如 比较不同方案治疗白血病的缓解率曲线，以了解哪种方案较优。,华西公共卫生学院卫生统计教研室,72,3. 生存过程的影响因素分析,例如 为了改善乳腺癌手术后的预后，应该了解影响病人预后的主要因素，从而在手术前、后加以预防和干预。,华西公共卫生学院卫生统计教研室,73,两样本生存率曲线的比较,华西公共卫生学院卫生统计教研室,74,例15.322例肺癌患者经随机化分配到放疗组和放化疗联合组，从缓解出院日开始随

22、访，随访时间(月)如下，试比较放疗组和放化疗联合组患者的生存曲线有无差别。 甲：放疗组1，2，3，5，6，9+，11，13，16，26，37+ 乙：放化疗联合组10，11+，14，18，22，22，26，32， 38， 40+，42+,华西公共卫生学院卫生统计教研室,75,图15.3放疗组和放化疗联合组患者的生存曲线,华西公共卫生学院卫生统计教研室,76,建立检验假设和确定检验水准 H0：放疗组与放化疗联合组患者的生存曲线位 置相同 H1：放疗组与放化疗联合组患者的生存曲线位 置不同 0.05,华西公共卫生学院卫生统计教研室,77,华西公共卫生学院卫生统计教研室,78,2. 计算统计量 值,将两组患者按生存时间统一从小到大排序，并标明组别、死亡数。 生存时间