Stata在Meta分析中的应用课件

1、Stata在Meta分析中的应用席波, 博士,讲师山东大学公共卫生学院xibo201020世纪70年代，Cochrane首次提出循证医学的思想，并组织了激素治疗早产孕妇的系统评价，挽救了成千上万的早产儿。自此在短短数十年年时间，循证医学的种子不但遍植全球，而且在20世纪末也开始在我国生根发芽。如今，这一理念业已被学院派的医生所接受，而且还在一定程度上改变了我国的医疗实践。按照CSCO候任主委、广东省人民医院副院长吴一龙教授的观点，循证医学的精髓可以概括为3点：一是证据的分级；二是建立起一套基于数理统计方法获取证据的研究体系；三是最佳证据、医生经验和患者需求三方面的紧密结合。就证据分级而言，按照

2、特定病种的特定疗法收集所有质量可靠的随机对照试验(randomizedcontrolledtrial，RCT)后，根据meta分析而得出的结论无疑是级别较高的证据。近年来，有关Meta分析的论文越来越多的见诸于国内外学术期刊。丁香园论文与基金频道与PubMed平台的小编收集了近年来国内外Meta分析的论文数据，制成图表，以利大家了解这方面的趋势变化情况。图1：Meta分析的世界“强”国 图2：近年来全球Meta分析文章的数量变化趋势 图3：近年来中国作者Meta分析文章的数量变化趋势 如图1、2、3所示，自打Meta分析进入学术界视野以来，被PubMed检索的Meta分析论文已经接近6万篇，呈

3、逐年快速增长的态势。发表Meta分析论文最多的三个国家分别是美国、英国和加拿大。中国紧随三强位居第四。从增长曲线来看，中国Meta分析论文虽然起步较晚，但增长迅速。相对全球Meta分析论文数量的变化趋势，中国Meta分析论文年度发表数量在2006年之前长期徘徊在100篇以下，但之后数量迅猛增长，增长曲线的斜率更是如火箭发射般陡峭。在发表Meta分析论文最多的十大期刊中，以循证医学为主要内容Cochrane Database Syst Rev杂志当仁不让，共发表了4491篇，雄踞榜首。一些老牌的综合性医学期刊，如BMJ和Lancet也刊登了大量的Meta分析论文。除此之外，挤进 10强榜单的的杂

4、志还有一个特点，即内容相对较为侧重内科，外科相当罕见。这可能与外科手术多有个人烙印有关。不过在现代医学中，内科与外科的界限已经越 来越模糊。外科手术越来越趋于微创，而内镜和介入等理念的的引入也使内科医生有了更多手术操作的机会。相信假以时日，这一特点或被逆转也未可知。图 5列出了发表中国作者Meta分析论文最多的10本期刊，前三名分别是Mol Biol Rep、Plos One和Cochrane Database Syst Rev，这其中还有三家中国杂志中华流行病学杂志、中华医学杂志和中国肺癌杂志，分别名列第6、7和10名。图6：中国Meta分析论文作者的地域分布图7：中国发表Meta分析论文的

5、十大“牛”人 推荐2个循证医学相关网站丁香园：/bbs/公卫论坛：/Xi B, Takeuchi F, Chandak GR, Kato N, Pan HW; AGEN-T2D Consortium, Zhou DH, Pan HY, Mi J. Common polymorphism near the MC4R gene is associated with type 2 diabetes: data from a meta-analysis of 123,373 individuals. Diabetologia. 2012 Oct;55(10):2660-6. (IF=6.8)Xi B,

6、 Chandak GR, Shen Y, Wang Q, Zhou D. Association between Common Polymorphism near the MC4R Gene and Obesity Risk: A Systematic Review and Meta-Analysis. PLoS One. 2012;7(9):e45731. (IF=4.0) Xi B, Zeng T, Liu L, Liang Y, Liu W, Hu Y, Li J.Association between polymorphisms of the renin-angiotensin sys

