用 stata进行单个率meta分析程序总结感谢版主对我的方法进行验证，这里整理一下方面大家研究谷歌的程序 （标红部分， 分批录入 stata12.0.可得到结果 ）clear input study cases total 1 20 1000 2 40 5000 3 30 1500 4 25 3300 end gen p = . gen se = ./ get proportions and std errors forv i =1(1)4 cii totali casesi 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 4 页 - - - - - - - - - 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 4 页 - - - - - - - - - 1 76 451 2 86 202 3 24 97 4 401 2502 end gen p = . gen se = .forv i =1(1)4 cii totali casesi 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 4 页 - - - - - - - - - 我自已编的程序结果见贴子中的图片：录入格式， 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 4 页 - - - - - - - - - 。
