ggcoxadjustedcurves是计算cox调整下的生存曲线,不同于单纯生存曲线。目前运行代码如下:
library("survival")
library("survminer")
lung$sex <- ifelse(lung$sex == 1,"Male", "Female")
fit <- coxph(Surv(time, status) ~sex + ph.ecog + age,data = lung)
ggcoxadjustedcurves(fit, data=lung,variable=lung$sex)
但是出现: