用stata进行单个率meta分析程序总结
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用stata进行单个率meta分析程序总结
感谢版主对我的方法进行验证，这里整理一下方面大家研究
谷歌的程序（标红部分，分批录入stata12.0.可得到结果。）clear
input study cases total
1 20 1000
2 40 5000
3 30 1500
4 2
5 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
1 76 451
2 86 202
3 2
4 97
4 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