Meta分析原理与实践课件

1、Meta分析实践原理vs.首都医科大学 陈卉Meta分析与系统评价Meta分析基本过程Meta分析实例一二三2 系统评价Meta分析Meta分析与系统评价 4 系统评价（systematic review，SR）全面收集针对某一具体问题的所有相关研究，并对这些独立研究进行严格的质量评价后，将它们的结果联合起来进行综合分析（合成）和评价5个阶段提出要评价的问题检索和筛选研究评价独立研究的质量综合分析数据和表达结果解释结果1.1 系统评价Meta分析足够的研究数量满足一定的相似性 5 Meta分析（meta-analysis）Meta分析也称为荟萃分析，狭义地讲，是指根据研究的试验设计类型以及反映

2、待评价问题的统计学指标，将多个独立研究的指标进行合成，并对结果进行假设检验通过Meta分析可以增大样本量，提高检验效能，增强结论的可靠性和客观性1.2 Meta分析例1 为评价肾移植术后CNI撤除后以TOR-I替代治疗的疗效，研究者收集了8个随机对照试验，观察试验组与对照组（持续CNI治疗）患者术后1年发生急性排斥反应的情况研究编号作者年份试验组ARR/所有对照组ARR/所有1Watson20050 / 190 / 192Stallone20032 / 202 / 203Baboolal20038 / 423 / 454Grinyo200410 / 444 / 435Schena20097 /

3、 5554 / 2756Johnson200121 / 2159 / 2157Lebranchu200917 / 9510 / 978Gonwa200222 / 10018 / 97资料来源：吴斌等. 中国循证医学杂志，2010，10(1)：33-39.文献实例 7 例2 为评价降压类药物对肾移植患者的影响，研究者收集了11个随机对照试验，观察试验组（CCB）与对照组（安慰剂）患者服药2周后的肾小球滤过率文献实例 8 研究编号作者年份试验组 例数 均数 标准差对照组 例数 均数 标准差1Gossmann20023154.921.52948.318.72Guerin198912341512331

4、03Harpers19964243.912.33639.412.64Kuypers20045350.116.2534215.45Madsen19982749182440166Rump2000244719.6264930.67Van der Schaff19952054202049168Venkat-Raman19992349202347159Wilkie199315641915652110Wilkie1993785277893211Wilkie1994176425176425资料来源：Cross NB, et al. Transplantation 2009; 88: 718 10 二分类资料

5、常用效应量RD：风险差或率差（risk difference或ratio difference）RR：风险比或相对危险度（risk ratio或relative risk）OR：比值比（odds ratio）选择依据队列研究或RCT：首选RR，两组基线水平一致时可选RD病例对照研究：只能计算OR，试验组和对照组的事件发生率均较低（通常低于20%）时可选RR需要知道每个研究中试验组与对照组事件发生及未发生的例数2.1 确定效应量 11 连续型资料加权均数差（weighted mean difference，WMD）以研究原有的度量真实反映了试验效应，在实际分析时易于理解和使用标准化均数差（sta

6、ndardized mean difference，SMD）度量单位不同和/或均数差相差较大时选用没有度量单位，只能说明统计学意义，难以结合临床意义进行结果的解释需要知道每个独立研究试验组和对照组的样本例数、原始均数及标准差2.1 确定效应量 12 异质性（heterogeneity）指独立研究的资料之间存在变异资料间异质性的程度决定了Meta分析的方法异质性检验方法计算Q统计量并进行卡方检验计算I2统计量并直接判断2.2 异质性分析 14 I2（I-squared）统计量法I2 25%：异质性较小25%50%：异质性较大两种方法的选择Q统计量法检验效能较低，检验结果不可靠I2统计量校正了研究

7、数量k对Q值的影响，判断异质性的结果更为稳健可靠2.2 异质性分析可接受 15 研究间存在较明显的异质性时的处理方法2.2 异质性分析异质性分析固定效应模型异质性原因亚组分析Meta回归随机效应模型无有放弃Meta分析无有 17 判断合并效应量的统计学意义Z检验法ZZ0.05则P0.05，合并效应量有统计学意义根据效应量的95%置信区间判断对于WMD、SMD或RD：95%CI不包含0，则P0.05对于OR或RR： 95%CI不包含1，则P0.052.3 合并效应量 18 用森林图（forest plot）展示分析过程及结果2.3 合并效应量异质性检验I2=0.0% 0.1095%CI不包含1合

