meta分析专题知识讲座

1、meta分析专题知识讲座meta分析专题知识讲座第1页二、检索资料资料类型：类数据：公开发表文件或汇报中方法学与结果个别数据。类数据：公开发表文件或汇报中表格数据。类数据：图形信息。类数据：工作汇报和上报资料。meta分析专题知识讲座第2页资料检索策略：1、先进行预检索，大致确定检索范围2、检索时进行必要限定，如：研究对象、语种、出版年限、出版类型等3、确保较高查全率4、计算机检索与手工检索相结合，并重视参考文件。5、注意经过其它渠道搜集，如：会议论文、专著内章节等。meta分析专题知识讲座第3页三、拟订纳入标准和排除标准纳入标准：如患者疾病类型、年纪、性别、种族、病情严重程度。描述性研究还是

2、分析性研究？病例对照研究还是队列研究？采取随机、盲法、干预办法等？剔除标准：如不一样杂志或媒体上发表同一结果不能重复选择。纳入标准和排除标准不宜过严或过宽meta分析专题知识讲座第4页四：纳入研究质量评价评价内在真实性、外在真实性和可能影响结果解释原因。内在真实性：考查独立研究是否存在偏倚及其影响程度，包含患者是否随机分配、试验组对照组可比性、是否采取盲法判断。当前真实性评价尚不完善1993年Guyatt提出了一个简便有效评价临床研究资料真实性和用途方法meta分析专题知识讲座第5页真实性评价表研究结果正确吗？（1）研究对象是否随机分配到治疗组和对照组？（2）进入试验患者对结论有贡献吗？ 随访

3、完成情况 不依从者仍在原来组中分析吗？（3）使用盲法了吗？（患者、医生、研究者）（4）试验开始时组间可比吗？（5）除了干预方法，对各组处理相同吗？ 结果是多少？meta分析专题知识讲座第6页（1）疗效是多少？依据资料性质计算疗效大小计数资料设对照组率为X，对照组率为Y则绝对危险降低许为X-Y相对危险降低值为（X-Y ）/X计量资料设对照组试验前后某指标改变量为试验组试验前后某指标改变量为Y2- Y1则绝对改变量为（X2-X1 ）-（ Y2- Y1）（2）疗效精度假设检验P值大小或95%可信限meta分析专题知识讲座第7页这个结果对治疗我患者有帮助吗？（1）结果能够用于我患者吗？（2）全部主要临

4、床结局均被考虑了吗？（3）治疗收益与潜在危险和费用比较meta分析专题知识讲座第8页五 、提取纳入文件信息包含基础信息、研究特征、结果测量等。所取信息必须可靠。对缺乏原始数据文件，在与原著者联络后依然没有数据，应该删除。六、资料统计学处理制订统计分析方案选择适当效应量指标：OR、RR、P值等纳入研究异质性检验模型选择及统计学分析meta分析专题知识讲座第9页七 敏感性分析采取不一样模型或排除某项研究，以考查Meta分析结论稳定性。八、形成结果汇报背景、目标、检索策略、选择标准、资料搜集方法、主要结果和评价者结论meta分析专题知识讲座第10页Meta分析统计学过程一、效应量确实定常见有率差、均

5、数之差、标准化差值（standardized mean difference, SMD）相对危险度、优势比、绝对危险降低率（ARR）、需要治疗患者数（NNT）。当缺乏上述效应量时，可酌情考虑用假设检验P值，或检验统计量作为效应量。meta分析专题知识讲座第11页1、分类资料效应量（1）病例对照研究效应量： OR=ad/bc对数方差：Var（ln oR）=1/a+1/b+1/c+1/d95%CI：expln oR1.96 Var（ln oR）（2）随访研究资料效应量1）相对危险度（RR）、归因危险度（AR）、率差（RD）Var（ln RR ）=（b/an1）+（d/cn2）meta分析专题知识讲

6、座第12页95%CI：expln RR1.96 Var（ln RR）2）归因危险度（ARR）Var（ ARR ）=（ab/n31）+（cd/n32）95%CI：ARR1.96 Var（ RR）ARR倒数称为需要治疗患者数（number need to treat NNT），即需要治疗多少名患者才能预防一次不利结局出现，是临床决议常见指标meta分析专题知识讲座第13页2、定量资料标准化均数之差组别 例数 治疗前后差值 标准差治疗组 n1 X1 s1 对照组 n2 X2 s2 累计 n1+ n2 x s 均数标准化差值定义 D=（ X1- X2）/ s合并标准差（s）=（ n1- 1） s21

