基于Stata的生存资料的meta分析代码

Meta分析一般步骤：

1.异质性：得到汇总结果以及各研究间数据差异

         I2：异质性大小（常见I2大于等于50%选择随机效应模型）

         亚组：二分类变量

         回归：连续性变量------纳入文献需要10篇

2.敏感性分析：评价结果稳定性

3.发表偏倚：评价meta分析最常见的系统误差

         识别：漏斗图、Begg检验、Egger检验------纳入文献需要>5篇或>8篇

         发表偏倚处理(剪补法、失安全系数法)---------存在发表偏倚时才进行此项分析

具体实施：

1.代码

gen lnHR = ln(hr)                            #数据转换
gen selnHR = (ln(ul) - ln(ll))/3.92

metan lnHR selnHR, label(namevar = author,yearvar = year) Random texts(120) eform effect(HR) boxsca(0.5) boxsha(1)             #异质性
metan lnHR selnHR, label(namevar = author,yearvar = year) Random texts(120) eform effect(HR) boxsca(0.5) boxsha(1) by()        #亚组（在by里面填写分组依据）

metaninf lnHR selnHR, label(namevar = author,yearvar = year) Random texts(120) eform effect(HR) boxsca(0.5) boxsha(1)             #敏感性

metafunnel lnHR selnHR                       #发表偏倚
metabias lnHR selnHR, graph(begg)            #begg检验
metabias lnHR selnHR, graph(egger)           #egger检验
metatrim lnHR selnHR, funnel                 #剪补法

metareg lnHR age, wsse(selnHR) bsest(reml)   #回归（大于10篇用）

2.剪补法的阐述

         我们采用漏斗图直观的展示发表偏倚，采用Egger和 Begg对漏斗图进行分析，分析结果显示Egger和Begg检验提示存在偏倚(P=0.000)。因此，考虑研究间存在发表偏倚，对漏斗图采用剪补法进一步分析，模型中加入xxx项研究后使得漏斗图对称，此时合并效应量为xxx(xxx,xxx)。