医学统计学生存分析

1、MetaMeta分析分析 第二部分第二部分 生存分析生存分析 第一部分第一部分 第一部分 在医学研究中，常常用在医学研究中，常常用追踪的方式追踪的方式来研来研 究事物发展的规律。如：了解某究事物发展的规律。如：了解某药物的疗药物的疗 效效，了解，了解手术的存活时间手术的存活时间，了解某，了解某医疗仪医疗仪 器设备使用寿命器设备使用寿命等等。等等。 生存资料的特点生存资料的特点 (1)(1)包含有包含有结局结局和和时间时间两个方面的信息。两个方面的信息。 (2)(2)结局资料结局资料一般为一般为二分类资料二分类资料，即结局是两对立，即结局是两对立 相互排斥的事件，如生存与死亡、有效与无效等。相互

2、排斥的事件，如生存与死亡、有效与无效等。 (3)(3)生存资料一般需经过生存资料一般需经过前瞻性随访观察前瞻性随访观察才能获得才能获得 ，随访往往从某统一的时点开始，到某规定的时点，随访往往从某统一的时点开始，到某规定的时点 结束，所以生存资料有时也称为随访资料。结束，所以生存资料有时也称为随访资料。 (4)(4)由于失访等原因使一些研究对象的生存时间难由于失访等原因使一些研究对象的生存时间难 判断，导致部分生存时间判断，导致部分生存时间数据不完整数据不完整。 生存资料的数据类型生存资料的数据类型 按观察对象生存时间的完整性可分为两种类型。按观察对象生存时间的完整性可分为两种类型。 (1)(1

3、)完全数据：完全数据：指观察起点到发生结局事件的时间指观察起点到发生结局事件的时间 明确、完整的资料，即具有明确完整生存时间的明确、完整的资料，即具有明确完整生存时间的 数据。数据。 (2)(2)截尾数据：截尾数据：指由于其他因素指由于其他因素( (非研究因素非研究因素) )导致导致 观察对象的生存时间难以明确判断，这种生存时观察对象的生存时间难以明确判断，这种生存时 间数据称为截尾数据。也称为截尾值、删失值或间数据称为截尾数据。也称为截尾值、删失值或 终检值。终检值。 不论截尾数据的产生原因为何，截尾生存时间的不论截尾数据的产生原因为何，截尾生存时间的 计算均为观察起点至截尾点所经历的时间，

4、而且计算均为观察起点至截尾点所经历的时间，而且 一般情况下其准确的生存时间要长于截尾时间，一般情况下其准确的生存时间要长于截尾时间， 常在此类数据的右上角标记常在此类数据的右上角标记“+ +”。 生存分析基本概念生存分析基本概念 u起始事件起始事件是反映生存时间起始特征的事件，如疾是反映生存时间起始特征的事件，如疾 病确诊、某种疾病治疗开始、接触毒物等。病确诊、某种疾病治疗开始、接触毒物等。 u在生存分析随防研究过程中，一部分研究对象可在生存分析随防研究过程中，一部分研究对象可 观察到死亡，可以得到准确的生存时间，它提供观察到死亡，可以得到准确的生存时间，它提供 的信息是完全；这种事件称为的信

5、息是完全；这种事件称为失效事件，失效事件，也称之也称之 为为死亡事件、终点事件。死亡事件、终点事件。 生存时间（生存时间（survival time)survival time)是指任何两个有联系事件之间的是指任何两个有联系事件之间的 时间间隔，常用时间间隔，常用t t表示。狭义的生存时间指患某疾病的病表示。狭义的生存时间指患某疾病的病 人从发病到死亡所经历的时间跨度，广义的生存时间定义人从发病到死亡所经历的时间跨度，广义的生存时间定义 为从某种起始事件到终点事件所经历的时间跨度。为从某种起始事件到终点事件所经历的时间跨度。 如：急性白血病病人从治疗开始到复发为止之间的缓解期，如：急性白血病病

6、人从治疗开始到复发为止之间的缓解期， 冠心病病人两次发作之间的时间间隔，戒烟开始到重新吸烟冠心病病人两次发作之间的时间间隔，戒烟开始到重新吸烟 之间的时间间隔，接触危险因素到发病的时间间隔等。之间的时间间隔，接触危险因素到发病的时间间隔等。 生存分析中最基本的问题就是计算生存时间，要明确规定事生存分析中最基本的问题就是计算生存时间，要明确规定事 件的起点、终点及时间的测度单位，否则就无法分析比较。件的起点、终点及时间的测度单位，否则就无法分析比较。 中位生存时间中位生存时间是指寿命中位数，表示有且只有是指寿命中位数，表示有且只有 50%50%的观察对象还可以活这么长时间。由于截尾的观察对象还可

