【问题标题】:Retrieve data from a summary in R从 R 中的摘要中检索数据
【发布时间】: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 &lt;- metaprop(4:1, c(10, 20, 30, 40)); str(summary(m1))
  • 我已经修改了问题和数据,现在可以工作了!谢谢!!

标签: r summary


【解决方案1】:

如果没有可重现的示例很难,但我建议不要使用 $ 作为提取运算符,而是尝试使用双括号 ([[]]) 拉出特定对象,然后使用单括号 (@ 987654323@) 将特定项目从该对象中拉出。

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2020-10-13
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2021-11-09
    • 1970-01-01
    相关资源
    最近更新 更多