8、并效应量RR=1.6495%CI: 1.192.27 19 发表性偏倚（publication bias）有统计学意义的研究结果（阳性结果）比无统计学意义的结果（阴性结果）更容易被投稿和发表产生来源设计方案的选择研究者的主观愿望资助者的利益识别方法漏斗图法Begg秩相关法Egger线性回归法2.4 识别发表性偏倚结合使用 20 漏斗图（funnel plot）以每个研究的效应量（或效应量的对数）为横轴、样本量（或效应量的标准误）为纵轴绘制的散点图如果漏斗图中的散点表现为明显不对称则认为存在发表性偏倚漏斗图方法是一种定性的图示法，研究个数较少时不宜使用2.4 识别发表性偏倚 21 Begg秩相关

9、法对效应量标准化后用秩相关分析来检验效应量是否与其标准误存在相关性假设检验得到P0.05时可以认为不存在发表性偏倚当纳入的独立研究的数量较少时，检验效能相对较低以每个研究的效应量（或效应量的对数）为纵轴、效应量的标准误为横轴绘制的散点图称为Begg漏斗图2.4 识别发表性偏倚 22 Egger直线回归法拟合标准正态离差对效应估计值的精确度（precision）的线性方程直线的截距代表不对称的程度，若假设检验P0.05则可以认为不存在发表性偏倚以每个研究的精确度为横轴、效应量的标准误为纵轴绘制的散点图称为Egger图， 截距的95%置信区间包含0 表示不存在发表性偏倚2.4 识别发表性偏倚 24

10、 纳入的研究及基本数据3.1 二分类资料Meta分析研究编号作者年份试验组ARR/所有对照组ARR/所有1Watson20050/190/192Stallone20032/202/203Baboolal20038/423/454Grinyo200410/444/435Schena20097/5554/2756Johnson200121/2159/2157Lebranchu200917/9510/978Gonwa200222/10018/97资料来源：吴斌等. 中国循证医学杂志，2010，10(1)：33-39. 25 在Stata中输入数据二分类资料输入试验组和对照组发生事件及未发生事件的例数

11、3.1 二分类资料Meta分析合并效应量并生成森林图采用M-H法的固定效应模型合并效应量RR 合并RR=1.645 (95%CI: 1.1902.274)，P=0.003试验组发生ARR的风险是对照组的1.645倍，CNI撤除后TOR-I替代治疗会使发生ARR的风险增加64.5% 27 3.1 二分类资料Meta分析0.0031.645 28 识别发表性偏倚Begg检验结果，P=0.5480.05，不存在发表性偏倚Egger检验结果，P=0.6020.05，不存在发表性偏倚3.1 二分类资料Meta分析0.5480.60295%CI包含0 29 纳入的研究及基本数据3.2 连续型资料Meta分

12、析研究编号作者年份试验组 例数 均数 标准差对照组 例数 均数 标准差1Gossmann20023154.921.52948.318.72Guerin19891234151233103Harpers19964243.912.33639.412.64Kuypers20045350.116.2534215.45Madsen19982749182440166Rump2000244719.6264930.67Van der Schaff19952054202049168Venkat-Raman19992349202347159Wilkie199315641915652110Wilkie19937852

13、77893211Wilkie1994176425176425资料来源：Cross NB, et al. Transplantation 2009; 88: 718 30 在Stata中输入数据对于连续型资料分别输入试验组和对照组的例数、均数和标准差3.2 连续型资料Meta分析 31 确定效应量并进行异质性检验纳入的研究均为RCT，数据为连续型资料，各研究的GFR均值相差不大，故选加权均数差WMD为效应量异质性检验统计量Q=6.48，P=0.7740.10，I 2=0.0% 50%，纳入的研究不存在异质性3.2 连续型资料Meta分析6.480.7440.0%合并效应量并生成森林图采用IV法的

14、固定效应模型合并效应量WMD 合并WMD=4.774 (95%CI: 1.9217.628)，P=0.001试验组比对照组的GFR高4.8个百分点 32 3.2 连续型资料Meta分析4.7740.0011.9217.628 33 合并效应量并生成森林图3.2 连续型资料Meta分析异质性检验I2=0.0% 0.1095%CI不包含0合并效应量WMD=4.7795%CI: 1.927.63识别发表性偏倚Begg检验结果，P=0.0620.05Egger检验结果，P=0.0390.05综合考虑，认为存在发表性偏倚 34 3.2 连续型资料Meta分析0.0620.039 35 识别发表性偏倚3.2 连续型资料Meta分析 36 亚组Meta分析按不同的临床特征、研究质量等分组后进行Meta分析例1中同时对术后2年和4年急性排斥反应做Meta分析