系統(tǒng)評價中的常用統(tǒng)計分析方法(2013.1.15)

1、1 1系統(tǒng)評價中的常用統(tǒng)計分析方法及系統(tǒng)評價中的常用統(tǒng)計分析方法及statastata軟件軟件2 2內容內容1、meta分析的用途2、meta分析的基本步驟3、meta分析中常用的統(tǒng)計指標4、meta分析中常用統(tǒng)計方法5、meta分析舉例及結果解釋6、meta分析結果的穩(wěn)定性和可靠性分析7、與stata有關的meta分析文獻及書籍介紹3 3MetaMeta分析簡介分析簡介 Meta分析（meta analysis）由beecher在1955年最先提出，并由Glass在1977年首次命名。目前，meta分析已經(jīng)在教育學、心理學、社會科學等領域得到越來越多廣泛的應用，在循證醫(yī)學和循證衛(wèi)生管理中也發(fā)

2、揮越來越重要的作用。4 4MetaMeta分析的用途分析的用途 在世界范圍內，對同一研究目的的研究可能有幾個，幾十個，甚至上百個學者在不同地區(qū)、不同年代進行研究并報告結果，但各學者在研究設計、對象選擇、樣本含量、指標選擇、統(tǒng)計方法等方面不完全相同，導致研究結果并不完全一致，對這些結果進行綜合評價和取舍是十分困難的，而meta分析正是對這些結果進行定量綜合定量綜合的適宜統(tǒng)計方法。5 5MetaMeta分析的用途分析的用途 通過meta分析可以達到以下目的：1、提高統(tǒng)計學檢驗效能 在進行假設檢驗時，能否得到“有統(tǒng)計學意義”的結果與樣本含量存在一定的關系，meta分析綜合多個同類研究結果，總樣本量增

3、大，因而可提高統(tǒng)計檢驗的效能。6 6MetaMeta分析的用途分析的用途 2、解決單個研究間的矛盾，評價結果的一致性 由于各種研究在設計、對象的選擇、樣本含量、試驗條件等方面不同，其研究結果的質量存在很大差異，一般綜述方法很難對研究結果進行取舍，而meta分析可以估計各個研究可能存在的偏倚以及異質性的來源，采用統(tǒng)計方法對各個研究結果進行定量綜合評價。7 7MetaMeta分析的用途分析的用途3、改善對效應量的估計 多個同類研究的結果可能在程度和方向上存在差異，有時甚至得到相互矛盾的研究結果。用meta分析的綜合研究結果可以估計各個研究效應量的平均水平，從而可得到一個供選擇的明確結論。8 8Me

4、taMeta分析的用途分析的用途4、解決既往單個研究未明確的新問題 meta分析可以探索單個研究未闡明的某些問題，發(fā)現(xiàn)既往研究存在的缺陷，繼而提出新的研究問題的思路。9 9MetaMeta分析的臨床應用分析的臨床應用l干預措施的評價 估計預防、治療、康復的效果和風險l診斷性試驗 提供更為可靠的診斷性試驗精確性計算l流行病學 提供更為可靠的病因/危險因素估計10101111MetaMeta分析的基本步驟分析的基本步驟1、提出問題，制定研究計劃 問題的提出需要系統(tǒng)復習大量文獻，meta分析所研究的問題一般可來自生物醫(yī)學研究領域中不確定或有爭議的問題。與其他研究一樣，meta分析課題的研究計劃包括研

5、究目的、現(xiàn)狀、意義、方法、數(shù)據(jù)收集與分析、結果解釋、報告撰寫等。1212MetaMeta分析的基本步驟分析的基本步驟2、檢索相關文獻 一般從欲研究的問題出手，確定相應的檢索詞及其之間的搭配關系，制定檢索策略和檢索范圍。對檢索結果要進行查全、查準與否的評價，這是至關重要的，否則會影響meta分析結論的可靠性和真實性。1313MetaMeta分析的基本步驟分析的基本步驟3、篩選納入文獻 以明確的納入和排除標準從檢索出的文獻中篩選合乎要求的文獻。在制定文獻納入和排除標準時，要考慮研究對象、設計類型、處理因素、結局效應、樣本大小、觀察年限、文獻發(fā)表時間和語種等方面的問題。1414MetaMeta分析的

