如何解决从嵌套的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
的好处”
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。