数据和库:
test <- tibble(start=c(1,2,5,2,1,7,3,4,8,8),
age=c(2,3,6,7,8,9,9,9,14,17),
event=c(1,1,0,1,1,1,1,0,0,0),
x=c(1,0,0,1,0,1,1,1,0,0),
sex=c(0,0,0,0,0,1,1,1,1,1))
library(tidyverse)
library(broom)
library(survival)
我想嵌套几个分组的tibble,创建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, event) ~ x + sex + age, data = test)),
tidied = map(fit, tidy),
glanced = map(fit, glance)
)
unnest(c(tidied, glanced), names_repair = "universal" )
Ben在评论中提供的答案:
"使用CCD_ 3在您的CCD_;