我有一个与ggsurvplot的动态标签和标题有关的一般问题。
使用来自包purrr的map
在嵌套列表中使用时,是否可以为ggsurvplot函数提供动态图例标签和图例标题?
我使用的数据集是NCCTG肺癌数据lung{survival}
,经过略微调整的数据包括ECOG得分低于3
下面的代码为每个预后因素(性别和ECOG评分,在变量名中定义)动态创建生存图。但是,我无法弄清楚如何向title
函数中的参数添加动态legend.labs
,legend.title
和ggsurvplot
,因为它们需要作为字符向量提供应该随每个变量名(性别和ph.ecog)而变化。
任何帮助将不胜感激。
谢谢!
HTJ
varnames <- c("sex", "ph.ecog")
lung <- lung %>% filter(ph.ecog<3)
res.tabs <- lung %>%
gather(key=VAR, value, varnames) %>%
nest( data = -VAR) %>% mutate(
fit = map(data, ~ coxph(Surv(time, status) ~ value, data = .x)),
survfit = map(data, ~ survfit(Surv(time, status) ~ value , data = .x, conf.type ="log")),
plots = ggsurvplot(survfit, data, conf.int = TRUE, title = "???", legend.labs = "???", legend.title = "???"))