Stata在Meta分析中的应用.pdf
《Stata在Meta分析中的应用.pdf》由会员分享,可在线阅读,更多相关《Stata在Meta分析中的应用.pdf(8页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2 0 0 7 年1 2 月第7 卷第6 期循证医学T h eJ o u r n a o fE v i d e n c e B a s e dM e d i c i n eD e e 2 0 0 7V 0 1 7N O 6循证医学中的医学统计学问题S t a t a 在M e t a 分析中的应用柏建岭1,钟文昭2郑明华3,陈峰(1 南京医科大学公共卫生学院流行病与卫生统计学系,南京2 1 0 0 2 9;2 广东省人民医院肿瘤中心,广州5 1 0 0 8 0;3 温州医学院附属第一医院感染内科,浙江温州3 2 5 0 0 0)摘要 S t a t a 是一个功能强大而又小巧玲珑的统计分析软件
2、,其M e t a 分析功能十分全面。本文结合3 个实例介绍rS t a t a 软件对二分类变量、连续性变量和生存分析等所涉及的S t a t a 命令及输出结果。和R e v i e wM a n a g e r 相比,S t a t a 可更个性化、更高效专业地完成M e t a 分析,且更易于为高影响力杂志所认同。关键词 统计学;软件;S t a t a;M e t a 分析 中图分类号 T P 3 1 9 文献标识码 A 文章编号1 6 7 1 5 1 4 4(2 0 0 7)0 6 0 3 6 3 0 6M e t a A n a l y s i si nS t a t aB A
3、IJ i a n l i n 9 1,Z H O N GW e n z h a 0 2,Z H E N GM i n g h u a 3,C H E NF e n 9 1(J D e p a r t m e n to fE p i d e m i o l o g ya n dB i a s t a t i s t i c s,S c h o o lo fP u b l i cH e a l t h,N a n j i n gM e d i c a lU n i v e r s i t y,N a n j i n g2 1 0 0 2 9,C h i n a;2 C a n c e rC e n
4、 t e r,G u a n g d o n gP r o v i n c i a lP e o p l e sH a s p i t a l,G u a n g z h o u5 1 0 0 8 0,C h i n a;3 D e p a r t m e n to fI n f e c t i o na n dL i v e rD i s e a s e s,T h eF i r s tA f f i l i a t e dH o s p i t a lo fW e n z h o uM e d i c a lC o l l e g e,Z h e j i a n gW e n z h o
5、u3 2 5 0 0 0,C h i n a)A b s t r a c t:S t a t ai sas m a l lb u tp o w e r f u ls o f t w a r e,w i t hc o m p r e h e n s i v ef u n c t i o nf o rm e t a a n a l y s i s T h i sm a n u s c r i p t,t o g e t h e rw i t ht h r e ec a s e s,i l l u s t r a t e dS t a t ap r o c e d u r e sa n dc o r
6、 r e s p o n d i n gr e s u l t sf o rp r o c e s s i n gd i c h o t o m o u sv a r i a b l e s,c o n t i n u o u sv a r i a b l e sa n ds u r v i v a la n a l y s i si nm e t a a n a l y s i s C o m p a r e dw i t hR e v i e wM a n a g e r,a n o t h e rp o p u l a rs o f t w a r ef o rm e t a a n a
7、 l y s i s,S t a t ai sm o r ep e r s o n a l i z e d,m o r ep r o f e s s i o n a lw i t he q u i v a l e n te f f i c a c y,a n dm o r ee a s i l yr e c o g n i z e db yh i g h i m p a c tm a g a z i n e s K e yw o r d s:s t a t i s t i c s;s o f t w a r e;S t a t a;m e t aa n a l y s i s1前言S t a
8、t a 是一个功能强大而又小巧玲珑的统计分析软件,最初由美国计算机资源中心(C o m p u t e rR e s o u r c eC e n t e r)研制,现为S t a t a 公司的产品。从1 9 8 5 年1 0 版问世以来。已连续推出1 0 个主要版本,并从4 0 版起进入W i n d o w s 时代。通过不断的更新和扩充,软件功能已日趋完善。它操作灵活、简单、易用,同时具有数据管理软件、统计分析软件、作者简介 柏建岭(1 9 7 8 一),男,江苏盐城人,生物统计学在读博士。通讯作者 陈峰,T e l:0 2 5 8 6 8 6 2 7 5 4,E-m a i l:d
9、r e h e n f e n g 1 6 3 C O B绘图软件、矩阵计算软件和程序语言的特点,在许多方面别具一格,和S A S、S P S S 一起并称为新的三大权威统计软件。S t a t a 的许多高级统计模块均是程序文件(A D O 文件),S t a t a 公司在这方面持开放的态度,允许用户自行修改、添加和发布A D O 文件,用户可随时到S t a t a 网站或者其他个人网站上寻找并下载所需的程序包安装后使用。这一特点使得全球的统计学家均乐于在S t a t a 上首先实现所研究的最新计算方法,并对外免费提供下载。从而使S t a t a 始终处于统计分析方法发展的最前沿,用
10、户几乎总是能很快找到最新统计计算方法的S t a t a 程序版本。这也使S t a t a 自身成为几大统计软件中升级最多、最频繁的一个。万方数据3 6 4循证医学2 0 0 7 年第7 卷第6 期R e v i e wM a n a g e r(简称R e v M a n)是国际C o c h r a n e协作网为系统评价(s y s t e m a t i cr e v i e w)I 作者所提供的专用软件,是C o c h r a n e 系统评价的一体化、标准化软件。从计算机软件的角度来看,它主要包括了C o c h r a n e 系统评价的文字处理与M e t a 分析两大功能
11、。是目前M e t a 分析专用软件中较成熟的软件之一。R e v M a n 虽然简单易学但只提供了对计数资料如有效率的优势比(o d d sr a t i o)、相对危险度(r e l a t i v er i s k)和差率(r i s kd i f f e r e n c e)的合并方法和计量资料如呈正态分布的资料(血压、身高等)的合并方法。同时,R e v M a n 缺乏诊断性试验的M e t a分析和M e t a 回归的功能。相对而言,S t a t a 的M e t a分析功能更全面和强大。该软件除了可以完成二分类变量和连续性变量的M e t a 分析,也可以进行M e t
12、a 回归分析、累积M e t a 分析、单个研究影响分析、诊断试验的M e t a 分析、剂量反应关系M e t a 分析、生存分析资料合并等几乎所有M e t a 分析方法;还可以对发表偏倚进行B e g g S 检验和E g g e r s 检验。在S t a t a 软件中,可以绘制M e t a 分析的相关图形,如森林图(f o r e s tp l o t)、漏斗图(f u n n e lp l o t)和L A b b e 图。国外文献中大部分M e t a 分析文章都是通过S t a t a 的m e t a a d o 模块完成。本文拟通过3 个实例介绍S t a t a 用于
13、M e t a 分析的方法。图1S t a t a9 2 启动后的界面2S t a t a 的界面功能和基本语法2 1S t a t a 的界面功能单击S t a t a 软件图标进入S t a t a,图1 为S t a t a9 2启动后的界面。除了W i n d o w s 版本软件通用的菜单栏、工具栏、状态栏等外,S t a t a 的默认启动界面主要由四个窗口构成,分述如下。2 1 1 结果窗口(S t a t aR e s u h s)位于界面右上部,软件运行中的所有信息,如所执行的命令、执行结果和出错信息等均在该处体现。窗口中使用不同的颜色区分不同的文本,如默认情况下白色表示命令
14、,红色表示错误信息。绿色和黄色为结果输出和注释。2 1 2 命令窗口(S t a t aC o m m a n d)位于结果窗口右下部。相当于D O S 中的命令行,此处用于键入需要执行的命令,回车后即开始执行,相应的结果则会在结果窗E l 中显示。2 1 3 命令回顾窗口(R e v i e w)位于界面左上部,所有执行过的命令会依次在该窗口中列出,单击后命令即被自动拷贝到命令窗口中:如果需要重复执行,用鼠标在“R e v i e w”窗口中双击相应的命令即可。2 1 4 变量名窗口(V a r i a b l e s)位于界面左下部列出当前数据集中的所有变量名称。2 2S t a t a
15、命令的基本语法S t a t a 命令的基本语法格式如下:特殊选项 关键词命令参数,命令选项中括号表示其中的内容不一定总是出现,下面对语句中的各元素分别加以解释。2 2 1 特殊选项是一些在大部分命令中通用的选项。由于执行的功能比较特殊,因此将他们提前,并使用空格和命令正文分隔。特殊选项中最常用的有分组执行相同语句的“b v”命令,按指定的条件重复执行的“f o r”命令等。2 1 2 关键词相当于一句话的主语,指明了所执行的是哪一条S t a t a 命令,关键词在一条命令中必须出现。大多数命令的关键词都是采用相关的英文单词简单易记,并且在S t a t a 中还允许对关键词进行缩写(每个命
16、令不同,无特殊规律),方便了使用。2 2 3 命令参数相当于一句话的谓语和宾语用于指明相应的命令在执行时需要使用的变量、参数等。大多数S t a t a 命令都需要指定参数,但也有例外,此时系统会自动按照缺省方式执行。比如“d e s c r i b e”命令,如果不指定任何参数,则系统会默认对当前使用的数据集中的所有变量进行描述。2 2 4 命令选项相当于一句话中的定语、状语、补语等修饰成分,用于对相应的命令进行限制或更精确的指定,万方数据柏建岭,等S t a t a 在M e t a 分析中的应用3 6 5在命令中不一定出现。3M e t a 分析中常用的命令表1 中列出了M e t a
17、分析时在S t a t a 中进行数据操作时常用的命令关键词 2 。表1S t a t a 中进行m e t a 分析时常用的命令关键词命令关键词执行的操作m e t a nm e t a c u n lm e t a i n fm e t a b i a sm e t a r e g固定效应模型和随机效应模型M e t a 分析累积M e t a 分析M e t a 分析中单个研究的影响M e t a 分析中检验发表偏倚M e t a 回归囿于篇幅,仅详细介绍M e t a 分析中最常用的m e t a n 和m e t a b i a s 命令。3 1m e t a n 命令m e t a
18、 n 是m e t aa n a l y s i s 命令,可以合并二分类及连续数据;而且对效应量及相应标准误或可信区间进行M e t a 分析;并能得到1 2 统计量进行异质性检验,同时给出森林图(f o r e s tp l o t)。其命令行格式为:m e t a n 变量名,选择项二分类资料m e t a 命令后跟四个变量。依次为试验组发生事件(如死亡)数、未发生事件(如未死亡)数、对照组发生事件(如死亡)数和未发生事件(如未死亡)数。连续性资料m e t a 命令后跟六个变量,分别是试验组的样本量、均数、标准差,对照组的样本量、均数、标准差,如m e t a nn lm ls d l
19、n 2m 2s d 2。已知效应量和标准误或效应量和可信区间上下限以O R 为例m e t a n 命令后跟两个变量(效应量的对数值和其标准误)如m e t a nl o g o rs e l o g o r,e f o r me f f e c t(o d d sr a t i o)或三个变量(效应量的对数值,可信区间上限的对数值,可信区间下限的对数值)如m e t a nl o g o rl o g u ll o g l l,e f o r me f f e c t(o d d sr a t i o)。3 1 1 用于二分类资料(四变量)的选项有:r r 串合并相对危险度(r i s kr
20、 a t i o s),为默认选项o r 水合并优势比(o d d sr a t i o s)r d 术合并率差(r i s kd i f f e r e n c e s)f i x e d 冰指定固定效应模型(M a n t e la n dH a e n s z e lm e t h o d),为默认选项f i x e d i 术指定固定效应模型(i n v e r s ev a r i a n c em e t h o d)p e t o 半指定P e t o 法合并优势比(o d d sr a t i o s)r a n d o m 木指定随机效应模型(D e r S i m o n
21、i a n&L a i r dm e t h o d),其异质性估计采用M a n t e l H a e n s z e lm o d e lr a n d o m i:|:指定随机效应模型(D e r S i m o n i a n&L a i r dm e t h o d),其异质性估计采用i n v e r s ev a r i a n c ef i x e de f f e c tm o d e l3 1 2 用于连续性资料(六变量)的选项有:c o h e n 术C o h e n 法合并标准差,为默认选项h e d g e s 半H e d g e 法合并标准差g l a s s
22、 木G l a s s 法合并标准差n o s t a n d a r d 术合并加权均数之差(非标化),否则合计标化均数之差f i x e d 术指定固定效应模型(i n v e r s ev a r i a n c em e t h o d),为默认选项r a n d o m 术指定随机效应模型(D e r S i m o n i a n&L a i r dm e t h o d)3 1 3 用于已知效应量和标准误或效应量和可信区间上下限选项有:e f r o m 术效应量及其可信区间指数化(仅用于合并R R 和O R)3 1 4 共用选项有:b y(b y v a r):指定分组变量亚组
23、分析时常用l a b e l(n a m e v a r=n a m e v a r ,y e a r v a r=y e a r v a r )水指定每个研究的标签b o x s c a()术指定森林图中方块的显示比例缺省为13 2m e t a b i a s 命令m e t a b i a s 命令用于B e g ga n dM a z u m d a r 等提出的秩相关检验及E g g e r 等提出的线性回归法用于检测漏斗图的对称性。命令行格式为:m e t a b i a s 变量名,选择项变量有三种形式:效应量及其标准误如m e t a b i a sl o g r rs e l
24、 o g r r,g r a p h(b e g g)。效应量及其方差如m e t a b i a sl o g r rv a r l o g r ri fs i t e=3,v a rg r a p h(e g g e r)。R R(O R)及其可信区间上下限如1 1 3 e t a b i a sI T1 1u l,c ib y(s i t e)。选项有:b y(b y v a r)术指定分组变量g r a p h(b e g g):l:画b e g g 漏斗图g r a p h(e g g e r):l:画e g g e r 不对称图4 实例分析例1 3 为了探讨用A s p i r i
25、 n 预防心肌梗塞(m y o c a r d i a li n f a r c t i o n,M I)后死亡的发生,美国在1 9 7 6 1 9 8 8 年问进行了7 个关于A s p i r i n 预防M I 后死亡 万方数据3 6 6循证医学2 0 0 7 年第7 卷第6 期的研究,其结果见表2,其中6 项研究的结果表明A s p i r i n 组与安慰剂组的M I 后死亡率的差别无统计学意义,只有1 项结果表明A s p i r i n 预防M I 后死亡有效并且差别有统计学意义。现根据表2 所提供的资料作M e t a 分析。表2A s p i r i n 预防心肌梗塞死亡的临
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Stata Meta 分析 中的 应用
限制150内