6、基本步驟分析的基本步驟4、提取納入文獻的數(shù)據(jù)信息 被meta分析采用的數(shù)據(jù)信息一般包括基本信息、研究特征、結果測量等內容、確定和選擇需要分析和評價的效應變量。1515MetaMeta分析的基本步驟分析的基本步驟5、納入文獻的質量評價 主要考察各個研究的質量間是否存在差異。質量高低可用權重來表示，也可用量表或評分系統(tǒng)來評價。1616MetaMeta分析的基本步驟分析的基本步驟6、數(shù)據(jù)的統(tǒng)計學處理 主要包括：明確資料類型、選擇恰當?shù)男笜耍贿M行同質性檢驗、選擇適合的統(tǒng)計分析模型；效應合并的參數(shù)估計與假設檢驗；效應合并值參數(shù)估計的圖示。1717MetaMeta分析的基本步驟分析的基本步驟7、敏感性

7、分析 目的是了解分析結論的穩(wěn)健性。主要通過以下方式來考察meta分析的結論是否因為采取不同的分析措施而發(fā)生根本性的變化。（1）選擇不同統(tǒng)計模型時，效應合并值點估計與區(qū)間估計的差異；（2）排除質量較差的文獻前后，結論的差異；（3）對文獻進行分層前后，結論的差異；（4）改變納入和排除標準，結論的差異。1818MetaMeta分析的基本步驟分析的基本步驟8、結果的分析與討論 包括異質性及其對效應合并值的影響；幾種設計類型的亞組分析；各種差異的識別與控制；meta分析結果的實際意義。1919常用數(shù)據(jù)類型常用數(shù)據(jù)類型計量（定量）數(shù)據(jù)，如身高，體重，肺活量等 二分類，如生存/死亡、有效/無效等分類數(shù)據(jù) 有

8、序，如-,+,+,+等 多分類 無序，如血型，職業(yè)等2020MeatMeat分析中常用的效應指標分析中常用的效應指標計數(shù)資料，常用的效應測量指標為：危險度（risk），又稱率。危險比率（risk ratio,RR），又稱為相對危險度（relative risk，RR）比值比（odds ratio,OR）危險差（risk difference，RD） 前兩個為相對效應指標，即試驗治療與對照治療何者為優(yōu)；后一指標為絕對效應。相對危險減少率（relative risk reduction，RRR）需處理病人數(shù)（number needed to treat，NNT）2121危險度（危險度（riskri

9、sk）危險度 24人沿一斜坡滑雪，6人滑到，則滑倒的危險度為：6人滑到/24人=0.25；Risk=發(fā)生某一事件的人數(shù)/觀察的總人數(shù)2222比值比值(odds)(odds)24人沿一斜坡滑雪，6人滑到，則滑倒的比值為：6人滑到/18人未滑倒=0.33Odds=發(fā)生某一事件的人數(shù)/未發(fā)生某事件的人數(shù)2323組間比較組間比較groupCommon cold合計YesNo治療組121（a）160(b)281對照組160(c)121(d)281合計281281562Kerugol 2006:incident of common cold2424危險度與比值的比較危險度與比值的比較危險度與比值的差異-K

10、erugol試驗的對照組中，281名人中有160人發(fā)病，則發(fā)病的概率為：-Risk=160/281=0.57 odds=160/121=1.32-Kerugol試驗的干預組中，281名人中有121人發(fā)病，則發(fā)病的概率為：-Risk=121/281=0.43 odds=121/160=0.76groupCommon cold合計YesNo治療組121160281對照組160121281合計2812815622525危險比率（危險比率（risk ratio,RRrisk ratio,RR）= =相對危險度（相對危險度（relative risk,RRrelative risk,RR）干預組的危險度

