Meta分析的常见问题
- 格式:doc
- 大小:224.50 KB
- 文档页数:4
【前言】
我也是一名Meta分析的初学者,仅在此分享我学习Stata10实现meta分析的详细步骤。
本系列文章主要和大家一同学习如何用Stata10实现meta分析,重点在于数据的要求和详细的菜单操作,而尽量避免探讨meta分析的原理,相关原理请参阅相关书籍及文献。
由于小生也是初学者,所以出现的错误敬请广大老师指明斧正,谢谢~
【编外】
安装好了stata,也下载到了meta分析模块,但是菜单却没有怎么办呢?
这是很多像我一样的新手都会遇到的问题,其实很简单,只需要用笔记本打开profile.do文件,拷贝下面的代码粘贴。
-------------------------------------------------------------------------------
if _caller() >= 8 {
window menu clear
window menu append submenu "stUser" "&Meta-Analysis"
window menu append item "Meta-Analysis" "Of Binary and Continuous (meta&n)" "db metan" window menu append item "Meta-Analysis" "Of Effects (&meta)" "db meta"
window menu append item "Meta-Analysis" "Of p-values (meta&p)" "db metap"
window menu append item "Meta-Analysis" "Cumulative (meta&cum)" "db metacum"
window menu append item "Meta-Analysis" "Regression (meta®" "db metareg"
window menu append item "Meta-Analysis" "Funnel Graph, metan-based (f&unnel)" "db funnel" window menu append item "Meta-Analysis" "Funnel Graph, vars reordered (funnel&2)" "db funnel2"
window menu append item "Meta-Analysis" "L'abbe Graph, metan-based (&labbe)" "db labbe" window menu append item "Meta-Analysis" "NNT, metan-based (metann&t)" "db metannt" window menu append item "Meta-Analysis" "Influence Analysis, metan-based (metan&inf)" "db metaninf"
window menu append item "Meta-Analysis" "Influence Analysis, meta-based (metain&f)" "db metainf"
window menu append item "Meta-Analysis" "Galbraith Plot for Heterogeneity (&galbr)" "db galbr"
window menu append item "Meta-Analysis" "Publication Bias (meta&bias)" "db metabias" window menu append item "Meta-Analysis" "Trim and Fill Analysis (metatrim)" "db met&atrim" window menu refresh
}
-------------------------------------------------------------------------------
若还有安装方面的疑问,请参考本版块置顶帖及精华帖。
【正文】
【连载1:小试牛刀Stata10实现meta分析的P值合并】
很多关于meta分析的国内著作或者相关文献,都是从P值合并入手,因为这个的操作最简单明了。
国内最常采用的一个研究是1986年至1988年国产雷尼替丁治疗消化性溃疡的四个RCT试验。
表格如下:
表国产雷尼替丁治疗消化性溃疡的疗效
----------------------------------------------
编号年份 P值(双侧)
----------------------------------------------
1 1988 1.0000
2 1987 1.0000
3 1986 0.0237
4 1986 0.6399
----------------------------------------------
但是要是想用Stata10实现P值合并,需要对输入的数据进行处理,需要将数据处理如下:
---------------
单侧P值
---------------
0.5000
0.5000
0.0118
0.3199
----------------
. metan tdeath tnodeath cdeath cnodeath, or chi2 label(namevar=trialid)
. metan tdeath tnodeath cdeath cnodeath, rd random
. metan tdeath tnodeath cdeath cnodeath, olevel(99) ilevel(95) sortby(year) label(namevar=trialid, yearid=year) nostats nowt
. local ovratio=$S_1
. funnel , sample ysqrt xlabel(0.1,0.5,1,5,10) ylabel(0,500,1000) overall(`ovratio')
. labbe tdeath tnodeath cdeath cnodeath, xlabel(0,0.25,0.5,0.75,1) ylabel(0,0.25,0.5,0.75,1)
. metan n1 m1 sd1 n2 m2 sd2, nostandard
. metan n1 m1 sd1 n2 m2 sd2, random xla(-2,-1,0,1,2) force saving(metagrph)。