医学统计学(11)•Meta分析•第二部分•生存分析•第一部分•生存分析•第一部分在医学研究中,常常用追踪的方式来研究事物发展的规律。如:了解某药物的疗效,了解手术的存活时间,了解某医疗仪器设备使用寿命等等。生存资料的特点(1)包含有结局和时间两个方面的信息。(2)结局资料一般为二分类资料,即结局是两对立相互排斥的事件,如生存与死亡、有效与无效等。(3)生存资料一般需经过前瞻性随访观察才能获得,随访往往从某统一的时点开始,到某规定的时点结束,所以生存资料有时也称为随访资料。(4)由于失访等原因使一些研究对象的生存时间难判断,导致部分生存时间数据不完整。生存资料的数据类型•按观察对象生存时间的完整性可分为两种类型。(1)完全数据:指观察起点到发生结局事件的时间明确、完整的资料,即具有明确完整生存时间的数据。(2)截尾数据:指由于其他因素(非研究因素)导致观察对象的生存时间难以明确判断,这种生存时间数据称为截尾数据。也称为截尾值、删失值或终检值。不论截尾数据的产生原因为何,截尾生存时间的计算均为观察起点至截尾点所经历的时间,而且一般情况下其准确的生存时间要长于截尾时间,常在此类数据的右上角标记“+”。生存分析基本概念起始事件是反映生存时间起始特征的事件,如疾病确诊、某种疾病治疗开始、接触毒物等。在生存分析随防研究过程中,一部分研究对象可观察到死亡,可以得到准确的生存时间,它提供的信息是完全;这种事件称为失效事件,也称之为死亡事件、终点事件。生存时间(survivaltime)是指任何两个有联系事件之间的时间间隔,常用t表示。狭义的生存时间指患某疾病的病人从发病到死亡所经历的时间跨度,广义的生存时间定义为从某种起始事件到终点事件所经历的时间跨度。如:急性白血病病人从治疗开始到复发为止之间的缓解期,冠心病病人两次发作之间的时间间隔,戒烟开始到重新吸烟之间的时间间隔,接触危险因素到发病的时间间隔等。生存分析中最基本的问题就是计算生存时间,要明确规定事件的起点、终点及时间的测度单位,否则就无法分析比较。中位生存时间是指寿命中位数,表示有且只有50%的观察对象还可以活这么长时间。由于截尾数据的存在,中位生存期的计算不同于普通的中位数,它可以利用生存函数公式或生存曲线图,令生存率为50%时,推算出生存时间。生存函数生存概率又称为生存率或生存函数,它表示一个病人的生存时间长于时间t的概率,用S(t)表示。以时间t为横坐标,S(t)为纵坐标所作的曲线称为生存率曲线,它是一条下降的曲线,下降的坡度越陡,表示生存率越低或生存时间越短,其斜率表示死亡速率。死亡函数表示死亡速率的大小。如以t为横坐,f(t)为纵坐标作出的曲线称为密度曲线,由曲线上可看出不同时间的死亡速率及死亡高峰时间。纵坐标越大,其死亡速率越高,如曲线呈现单调下降,则死亡速率越来越小,如呈现峰值,则为死亡高峰。生存分析的主要内容:1.描述生存过程(估计生存函数)2.比较生存过程(比较生存函数)3.影响生存时间的因素分析SPSS中的菜单位置生存率的估计方法有参数法和非参数法。常用非参数法,非参数法主要有二个,即乘积极限法与寿命表法。乘积极限法,又称Kaplan-Meier法,主要用于观察例数较少(n50)而未分组的生存资料;寿命表法适用于观察例数较多而分组的资料,不同的分组寿命表法的计算结果亦会不同,当分组资料中每一个分组区间中最多只有1个观察值时,寿命表法的计算结果与乘积极限法完全相同。描述生存过程•乘积极限法:Kaplan-Meier•寿命表法:寿命表•因素分析:Cox回归【例1】某治疗方案治疗晚期肝癌患者7例,随访时间(月)记录如下。试分析其生存过程。=╳SPSS软件操作•第一步:建立变量。•第二步:输入原始数据•第三步:生存分析(1)•第三步:生存分析(2)•第三步:生存分析(3)•第四步:结果解读(1)结果解读:生存表基本描述•第四步:结果解读(2)结果解读:生存表统计描述中位生存时间及95%CI•第四步:结果解读(3)结果解读:生存函数图【例2】某医院对100例胰腺癌切除术后的患者进行随访,得资料如下。试分析其生存过程。•(7)=(6)-(5)/2•(8)=(4)/(7)•(9)=1-(8)SPSS软件操作•第一步:建立变量。•第二步:输入原始数据•第三步:加权个案•第四步:生存分析(1)•第四步:生存分析(2)•第三步:生存分析(3)•第四步:结果解读(1)结果解读:寿命表描述•第四步:结果解读(2)结果解读:生存函数生存曲线的比较•对数秩检验:其基本思想是,在假定无效假设(两总体生存曲线相同)成立的前提下,可根据不同日期两种处理的期初人数和死亡人数,计算各种处理在各个时期的理论死亡数。若无效假设成立,则实际死亡数与理论死亡数不会相差太大。•Gehan比分检验:其基本思想是,在假定无效假设(两总体生存曲线相同)成立的前提下,则两样本来自分布相同的总体,两样本的Gehan比分合计V值应为0,若V值偏离0太远,则无效假设成立的可能性就很小。【例3】某医生将20例某恶性肿瘤患者随机分为两组,一组采用纯中药治疗(甲组),另一组采用中西医结合治疗(乙组),观察终点为死亡,从缓解出院日开始随访。试分析两种治疗方案的效果有无差别?•(6)/(11)=总死亡率•甲组理论死亡数=(9)*总死亡率•乙组理论死亡数=(10)*总死亡率SPSS软件操作•第一步:建立变量。•第二步:输入原始数据•第四步:生存分析(1)•第四步:生存分析(2)•第三步:生存分析(3)•Log-rank•第三步:生存分析(4)•第四步:结果解读(1)结果解读:基本描述•第四步:结果解读(2)•第四步:结果解读(3)结果解读:生存表统计描述中位生存时间及95%CI•第四步:结果解读(4)结果解读:组间比较的LogRank检验•第四步:结果解读(5)结果解读:生存曲线【例4】20例肺癌患者随机分为放疗组和放疗加中药联合治疗组,从缓解出院日开始随访,随访时间如下。试分析两种治疗方案的效果有无差别?组别随访时间放疗组151719+232835+46525560放+中1926384048+5358+626875•R1:肯定比他生存时间短的病例数•R2:肯定比他生存时间长的病例数•相当于U检验,与1.96比较SPSS软件操作•第一步:建立变量。•第二步:输入原始数据•第三步:生存分析(1)•第三步:生存分析(2)•第三步:生存分析(3)•第四步:结果解读(1)结果解读:寿命表•第四步:结果解读(2)结果解读:中位生存时间•第四步:结果解读(3)结果解读:生存曲线•第四步:结果解读(4)结果解读:组间比较Z=1.706,p=0.191COX比例风险回归模型在医学研究中,观察对象生存时间往往受多种因素的影响,如胃癌手术后的生存时间,除了与治疗方案有关外,还可能与患者年龄、体质、病情轻重、病理类型、用药等情况有关。统计学上将这些因素称为协变量。由于生存时间资料常存在截尾值,生存时间t往往不能满足正态分布和方差齐性的要求,一般不适宜用参数方法(如多元线性回归等)来分析生存时间与各协变量之间的关系。为解决这类问题,英国生物统计学家D.RCox于1972年提出比例风险回归模型用于分析带有协变量的生存时间资料。【例5】50例某恶性肿瘤的生存时间及有关预后因素资料如表。X1为患者年龄、X2为病理类型(0=低分化,1=高分化),X3为有无淋巴结浸润(0=无,1=有),X4为缓解出院后的巩固治疗(0=无,1=有),T为患者的生存时间,S为患者的随访结局(0=截尾,1=死亡)。试进行生存分析。SPSS软件操作•第一步:建立变量。•第二步:输入原始数据•第三步:生存分析(1)•第三步:生存分析(2)•第三步:生存分析(3)•第三步:生存分析(4)•第四步:结果解读(1)结果解读:有意义的因素,相对危险度RR及其95%CI。•X3为有无淋巴结浸润(0=无,1=有),RR1,危险因素•X4为缓解出院后的巩固治疗(0=无,1=有),RR1,保护因素•第四步:结果解读(2)结果解读:生存函数,具体到个案的生存函数估计值。PH(比例风险)假定判断满足前提条件才可以进行Cox回归是否满足前提条件,看以因素为分类依据的生存曲线是否不交叉,若有交叉则表示不满足前提条件•第五步:PH假定判断(1)•第五步:PH假定判断(2)•第六步:结果解读(1)•第六步:结果解读(2)同理将X3改为X4•重复相同的步骤•第五步:PH(比例风险)假定判断(1)•第五步:PH假定判断(2)•第六步:结果解读(3)•第六步:结果解读(4)•Meta分析•第二部分研究时有无设计干预因素?实验性研究观察性研究是否随机有无对照组分析性研究描述性研究时间方向队列研究非随机对照试验随机对照试验病例对照研究横断面研究有无有无否是暴露→结局结局→暴露Meta分析•系统综述和meta分析的基本概念•系统综述的基本步骤介绍•Meta分析•软件操作-----RevMan对某一具体的临床问题,系统、全面地收集所有已发表或未发表的相关的临床研究文章,用统一、科学的评价标准筛选出合格的研究质量评价,应用统计学方法定量综合/描述性方法进行定性综合,得出可靠的结论,并随着新的临床研究结果的出现及时作出更新。系统综述Meta分析•Ameta-analysisisatwo-stageprocess.–提取单个研究的数据,并估计其进行点估计和可信区间;–决定是否合适将结果汇总,若是,计算其汇总值。•Meta分析不仅是简单将单个研究的数据累加90系统综述的特征:最佳证据•规范的临床问题•全面、完整的资料•对原始研究的质量评价,纳入合格的研究•统计学综合(meta分析)/描述性综合(偏倚)•进行敏感性分析和亚组分析•规范的评价报告(明确的方法学可允许重复)•随着新的研究结果出现进行及时更新。••综述•Meta分析•SR•From:CriticalAppraisalSkillsProgramme(CASP),Oxford.Asystematicreviewmayhaveastatisticalcombinationofstudies(ameta-analysis)butitdoesnothaveto.92传统综述和SR比较研究问题:涉及范围泛常集中于某一问题文献来源:不全面明确,常为多渠道检索方法:常未说明有明确检索策略文献选择:有潜在偏倚有明确入选/排除标准文献评价:方法不统一有严格评价方法结果合成:定性定量/定性结论推断:有时遵循研究依据大多遵循研究依据结果更新:不定期更新依据新试验定期更新传统综述SR93为什么进行系统评价/meta分析•信息太多、医生太忙、时间太少•减少偏倚和随机误差•提高统计效度•探索研究结果的不一致性•临床决策依据•为以后的知识更新提供一个系统简易的方法为临床实践提供一个好的证据临床医生需要获得明确的信息系统综述的类型94•定量SR/定性SR•干预研究SR•诊断试验准确性SR•频率meta分析•生存资料meta分析•Networkmeta分析•。。。。。。•近年来全球Meta分析文章的数量变化趋势•近年来中国作者Meta分析文章的数量变化趋势•发表Meta分析论文最多的十大期刊•最青睐中国meta分析论文的十大杂志Meta分析•系统综述和meta分析的基本概念•系统综述的基本步骤介绍•Meta分析•软件操作-----RevManSR步骤(2008CochraneSR手册)1.建立一个规范化的问题2.制定纳入研究的标准3.检索研究4.筛选研究和收集数据5.研究质量的严格评阅6.分析数据并在可能的情况下进行meta-分析7.解决报告的偏倚8.陈述结果并制作结果摘要表格9.解释结果,得出结论100PICO_S始终贯穿SRSR1:建立一个规范化问题•PICO的含义•临床问题vs基础问题•PICO的概念延伸•PICO_S始终贯穿SRPopulationInterventionComparatorOutcomesSt