Meta分析原理與實(shí)踐

Meta分析實(shí)踐原理vs.首都醫(yī)科大學(xué)陳卉Meta分析與系統(tǒng)評(píng)價(jià)Meta分析基本過程Meta分析實(shí)例一二三2

系統(tǒng)評(píng)價(jià)Meta分析Meta分析與系統(tǒng)評(píng)價(jià)

4

系統(tǒng)評(píng)價(jià)（systematicreview，SR）全面收集針對(duì)某一具體問題的所有相關(guān)研究，并對(duì)這些獨(dú)立研究進(jìn)行嚴(yán)格的質(zhì)量評(píng)價(jià)后，將它們的結(jié)果聯(lián)合起來(lái)進(jìn)行綜合分析（合成）和評(píng)價(jià)5個(gè)階段提出要評(píng)價(jià)的問題檢索和篩選研究評(píng)價(jià)獨(dú)立研究的質(zhì)量綜合分析數(shù)據(jù)和表達(dá)結(jié)果解釋結(jié)果1.1系統(tǒng)評(píng)價(jià)Meta分析足夠的研究數(shù)量滿足一定的相似性

5

Meta分析（meta-analysis）Meta分析也稱為薈萃分析，狹義地講，是指根據(jù)研究的試驗(yàn)設(shè)計(jì)類型以及反映待評(píng)價(jià)問題的統(tǒng)計(jì)學(xué)指標(biāo)，將多個(gè)獨(dú)立研究的指標(biāo)進(jìn)行合成，并對(duì)結(jié)果進(jìn)行假設(shè)檢驗(yàn)通過Meta分析可以增大樣本量，提高檢驗(yàn)效能，增強(qiáng)結(jié)論的可靠性和客觀性1.2Meta分析確定效應(yīng)量異質(zhì)性分析合并效應(yīng)量識(shí)別發(fā)表性偏倚Meta分析基本過程例1為評(píng)價(jià)腎移植術(shù)后CNI撤除后以TOR-I替代治療的療效，研究者收集了8個(gè)隨機(jī)對(duì)照試驗(yàn)，觀察試驗(yàn)組與對(duì)照組（持續(xù)CNI治療）患者術(shù)后1年發(fā)生急性排斥反應(yīng)的情況研究編號(hào)作者年份試驗(yàn)組ARR/所有對(duì)照組ARR/所有1Watson20050/190/192Stallone20032/202/203Baboolal20038/423/454Grinyo200410/444/435Schena20097/5554/2756Johnson200121/2159/2157Lebranchu200917/9510/978Gonwa200222/10018/97資料來(lái)源：吳斌等.中國(guó)循證醫(yī)學(xué)雜志，2010，10(1)：33-39.文獻(xiàn)實(shí)例

7

例2為評(píng)價(jià)降壓類藥物對(duì)腎移植患者的影響，研究者收集了11個(gè)隨機(jī)對(duì)照試驗(yàn)，觀察試驗(yàn)組（CCB）與對(duì)照組（安慰劑）患者服藥2周后的腎小球?yàn)V過率文獻(xiàn)實(shí)例

8

研究編號(hào)作者年份試驗(yàn)組例數(shù)均數(shù)標(biāo)準(zhǔn)差對(duì)照組例數(shù)均數(shù)標(biāo)準(zhǔn)差1Gossmann20023154.921.52948.318.72Guerin19891234151233103Harpers19964243.912.33639.412.64Kuypers20045350.116.2534215.45Madsen19982749182440166Rump2000244719.6264930.67VanderSchaff19952054202049168Venkat-Raman19992349202347159Wilkie199315641915652110Wilkie1993785277893211Wilkie1994176425176425資料來(lái)源：CrossNB,etal.

Transplantation2009;88:7–18

9

效應(yīng)量（effectsize，ES）Meta分析提供給系統(tǒng)評(píng)價(jià)的證據(jù)指標(biāo)確定依據(jù)試驗(yàn)類型隨機(jī)對(duì)照試驗(yàn)（randomcontrolledtrial，RCT）隊(duì)列研究（cohortstudy）病例對(duì)照研究（case-controlstudy）資料類型二分類資料（binarydata）：觀察“是否發(fā)生急性排斥反應(yīng)”連續(xù)型資料（continuousdata）：觀察“腎小球?yàn)V過率”2.1確定效應(yīng)量

10

