搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  meta dsl.pdf
资料下载链接地址: https://bbs.pinggu.org/a-1564321.html
附件大小:
183.4 KB   举报本内容
原题如下:
某学者拟分析鼻咽癌与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 (lower95% 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.795% 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过程哪里错了呢?请教高手!谢谢!



    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2025-12-30 10:13