用stata进行单个率meta分析程序总结_meta分析stata实例

用stata进行单个率meta分析程序总结由刀豆文库小编整理，希望给你工作、学习、生活带来方便，猜你可能喜欢“meta分析stata实例”。

用stata进行单个率meta分析程序总结

感谢版主对我的方法进行验证，这里整理一下方面大家研究

谷歌的程序（标红部分，分批录入stata12.0.可得到结果。）

clear

input study cases total 20 1000 40 5000 30 1500 25 3300 end

gen p =.gen se =.// get proportions and std errors forv i =1(1)4 {

cii total[`i'] cases[`i']

qui replace p = r(mean)in `i'

qui replace se = r(se)in `i' }

// get the inverse variance-weighted proportion

// use the official Stata-vwls-command

gen cons =1

vwls p cons, sd(se)

// use the user written-metan-command

// for fixed-effects meta-analysis

metan p se, nograph fixed

// for random-effects meta-analysis

metan p se, nograph random

我的数据，用谷歌方法运行的命令：

clear

input study cases total

451 2 86 202 3 24 97 401 2502 end

gen p =.gen se =.forv i =1(1)4 {

cii total[`i'] cases[`i']

qui replace p = r(mean)in `i'

qui replace se = r(se)in `i' } gen cons =1

vwls p cons, sd(se)metan p se, nograph fixed metan p se, nograph random

我自已编的程序结果见贴子中的图片：

录入格式，r n clear input study r n 1 0.831 154 2 0.828 134 3 0.88 100 end

generate ser=sqrt(r*(1-r)/n)metan r ser, fixed label(namevar=study)metan r ser, random label(namevar=study)metafunnel r ser