第二十二章META分析实习指导

第二十二章Meta分析［教学要求］了解：Meta分析的偏倚和应用。熟悉：Meta分析的统计方法。掌握：Meta分析的基本含义。［重点难点］第一节Meta分析概述一、Meta分析的含义对同一研究目的或项目的许多研究结果，由于各学者在研究设计、对象选择、样本含量、指标选择、统计方法等方面不完全相同，导致研究结果并不完全一致，对形成最后的综合结论产生重要影响。英国教育心理学家GeorgeGlass给meta分析的定义是"以综合研究结果为目的而对大量单项研究结果进行统计分析”。Meta分析实质上就是汇总相同研究目的的多项研究结果并分析评价其合并效应量的一系列过程，即通过综合多项研究结果而提供一个量化的平均效果或联系来回答问题。二、Meta分析的基本步骤提出问题，制定研究计划可通过查阅和分析大量文献而提出meta分析所研究的问题。Meta分析课题的研究计划包括研究目的、现状、意义、方法、数据收集与分析、结果解释、报告撰写等。检索相关文献一般从研究问题入手，确定相应的检索词及其之间的搭配关系，制定检索策略和检索范围。筛选纳入文献以明确的纳入和剔除标准从检索出的文献中筛选合乎要求的文献。提取纳入文献的数据信息被meta分析采用的数据信息一般包括基本信息、研究特征、结果测量等内容，确定和选择需要分析和评价的效应变量。纳入文献的质量评价主要考察各个研究是否存在偏倚（如选择偏倚、失访偏倚、发表偏倚等）及其影响程度。数据的统计学处理包括：明确资料类型、选择恰当的效应指标；进行同质性检验、选择适合的统计分析模型；效应合并值的参数估计与假设检验；效应合并值参数估计的图示。敏感性分析就不同统计模型的选择、质量较差的文献的剔除前后、文献进行分层的前后和纳入和剔除标准的改变前后进行比较，以了解meta分析结论的稳定性结果的分析与讨论包括异质性及其对效应合并值的影响、几种设计类型的亚组分析、各种偏倚的识别与控制和meta分析结果的实际意义。第二节Meta分析的统计方法Meta分析的统计方法包括固定效应模型和随机效应模型。固定效应模型假设各独立研究来自同一总体的样本，各研究的效应值只是总体参数的一次实现；各研究间的差异只由抽样误差引起，不同研究间的变异很小；一项个体研究内部的方差与另一项个体研究内部的方差相等（方差齐性）。随机效应模型假定各项研究来自不同的总体；研究间的变异较大，既包括各项研究内部的方差，也包括个体研究间的方差；每项研究有其相应的总体参数，meta分析的合并效应值是多个不同总体参数的加权平均。在选择meta分析的统计模型时，首先要对各研究作同质性检验（homogeneitytest），若检验结果不拒绝零假设，即各研究间的差异没有统计学意义，可采用固定效应模型，其主要统计方法包括Mantel-Haenszel法、Peto法和GeneralVariance-Based法；若拒绝零假设，则认为研究间存在异质性，应采用随机效应模型，其统计方法主要是DerSimonianandLaird法。利用meta分析对纳入研究的文献数据进行统计处理时，首先要明确资料的类型及结局变量，然后对待合并的多个研究进行同质性检验，选择适宜的统计分析模型。一、连续型变量资料的meta分析方法对连续型变量资料进行meta分析可选择标准化均数差作为效应变量。根据研究间的齐性情况，选择不同的统计分析模型。若研究间的齐性情况好，则选择固定效应模型；否则，选择随机效应模型，具体统计方法是DerSimonianandLaird法。二、分类变量资料的meta分析方法对于分类变量资料，主要讨论四格表资料的meta分析。能够形成四格表资料的研究方法最常见的有临床随机化试验、病例-对照研究、队列研究和诊断试验评价等。这些研究的数据基本格式见教材表22-4。分析的效应指标可以是OR、RD等。通过齐性检验，若研究间的齐性情况好，则选择固定效应模型，具体统计方法有Mantel-Haenszel法（简称M-H法）、Peto法、Fleiss法以及Generalvariance-based法；否则，选择随机效应模型，具体统计方法是DerSimonianandLaird法。值得注意的是，有时各项研究并不提供四格表的原始资料，而只有RR和95%CI（RRl，RR/，此时可按照教材的公式（22-35）、（22-36）和（22-37）估计方差。第三节Meta分析的偏倚与应用一、发表偏倚的识别与控制在meta分析中可能出现的偏倚主要包括抽样偏倚、选择偏倚和研究内偏倚。抽样偏倚：查找有关文献时产生的偏倚，包括发表偏倚、查找偏倚、索引偏倚、引文偏倚和语种偏倚等。选择偏倚：根据文献的纳入和剔除标准选择符合meta分析的文献时产生的偏倚，包括纳入标准偏倚和选择者偏倚等。研究内偏倚：在资料提取时产生的偏倚。包括提取者偏倚、研究质量评分偏倚和报告偏倚。发表偏倚是meta分析中最常见的偏倚，是指具有统计学意义的研究结果较无统计学意义或无效的结果被报告和发表的可能性更大。识别和控制发表偏倚的方法包括漏斗图分析、线性回归法、秩相关检验法和失安全系数法。二、Meta分析的用途应用范围：病因学研究中因果联系的强度与特异性；各种干预措施效果的程度、特异性以及卫生经济学问题；卫生策略的效果评价。应用目的：提高统计学检验效能；解决单项研究间的矛盾，评价结果的一致性；改善对效应量的估计；解决既往单项研究尚未明确的新问题。

