meta分析中結(jié)果有些會產(chǎn)生偏倚，偏倚的識別方法有: 漏斗圖法、 Begg秩相關(guān)法、 Egger 回歸法。我們今天來介紹使用stata繪制meta分析漏斗圖，雖然講得有點(diǎn)遲了，但是還要講一講，也許有些人需要呢。
首先要進(jìn)行stata的meta分析模塊進(jìn)行安裝，不懂的看看我的上一篇文章《Stata繪制(meta分析)森林圖（1）》，安裝好以后我們先要導(dǎo)入數(shù)據(jù)（回復(fù)：森林圖數(shù)據(jù)1，可以獲得數(shù)據(jù)），我就直接黏貼進(jìn)去了
一個很簡單的數(shù)據(jù)，stud表示研究名稱，a為實(shí)驗(yàn)組陽性人數(shù)，b為實(shí)驗(yàn)組總?cè)藬?shù)，c為對照組陽性人數(shù)，d為對照組總?cè)藬?shù)
因?yàn)槁┒穲D是對森林圖進(jìn)行的檢驗(yàn)，所以我們要先進(jìn)行森林圖繪制，具體可以看文章《Stata繪制(meta分析)森林圖（1）》，我這里就直接代碼了
metan a b c d, label(namevar=study) fixed rr
生成圖后也生成了部分?jǐn)?shù)據(jù)，我們在數(shù)據(jù)框看一下，我們要對_ES這個數(shù)據(jù)（下圖）轉(zhuǎn)換一下，取它的對數(shù)，等下用它來作圖
gen logrr=log( _ES )
生成了logrr這個變量就可以作圖，可以界面操作或者代碼,界面操作：
代碼操作就是：
metabias logrr _selogES, graph(begg)
有時候stata會報錯，顯示command metabias is unrecognized
這是找不到metabias這個函數(shù)的意思，我們可以把這個安裝模塊的命令復(fù)制到stata重新安裝一遍就可以了（公眾號回復(fù)：代碼，可以獲得） ，
我們繪制好的漏斗圖是黑色的，想要變成白底的就直接把它復(fù)制黏貼進(jìn)word就可以了
除了漏斗圖，stata還給出了Begg’s Test和Egger’s test的結(jié)果，這也是很重的結(jié)果，P大于0.05表示沒有偏倚
Egger圖的繪制方法也是類似的
代碼是：
metabias logrr _selogES, graph(egger)
想具體了解偏倚怎么判斷的可以看這篇文獻(xiàn)，我覺得寫得很好。
本文結(jié)束。
本站僅提供存儲服務(wù)，所有內(nèi)容均由用戶發(fā)布，如發(fā)現(xiàn)有害或侵權(quán)內(nèi)容，請
點(diǎn)擊舉報。
