田金徽 李 倫 楊克虎,3

·循證醫(yī)學(xué)方法學(xué)·

頻率統(tǒng)計(jì)方法網(wǎng)狀Meta分析在STATA軟件中的實(shí)現(xiàn)

田金徽1,2李 倫2,3楊克虎1,2,3

2002年，Lumley[1]提出了網(wǎng)狀Meta分析，在其后的10多年時(shí)間內(nèi)，網(wǎng)狀Meta分析迅速得到臨床科研工作者青睞，并得到國(guó)際衛(wèi)生技術(shù)評(píng)估組織的認(rèn)可[2,3]。網(wǎng)狀Meta分析可以采用頻率統(tǒng)計(jì)學(xué)方法和貝葉斯方法來實(shí)現(xiàn)。頻率統(tǒng)計(jì)學(xué)方法通過統(tǒng)計(jì)樣本，基于假設(shè)檢驗(yàn)和推斷結(jié)論；貝葉斯方法基于貝葉斯原理，根據(jù)先驗(yàn)概率推斷后驗(yàn)概率。貝葉斯方法計(jì)算能力靈活且強(qiáng)大， 但需要較高的統(tǒng)計(jì)學(xué)知識(shí)；頻率統(tǒng)計(jì)學(xué)方法簡(jiǎn)單易懂。有研究顯示頻率統(tǒng)計(jì)學(xué)方法和貝葉斯方法的網(wǎng)狀Meta分析結(jié)果同樣可靠[4]。本文以1個(gè)實(shí)例系統(tǒng)介紹如何利用頻率統(tǒng)計(jì)學(xué)方法在Stata 軟件中實(shí)現(xiàn)網(wǎng)狀Meta分析的數(shù)據(jù)分析，并呈現(xiàn)證據(jù)關(guān)系圖、證據(jù)貢獻(xiàn)圖、不一致性檢驗(yàn)圖、發(fā)表偏倚檢驗(yàn)圖、網(wǎng)狀Meta分析結(jié)果圖和排序圖，以期更好地開展網(wǎng)狀Meta分析。

1 數(shù)據(jù)來源

以Li等[5]最近發(fā)表在PLoSOne的網(wǎng)狀Meta分析中單孔腹腔鏡膽囊切除術(shù)(A)、三孔腹腔鏡膽囊切除術(shù)(B)、四孔腹腔鏡膽囊切除術(shù)(C)和四孔小切口腹腔鏡膽囊切除術(shù)(D)術(shù)后麻醉劑應(yīng)用情況的數(shù)據(jù)作為本文數(shù)據(jù)來源，4種術(shù)式結(jié)局的原始數(shù)據(jù)如表1所示。

表1 原始數(shù)據(jù)

注n/N:各組術(shù)后麻醉劑應(yīng)用的病例數(shù)/接受治療措施總病例數(shù)

2.2 網(wǎng)狀Meta分析的常用命令

2.2.1 mvmeta命令(方法1) 由Chaimani等[6]提出，主要基于White開發(fā)的mvmeta程序包[7,8]。需將數(shù)據(jù)整理為直接比較格式和網(wǎng)狀Meta分析格式。直接比較格式將所有的研究拆分成兩兩比較的形式，如存在3臂研究拆分為3組數(shù)據(jù)，4臂研究拆分為6組數(shù)據(jù)。表2顯示了表1原始數(shù)據(jù)經(jīng)整理后的直接比較數(shù)據(jù)。

表2 直接比較數(shù)據(jù)

注n/N:各組術(shù)后麻醉劑應(yīng)用的病例數(shù)/接受治療措施總病例數(shù)