11、為=121/281=0.43對照組的危險度為=160/281=0.57groupCommon cold合計YesNo治療組121160281對照組160121281合計281281562危險比率=121/281=0.43=0.76 160/281 0.57RR=1，表明干預組與對照組無差異2626比值比（比值比（odds ratio,ORodds ratio,OR）干預組的比值比為=121/160=0.76對照組的比值比為=160/121=1.32groupCommon cold合計YesNo治療組121160281對照組160121281合計281281562比值比=121/160=0.76

12、=0.58 160/121 1.32OR=1，表明干預組與對照組無差異2727效應指標的含義效應指標的含義率，即危險度，用來說明某事發(fā)生可能性大小的指標；RR:相對危險度，隊列研究中常用指標。說明暴露組的發(fā)生率是非暴露組的多少倍。RR越大，表明暴露與結局關聯(lián)的強度越大。相對危險度與關聯(lián)的強度見下表2828OROROR即比值比，也叫優(yōu)勢比。常用在病例-對照研究中。OR 值的意義與相對危險度相似，即暴露組的危險性為非暴露組的多少倍！在慢性病中，或疾率小于5%時，ORRR2929RDRDRD，即（rate difference）,又稱特異危險度、歸因危險度或超額危險度，或絕對危險度減少率（ARR）。

13、意義：對照組發(fā)病率與暴露組發(fā)病率相差的絕對值！如上例，治療組的發(fā)病率為0.43，對照組的發(fā)病率為0.57，則RD=0.14，表示治療組的發(fā)病率比對照組低14%。3030RRRRRRRRR，相對危險減少率（relative risk reduction），反應實驗組與對照組某病發(fā)生率的相對變化量。如上例，治療組的發(fā)病率為0.43，對照組的發(fā)病率為0.57，則RRR=1-0.76=0.24，表示治療可減少24%的發(fā)病率，或由于治療，使得24%的免于發(fā)病。3131NNTNNTNNT，需處理病人數(shù)（number needed to treat）,意義為：在一定的時間內，對患者采用某種防治措施后，得到一

14、例有利結果需要防治的病例數(shù)。如上例，治療組的發(fā)病率為0.43，對照組的發(fā)病率為0.57，則ARR=0.14，NNT=1/0.14=7.14,含義為:只需治療7人，即可治愈好一人。可見，NNT越小，其防治效果越好，臨床意義越大。3232四、Meta 分析中的統(tǒng)計相關內容NNH . ：對患者采用某種防治措施處理，得到一：對患者采用某種防治措施處理，得到一例副作用需要處理的病例數(shù)例副作用需要處理的病例數(shù)（the nubmer of needed to the nubmer of needed to harm one more patients form the therapy, harm one m

15、ore patients form the therapy, NNHNNH） .計算公式：計算公式： NNH=1/ NNH=1/ | |P1P1g g P2P2g g| | =1/ARI =1/ARI .由公式可見：由公式可見：NNHNNH值越小，該防治措施引起的副作用值越小，該防治措施引起的副作用越大，其值越大越好。越大，其值越大越好。 .例如：某治療措施引起的副反應發(fā)生率例如：某治療措施引起的副反應發(fā)生率64%64%，而在對，而在對照組中出現(xiàn)類似的副反應發(fā)生率為照組中出現(xiàn)類似的副反應發(fā)生率為37%37%，計算其，計算其ARI= ARI= |37%|37%64%|=27% 64%|=27%

