r-从嵌套的髋关节模型中获得置信区间并用扫帚进行exp



数据和库:

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_;

最新更新