網(wǎng)狀Meta分析格式的數(shù)據(jù)呈現(xiàn)2個(gè)變量y和S。y指相對(duì)于共同對(duì)照結(jié)果的對(duì)數(shù)。在頻率統(tǒng)計(jì)法網(wǎng)狀Meta分析中，需要設(shè)立一個(gè)共同對(duì)照，然后基于共同對(duì)照計(jì)算其他干預(yù)措施的效果。在本文實(shí)例中，由于四孔腹腔鏡膽囊切除術(shù)為傳統(tǒng)腹腔鏡膽囊切除術(shù)，設(shè)為共同對(duì)照。但對(duì)于不存在共同對(duì)照的研究，需要賦予一個(gè)共同對(duì)照，使其總樣本數(shù)和反應(yīng)數(shù)較小，如分別為0.01和0.001，進(jìn)而計(jì)算其他所有干預(yù)措施的結(jié)果。S為方差(VAR)，是標(biāo)準(zhǔn)誤差(對(duì)數(shù)標(biāo)準(zhǔn)誤差)平方。對(duì)于雙臂研究，VAR為相對(duì)于共同對(duì)照結(jié)果對(duì)數(shù)的方差；對(duì)于3或4臂研究，除了共同對(duì)照之外的干預(yù)措施間還存在協(xié)方差。以3臂研究為例，分別存在3個(gè)干預(yù)措施A、B和C，VarAB、VarAC、VarBC分別表示3者兩兩比較時(shí)的方差，3者之間的關(guān)系如下：VarBC=(VarAB+VarAC- VarBC)，即VarBC=(VarAB+VarAC)/2。表2結(jié)果可以進(jìn)一步整理為表3的形式。

之后根據(jù)mvmeta程序包的相關(guān)命令(表4)，可以呈現(xiàn)網(wǎng)狀Meta分析所有結(jié)果的圖形表示方法，如證據(jù)關(guān)系圖、證據(jù)貢獻(xiàn)圖、不一致性檢驗(yàn)圖、發(fā)表偏倚檢驗(yàn)圖、網(wǎng)狀Meta分析結(jié)果圖和網(wǎng)狀Meta分析結(jié)果排序圖。該程序包的網(wǎng)狀Meta分析結(jié)果為一致性模型下的分析結(jié)果，同時(shí)提供不一致性因子對(duì)網(wǎng)狀Meta分析的不一致性進(jìn)行分析，彌補(bǔ)了該程序包不能提供不一致性模型結(jié)果的局限。本文實(shí)例證據(jù)貢獻(xiàn)圖結(jié)果顯示，如在A與B網(wǎng)狀Meta分析中，A與B、C和D直接比較的權(quán)重分別為57.9%、11.3%和6.9%；B與C、D間接比較的權(quán)重分別為33.9%和3.4%。結(jié)果排序顯示，對(duì)于術(shù)后麻醉劑應(yīng)用，A、B、C和D術(shù)式成為最佳干預(yù)措施的概率分別為0.64、0.10、0.02和0.24。

表3 網(wǎng)狀Meta分析格式的數(shù)據(jù)

注 yA、yB、yD分別為單孔腹腔鏡膽囊切除術(shù)、三孔腹腔鏡膽囊切除術(shù)、四孔小切口腹腔鏡膽囊切除術(shù)相對(duì)于四孔腹腔鏡膽囊切除術(shù)術(shù)后麻醉劑的應(yīng)用結(jié)果OR的對(duì)數(shù)，SAA、SBB、SDD分別為對(duì)應(yīng)的方差，SAB、SAD、SBD分別為多臂研究中協(xié)方差。本研究賦予共同對(duì)照時(shí)，反應(yīng)數(shù)賦予0.000 331，總樣本數(shù)賦予0.001

表4 mvmeta主要命令和功能

2.2.2 network命令(方法2) 為White等[9]開發(fā)，也是基于mvmeta進(jìn)行網(wǎng)狀Meta分析，該程序包置入Dias等[10]點(diǎn)分法，可對(duì)網(wǎng)狀Meta分析的不一致性進(jìn)行分析。可以一步實(shí)現(xiàn)網(wǎng)狀Meta分析數(shù)據(jù)的準(zhǔn)備。可以得到一致性和不一致性模型的網(wǎng)狀Meta分析結(jié)果，以森林圖的形式呈現(xiàn)，但network命令無法制作證據(jù)貢獻(xiàn)圖，提供的森林圖不顯示比較結(jié)果。本文實(shí)例排序結(jié)果顯示，對(duì)于術(shù)后麻醉劑應(yīng)用，A、B、C和D術(shù)式成為最佳干預(yù)措施的概率分別為0.66、0.10、0.02和0.21，與方法1結(jié)果基本一致。

2.3.3 Indirect命令(方法3) 由Miladinovic等[11]提出的

表5 network命令和功能

