从嵌套的coxph模型中获取带有扫帚的置信区间和exp

问题描述

数据和库：

test <- tibble(start=c(1,2,5,1,7,3,4,8,8),age=c(2,6,9,14,17),event=c(1,0),x=c(1,sex=c(0,1))

library(tidyverse)
library(broom)
library(survival)

我想嵌套几个分组的小对象并创建coxph对象，并以整洁的眼神​​（从扫帚包中）提取和嵌套数据。在整洁的输出中，我还希望对数据取幂并置信区间。这有效：

coxph_obj <- (coxph(Surv(start,event) ~ x + sex + age,test)) 
tidy(coxph_obj,exponentiate = TRUE,conf.int = TRUE)

但是，我不知道如何让exponentiate = TRUE,conf.int = TRUE在下面的tidied = map(fit,tidy)中工作：

test %>%
  nest(data = -sex) %>% 
  mutate(
    fit = map(data,~ coxph(Surv(start,data = test)),tidied = map(fit,tidy),glanced = map(fit,glance)
  )

 unnest(c(tidied,glanced),names_repair = "universal" )

解决方法

Ben在评论中提供的答案：

“使用tidied = map(fit,tidy,exponentiate = TRUE,conf.int = TRUE)会给您带来mutate的好处”