7、tem genes and breast cancer risk: a meta-analysis. Breast Cancer Res Treat. 2011 Nov;130(2):561-8. (IF=4.4)系统综述与Meta分析的关系 Meta分析是系统综述的一种类型 Meta分析是系统综述进行定量合成的统计方法Meta分析步骤选题-Novelty检索文献- (Pubmed and Embase (or Scopus)根据入选标准选择合格的研究复习每个研究并进行质量评估- (Jadad评分)提取信息，填写摘录表，建立数据库绘制森林图-(Stata/Revman)异质性检验(如果有异质性

8、，采用Meta回归探讨来源)亚组分析敏感性分析检验发表偏倚例1 为了探讨用Aspirin预防心肌梗塞(MI)后死亡的发生，美国在1976-1988年间进行了7个关于Aspirin预防MI后死亡的研究，其结果见表格，其中6项研究的结果表明Aspirin组与安慰剂组的MI后死亡率的差别无统计学意义，仅1项结果表明Aspirin预防 MI后死亡有效。Forest plot命令: metan death1 live1 death2 live2, or label (namevar=study, yearvar=year) metan death1 live1 death2 live2, or labe

9、l (namevar=study, yearvar=year) xlab(0.45,1,2.19) metan death1 live1 death2 live2, or label (namevar=study, yearvar=year) xlab(0.45,1,2.19) textsize(150)Sensitivity analysis命令： generate logor=log(_ES) generate selogor=_selogES metainf logor selogor, eform id(study) random printFirst authorPublicatio

10、n yearExperimental treatmentControl treatmentNoMeanSDNoMeanSDJasni1968916.28.7938.112.8Dick19702417.682440.713Lee19732130.516.52141.419.8Berry19741213111223.711.1Lee19741814.612.41826.415.1Stenberg1992216.31.72111.12.5Geital19952010.84.72016.37.7Forest plot 命令： metan n1 mean1 sd1 n2 mean2 sd2, label

11、 (namevar=study, yearvar=year) random textsize(180)Sensitivity analysis命令： metainf _ES _seES, eform id(study) random printPublication bias命令： metabias _ES _seES, graph(begg)Source of heterogeneity命令： metareg _ES year n1 n2, wsse (_seES) metan n1 mean1 sd1 n2 mean2 sd2,label (namevar=study, yearvar=y

12、ear) random textsize(160) by(group_n1)例3 探讨MC4R基因多态性与2型糖尿病的关系Frayling TM, Timpson NJ, Weedon MN, et al. A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science 2007; 316:889894 FTO 肥胖 2型糖尿病 Adjusted for BMIOR=1.03, 95%CI 0.96-1.10,

13、p=0.44OR=1.15, 95%CI 1.09-1.23, p=910-6OR=1.32, 95%CI 1.26-1.39, p=310-26Li H, Kilpelinen TO, Liu C,et al. Association of genetic variation in FTO with risk of obesity and type 2 diabetes with data from 96,551 East and South Asians. Diabetologia. 2012 Apr;55(4):981-95. gene lnor=ln(or) gene lnci1=ln

14、(ci1) gene lnci2=ln(ci2) metan lnor lnci1 lnci2, eform label (namevar=study) textsize(160) xlab(0.28,1,3.13) gene logor=log(or) gene selogor=(log(ci2)-log(ci1)/3.92 metainf logor selogor, id (study) eform printmetabias logor selogor, graph(begg)metan lnor lnci1 lnci2,eform label (namevar=study) text

15、size(130) xlab(0.27,1,3.13) by (ethnicity)第一讲，Stata12.0软件的下载安装及meta分析模块的正确安装 /detail/4115006第二讲：Stata常用meta分析命令及二分类变量森林图的绘制 /detail/4115032第三讲：连续性变量森林图的绘制及如何在stata中实现亚组分析/detail/4116805 第四讲：Stata检测发表偏倚及绘制漏斗图/detail/4117081第五讲：Stata进行敏感性分析及结果解读/detail/4150975第六讲：Stata进行累计meta分析及结果解读/detail/4150978第七讲：Stata进行meta回归分析与回归结果的解读/deta