Stata在Meta分析中的应用(LHF).pptx

Stata在Meta分析中的应用Meta分析在Stata中进行数据操作时常用到的命令关键词表1Meta分析常用命令关键词命令关键词执行的命令metan固定效应模型与随机效应模型meta分析metafunnel漏斗图的绘制metabiasMeta分析中检验发表偏倚midas用于诊断性meta分析metainfMeta分析中单个研究的影响metaregMeta回归分析结果窗口以Stata12.0版本为例命令窗口以Stata12.0版本为例命令回顾窗口以Stata12.0版本为例变量名窗口以Stata12.0版本为例一、Stata实现Meta分析的基本命令不同的文献可能选择不同的效应量，如二分类资料的RR、OR或连续型变量的SMD（合并标准化均数差）等；Stata提供了统一的的命令metan，metan命令可以合并二分类及连续数据，对效应量及相应标准误或可信区间进行Meta分析，并能得到P统计量进行异质性检验，同时给出森林图。其命令格式为：metan变量名,[选择项]变量名是必须的，选择项是可供选择的，依具体情况而定。变量名（提取的资料类型）二分类资料metan命令后跟四个变量。依次为试验组发生事件(如死亡)数、未发生事件(如未死亡)数、对照组发生事件(如死亡)数和未发生事件(如未死亡)数，如:metandeath1live1death2live2连续性资料metan命令后跟六个变量，分别是试验组的样本量、均数、标准差，对照组的样本量、均数、标准差，如:metann1m1sd1n2m2sd2变量名metan变量名,[选择项]效应值及其标准误，如对数化OR及log(OR)的标准差、SMD及其标准差效应量及其方差效应量及其95%置信区间(生存分析）通常我们可以先利用metan进行计算，Stata会自动生成相关的效应量，如_ES、_selogES，我们可以直接利用这些效应量。选择项metan变量名,[选择项]用于二分类资料的(四变量)的选项有：rr合并相对危险度(riskratio),为默认选项or合并优势比(oddsratios)fixed指定固定效应模型peto指定Peto法合并优势比random指定随机效应模型选择项用于连续性资料（六变量）的选项有：cohen用Cohen法合并标准差，为默认选项hedges用Hedge法合并标准差nostandard合并加权均数之差(非标化)，否则合计标化均数之差fixed指定固定效应模型，为默认选项random指定随机效应模型表2阿司匹林预防心肌梗塞死亡的7项临床试验结果StudynamePublicationyearAspiringroup?Placebogroupdeathslives?deathslivesMRC-1197449566?67557CDP197644714?64707MRC-21979102730?126724GASP197932285?38271PARIS198085725?52354AMI2192038ISIS-2198815707017体操作步骤第一步，在Stata中录入数据。在Stata命令窗口中键入edit或在StataData菜单栏中单击Dataeditor，即可进入Stata的数据编辑器，按列输入数据，见下图。第二步，进行Meta分析。metandeath1live1death2live2,or连续型数据：有关短程小剂量强的松vs.安慰剂或非甾体抗炎药治疗类风湿性关节炎的7个临床随机对照试验(RCTs),观察类风湿性关节炎患者的关节压痛指数具体操作步骤第一步，在Stata中录入数据。其中变量study表示纳入的研究,year表示年份,n1表示试验组的样本含量,mean1表示试验组的均数,sd1表示试验组的标准差,n2表示对照组的样本含量,mean2表示对照组的均数,sd2表示对照组的标准差，见下图。第二步，进行Meta分析。Stata命令如下：metann1mean1sd1n2mean2sd2,randomcohen漏斗图的绘制与不对称性检验漏斗图、线性回归法（Egger法）、秩相关检验（Begg法）metafunnel（漏斗图）funnel（漏斗图）metabias（Egger法、Begg法）如果不存在发表偏倚，理论上所纳入的每个独立的研究效应值的点估计值围绕研究真实值对称分布而呈现一倒置的漏斗。metafunnelfunnel以阿司匹林为例：genlogor=log(_ES)rename_selogESselogorme