二分類資料常用效應(yīng)量RD：風(fēng)險(xiǎn)差或率差（riskdifference或ratiodifference）RR：風(fēng)險(xiǎn)比或相對(duì)危險(xiǎn)度（riskratio或relativerisk）OR：比值比（oddsratio）選擇依據(jù)隊(duì)列研究或RCT：首選RR，兩組基線水平一致時(shí)可選RD病例對(duì)照研究：只能計(jì)算OR，試驗(yàn)組和對(duì)照組的事件發(fā)生率均較低（通常低于20%）時(shí)可選RR需要知道每個(gè)研究中試驗(yàn)組與對(duì)照組事件發(fā)生及未發(fā)生的例數(shù)2.1確定效應(yīng)量

11

連續(xù)型資料加權(quán)均數(shù)差（weightedmeandifference，WMD）以研究原有的度量真實(shí)反映了試驗(yàn)效應(yīng)，在實(shí)際分析時(shí)易于理解和使用標(biāo)準(zhǔn)化均數(shù)差（standardizedmeandifference，SMD）度量單位不同和/或均數(shù)差相差較大時(shí)選用沒有度量單位，只能說明統(tǒng)計(jì)學(xué)意義，難以結(jié)合臨床意義進(jìn)行結(jié)果的解釋需要知道每個(gè)獨(dú)立研究試驗(yàn)組和對(duì)照組的樣本例數(shù)、原始均數(shù)及標(biāo)準(zhǔn)差2.1確定效應(yīng)量

12

異質(zhì)性（heterogeneity）指獨(dú)立研究的資料之間存在變異資料間異質(zhì)性的程度決定了Meta分析的方法異質(zhì)性檢驗(yàn)方法計(jì)算Q統(tǒng)計(jì)量并進(jìn)行卡方檢驗(yàn)計(jì)算I2統(tǒng)計(jì)量并直接判斷2.2異質(zhì)性分析

13

Q統(tǒng)計(jì)量檢驗(yàn)法θi是第i個(gè)研究的效應(yīng)量，是所有研究的效應(yīng)量的匯總，wi是第i個(gè)研究的權(quán)重，共有k個(gè)研究Q統(tǒng)計(jì)量服從自由度為k-1的χ2分布若Q統(tǒng)計(jì)量對(duì)應(yīng)的P值>0.10，則認(rèn)為多個(gè)研究間不具有異質(zhì)性，否則認(rèn)為多個(gè)研究之間存在異質(zhì)性2.2異質(zhì)性分析

14

I2（I-squared）統(tǒng)計(jì)量法I2≤25%：異質(zhì)性較小25%<I2≤50%：異質(zhì)性中等I2>50%：異質(zhì)性較大兩種方法的選擇Q統(tǒng)計(jì)量法檢驗(yàn)效能較低，檢驗(yàn)結(jié)果不可靠I2統(tǒng)計(jì)量校正了研究數(shù)量k對(duì)Q值的影響，判斷異質(zhì)性的結(jié)果更為穩(wěn)健可靠2.2異質(zhì)性分析可接受

15

研究間存在較明顯的異質(zhì)性時(shí)的處理方法2.2異質(zhì)性分析異質(zhì)性分析固定效應(yīng)模型異質(zhì)性原因亞組分析Meta回歸隨機(jī)效應(yīng)模型無(wú)有放棄Meta分析無(wú)有

16

合并效應(yīng)量的過程計(jì)算各研究的效應(yīng)量及標(biāo)準(zhǔn)誤為每個(gè)研究賦予權(quán)重對(duì)各研究的效應(yīng)量加權(quán)合并涉及方法的選擇2.3合并效應(yīng)量固定效應(yīng)模型隨機(jī)效應(yīng)模型二分類資料M-H法M-H法IV法IV法Peto法*D-L法連續(xù)型資料IV法IV法

D-L法*注：僅用于計(jì)算OR值M-H法（Mantel-Haenszel）Peto法（M-H的改進(jìn)方法）IV法（InverseVariance）D-L法（DerSimonian-Laird）

17

判斷合并效應(yīng)量的統(tǒng)計(jì)學(xué)意義Z檢驗(yàn)法Z>Z0.05則P<0.05，合并效應(yīng)量有統(tǒng)計(jì)學(xué)意義根據(jù)效應(yīng)量的95%置信區(qū)間判斷對(duì)于WMD、SMD或RD：95%CI不包含0，則P<0.05對(duì)于OR或RR：95%CI不包含1，則P<0.052.3合并效應(yīng)量

