从类“summary.coxph”的对象中提取摘要效果答案

【问题标题】：Extract summary effect from object of class "summary.coxph"从类“summary.coxph”的对象中提取摘要效果
【发布时间】：2016-04-25 06:06:15
【问题描述】：

我有“summary.coxph”类的对象

Fixed-effects meta-analysis
Call: meta.summaries(d = A[[i]]$hr, se = A[[i]]$se, method = c("fixed"), 
    names = A[[i]]$metabolite, conf.level = 0.95)
Summary effect=1.09   95% CI (1.05, 1.13)
Estimated heterogeneity variance: 0  p= 0.398

我需要提取Summary effect、95%CI 和p。但我所有的尝试都失败了。你能帮帮我吗？

我的数据如下所示：

$             :List of 14
  ..$ effects        : num [1:3] 1.07 1.07 1.14
  ..$ stderrs        : num [1:3] 0.0298 0.0421 0.0436
  ..$ summary        : num 1.09
  ..$ se.summary     : num 0.0213
  ..$ test           : num [1:2] 51.2 0
  ..$ het            : num [1:3] 1.841 2 0.398
  ..$ call           : language meta.summaries(d = A[[i]]$hr, se = A[[i]]$se, method = c("fixed"), names = A[[i]]$metabolite, conf.level = 0.95)
  ..$ names          : chr [1:3] "XXL.VLDL.TG" "XXL.VLDL.TG" "XXL.VLDL.TG"
  ..$ tau2           : num 0
  ..$ variance.method: chr "fixed"
  ..$ weights        : num [1:3] 1124 563 525
  ..$ weight.method  : chr "fixed"
  ..$ conf.level     : num 0.95
  ..$ logscale       : logi FALSE
  ..- attr(*, "class")= chr "meta.summaries"
 - attr(*, "class")= chr "summary.coxph"

感谢您的帮助。

【问题讨论】：

看rmeta::print.meta.summaries的最后几行

标签： r

【解决方案1】：

根据@rawr 的评论，你会发现：

给定一个summary.coxph对象x，使用x$het[3]作为p值， x$summary 用于总结效果，以及 x$summary + qnorm(c(0.025, 0.975)) * x$se.summary 表示 95% CI。

【讨论】：