16、， NNH=1/27%=4NNH=1/27%=4。即該治療措施每處理。即該治療措施每處理4 4個病例，就會出現(xiàn)個病例，就會出現(xiàn)1 1例副反應。例副反應。3333MeatMeat分析中常用的效應指標分析中常用的效應指標計量資料計量資料，常用的效應測量指標為： 權重的均數(shù)差(weighted mean difference, WMD)當采用同樣測量方法測量同一個指標時； 標準化的均數(shù)差(standardisedmean difference, SMD): 對同一治療效應采用不同的測量方法或單位，如測定疼痛采用不同的量表測量。均數(shù)差的計算均數(shù)差的計算Mean difference干預組均數(shù)-對照組均

17、數(shù)3434權重的均數(shù)差權重的均數(shù)差（weighted Mean differenceweighted Mean difference，WMDWMD）什么情況下使用：每一試驗中以標準方式對連續(xù)變量進行測量時使用 Hbmg/dl;BPmmHg以均數(shù)和標準差表示匯總的結果以自然單位表示，如mmHg臨床醫(yī)生容易理解權重的均數(shù)差(WMD)=0，表明治療組與對照組無差異！3535標準化的均數(shù)差標準化的均數(shù)差standard Mean differencestandard Mean difference（SMDSMD）什么情況下使用：每一試驗中以不同的測量單位對同一結局指標進行測量時使用如如痛的的重程度，采

18、用不同量表抑程度的測量：0-20，0-100不可能進行直接比較以均值和標準差表示(mean value)(standard deviation, SD)匯總的結果不標注單位3636MeatMeat分析中常用的效應指標分析中常用的效應指標可信區(qū)間，可信區(qū)間，CI，confidence inteval使用效應量來估計兩組間效應值的差異，如OR=0.736；從統(tǒng)計學看，這種效應差異可能產(chǎn)生于機遇的作用，如何保證這種差異是由于治療效應的作用而不是隨機作用？衡量的指標為可信區(qū)間可信區(qū)間。3737可信區(qū)間可信區(qū)間CICI一項研究提供效應大小的估計值可信區(qū)間是指真實值的測量值所分布的范圍95%CI的含義對測

19、量值的分布范圍，人們有95%的信心認為樣本（人群）測量真值是在此范圍。如果可信區(qū)間較窄，說明估計的效應值更精確；如果可信區(qū)間較寬，說明效應的估計值的精確度較差。38383939權重（加權）權重（加權）weighweighl在meta-分析中, , 簡單地進行效應值的均數(shù)計算會造成誤導l每一個研究被賦予一定的權重 大樣本研究由于效應估計值更為精確( (更可靠) )，在匯總分析中具有較大的權重 小樣本研究由于估計值不夠精確( (欠可靠) )，賦予的權重較小Meta分析的合并效應量為各獨立研究的加權平均數(shù)加權平均數(shù)4040權重的計算權重的計算各研究對合并效應量的貢獻大小不等，根據(jù)事件發(fā)生率和樣本含量

20、分配權重，然后根據(jù)各研究的權重計算出合并效應量平均權重的計算公式為：Ti：研究i的療效估計值Wi：該研究在全部研究中所分配的權重。4141二分類變量二分類變量metameta分析常用的統(tǒng)計方法分析常用的統(tǒng)計方法倒方差法（woolfs method）M-H法（Mantel-Haenszel）Peto法4242二分類變量效應量方差的計算公式二分類變量效應量方差的計算公式4343倒方差法倒方差法基本思想：根據(jù)標準誤（SE）計算方差的倒數(shù)而獲得各研究的權重。通過此法，具有較小標準誤的大規(guī)模研究獲得的權重大于具有較大標準誤而規(guī)模較小的研究，這種分配權重的方法減少了合并效應量的不確定性。使用固定效應模型分

21、配權重。平均權重計算公式為：Ti：研究i的療效估計值Si：該估計值的標準誤。4444M-HM-H法法權重中反應了標準誤的變化，適合小樣本含量和事件發(fā)生較多的研究。4545PetoPeto法法適用于大型研究的小效應量的合并分析，包括生存資料、實際事件數(shù)（O）和理論數(shù)（E）差值的比較。步驟為：1、計算各研究的實際數(shù)（O，Observed enents）和理論數(shù)（E,expected enents ）。2、計算各研究（O-E）及其總和。3、計算各研究的精確方差V及其總和（V：研究權重）；4、比較實際差值與其方差：4646倒方差法、倒方差法、M-HM-H法和法和PetoPeto法權重差異法權重差異三種