7、以活这么长时间。由于截尾 数据的存在，中位生存期的计算不同于普通的中数据的存在，中位生存期的计算不同于普通的中 位数，它可以利用生存函数公式或生存曲线图，位数，它可以利用生存函数公式或生存曲线图， 令生存率为令生存率为50%50%时，推算出生存时间。时，推算出生存时间。 生存函数生存函数 生存概率又称为生存率或生存函数，它表示一个病人的生存生存概率又称为生存率或生存函数，它表示一个病人的生存 时间长于时间时间长于时间t t的概率，用的概率，用S(t)S(t)表示。表示。 以时间以时间t t为横坐标，为横坐标，S(t)S(t)为纵坐标所作的曲线称为生存率曲线为纵坐标所作的曲线称为生存率曲线, ,

8、 它是一条下降的曲线，下降的坡度越陡，表示生存率越低或它是一条下降的曲线，下降的坡度越陡，表示生存率越低或 生存时间越短，其斜率表示死亡速率。生存时间越短，其斜率表示死亡速率。 死亡函数死亡函数 表示死亡速率的大小。如以表示死亡速率的大小。如以t t为横坐，为横坐，f(t) f(t) 为纵坐标作出的曲线称为密为纵坐标作出的曲线称为密 度曲线，由曲线上可看出不同时间的死亡速率及死亡高峰时间。纵坐度曲线，由曲线上可看出不同时间的死亡速率及死亡高峰时间。纵坐 标越大，其死亡速率越高，如曲线呈现单调下降，则死亡速率越来越标越大，其死亡速率越高，如曲线呈现单调下降，则死亡速率越来越 小，如呈现峰值，则为

9、死亡高峰。小，如呈现峰值，则为死亡高峰。 p生存分析的主要内容：生存分析的主要内容： 1. 1.描述生存过程（估计生存函数）描述生存过程（估计生存函数） 2.2.比较生存过程（比较生存函数）比较生存过程（比较生存函数） 3.3.影响生存时间的因素分析影响生存时间的因素分析 SPSS中的中的 菜单位置菜单位置 生存率的估计方法有参数法和非参数法。常用非生存率的估计方法有参数法和非参数法。常用非 参数法，非参数法主要有二个，即乘积极限法与参数法，非参数法主要有二个，即乘积极限法与 寿命表法。寿命表法。 乘积极限法，又称乘积极限法，又称Kaplan-MeierKaplan-Meier法，法，主要用于

10、观主要用于观 察例数较少（察例数较少（n50n1,RR1,危险因素危险因素 X4X4为缓解出院后的巩固治疗（为缓解出院后的巩固治疗（0=0=无，无，1=1=有），有），RR1,RR6 weeks) High risk.12 weeks: 31/110 missing from intervention group; 18/113 missing from control group. Reasons differ across groups. Selective reporting (reporting bias) High risk.Three rating scales for cogni

11、tion listed in Methods, but only one (with statistically significant results) is reported. 漏斗图 Funnel plot SR8:陈述结果 森林图森林图 结果总结表结果总结表 SRSR目的是为卫生决策服务，所以同时要用目的是为卫生决策服务，所以同时要用 科普语言科普语言 SR9SR9：解释结果，得出结论：解释结果，得出结论 1. 1.考虑研究的局限性：如发表偏倚考虑研究的局限性：如发表偏倚 2.2.考虑证据强度考虑证据强度 效应大小，变异效应大小，变异 比较其它综述比较其它综述 3.3.考虑可应用性考虑

12、可应用性 外延性：结果还适合于哪些患者外延性：结果还适合于哪些患者 治疗描述：若应用这个治疗，需详细描述治疗描述：若应用这个治疗，需详细描述 4. 4. 考虑考虑benefit/harmbenefit/harm 5.5.考虑经济学评价考虑经济学评价 可支付性可支付性 相对于其它治疗相对于其它治疗 6.6.考虑进一步研究考虑进一步研究 是否需要进一步研究来明确答案是否需要进一步研究来明确答案 得出结论得出结论 综述的结果对临床实践的影响综述的结果对临床实践的影响 综述的结果对进一步研究的影响综述的结果对进一步研究的影响 MetaMeta分析分析 系统综述和系统综述和metameta分析的基本概念

