以下是柳叶刀兄回答一个关于如何进行egger偏倚计算的帖子。原贴见下面地址：在这个帖子中，柳叶刀兄给出了如何在stata中进行偏倚计算的方法，小弟还有几处不明白，恳请柳叶刀兄及各位老师指教，谢谢这很简单，用stata软件分析即可，过程如下（这是直接的方法，当然也可以先算出OR的标准误）：1.打开stata的数据编辑表，输入各个II/(ID十DD)OR值，然后同样的输入各OR值的95%CI的下限和上限。变量自定，我这里用low代表OR的下限，up代表上限制。2. 在stata命令窗口输入：metabias OR low up, ci graph(egger)命令。这里括号内的egger表示egger's作图法。即可以输出结果如下：请问各位老师：1 在metabias命令中的ci代表什么意思？2 如何求OR的标准误（用哪个公式）？3 在连续变量中，是不是同样输入mean differences， low ，up三列，即可在stata中进行偏倚的计算？恳请各位老师指教，谢谢！！太复杂了，呵呵，我不懂，帮你顶起来，呵呵，期待回答！谢谢甫寸兄！第二个问题我已经知道了：gen selogor=sqrt(1/a+1/b+1/c+1/d)还请各位老师继续指教其它问题，呵呵ci indicates that exp(theta) and its confidence interval, (ll, ul), were supplied on the command line instead of theta and se_theta. Option varshould not be specified when option ci is specified.如：metabias rr ll ul, ci by(site)第三个问题：可能不行。metabais是基于meta命令的，而后者只能处理分类数据；不正之处，请高手指正。谢谢张主任的回复！也就是说ci的输入，代表前面命令中输入的是疗效已经其上下限，而不需要输入疗效的对数及对数的标准误？exp（theta）是不是疗效的意思？Option var should not be specified when option ci is specified.这句不懂，能否麻烦您解释一下，谢谢对于第三个问题，我看到以前stata兄发的一个帖子，里面提供了如何运用stata进行meta分析，祥见?bid=96&id=6488335&sty=3&keywords=stata在stata兄的讲稿中，好像对于连续变量，也可采用metan命令进行meta分析的。紫竹林之舞 wrote:谢谢张主任的回复！也就是说ci的输入，代表前面命令中输入的是疗效已经其上下限，而不需要输入疗效的对数及对数的标准误？exp（theta）是不是疗效的意思？Option var should not be specified when option ci is specified.这句不懂，能否麻烦您解释一下，谢谢对于第三个问题，我看到以前stata兄发的一个帖子，里面提供了如何运用stata进行meta分析，祥见?bid=96&id=6488335&sty=3&keywords=stata在stata兄的讲稿中，好像对于连续变量，也可采用metan命令进行meta分析的。1.exp（theta）:the risk (or odds) ratio;（也就是指效应指标如：or/rr/） ll : the lower limit of the risk ratio's confidence interval ul : the upper limit of the risk ratio's confidence interval2.如果指定ci开关，则不需要var开关。3.stata兄的文章metabias分析的例子是meta-analysis in stata中的例子，还是分类数据。请再仔细看一下。4.可仔细看看metabias命令的帮助。5.metan命令可以用于连续数据的meta，而meta命令则不可以，而metabias是基于后者的，所以我推测可能不能用于连续数据的偏倚估计。请注意，stata中用于meta分析的基本命令有两个：metan和meta，前者应用范围较后者广。紫竹林之舞兄好像没认真看帖子啊，呵呵谢谢张主任的回复！！我以前并不知道所有命令基于meta和metan两个根本。我还以为meta和metan就是单纯的两个命令呢，呵呵又看了一下stata兄的讲稿，metan命令既可以用作连续变量的meta分析，又可以用作分类变量的meta分析。但metabias只做了分类变量的偏倚检测，呵呵stata兄所作的这个metabias的结果，是不是同时采用Begg法（漏斗图）和Egger法（线形回归）对偏倚进行检测？那连续变量的偏倚检测又用什么命令呢？ （缩略图，点击图片链接看原图）metabias就是采用Begg法和Egger法。不过方法不一样：Begg法是采用秩相关检验，而Egger法是采用线性回归法。紫竹林之舞 wrote:那连续变量的偏倚检测又用什么命令呢？个人认为要用到metan、funnle命令。运行metan后，立即执行funnle命令，试试看。谢谢张主任的回复！！我自己又搜索了一下，找到一个ppt课件,详细信息如下：名称：Meta-analysis in Stata Day 2 morning Publication bias and heterogeneity作者：Darren Greenwood Yu-Kang Tu单位：Biostatistics Unit, Centre for Epidemiology & Biostatistics, University of Leeds里面涉及了很多stata中的命令，发上来分享一下，呵呵 metainstata.rar (64.75k)继续请教张主任，在这个图中，如何解读这些结果？是不是在Begg法中，主要看第二个Pr > |z| = 这个指标？如果大于0.05，代表无偏倚？Egger法中主要看bias一行的Pr > |z|指标，如果大于0.05的话，代表无偏倚？ screen.width-333)this.width=screen.width-333" width=640 height=422 title="Click to view full 97817636_snap.jpg (640 X 422)" border=0 align=absmiddle>又搜索到一篇专门介绍stata中关于meta分析命令的文章，发上来一同分享，呵呵 STATA_Metaanalysis_commands_V6_March2004.rar (154.38k)紫竹林之舞 wrote:继续请教张主任，在这个图中，如何解读这些结果？是不是在Begg法中，主要看第二个Pr > |z| = 这个指标？如果大于0.05，代表无偏倚？Egger法中主要看bias一行的Pr > |z|指标，如果大于0.05的话，代表无偏倚？建议用Egger法，如图中p=0.007，说明有偏倚。支持,不知道怎么进行funnle命令,看说明也不明白,请帮助下载了funnel模块。funnel具体该怎么用俺还不知道就简单用以下命令测试了一下：funnel logor selogor发现，采用funnel命令只能生成一个漏斗图。但具体的的偏倚检测，并没有出现。同时不知道如何采取Begg法或者Egger法。请各位老师指教，谢谢！！紫竹林之舞 wrote:下载了funnel模块。funnel具体该怎么用俺还不知道就简单用以下命令测试了一下：funnel logor selogor发现，采用funnel命令只能生成一个漏斗图。但具体的的偏倚检测，并没有出现。同时不知道如何采取Begg法或者Egger法。请各位老师指教，谢谢！！funnel命令 当然只生成一个漏斗图，它是用漏斗图检验偏倚的，不是采用Begg法或者Egger法。检验发展偏倚，可以用漏斗图检验，也可以用统计学方法检验如Begg法或者Egger法。还请紫兄查阅相关文献，呵呵非常感谢张主任每次都及时的给俺帮助！！！这么看来，连续变量不可能采用基于meta的metabias命令来进行偏倚的检测。那有没有基于metan的相关命令，可以对连续变量进行偏倚检测呢（采用Begg法或Egger法）？紫竹林之舞 wrote:非常感谢张主任每次都及时的给俺帮助！！！这么看来，连续变量不可能采用基于meta的metabias命令来进行偏倚的检测。那有没有基于metan的相关命令，可以对连续变量进行偏倚检测呢（采用Begg法或Egger法）？funnel命令就是基于metan检验偏倚的命令，另一个是labbe(是检验分类数据的），funnel是采用漏斗图方法检验偏倚的，通过观察漏斗图是否对称而得出结论，并不是通过统计学方法如Begg法或Egger法而得出来的。metan命令执行结束后，不要执行其他命令，马上键入funnel命令，就可以得到漏斗图。funnel命令是仅仅靠视觉判断漏斗图是否对称来检验偏倚。那有没有命令是采用Begg法或Egger法对连续变量进行偏倚检测呢？（就像metabias可以用Begg法和Egger法对分类变量的偏倚进行检测一样）请继续指教，谢谢！紫竹林之舞 wrote:funnel命令是仅仅靠视觉判断漏斗图是否对称来检验偏倚。那有没有命令是采用Begg法或Egger法对连续变量进行偏倚检测呢？（就像metabias可以用Begg法和Egger法对分类变量的偏倚进行检测一样）请继续指教，谢谢！没有明白了，又学了很多东西，呵呵谢谢指教！！zhangdog wrote:第三个问题：可能不行。metabais是基于meta命令的，而后者只能处理分类数据；不正之处，请高手指正。那metareg应该也是基于meta命令的吧，那也不能应用于连续变量的回归分析了？突然发现metareg命令好像可以用于连续变量。请各位战友到俺新开的帖子指教，谢谢！！！>