Indirect命令可以同時(shí)分析和呈現(xiàn)直接比較和間接比較Meta分析結(jié)果，僅用于簡(jiǎn)單的網(wǎng)絡(luò)，具體命令為indirect logOR selogOR study order, random eff(OR) eform trta(t1) trtb(t2)，其中random/fixed是效應(yīng)模型，eff()是結(jié)果效應(yīng)的名稱，trta()和trtb()分別是所要分析的干預(yù)措施。

本文實(shí)例中，如當(dāng)A和B的效果可以通過C和D得到間接比較結(jié)果的時(shí)候，該命令就不可以同時(shí)對(duì)多個(gè)路徑的間接比較結(jié)果進(jìn)行分析。且在采用該程序包的時(shí)候，必須嚴(yán)格地進(jìn)行編號(hào)，即先運(yùn)行哪一組，根據(jù)兩組直接比較的結(jié)果計(jì)算間接比較的結(jié)果，編號(hào)不同該程序包可能不產(chǎn)生結(jié)果或者產(chǎn)生錯(cuò)誤的結(jié)果。

3 總結(jié)

如采取頻率統(tǒng)計(jì)法進(jìn)行網(wǎng)狀Meta分析，上述3個(gè)程序包或者命令單一應(yīng)用都不足以制作完美的網(wǎng)狀Meta分析。如方法1和方法2均可實(shí)現(xiàn)干預(yù)措施相對(duì)于共同對(duì)照的治療效果，不能夠產(chǎn)生不同干預(yù)措施之間的效果差異，方法3只適用于不存在直接比較的情況下。因此，可以利用方法工程序包制作數(shù)據(jù)，進(jìn)行不一致性的檢驗(yàn)、結(jié)果排序；利用方法1的程序包，采用圖的形式展示證據(jù)關(guān)系圖、網(wǎng)狀Meta分析結(jié)果、不一致性、結(jié)果排序以及證據(jù)貢獻(xiàn)圖；利用方法3呈現(xiàn)間接比較的結(jié)果。

[1]Lumley T. Network meta-analysis for indirect treatment comparisons Stat Med,2002, 21(16): 2313-2324

[2]Lee AW. Review of mixed treatment comparisons in published systematic reviews shows marked increase since 2009. J Clin Epidemiol,2014, 67(2):138-143

[3]Li L, Tian J, Tian H, et al. Network meta-analyses could be improved by searching more sources and by involving a librarian. J Clin Epidemiol,2014, 67(9):1001-1007

[4]Carlin BP, Hong H, Shamliyan TA, et al. Case Study Comparing Bayesian and Frequentist Approaches for Multiple Treatment Comparisons Agency for Healthcare Research and Quality (US) ,2013

[5]Li L, Tian J, Tian H, et al. The Efficacy and safety of different kinds of laparoscopic cholecystectomy: a network meta analysis of 43 randomized controlled trials. PLoS One, 2014, 9(2): e90313

[6]Chaimani A, Higgins JP, Mavridis D, et al. Graphical tools for network meta-analysis in STATA. PLoS One, 2013,8(10): e76654

[7]White IR. Multivariate random-effects meta-regression: updates to mvmeta. Stata Journal,2011,11(2): 255-270

[8]White I. MVMETA: Stata module to perform multivariate random-effects meta-analysis. https://ideas.repec.org/c/boc/bocode/s456970.html

[9]White IR, Barrett JK, Jackson D, et al. Consistency and inconsistency in network meta-analysis: model estimation using multivariate meta-regression. Research Synthesis Methods,2012, 3(2): 111-125

[10]Dias S, Welton NJ, Caldwell DM, et al. Checking consistency in mixed treatment comparison meta-analysis. Stat Med,2010, 29(7-8): 932-944

[11]Miladinovic B, Hozo I, Chaimani A, et al. Indirect treatment comparison. Stata Journal,2014,14(1):76-86

(本文編輯：丁俊杰)

10.3969/j.issn.1673-5501.2014.06.014

1 蘭州大學(xué)基礎(chǔ)醫(yī)學(xué)院循證醫(yī)學(xué)中心 蘭州，730000；2 甘肅省循證醫(yī)學(xué)與臨床轉(zhuǎn)化重點(diǎn)實(shí)驗(yàn)室 蘭州，730000;3 蘭州大學(xué)第一臨床醫(yī)學(xué)院 蘭州，730000

楊克虎，E-mail：yangkehuebm2006@163.com

2014-09-03

2014-11-15)