全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5308 6
2014-05-30
原题如下:
某学者拟分析鼻咽癌与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过程哪里错了呢?请教高手!谢谢!

二维码

扫码加我 拉你入群

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

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

全部回复
2014-5-30 05:23:52
看来,也许教材上的有误。我用本论坛上的MedCalc给你算了一下,结果与R的差不多。
附件列表
1.jpg

原图尺寸 31.14 KB

1.jpg

2.jpg

原图尺寸 203.67 KB

2.jpg

二维码

扫码加我 拉你入群

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

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

2014-5-30 06:16:56
书上的结果也许是不对的。
用Comprehensive Meta Analysis V2加以分析,结果与R的结果一致。
附件列表
1.jpg

原图尺寸 99.97 KB

1.jpg

2.jpg

原图尺寸 108.56 KB

2.jpg

二维码

扫码加我 拉你入群

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

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

2014-5-30 07:59:34
哈哈。。。闲来无事,用Stata做了个。。供参考! QQ图片20140530075741.jpg
附件列表
QQ图片20140530075741.jpg

原图尺寸 50.05 KB

QQ图片20140530075741.jpg

二维码

扫码加我 拉你入群

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

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

2014-5-30 09:26:34
感谢各位的积极试验和回复啊!

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

meta dsl.pdf

大小:183.4 KB

 马上下载

如果要是论坛收论坛币了,请跟我说,我还给您

二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群