18

用森林圖（forestplot）展示分析過程及結(jié)果2.3合并效應(yīng)量異質(zhì)性檢驗(yàn)I2=0.0%<50%P=0.576>0.1095%CI不包含1合并效應(yīng)量RR=1.6495%CI:1.19﹣2.27

19

發(fā)表性偏倚（publicationbias）有統(tǒng)計(jì)學(xué)意義的研究結(jié)果（陽(yáng)性結(jié)果）比無(wú)統(tǒng)計(jì)學(xué)意義的結(jié)果（陰性結(jié)果）更容易被投稿和發(fā)表產(chǎn)生來(lái)源設(shè)計(jì)方案的選擇研究者的主觀愿望資助者的利益識(shí)別方法漏斗圖法Begg秩相關(guān)法Egger線性回歸法2.4識(shí)別發(fā)表性偏倚結(jié)合使用

20

漏斗圖（funnelplot）以每個(gè)研究的效應(yīng)量（或效應(yīng)量的對(duì)數(shù)）為橫軸、樣本量（或效應(yīng)量的標(biāo)準(zhǔn)誤）為縱軸繪制的散點(diǎn)圖如果漏斗圖中的散點(diǎn)表現(xiàn)為明顯不對(duì)稱則認(rèn)為存在發(fā)表性偏倚漏斗圖方法是一種定性的圖示法，研究個(gè)數(shù)較少時(shí)不宜使用2.4識(shí)別發(fā)表性偏倚

21

Begg秩相關(guān)法對(duì)效應(yīng)量標(biāo)準(zhǔn)化后用秩相關(guān)分析來(lái)檢驗(yàn)效應(yīng)量是否與其標(biāo)準(zhǔn)誤存在相關(guān)性假設(shè)檢驗(yàn)得到P>0.05時(shí)可以認(rèn)為不存在發(fā)表性偏倚當(dāng)納入的獨(dú)立研究的數(shù)量較少時(shí)，檢驗(yàn)效能相對(duì)較低以每個(gè)研究的效應(yīng)量（或效應(yīng)量的對(duì)數(shù)）為縱軸、效應(yīng)量的標(biāo)準(zhǔn)誤為橫軸繪制的散點(diǎn)圖稱為Begg漏斗圖2.4識(shí)別發(fā)表性偏倚

22

Egger直線回歸法擬合標(biāo)準(zhǔn)正態(tài)離差對(duì)效應(yīng)估計(jì)值的精確度（precision）的線性方程直線的截距代表不對(duì)稱的程度，若假設(shè)檢驗(yàn)P>0.05則可以認(rèn)為不存在發(fā)表性偏倚以每個(gè)研究的精確度為橫軸、效應(yīng)量的標(biāo)準(zhǔn)誤為縱軸繪制的散點(diǎn)圖稱為Egger圖，

截距的95%置信區(qū)間包含0表示不存在發(fā)表性偏倚2.4識(shí)別發(fā)表性偏倚二分類資料Meta分析連續(xù)型資料Meta分析Meta分析實(shí)例

24

納入的研究及基本數(shù)據(jù)3.1二分類資料Meta分析研究編號(hào)作者年份試驗(yàn)組ARR/所有對(duì)照組ARR/所有1Watson20050/190/192Stallone20032/202/203Baboolal20038/423/454Grinyo200410/444/435Schena20097/5554/2756Johnson200121/2159/2157Lebranchu200917/9510/978Gonwa200222/10018/97資料來(lái)源：吳斌等.中國(guó)循證醫(yī)學(xué)雜志，2010，10(1)：33-39.

25

在Stata中輸入數(shù)據(jù)二分類資料輸入試驗(yàn)組和對(duì)照組發(fā)生事件及未發(fā)生事件的例數(shù)3.1二分類資料Meta分析

26

確定效應(yīng)量并進(jìn)行異質(zhì)性檢驗(yàn)納入的研究均為RCT，數(shù)據(jù)為二分類資料，故選用RR異質(zhì)性檢驗(yàn)統(tǒng)計(jì)量Q=4.76，P=0.576>0.10，I2=0.0%<50%，可以認(rèn)為納入的研究不存在異質(zhì)性3.1二分類資料Meta分析合并效應(yīng)量并生成森林圖采用M-H法的固定效應(yīng)模型合并效應(yīng)量RR合并RR=1.645

(95%CI:1.190﹣2.2