【发布时间】:2018-01-01 19:20:39
【问题描述】:
我正在使用包 meta 在 R 中进行元分析。
我正在处理的数据是这样的:
dat1 <- data.frame(E1 = c(1, 1, 34, 24),
N1 = c(45, 34, 47, 34),
E2 = c(23, 12, 12, 56),
N2 = c(34, 45, 36, 40),
Group = c('A', 'A', 'B', 'B'))
dat1
E1 N1 E2 N2 Group
1 45 23 34 A
1 34 12 45 A
34 47 12 36 B
24 34 56 40 B
使用的代码是:
library(meta)
metapr <- metaprop(event = E1, n = N1, data = dat1, byvar = Group, comb.fixed = FALSE)
ss <- summary(metapr)
我需要使用摘要中的一些数据创建一个表,但我无法使用以下语法检索它:
ss$......
我需要以下值:
y <- c('Variable', 'Number studies','Pooled proportion', '95% CI', 'I^2%',
'p-value') #for the study
x1 <- c('Variable', 'Number studies','Pooled proportion', '95% CI', 'I^2%')
#for each subgroup
x2 <- c('Variable', 'Number studies','Pooled proportion', '95% CI', 'I^2%')
#for each subgroup
结果是:
y <- c('E1', 4, 0.2534, '[0.0501; 0.6861]',91.8, 0.0001)
x1 <- c('E1A', 2,0.0256, '[0.0064; 0.0965]', 0.0)
x2 <- c('E!B', 2,0.7160, '[0.6086; 0.8034]', 0.0)
谢谢!!
【问题讨论】:
-
如果您向reproducible example 提供示例输入数据和所需输出,则更容易为您提供帮助。这样我们就可以运行和测试代码了。
-
E1 = 1ER 和 N1 = 1N 吗?
-
是的,对不起。我已经修改了。
-
您的代码包含一个或多个错误,因此无法重现该示例。不过,这段代码可能会给你一些关于你感兴趣的对象结构的想法:
library(meta);m1 <- metaprop(4:1, c(10, 20, 30, 40));str(summary(m1)) -
我已经修改了问题和数据,现在可以工作了!谢谢!!