用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