22、方法由于在分配權重方法上的差異，各研究的權重可出現(xiàn)很大的差異，但對合并效應量的影響可能不會太大。OR、RR、RD的選擇一般原則如下：小事件選擇Peto OR或M-H OR，大事件選擇RR。但由于OR和RR是相對效應量，往往不能反映實際效應量的大小。如兩種不同藥物做了兩個治療的研究，每個研究的治療組和對照組各100人，研究1治療組治愈40人，對照組治愈20人，則OR=2；研究2的治療組治愈數(shù)為20例，對照組為10例，RR=2，僅從RR看療效相同。但使用RD來表示時，前者的RD為20%，后者為10%，療效的優(yōu)劣一目了然。4747連續(xù)性變量連續(xù)性變量metameta分析常用的統(tǒng)計方法分析常用的統(tǒng)計方

23、法1、均數(shù)差：MD=2、權重的計算：與倒方差法權重計算公式同，即3、按照公式合并效應量。iixx21加權均數(shù)差加權均數(shù)差與標準化均數(shù)差標準化均數(shù)差加權均數(shù)差標準化均數(shù)差1、均數(shù)差的計算和權重的計算，同上。2、均數(shù)差標準化：SMD=MD/SD，SD的計算見下。3、按照公式合并效應量，公式同上。iiiiwseww1)(4848SMDSMD的計算的計算Standardized mean difference solution：express in standard units:effect size(ES) SMD=MD/SD=ESgroup樣本量均數(shù)標準差干預組n1S1對照組n2S2221*) 1