13、分析的基本概念 系统综述的基本步骤介绍系统综述的基本步骤介绍 MetaMeta分析分析 软件操作软件操作-Rev Man-Rev Man MetaMeta分析软件：分析软件： RevMan 操作简易操作简易 系统综述和系统综述和meta分析分析 Stata 功能强大功能强大 发表偏倚定量及图片编辑发表偏倚定量及图片编辑 Comprehensive meta analysis 个性化绘图个性化绘图 Meta-dics R SAS SPSS MetaMeta分析涉及的几个统计处理分析涉及的几个统计处理 统计处理流程：统计处理流程： a异质性检验（齐性检验）。异质性检验（齐性检验）。 b统计合并效应

14、量并进行统计推断。统计合并效应量并进行统计推断。 c采用采用“倒漏斗图倒漏斗图”了解潜在的发表了解潜在的发表 偏倚。偏倚。 异质性检验（齐性检验）异质性检验（齐性检验） MetaMeta分析效应合并时的变异分解：分析效应合并时的变异分解： 1 1）研究内变异：各个研究的样本量不同，抽样误）研究内变异：各个研究的样本量不同，抽样误 差不同；差不同； 2 2）研究间变异：各个研究在设计、对象、偏倚控）研究间变异：各个研究在设计、对象、偏倚控 制等方面存在差异。制等方面存在差异。 异质性检验基本思想：异质性检验基本思想： 假如研究资料间的真实效应量一致，那么假如研究资料间的真实效应量一致，那么 实际

15、效应量间的差异可认为是由误差引起的。实际效应量间的差异可认为是由误差引起的。 异质性检验（齐性检验）异质性检验（齐性检验） 识别异质性：识别异质性：观察几个独立研究的差异情观察几个独立研究的差异情 况，是否合并？况，是否合并？ Q检验检验 p值值0.1表示有异质性表示有异质性 合并效应合并效应 RevMen软件计算结果：软件计算结果： 总体效应检验，总体效应检验，p0.05p0.05有意义有意义 计量资料：均数之差计量资料：均数之差 氟对女童氟对女童 I I I I 掌骨皮质厚度的影响掌骨皮质厚度的影响 标准化均数之差: di =( 1i - 2i) /Sci , i = 1, 2, 3 ,k

16、 权重系数: wi = n 1i + n 2i x x x 合并效应合并效应 合并效应合并效应 亚组分析亚组分析 Subgroup AnalysisSubgroup Analysis 根据性别、年龄、种族等进行分组分析根据性别、年龄、种族等进行分组分析 发表偏倚发表偏倚 Publication BiasPublication Bias 发表偏倚发表偏倚是指“统计学上有意义”的阳性研究 结果较“统计学上没有意义”的阴性研究结果或 无效的(null)的研究结果更容易被发表。 发表偏倚对meta分析结果的真实性和可靠性有 很大的影响，特别是当入选meta分析的研究主要 是以小样本研究小样本研究为主时

17、，发表偏倚常使meta分析 的效应合并值被高估效应合并值被高估，甚至使结论逆转，产生误 导，即本来没有统计学意义的meta分析结果变为 有统计学意义的结果。 发表偏倚发表偏倚 Publication BiasPublication Bias 发表偏倚发表偏倚 Publication BiasPublication Bias 绘制漏斗图分析绘制漏斗图分析 P|t| 95% Conf. Interval -+- slope | .4270146 .286694 1.49 0.197 -.3099558 1.163985 bias | .0547229 .947982 0.06 0.956 -2.3

18、82143 2.491589 MetaMeta分析分析 系统综述和系统综述和metameta分析的基本概念分析的基本概念 系统综述的基本步骤介绍系统综述的基本步骤介绍 MetaMeta分析分析 软件操作软件操作-Rev Man-Rev Man RevMan RevMan软件是软件是CochraneCochrane国际协作网制作和保存国际协作网制作和保存 CochraneCochrane系统评价的一个程序，由北欧系统评价的一个程序，由北欧CochraneCochrane中心中心 制作和更新，是制作和更新，是CochraneCochrane系统评价的一体化、标准化系统评价的一体化、标准化 软件。软

