ggfortify不支持幸存多个协变量？

Question

通过添加＆＃34; age＆＃34;

，可以从文档示例中修改以下示例

d.coxph <- (survfit(Surv(time, status) ~ sex+age, data = lung))
autoplot(d.coxph)

我会收到以下错误：

  

levels<-中的错误（*tmp*，值= if（nl == nL）as.character（标签）其他paste0（标签，：     因子水平[41]重复

     

输入一个帧编号，或输入0以退出

     

1：autoplot（d.coxph）＆gt; levels<-（*tmp*中的错误，值= if（nl == nL）as.character（标签），否则为paste0（标签，：     因子水平[41]重复

     

输入一个帧编号，或输入0以退出

     

1：autoplot（d.coxph）   2：autoplot.survfit（d.coxph）   3：fortify（object，surv.connect = surv.connect，fun = fun）   4：fortify.survfit（object，surv.connect = surv.connect，fun = fun）   5：因子（rep（groupIDs，model $ strata），levels = groupIDs）

     

2：autoplot.survfit（d.coxph）   3：fortify（object，surv.connect = surv.connect，fun = fun）   4：fortify.survfit（object，surv.connect = surv.connect，fun = fun）   5：因子（rep（groupIDs，model $ strata），levels = groupIDs）

Answer 1

一种可能的解决方案是将连续变量age分为几类，并考虑幸存公式中性别与年龄之间的相互作用：

library(survival)
data(lung)
lung$age2cat <- cut(lung$age,breaks=2)
lung$sex <- factor(lung$sex, labels=c("F","M"))
d.coxph <- survfit(Surv(time, status) ~ interaction(sex,age2cat), data = lung)

autoplot(d.coxph, conf.int=F, surv.size=1)