软件meta分析命令-2018朱教授原创

1、R软件Meta分析命令集合此命令由医药加朱教授编辑原创，付出大量劳动，请学员不要外传！1、计数资料Meta分析 (RR，OR， RD，Peto-OR)注释命令加载meta包library(meta)调取表格data-read.table(file.choose(),header=T,sep=,)浏览和编辑导入数据edit(data)执行计数资料meta分析metabin(case1,tot1,case0,tot0, data=data, sm=RR, method=MH, comb.fixed=T,comb.random=F, studlab=study)将meta结果命名为meta88文件m

2、eta88-metabin(case1,tot1,case0,tot0,data=data, method=MH,sm=RR,comb.fixed=T,comb.random=F,studlab=study)作森林图forest(meta88)做拉贝图labbe(meta88)做加尔布雷斯图radial(meta88)做博雅图baujat(meta88)敏感性分析metainf(meta88)敏感性分析制图forest(metainf(meta88)漏斗图funnel(meta88)剪补法分析并作图trimfill(meta88)funnel(trimfill(meta88)偏倚回归分析met

3、abias(meta88,method=linreg)回归分析metareg(meta88,TMN+AGE)注意：发表偏倚计算需要10个以上研究数目，否则报错。其他：exlude=1， byvar=，subset=, method=,title=, complab=, outclab=, method=MH,Inverse,GLMM,Peto。发生率低于1%为稀有事件，选取Peto-OR。2、计量资料Meta分析注释命令调取表格data-read.table(file.choose(),header=T,sep=,)编辑和浏览数据edit(data)做计量meta分析metacont(n1,m

4、1,sd1,n2,m2,sd2,data=data, sm=MD,comb.random=T,comb.fixed=F, studlab=study,byvar=race)命名为meta99meta99-metacont(n1,m1,sd1,n2,m2,sd2,data=data,sm=MD, comb.random=T,comb.fixed=F,studlab=study,byvar=race)作森林图forest(meta99)做拉贝图labbe(meta99)做加尔布雷斯图radial(meta99)做博雅图baujat(meta99)敏感性分析metainf(meta99)敏感性分析制

5、图forest(metainf(meta99)漏斗图funnel(meta99)剪补法分析并作图trimfill(meta99)funnel(trimfill(meta99)偏倚回归分析metabias(meta99,method=linreg)回归分析metareg(meta99,TMN+AGE)3、比率资料注释命令调取表格data-read.table(file.choose(),header=T,sep=,)编辑和浏览数据edit(data)做meta分析meta66=metaprop(events,total,studlab=study,sm=PFT,data=data) meta66

6、作森林图forest(meta66)做拉贝图labbe(meta66)做加尔布雷斯图radial(meta66)做博雅图baujat(meta66)敏感性分析metainf(meta66)敏感性分析制图forest(metainf(meta66)漏斗图funnel(meta66)剪补法分析并作图trimfill(meta66) funnel(trimfill(meta66)偏倚回归分析metabias(meta66,method=linreg)回归分析metareg(meta66,TMN+AGE)Available methods=Inverse variance；generalised li

7、near mixed model (GLMM)Logit transformation (sm=PLOGIT, default) Log transformation (sm=PLN) Freeman-Tukey Double arcsine transformation (sm=PFT) Arcsine transformation (sm=PAS) Raw, i.e. untransformed, proportions (sm=PRAW) 生存资料、泊松资料、相关系数和均数等（4 , 5, 6, 7, 8）数据类型命令data-read.table(file.choose(),heade

8、r=T,sep=,)edit(data)4、生存资料metagen(logHR, Se, sm=REML,data=data,comb.fixed=T,comb.random=F, studlab=study) (Log Hazard Ratios, Standard errors)Generic Inverse Variance Method:REML，ML，HS，SJ，HE，EB5、事件和人时-两组比较(密度、频次)metainc(event.e,time.e,event.c,time.c,sm=IRR,data=data, comb.fixed=F,comb.random=T,studl

9、ab=study,method=Cochran)泊松资料Note: Incidence Rate Ratio (sm=IRR), Incidence Rate Difference (sm=IRD) Method=MH, Inverse, Cochran, GLMM6、事件和人时-单组metarate(events,times,sm=IRLN,data=data, comb.fixed=T,comb.random=F, studlab=study,method=Inverse)泊松资料Note: Log transformation (sm=IRLN, default) Square root

10、 transformation (sm=IRS) Freeman-Tukey Double arcsine transformation (sm=IRFT) No transformation(sm=IR) Method= Inverse, GLMM7、相关系数 + 样本数metacor(cor, n, sm=zcor, data=data, comb.fixed=T,comb.random=F, studlab=study)sm=ZCOR or COR Method:Inverse Variance Weighting 8、均数和标准差metamean(n, mean, sd, sm=MRA

11、W,data=data,comb.fixed=T,comb.random=F, studlab=study)Raw, i.e. untransformed, means (sm=MRAW, default) Log transformed means (sm=MLN) Method:Inverse Variance methodinverse variance weightingNote: Calculation of fixed and random effects estimates for meta-analyses with continuous outcome data; inver

12、se variance weighting is used for pooling. Calculation of fixed effect and random effects estimates (risk ratio, odds ratio, risk difference, or arcsine difference) for meta-analyses with binary outcome data. Mantel-Haenszel, inverse variance, Peto method, and generalised linear mixed model (GLMM) a

13、re available for pooling. For GLMMs, the HYPERLINK :25095/library/metafor/html/rma.glmm.html rma.glmm function from R package metafor (Viechtbauer 2010) is called internally. Restricted maximum-likelihood estimator (method.tau=REML) Maximum-likelihood estimator (method.tau=ML) Hunter-Schmidt estimator (method.tau=HS) Sidik-Jonkman estimator (method.tau=SJ) Hedges estimator (method.tau=HE) Empirical Bayes estimator (method.tau=EB).Meta-analysis of incidence rates (metainc)Generic inverse variance meta-analysis (metagen) Meta-analysis of single correlations (metacor) Meta-analysis of single mea