19、件。 该软件的主要特点：该软件的主要特点： 可制作和保存可制作和保存CochraneCochrane系统评价的计划书和全文；系统评价的计划书和全文； 可对录入的数据进行可对录入的数据进行MetaMeta分析并将森林图的分析结果分析并将森林图的分析结果 以图表形式展示；以图表形式展示； 可对可对CochraneCochrane系统评价进行更新。系统评价进行更新。 /revman/download 也可以在也可以在OAOA上下载安装后使上下载安装后使 用用 账账 号号 密密 码码 1.1.熟悉界面（熟悉界面（1 1） 1.1.熟悉界面（熟悉界面（2 2）

20、 干预评价干预评价 诊断试验评诊断试验评 价价 方法学评价方法学评价 同类系统评价同类系统评价 1.1.熟悉界面（熟悉界面（3 3） 1.1.熟悉界面（熟悉界面（3 3） 菜单栏菜单栏 工具栏工具栏 文章结构文章结构 正文正文 摘要摘要 概要，语言概要，语言 更通俗，可更通俗，可 供免费下载供免费下载 背景和目的背景和目的 研究方法研究方法 研究结果研究结果 讨论讨论 2.2.撰写正文撰写正文-标题标题 2.2.撰写正文撰写正文-作者作者 双击双击 1.若已经注册，若已经注册， 则登录则登录Archie作作 者数据库调出信者数据库调出信 息；息； 2.还没注册，则还没注册，则 手动输入信息手动

21、输入信息 2.2.撰写正文撰写正文-作者作者 2.2.撰写正文撰写正文-时间时间 2.2.撰写正文撰写正文-摘要、方法、结果、讨论摘要、方法、结果、讨论 3.3.从从NENE数据库中导出文献数据库中导出文献 4.导入文献到Revman 4.导入文献到Revman 4.4.画流程图画流程图 4.4.画流程图画流程图 5.Meta5.Meta分析分析(1)(1) 5.Meta5.Meta分析分析(1)(1) 5.Meta5.Meta分析分析(2)(2) 二分类变量：二分类变量： RR/OR/ HR/RD 期望方差法期望方差法 一般倒方差法一般倒方差法 连续型变量：连续型变量： 加权均数加权均数 W

22、MD、标准、标准 化均差化均差SMD 5.Meta5.Meta分析分析(2)(2) 分类资料，小概率事件、分类资料，小概率事件、实验实验 组和对照组样本量严重不平衡组和对照组样本量严重不平衡 分类资料，较为稳健，应用广泛分类资料，较为稳健，应用广泛 连续型变量资料连续型变量资料 5.Meta5.Meta分析分析(3)(3) 5.Meta5.Meta分析分析(4)(4) 5.Meta5.Meta分析分析(5)(5) low or unclear risk, high risk 5.Meta5.Meta分析分析(5)(5) 5.Meta5.Meta分析分析(5)(5) 5.Meta5.Meta分析

23、分析(5)(5) 5.Meta5.Meta分析分析(5)(5) 5.Meta5.Meta分析分析(5)(5) 5.Meta5.Meta分析分析(5)(5) 谢谢 谢！谢！ 生存分析基本概念生存分析基本概念 u起始事件起始事件是反映生存时间起始特征的事件，如疾是反映生存时间起始特征的事件，如疾 病确诊、某种疾病治疗开始、接触毒物等。病确诊、某种疾病治疗开始、接触毒物等。 u在生存分析随防研究过程中，一部分研究对象可在生存分析随防研究过程中，一部分研究对象可 观察到死亡，可以得到准确的生存时间，它提供观察到死亡，可以得到准确的生存时间，它提供 的信息是完全；这种事件称为的信息是完全；这种事件称为失

24、效事件，失效事件，也称之也称之 为为死亡事件、终点事件。死亡事件、终点事件。 COXCOX比例风险回归模型比例风险回归模型 在医学研究中，观察对象生存时间往往受多种因在医学研究中，观察对象生存时间往往受多种因 素的影响，如胃癌手术后的生存时间，除了与治疗素的影响，如胃癌手术后的生存时间，除了与治疗 方案有关外，还可能与患者年龄、体质、病情轻重、方案有关外，还可能与患者年龄、体质、病情轻重、 病理类型、用药等情况有关。统计学上将这些因素病理类型、用药等情况有关。统计学上将这些因素 称为协变量。称为协变量。 由于生存时间资料常存在截尾值，生存时间由于生存时间资料常存在截尾值，生存时间t t往往往往