7、+ （n2 -1） s22 / （n1+ n2 -2）meta分析专题知识讲座第14页Var（ D ）=（ n1+ n2 / n1 n2 ）+ D2 /2 n1+2 n2 ）95%CI：D1.96 Var（ D）3、其它效应量假如纳入Meta分析原始研究没有给出结局变量，仅汇报了假设检验统计量或P值大小，此时可将这些统计量转化为效应尺度。如连续性资料给出了统计量Z，可将Z转化为效应尺度meta分析专题知识讲座第15页 效应尺度 = Z1/n1+ 1/n2 4、效应量描述常使用森林图（forest plot） ，森林图为多个原始文件效应量及其可信区间绘成，纵坐标为效应尺度，横坐标为原始文件编号，

8、按照一定次序，将各研究效应量及其可信区间绘制到图上。森林图可用于描述每个原始研究效应量分布及其特征，同时展开研究间结果差异情况。meta分析专题知识讲座第16页二、异质性检验基础思想假设研究资料间效应量一致，全部研究都来自同一总体，那么实际效应间差异可认为是由抽样误差造成，但若效应间差异过大，即这些变异不但仅是抽样引发，纳入Meta分析各个独立研究可能来自不一样总体。异质性检验常见加权检验或Peto法检验。meta分析专题知识讲座第17页加权检验基础步骤、建立检验水准、计算统计量值（）或（） 为第个研究权重，可用其方差倒数作为权重；为第个研究效应量，为全部研究平均效应量： 且服从于自由度为卡方

9、分布。meta分析专题知识讲座第18页例1：某学者欲研究急性心急梗塞患者心肺复苏后溶栓治疗效果情况，搜集了五个随机对照试验结果溶栓治疗组对照组 累计研究编号存活死亡存活死亡 20 13 272 268 573 2 3 4 7 13 27 3 55 12 133 103 303 4 31 11 63 5 2 157 5 83 49 47 86 265 累计 192 89 522 86 265meta分析专题知识讲座第19页权重法Q检验统计量计算过程研究 OR lnoR VarlnoR lnoR （lnoR）2 1 1.52 0.42 0.13 7.44 3.10 1.29 2 1.39 0.33

10、 0.80 1.25 0.41 0.14 3 3.55 1.27 0.12 8.42 10.67 13.51 4 2.33 0.84 0.16 6.32 5.33 4.50 5 3.1 1.13 0.07 15.30 17.31 19.58累计 38.73 36.82 39.02依据上述公式以lnoR 为效应量，以其方差为权重，计算Q检验统计量。Q=（） =39.02-（36.8）2/38.73=4.02自由度=5-1=4，P0.05, 异质性检验无统计学意义，研究效应间满足同质性。meta分析专题知识讲座第20页三、合并效应量假设检验合并效应量实际上多个研究效应量加权平均值，普通可分两步进行

11、预计，首先逐一计算每个研究效应量及其95%可信区间；然后依据资料类型与一致性检验结果，选择适当统计模型，预计合并效应量，必要时做假设检验。资料同质：固定效应模型（Fixed effect model）资料异质：随机效应模型（random effect model ）meta分析专题知识讲座第21页1、计数资料:常见四格表资料，主要起源于临床试验、病例对照和队列研究。效应量合并采取Peto法，也称改进Mantel-Haenszel法，分析步骤：（1）整理资料：每个研究对应一个四格表（2）计算第i个研究效应量OR值（3）计算异质性检验统计量Q（ Peto法）Q= Ei为第i个研究理论频数，Vi为第

12、i个研究方差Ei=n1im1i/Ti Vi=（n1im1in2im2i）/Ti2（Ti-1）meta分析专题知识讲座第22页（4）计算合并OR合并及其 95 % CI：OR合并=exp（ai-Ei）/ Vi95 % CI： Exp（ ln oR合并1.96/ Vi ）（5） OR合并假设检验：X2=（ai-Ei）2/ Vimeta分析专题知识讲座第23页以例1资料说明固定效应模型（ Peto法）预计过程及结果研究 ai Ei Vi （ai- Ei） （ai- Ei）2 oR 1 20 16.82 7.79 3.18 1.30 1.52 2 3 2.59 1.26 0.41 0.13 1.39

