网状Meta分析在Stata软件中的实现-曾宪涛

2013,13(11):13871391•1387•©2013•1,23131.4420002.4420003.430071WinBUGSRMetaStataRStataADORStataStataWinBUGSMetaMetaStataWinBUGSImplementationofNetworkMeta-AnalysisUsingStataSoftwareZENGXian-tao1,2,XUChang3,ZHANGChao1,LISheng31.CenterforEvidence-BasedMedicineandClinicalResearch,TaiheHospital,HubeiUniversityofMedicine,Shiyan442000,China;2.DepartmentofStomatology,TaiheHospital,HubeiUniversityofMedicine,Shiyan442000,China;;3.DepartmentofUrology,ZhongnanHospital,WuhanUniversity,Wuhan430071,ChinaAbstractTheWinBUGSsoftwarecanbecalledfromeitherR(providedR2WinBUGSasanRpackage)orStatasoftwarefornetworkmeta-analysis.UnlikeR,StatasoftwareneedstocreaterelevantADOscriptsatfirstwhichsimplifyoperationprocessgreatly.SimilarwithR,Statasoftwarealsoneedstoloadanotherpackagewhendrawingnetworkplots.Thisarticlebrieflyintroduceshowtoimplementnetworkmeta-analysisusingStatasoftwarebycallingWinBUGSsoft-ware.KeywordsNetworkmeta-analysis;Statasoftware;WinBUGSsoftware;Networkplot;BayesianstatisticsDOI:10.7507/1672-2531.2013023720122012A05020112011CZX012013EBM20130041984//MetaEmail:[e-mail]-analysis10Meta[1,2]RStataMetaStataMetaStataMetaStataWinBUGSJohnThompson[3]WinBUGS-fromstataStatamvmetaAnnaChaimani[4]RR2WinBUGSMeta[5]StataWinBUGSMeta1StataWinBUGSStataWinBUGSMetaMetaWinBUGS[6]FStata12.0WinBUGS1.4.3ADO[7]:\stata12.0\ado\baseStatahelpWinBUGShelpChinJEvid-basedMed2013,13(11):1387-1391•1388•CJEBM©2013EditorialBoardofChinJEvid-basedMed[8]2.1netmodel.txtWinBUGSWinBGUSStataWinBUGSWinBUGSMetaDoodle[8]4.5MetaWinBUGS2.2netlist.txtStataStataWinBUGSnetlist.txtlist(ns=64,nt=14)nsntnet.txtt[,1]t[,2]t[,3]r[,1]r[,2]r[,3]n[,1]n[,2]n[,3]na[]1277376831521501543121266786612412211831212557779121120119314NA40205NA161324NA214NA3952NA122125NA21113NA4848NA5355NA21213NA3746NA6062NA21214NA4149NA7275NA21214NA5656NA7984NA21214NA4549NA8278NA2ENDt[,1]t[,2]t[,3]t[,3]NAr[,]n[,]na[]na[]3ENDENDStatalogscript2.3Inits.txtWinBUGS“Inits.txt”list(d=c(NA,0,0,0,0,0,0,0,0,0,0,0,0,0),sd=1,mu=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0))list(d=c(NA,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1),sd=4,mu=c(-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3,-3))list(d=c(NA,2,2,2,2,2,2,2,2,2,2,2,2,2),sd=2,mu=c(-3,5,-1,-3,7,-3,-4,-3,-3,0,-3,-3,0,3,5,-3,-3,-1,-3,-7,-3,-3,5,-1,7,0,1,-4,5,0,-3,5,-1,-3,7,-3,-4,-3,-3,0,-3,-3,0,3,5,-3,-3,-1,-3,-7,-3,-3,5,-1,7,0,1,-4,5,0,1,-1,5,-4))2.4OutIndex.txt500500or[1,2]1500or[1,3]5001000or[1,4]10011500or[1,5]150120002013,13(11):13871391•1389•©2013•or[12,13]4500145500or[12,14]4550146000or[13,14]4600146500tau4650147000sd4700147500Stata2.5StataStatanetscript.txtStatascriptdisplay('log')check('f:\network\netmodel.txt')data('f:\network\netlist.txt')data('f:\network\net.txt')compile(1)inits(1,'f:\network\netInits.txt')gen.inits()update(500)set(or)set(tau)set(sd)update(500)coda(*,'f:\network\netout')quit()quit()log3MetaStataWinBUGSwbrun,script(f:\network\netscript.txt)WinBUGS(F:\WinBUGS14\WinBUGS14.exe)wbcoda,root(f:\network\netOut)clearsum3wbrunStataWinBUGSStataMetawbcodaStatasumWinBUGSscriptWinBUGS144.1STATAroutines[4]:\stata12.0\ado\basenetworkplot7Stata2AnnaChaimani[4]2t1t2r1r2n1n24.2networkplott1t2,[noweight]t1t2[]noweightStataMeta[9]4.3StataintervalplotifplotChinJEvid-basedMed2013,13(11):1387-1391•1390•CJEBM©2013EditorialBoardofChinJEvid-basedMed[10]StataWinBUGSMetaR[5]StataWinBUGSMetaStataMetaStataifplotStata2Stata1Meta2013,13(11):13871391•1391•©2013•StataprogramStataMetaWinBUGS1,,,.Meta:Meta.,2012,4(5):399-402.2LumleyT.Networkmeta-analysisforindirecttreatmentcompari-sons.StatMed,2002,21(16):2313-2324.3ThompsonJ,PamlerT,MorenoS.BayesianAnalysisinstatawithWinBUGS.TheStataJournal,2006,6(4):530-549.4ChaimaniA.UsingSTATAforNetworkMeta-Analysis[EB/OL].Availablefrom:(9):1137-1144.6,,,.Meta:MetaWinBUGS.,2012,4(5):395-398.7WinBUGSfromStataFileListing.Availablefrom:(9):1137-1144.9,,.Meta.,2013,13(11):1387-1391.10WhiteIR,BarrettJK,JacksonD,etal.Consistencyandinconsist-encyinnetworkmeta-analysis:modelestimationusingmultivariatemeta-regression.ResearchSynthesisMethods,2012,3(2):111-125.2013–07–042013–09–18