［案例讨论参考答案］见教材表22-8的资料。研究者直接将原始数据进行合并分析是不恰当的，因为每个作者在应用病例对照研究时研究设计、对象选择、样本含量、指标选择、统计方法等方面可能不完全相同，因此应该考虑各研究间可能存在的偏倚和异质性，通过综合分析得到恰当的结论。研究者虽然作了各研究间结果的异质性检验，但未作具体分析即采用了固定效应模型的统计分析方法，从异质性检验结果来看显然是不妥当的，异质性检验Q=26.41，P<0.05，说明3项研究结果间存在异质性，因此应该选择随机效应模型进行分析。根据异质性检验结果，计算过程和结果如下:1）计算各项研究的OR.:3项研究的OR值分别为10.22、1.86和43.5。2）用M-H法估计ORmh值ORMHZEorT.ORMHZEorT.TITTi58.826.62=8.883）计算异质性检验的Q值Q=Z虬Un(OR)-ln(ORJTimhi=0.03+10.40+2.01=12.444)计算校正因子D若QVk-1，则D=0；若QNk-1，则按下式计算D：Q-(k-1)]Z牛

D=i——fybc)2z(b乙〒-yIkTiJ本例Q=12.44Nk-1=3-1=2，应按照上式计算校正因子。_Q-(k-1)'~T112.44-G-1)]x6.62—ifybe¥yfb£TT"TkTfyE]2-yfbkTJi5)用DerSimonian-Laird法计算权重(w.)i=3.046.622-21.131W='D+(1/虬)Ti6)用DerSimonian-Laird法计算合并的OR值2wln(OR))1.70OR=expi-=exp()=8.16DLF|乙wI0.81i7)计算95%CI、expln(OR)土1.96/.JZw^=exp(ln(8.16)土fybc)2z(b乙〒-yIkTiJ本例Q=12.44Nk-1=3-1=2，应按照上式计算校正因子。_Q-(k-1)'~T112.44-G-1)]x6.62—ifybe¥yfb£TT"TkTfyE]2-yfbkTJi5)用DerSimonian-Laird法计算权重(w.)i=3.046.622-21.13与HBV感染之间存在联系。［电脑实验及结果解释］实验22-1连续型变量的固定效应模型程序22.1均数之差的Meta分析（固定效应模型）SAS程序及说明仃号程序说明01DATAa;建立数据集a；定义并输入变量，分别代表每项研究中药物组与02INPUTn1ix1is1in2ix2is2i;对照组的例数、均数和标准差；03spi=SQRT(((n1i-1)*s1i**2+(n2i-1)*s2i**2)/(n1i+n2i-2));计算每项研究的合并方差；04di=(x1i-x2i)/spi;计算每项研究的标准化均数差；05wi=2*(n1i+n2i)/(8+di**2);06wid=wi*di;wid2=wi*di**2;计算各研究的权重；07s_num=5;08CARDS;09241.270.70251.630.70331.290.96301.750.93441.282.41402.211.94351.321.85361.911.80431.221.70411.981.75;定义研究数；15DATAb;SETa;建立数据集b,读入数据集a中数据；

swi1+wi;swidl+wid;swid21+wid2;id=_N_;IF_N_=s_numTHENDO;swi=swi1;swid=swid1;swid2=swid21;END;PROCSORT;BYDESCENDINGid;DATAc;SETb;hbd=swid/swi;df=s_num-1;q=swid2-(swid**2)/swi;p1=1-PROBCHI(q,df);low=hbd-(1.96/(swi**0.5));up=hbd+(1.96/(swi**0.5));chisq=(swid**2)/swi;p=1-PROBCHI(chisq,df);PROCPRINT;VARn1ix1is1in2ix2is2idispiwiwidwid2;VARswiswidswid2hbdqp1lowupchisqp;RUN;运行结果：Output窗口：Obsnlixlislin2ix2is2idi计算效应合并值；计算同质性检验的自由度df;计算同质性检验统计量q；计算q对应的P值；计算效应合并值95%置信区间的下限；计算效应合并值计算效应合并值；计算同质性检验的自由度df;计算同质性检验统计量q；计算q对应的P值；计算效应合并值95%置信区间的下限；计算效应合并值95%置信区间的上限；计算效应合并值假设检验的统计量chisq；计算chisq对应的P值；调用PRINT过程输出结果；指定输出变量；运行程序；spiwiwidwid2swiTOC\o"1-5"\h\z351.321.85361.911.80-0.323321.8248117.5211-5.664931.83159.441.282.41402.211.94-0.422912.1990320.5408-8.686993.67385.331.290.96301.750.93-0.486330.9458615.2977-7.439773.61820.241.270.70251.630.70-0.514290.7000011.8580-6.098383.13631.Obsswidswid2hbdqp1lowupchisqp-36.925216.2416-0.430770.335500.98741-0.64246-0.2190715.9061.0031477................................实验22-2连续型变量的随机效应模型行号0102程序DATAa1;INPUTn1ix1is1in2ix2is2i;行号0102程序DATAa1;INPUTn1ix1is1in2ix2is2i;03spi=SQRT(((n1i-1)*s1i**2+(n2i-1)*s2i**2)/(n1i+n2i-2));建立数据集a1；定义并输入变量，分别代表每项研究中药物组与对照组的例数、均数和标准差；计算每项研究的合并方差；04di=(x1i-x2i)/spi;计算每项研究的标准化均数差；05wi=n1i+n2i;wid=wi*di;计算每项研究权重；06wid2=wi*di**2;07s_num=5;定义研究数；08CARDS;09135.04.7136.53.810304.91.7506.12.3113522.53.42524.910.7122012.51.472012.31.661386.50.7687.381.4114;15DATAa2;SETa1;建立数据集a2，i买入数据集a1中数据；16swi1+wi;swid1+wid;swid21+wid2;17id=_N_;18IF_N_=s_numTHENDO;19swi=swi1;swid=swid1;20swid2=swid21;21END;22PROCSORT;BYDESCENDINGid;23DATAa3;SETa2;24hbd=swid/swi;计算效应的加权均数；25sd2=(swid2/swi)-hbd**2;计算效应的加权方差；26se2=4*s_num/swi*(1+hbd**2/8);计算统计量的标准误；27chisq=s_num*sd2/se2;df=s_num-1;计算同质性检验统计量chisq和自由度df；28p=1-PROBCHI(chisq,df);计算chisq对应的P值；29IFsd2>se2THENsdel=SQRT(sd2-se2);计算随机效应模型的s；30ELSEsdel=0;31low=hbd-1.96*sdel;计算随机效应模型效应合并值95%置信区间下限;32up=hbd+1.96*sdel;计算随机效应模型效应合并值95%置信区间上限;33sdf=se2**0.5/s_num**0.5;计算固定效应模型的，-；d34flow=hbd-1.96*sdf;计算固定效应模型效应合并值95%置信区间下限;35fup=hbd+1.96*sdf;计算固定效应模型效应合并值95%置信区间上限;36PROCPRINT;调用PRINT过程输出结果；37VARn1ix1is1in2ix2is2idispiwiwidwid2;指定输出变量；38VARswiswidswid2hbdlowupchisqpflowfupsd2se2sdf;39RUN;运行结果：Output窗口：运行程序；Obsnlixlislin2ix2is2idispiwiwidwid2swiswidswid2186.50.7687.381.41-0.776951.1326316-12.43129.6585222-81.800946.0898

234520353013134.95.0.47.401.704.70201225245013.301..9010.6.1026.503660.127561.70-0.326147..30-0.572232.80-0.3509845678835879.09707.27376405.102460-19.568480-45.778226-9.12550.65096.3821.26.19563.2028...........Obs1-0.hbdlow36847-0.36847-0up.36847chisqp3.920570.41686-0.63379flow-0.10316fup0.071840sd2se20.0916190.13537sdf................实验22-3分类资料的固定效应模型程序22-3分类资料的Meta分析的固定效应模型SAS程序及说明（Peto法）行号程序说明01DATAa;02INPUTnainbincindi;u匕03ori=(nai*ndi)/(nbi*nci);04m1i=nai+nci;n1i=nai+nbi;05n2i=nci+ndi;m2i=nbi+ndi;06ti=nai+nbi+nci+ndi;07ei=(m1i*n1i)/ti;08s_num=5;建立数据集a；定义并输入每项研究中治疗组和对照组的发生数和未发生数变量；计算每项研究的。日；计算每项研究的行、列合计；计算每项研究的频数合计；计算每项研究中某事件发生数的期望值;定义研究数；09CARDS;49566675574471464707272903227701DATAa;02INPUTnainbincindi;u匕03ori=(nai*ndi)/(nbi*nci);04m1i=nai+nci;n1i=nai+nbi;05n2i=nci+ndi;m2i=nbi+ndi;06ti=nai+nbi+nci+ndi;07ei=(m1i*n1i)/ti;08s_num=5;151617;DATAb;SETa;czi=nai-ei;vi=m1i*m2i*n1i*n2i/((ti**2)*(ti-1));建立数据集上读入数据集3中数据；计算每项研究中某事件实际发生数与期望值之差；计算每项研究中某事件发生数的方差；计算同质性检验统计量q；计算q对应的P值；q=saevi-(sczi**2/svi);p1=1-PROBCHI(ABS(q),s_num-1);END;PROCSORT;BYDESCENDINGid;DATAc;SETb;lnor=sczi/svi;aei=czi/vi;sczi1+czi;svi1+vi;saei1+aei;saevi1+czi**2/vi;id=_N_;IF_N_=s_numTHENDO;sczi=sczi1;svi=svi1;saei=saei1;saevi=saevi1;建立数据集c计算同质性检验统计量q；计算q对应的P值；q=saevi-(sczi**2/svi);p1=1-PROBCHI(ABS(q),s_num-1);END;PROCSORT;BYDESCENDINGid;DATAc;SETb;lnor=sczi/svi;30orp=EXP(lnor);计算合并的OR；p31chisq=sczi**2/svi;计算ORp对应的为2值；32p=1-PROBCHI(chisq,1);计算X2对应的P值；33low=EXP((sczi-1.96*(svi**0.5))/svi);计算OR的95%置信区间下限;p34up=EXP((sczi+1.96*(svi**0.5))/svi);计算OR的95%置信区间上限;p35PROCPRINT;调用PRINT过程输出结果；36VARnainbincindiorieicziaeiviqp1指定输出变量；chisqporplowupsczisaeisvi;37RUN;运行程序；运行结果：Output窗口:Obsnainbincindiorieicziaeiviqp1185725523540.7981491.258-6.2582-0.2312827.05890.597480.9633421027301267240.80287112.780-10.7800-0.2186749.2971..327290322770.8059329.877-2.8770-0.2150413.3790..444714647070.6807653.541-9.5409-0.3800025.1075..549566675570.7197157.579-8.5787-0.3261326.3048..Obschisqporplowupsczisaeisvi110.2492.0013674360.763780.647620.90078-38.0348-1.37112141.1472345实验22-4分类资料的固定效应模型程序22-4分类资料的Meta分析固定效应模型SAS程序及说明(Generalvariance-based法)行号程序说明01DATAa;建立数据集a；02INPUTnainbincindi;定义并输入每项研究中治疗组和对照组的发生数和未发生数变量；03m1i=nai+nci;n1i=nai+nbi;计算每项研究的行、列合计；04n2i=nci+ndi;m2i=nbi+ndi;05ti=nai+nbi+nci+ndi;计算每项研究的频数合计；06cai=nai/n1i;coi=nci/n2i;计算危险度；07rdi=abs(cai-coi);计算危险度的差值；08vrdi=m1i*m2i/(n1i*n2i*ti);计算每项研究的方差；09wi=1/vrdi;wirdi=wi*rdi;计算每项研究的权重；10s_num=5;定义研究数；11CARDS;12182106826498613348127320461014490167834883415157998188876163461340312768

17;18DATAb;建立数据集b；IFN=1THENDO19UNTIL(LAST);20SETaEND=LAST;21swirdi+wirdi;swi+wi;svrdi+vrdi;22END;SETa;读入数据集a中数据；23rdg=swirdi/swi;计算合并的危险度差值；24ri2=rdi-rdg;计算每项研究危险度差值与合并危险度差值之差;25wri2=wi*(ri2**2);26vrdg=1/swi;27DATAc;建立数据集c；IFN=1THENDO28UNTIL(LAST);29SETbEND=LAST;30swri2+wri2;31END;SETb;32DATAd;SETc;建立数据集d,读入数据集c中数据；33q=swri2;计算同质性检验统计量q；34p=1-PROBCHI(q,s_num-1);计算q对应的P值；35low=rdg-1.96/SQRT(swi);计算合并危险度差值的95%置信区间下限；36up=rdg+1.96/SQRT(swi);计算合并危险度差值的95%置信区间上限；37PROCPRINT;调用PRINT过程输出结果；38VARnainbincindirdiwiwirdi;指定输出变量；39VARrdgqplowup;40RUN;运行程序；运行结果：Output窗口:ObsnainbincindirdiwiwirdirdgqPlowup118210682649860.0656004264.07279.7230.0594065.800850.214520.0452910.073521234812732046100.0359323091.13111.0700.0594065.800850.214520.0452910.073521349016783488340.0684014078.10278.9480.0594065.800850.214520.0452910.07352141579981888760.0407614217.86171.9240.0594065.800850.214520.0452910.073521534613403127680.0836693631.09303.8110.0594065.800850.214520.0452910.073521实验22-5分类资料的随机效应模型程序22-5分类资料的Meta分析随机效应模型程序及说明（DerSimonian-LairdV^彳亍号程序说明01OPTIONSNODATENONUMBER;指定输出结果中不包含日期和页码；02DATAr1;建立数据集r1；03INPUTaibicidi;定义并输入四格表变量；ori=ai*di/(bi*ci);ti=ai+bi+ci+di;x=(bi*ci)/ti;x2=x**2;x1=x*ori;CARDS;31234598133912053137951481057304312732855184917;DATAr2;IF_N_=1THENDOUNTIL(LAST);SETr1END=LAST;y+x;y1+x1;y2+x2;END;SETr1;ormh=y1/y;q1=x*(LOG(ori)-LOG(ormh))**2;DATAr3;IF_N_=1THENDOUNTIL(LAST);SETr2(KEEP=q1)END=LAST;q+q1;END;SETr2;k=6;IF(q>=(k-1))THENd=(q-k+1)*y/(y**2-y2);ELSEd=0;wi=1/(d+(1/x));zi=wi*LOG(ori);DATAr4;IF_N_=1THENDOUNTIL(LAST);SETr3(KEEP=wizi)END=LAST;swi+wi;szi+zi;END;SETr3;LABELl='95%CI下限'u='95%CI上限'；ordl=EXP(szi/swi);l=EXP(LOG(ordl)-1.96/SQRT(swi));u=EXP(LOG(ordl)+1.96/SQRT(swi));PROCPRINTDATA=r4;TITLE'表1';RUN;PROCPRINTDATA=r4(FIRSTOBS=1OBS=1)ROUNDLABEL;TITLE'表2';RUN;计算每项研究的。日和频数合计;建立数据集r2；计算M-H法估计的ORmH＜；计算统计量Q；建立数据集r3；定义研究的项数;定义校正因子。；计算权重；建立数据集r4；计算DerSimonian-Laird法的合并OR值；计算合并彼的95%置信区间下限；计算合并彼的95%置信区间上限；调用PRINT过程输出r4中数据；定义标题；调用PRINT过程输出r4中第1行观测值，保留2位有效位，用标签代替变量名；定义标题；运行程序；04050607080910111213141516171819202122232425262728293031323334353637383940414243运行结果：Output窗口：表1Obs1wi0.96949zi1.43704swi8.21582szi7.09181q10.47151q10.1146ai31bi2ci345di98ori4.40290ti476x1.4495820.673641.534888.215827.091811.6340410.1146133912059.761903120.8750030.967491.836068.215827.091811.4044310.114631379516.670891641.4451241.677851.554968.215827.091810.0008610.1146481057302.526321453.9310351.924190.611978.215827.091811.9806010.1146431273281.374431565.6153862.003160.116898.215827.091814.6231810.1146551849171.060091396.34532Obsx2x1yy1y2ormhkdluordl12.10136.3823519.661448.939992.20402.4891360.341621.196484.697322.3707120.76568.5416719.661448.939992.20402.4891360.341621.196484.697322.3707132.08849.6402419.661448.939992.20402.4891360.341621.196484.697322.37071415.45309.9310319.661448.939992.20402.4891360.341621.196484.697322.37071531.53257.7179519.661448.939992.20402.4891360.341621.196484.697322.37071640.26316.7266219.661448.939992.20402.4891360.341621.196484.697322.37071表2Obswiziswisziq1qaibicidioritixx210.971.448.227.090.4710.11312345984.44761.452.1Obsx1yy1y2ormhkd95%CI下限95%CI上限ordl16.3819.6648.9492.22.4960.341.24.72.37［思考与练习的参考答案］一、简答题1、Meta分析是以综合研究结果为目的而对大量单个研究结果进行统计分析的方法，meta分析实质上就是汇总相同研究目的的多个研究结果并分析评价其合并效应量的一系列过程，即通过综合多个研究结果而提供一个量化的平均效果或联系来回答问题。2、Me