24、2(*) 11(2221nnSnSnSD1x2x1S4949MetaMeta分析中的異質性（分析中的異質性（heterogeneityheterogeneity）系統(tǒng)綜述中納入的研究會有較大不同，研究中任何變異的情況都可被稱之為異質性。異質性總是存在的，但鑒定研究間的異質性十分重要。5050異質性的類型異質性的類型臨床異質性(概念上的異質性)：如對象特征、診斷、干預、對照、研究地點、評價結局等不同方法學異質性：研究設計與質量不同統(tǒng)計學上的異質性：不同試驗中觀察得到的效應，其變異性超過了機遇( (隨機誤差) )本身所致的變異性5151臨床異質性產(chǎn)生的原因臨床異質性產(chǎn)生的原因對象：研究間納入/ /

25、排除標準差異研究設計：隨機、盲法、樣本大小研究場所：入院偏倚治療措施：變異性，辨證論治，復雜干預評價結局：結局定義與測量不同平均效應值是否具有臨床意義?5252如何檢驗異質性如何檢驗異質性圖形審查圖形審查 如果每個研究結果的可信區(qū)間互不重疊，則研究間的差異很可能具有統(tǒng)計學上的顯著性。對同質性檢驗的統(tǒng)計學檢驗對同質性檢驗的統(tǒng)計學檢驗1 1、異質性的定性分析、異質性的定性分析 采用卡方檢驗和采用卡方檢驗和P P值。值?？ǚ街翟贑ochrane中又稱Q值。Q值相對于自由度（文獻數(shù)-1）越大，P值越小，則存在異質性的可能性越大；反之越小。該法不是很可靠的檢驗尤其是對于小樣本試驗或試驗數(shù)有限。2 2、異

26、質性的定量分析、異質性的定量分析I2 =(Q-df)/Q100%, Q=卡方值，df=n-1; I2 50% 較明顯的異質性， I2 75% 高度異質性25353如何處理研究間存在的顯著性異質性如何處理研究間存在的顯著性異質性 提醒系統(tǒng)綜述的讀者注意，并對匯總資料作出解釋時要慎重探討造成差異的可能原因，并力圖加以解釋 核對資料提取有無錯誤？ 研究設計的質量如何? 劑量有無差異? 對象特征方面有無差異?5454模型的選擇模型的選擇 依據(jù)異質性檢驗的結果： 如果異質性檢驗的結果P0.05，可認為各研究間不存在異質性，可采用固定效應模型，主要統(tǒng)計方法包括：Mantel-Haenszel法、Peto法

27、和General variance-Based法； 如果異質性檢驗的結果P0.05，可認為各研究間存在異質性，應采用隨機效應模型，其主要的統(tǒng)計方法有DerSimonianLaird法。5555StataStata軟件簡介軟件簡介Stata軟件由美國計算機資源中心研制，現(xiàn)為stata公司的產(chǎn)品。1985年推出。小巧玲瓏、操作簡單、統(tǒng)計方法齊全，完美實現(xiàn)了使用簡單和功能強大兩者的結合，受到初學者和高級用戶的歡迎。與SAS、SPSS一起被并稱為新的三大權威統(tǒng)計軟件，是WHO和國際臨床流行病學網(wǎng)推薦使用的統(tǒng)計分析軟件。5656MetaMeta分析舉例分析舉例定性變量的meta分析定量變量的meta分析

28、方法累積meta分析Meta回歸分析5757MetaMeta分析舉例分析舉例-定性資料定性資料Silagy and Ketteridge reported a systematic review of randomized controlled trials investigating the effects of physician advice on smoking cessation.數(shù)據(jù)如下：5858MetaMeta分析舉例分析舉例-定性資料定性資料a：干預組戒煙人數(shù)；r1:干預組總人數(shù)；c：對照組戒煙人數(shù)；r2:對照組總人數(shù)；5959StateState軟件操作步驟如下軟件操作步驟如

29、下1、命令方式：gen b=r1-a; /*計算出干預組未發(fā)生人數(shù)*/gen d=r2-c; /*計算出對照組未發(fā)生人數(shù)*/ metan a b c d, label(namevar=name, yearvar=year) fixed rr xlabel(0.1,0.2,0.5,0.75,1,1.25,2,5,10)6060StataStata軟件操作步驟如下軟件操作步驟如下2、菜單方式：2.1 分別計算干預組和對照組未發(fā)生人數(shù)。操作如下：DataCreate or Change data Create new variable,出現(xiàn)如下對話框，如圖填入即可。同理可計算d。6161StataS

30、tata軟件操作步驟如下（續(xù)）軟件操作步驟如下（續(xù)）2、菜單方式2.2 meta分析。操作如下：User菜單Meta-Analysis of Binary and Continuous（metan）,彈出如下對話框：6262StataStata軟件操作步驟如下（續(xù)）軟件操作步驟如下（續(xù)）菜單方式2.3 上圖graph opts選項卡設置如下：設置坐標軸刻度6363結果結果6464結果結果6565L LAbbeAbbe圖的制作過程圖的制作過程異質性評價工具異質性評價工具User菜單Meta-AnalysisLabbe，出現(xiàn)如下對話框并設置如下。66666767L LAbbeAbbe圖的作用圖的作

31、用判斷研究結果有無異質性異質性。LAbbe圖是以研究中的干預組事件發(fā)生率相對于對照組事件發(fā)生率作圖，若研究結果同質，則所有點呈直線分布，偏離線過遠，則表明研究結果不同質。圖中圓圈的大小表示效應量的大小。6868GalbraithGalbraith星狀圖星狀圖異質性評價工具異質性評價工具User菜單Meta-Analysisgalbraith plot for，出現(xiàn)如下對話框并設置如下。69697070GalbraithGalbraith星狀圖作用及說明星狀圖作用及說明判斷研究結果有無異質性異質性。原理：Galbraith星狀圖是以標準化估計值（如logOR/lnRR/SE相對于其標準誤的倒數(shù)做

32、圖，若散點圖斜率較為接近則說明研究間同質）。如上圖，有一個研究落在了95%CI線以外，其它研究的散點斜率差別不是很大，提示研究間有一定同質性，但可接受。7171異質性的統(tǒng)計量異質性的統(tǒng)計量HH7272異質性評價統(tǒng)計量異質性評價統(tǒng)計量HH的計算的計算是通過對統(tǒng)計量Q進行自由度（文獻數(shù)）的校正，H統(tǒng)計量：7373異質性統(tǒng)計量異質性統(tǒng)計量HH的意義的意義7474漏斗圖的制作過程漏斗圖的制作過程發(fā)表性偏倚評價工具發(fā)表性偏倚評價工具User菜單Meta-Analysisfunnel Graph,ver，出現(xiàn)如下對話框并設置如下。 metafunnel _ES _selogES, var xtitle(r

33、isk ration) ytitle(1/SE(Effect size)75757676發(fā)表偏倚的評價發(fā)表偏倚的評價EggerEgger、BeggBegg檢驗及檢驗及EggerEgger圖、圖、BeggBegg圖圖步驟：User菜單Meta-Analysispublication Bias，出現(xiàn)如下對話框并設置如下。7777發(fā)表偏倚檢驗結果發(fā)表偏倚檢驗結果7878本質為旋轉的漏斗圖本質為旋轉的漏斗圖79798080剪補法（剪補法（metatrimmetatrim）發(fā)表偏倚識別工具發(fā)表偏倚識別工具User菜單Meta-Analysistrim and Fill Analysis，出現(xiàn)如下對話框并

34、設置如下。81818282敏感性分析敏感性分析metainfmetainf敏感性分析，即通過比較納入或排除某些研究進行的meta分析結果，探討被去除的研究對合并效應的影響程度，從而判斷結果的穩(wěn)定性。且通過比較可了解異質性的來源。操作過程：User菜單Meta-Analysisinfluence analysis，metan-based，出現(xiàn)如下對話框并設置如下。8383848485858686NNTmetanntNNTmetannt8787。8888MetaMeta分析舉例分析舉例-定量資料定量資料例：比較氟含量對身體發(fā)育的影響，以女童掌II皮質厚度為指標，收集11個研究的數(shù)據(jù)如下，試進行me

35、ta分析。8989StateState軟件操作步驟如下軟件操作步驟如下操作如下：User菜單Meta-Analysis of Binary and Continuous（metan）,彈出如下對話框：9090結果結果91919292漏斗圖的制作漏斗圖的制作93939494累積累積metameta分析分析意義 累積meta分析（cumulative meta-analysis）指當一個新的研究出現(xiàn)時，將其納入到以前的meta分析中，并重新進行meta分析。 在實施時，一般按照臨床試驗的先后順序，將先實施的試驗進行meta分析，以后每出現(xiàn)一個就納入一個，這樣可提高效應的精確性。9595累積累積me

36、tameta分析舉例分析舉例（選自（選自Meta-analysis in stataMeta-analysis in stata中的累積中的累積metameta分析章節(jié)：鏈激酶治療心肌梗死）分析章節(jié)：鏈激酶治療心肌梗死）干預組對照組9696StataStata操作步驟操作步驟Step1：計算兩組未死亡人數(shù)，命令如下： gen nondeath1=pop1-deaths1gen nondeath0=pop0-deaths0step2：計算log(OR)值及其方差，命令如下： gen logor=log( deaths1/ nondeath1)/( deaths0/ nondeath0) gen

37、varlogor=1/ deaths1+1/ nondeath1+1/ deaths0+1/ nondeath0Step3：meta分析，命令為： meta logor varlogor,var eform graph(f) id( trialnam) xlab(0.1,0.5,1,2,10) ltr(0.1) rtr(10) cline xline(1) print b2(odds ratio) fmult(1.5)9797MetaMeta分析結果分析結果98989999上例的菜單操作方式上例的菜單操作方式操作如下：User菜單Meta-Analysis cumulative,彈出如下對話框

38、：100100上例的菜單操作方式上例的菜單操作方式101101累積累積MetaMeta分析分析命令如下： metacum deaths1 nondeath1 deaths0 nondeath0, label(namevar=trialnam, yearvar=year) fixed or xlabel(0.1,0.5,1,2,10)102102累積累積MetaMeta運行結果運行結果103103MetaMeta回歸回歸Meta回歸的意義 各研究的療效存在異質性時，可用meta回歸(meta-regression)對療效與研究特性的關系進行分析。 meta回歸是亞組分析的一種擴展，可進行多因素分

39、析，一般要求研究的文獻數(shù)量不少于10個。104104MetaMeta回歸舉例回歸舉例選自選自Meta-analysis in stataMeta-analysis in stata中的中的metameta回歸分析回歸分析章節(jié)：章節(jié)：BCGBCG預防結核）預防結核）105105MetaMeta回歸舉例回歸舉例上例變量的意義為106106分析過程分析過程gen logor=log( rt/( nt- rt)/( rc/( nc- rc) gen selogor=sqrt(1/rc)+(1/( nc- rc)+(1/ rt)+(1/(nt-rt) meta logor selogor,eform g

40、raph(r) cline xlab(0.5,1,1.5) xline(1) boxsh(4) b2(odds ratio-log scale)107107MetaMeta分析結果分析結果108108Sort lat此圖和上圖是一樣的，只不過是按照lat排序后顯示。109109MetaMeta回歸回歸命令：metareg logor lat, wsse(selogor) bsest(reml) noiter graph結果為：110110MetaMeta回歸分析結果回歸分析結果111111MetaMeta回歸舉例回歸舉例2 2（該例選自（該例選自“丁香園網(wǎng)站丁香園網(wǎng)站”）112112步驟步驟1

41、、先做單因素的meta回歸（itt，alloc,phd無意義）2、將有意義的因素再做多因素的meta回歸（見后）3、結果的解釋113113結果結果114114metameta分析結果的穩(wěn)定性和可靠性分析分析結果的穩(wěn)定性和可靠性分析穩(wěn)定性分析，有：符合方案集分析；意向性分析失安全數(shù)；最佳結果演示最差結果演示以上內容詳見循證醫(yī)學（楊克虎主編）115115metameta分析結果的穩(wěn)定性和可靠性分析分析結果的穩(wěn)定性和可靠性分析發(fā)表偏倚分析，有：漏斗圖法；Beggs法；Eggers法；剪補法；116116與與statastata有關的有關的metameta分析文獻解讀分析文獻解讀文獻1：stata在M

42、eta分析中的應用.柏建玲等，2007; 介紹了在進行meta分析時，stata軟件與RevMan軟件相比的優(yōu)勢；并介紹了meta分析的常用命令；通過列舉3個例子，分別介紹了二分類計數(shù)資料的meta分析、計量資料的meta分析和生存分析數(shù)據(jù)的meta分析；最后，將stata軟件和RevMan軟件的功能做了比較。117117與與statastata有關的有關的metameta分析文獻解讀分析文獻解讀( (續(xù)續(xù)) )文獻2：應用stata軟件的meta分析完成醫(yī)學科研定量評價.王靜等，2008年.介紹了metan,meta,metap,metacum(舉例),metareg（舉例）,敏感性分析命令metainf和metaninf;發(fā)表偏倚的識別命令metabias和metatrim的用法，并介紹了發(fā)表性偏倚的識別工具，如漏斗圖、Beggs檢驗和Eggers檢驗、剪補法如何解讀。118118與與statastata有關的有關的metameta分析文獻解讀分析文獻解讀( (續(xù)續(xù)) )文獻3：stata軟件在臨床試驗計量資料效應比較的Meta分析中的應用.莫傳偉等,2005.該文獻重