关于R的meta分析，请高高手帮忙分析下

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

原题如下：
某学者拟分析鼻咽癌与EB病毒感染的关系，资料结果见下表
研究        病例组（EB+）        对照组（EB+）        病例组（EB-）        对照组（EB-）
1                       31                          2                              320                         120
2                       20                          3                              72                           205
3.                      31                          3                              79                             51
4.                      62                         10                             57                             30
5.                      43                         12                             60                             32
6.                      55                         18                             49                             17


我用R写的命令如下
library(rmeta)
library(grid)
ai=c(31,20,31,62,43,55)
ci=c(2,3,3,10,12,18)
bi=c(320,72,79,57,60,49)
di=c(120,205,51,30,32,17)
toti=ai+ci
toci=bi+di
m235r=meta.DSL(toti,toci,ai,bi)
m235r
summary(m235r)

结果如下
> m235r
Random effects ( DerSimonian-Laird ) meta-analysis
Call: meta.DSL(ntrt = toti, nctrl = toci, ptrt = ai, pctrl = bi)
Summary OR= 3.7    95% CI ( 1.68, 8.14 )
Estimated random effects variance: 0.69
> summary(m235r)
Random effects ( DerSimonian-Laird ) meta-analysis
Call: meta.DSL(ntrt = toti, nctrl = toci, ptrt = ai, pctrl = bi)
------------------------------------
        OR (lower  95% upper)
[1,]  5.81    1.37      24.66
[2,] 18.98    5.48      65.78
[3,]  6.67    1.94      22.97
[4,]  3.26    1.46       7.27
[5,]  1.91    0.88       4.13
[6,]  1.06    0.49       2.28
------------------------------------
SummaryOR= 3.7  95% CI ( 1.68,8.14 )
Test for heterogeneity: X^2( 5 ) = 19.23 ( p-value 0.0017 )
Estimated random effects variance: 0.69
可是教材上的答案却是OR=3.07，95%CI（1.38，6.82）。我又根据书上的D-L算法一步步算了，与书上的结果一致，那么我这个R过程哪里错了呢？请教高手！谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：meta分析 Meta 高高手 ETA estimated EB病毒 鼻咽癌 资料

看来，也许教材上的有误。我用本论坛上的MedCalc给你算了一下，结果与R的差不多。

书上的结果也许是不对的。
用Comprehensive Meta Analysis V2加以分析，结果与R的结果一致。

哈哈。。。闲来无事，用Stata做了个。。供参考！

感谢各位的积极试验和回复啊！

问题是我是用教材上的原理一步一步演算的结果，跟教材答案一样，而跟我们的不一样。我把教材meta D+L原理过程图片附上来，大家帮忙看看是不是有漏洞或错误。再次谢谢大家了！

这个题目和卫生统计学,人民卫生出版社,第7版,第20章,例20-6一样.这个教材在计算校正因子时,使用的权重是wi=(bi*ci)/Ti,这个权重是Mantel-Haenszel法计算OR值时赋于的权重,我感觉应该使用inverse-variance法的权重计算方法,也就是wi=1/[(1/ai+1/bi+1/ci+1/di)].

jzdfenger 发表于 2015-12-17 23:21
这个题目和卫生统计学,人民卫生出版社,第7版,第20章,例20-6一样.这个教材在计算校正因子时,使用的权重是wi= ...

谢谢你的解答！