13、3 55 41.57 12.33 13.43 14.63 3.55 4 31 25.15 7.44 5.85 4.61 2.33 5 83 64.75 16.62 18.25 20.03 3.10 累计 192.0 150.88 45.43 41.12 40.70 2.47meta分析专题知识讲座第24页依据上表计算结果异质性Q检验Q=3.48 P0.05,异质性检验无统计学意义，说明研究间效应量满足同质性。OR合并= （ ai- Ei） / Vi= exp（0.9051）=2.47；95 % CI：1.853.33OR合并假设检验：X2合并=37.2，P0.05 采取固定效应模型。合并效应尺

14、度d合并=（ -25.22）/63.45=-0.397，其中95%CI为-0.64-0.14meta分析专题知识讲座第28页随机效应模型随机效应模型由Dersimonian和Larid于1986年提出，故称DL模型，其关键是对每个研究权重进行校正，即以研究内方差和研究间方差之和倒数作为权重纳入分析。1、计数资料（1）整理资料（2）计算第个研究效应量OR OR=ad/bcmeta分析专题知识讲座第29页（3）计算MH法比值比OR MH = MHOR/ MHMH加权系数MH=bc/T（4）进行异质性检验并计算D：= MHln（oR）-ln（ORMH）2 0 k-1D= -(k-1) MH/(MH)

15、2- (2MH) 若k-1meta分析专题知识讲座第30页(5)计算DL加权系数和方差SDL2=1/（D+1/MH） SDL2=1/ （6）计算DL法对数比值比ORDL及 Ln（ORDL）= Ln（OR）/ 95%CI= ORDLe1.96S2DL（6）合并效应量假设检验X2=Ln（ORDL）2/SDL2meta分析专题知识讲座第31页例3 ：5个随机对照试验研究结果 治疗组 对照组 研究 有 无 有 无 oR MH （oRMH） MH（ln oR-lnORMH）2 1 9 14 17 9 0.34 4.86 1.65 1.12 2 11 49 27 34 0.28 10.93 3.09 4.

16、84 3 6 21 22 13 0.17 7.45 1.26 1.39 4 79 60 74 60 1.07 16.26 17.36 7.15 5 8 28 16 17 0.30 6.49 1.97 2.29累计 46 25.34 25.80本例k=5，=25.80k-1，D=0.62， SDL2=0.15Ln（ORDL）= Ln（OR）/ = -1.03 ORDL=0.3695%CI为（0.17 ，0.77）meta分析专题知识讲座第32页2、计量资料（1）计算效应尺度d和合并标准差：同固定效应模式（2） d加权均数和加权方差预计值 d= d/ S2d=（ d2/ ）- d2式中权重=Ni

17、为各个研究样本含量之和。（3）计算异质校正因子S2e=（4k/ ）（1+ d2/8）（4）标准误计算当S2d S2e时，则校正标准误为 S= S2d- S2e当S2d S2e时，则校正标准误为 S = Se /kmeta分析专题知识讲座第33页（5）异质性检验 X2=k S2d/ S2e（6）计算平均效应尺度95%CI d1.96 S meta分析专题知识讲座第34页例4 慢性乙肝患者血清SIL-2R含量Meta分析结果 病例组 对照组 研究 n1i x1i s1i n2i x2i S2i sci di wi di wi d2i wi 1 40 232 63.2 20 70.5 32.4 55

18、.04 2.93 60 175.8 515.09 2 28 519 176 31 238 205 191.81 1.46 59 86.14 125.76 3 30 500.5 391.9 30 384 413.9 403.05 0.29 60 17.40 5.05 4 57 402.8 136.1 25 247.2 413.9 253.69 0.61 82 50.02 30.51累计 261 329.36 676.41 由上表可知d=329.36/261=1.262 S2d=676.41/261-1.262=0.999 S2e=0.074 异质性检验X2=k S2d/ S2e=54.00 P0.05 故采取随机效应模型因S2dS2e 所以S=0.962 平均效应尺度d95%CI为d 1.96 S 为（-0.62，3.15）meta分析专题知识讲座第35页Meta分析偏倚及其控制一 偏倚种类1、选择偏倚（1）发表偏倚：没有统计学意义或无效结论不发表文件（2）数据标引不