我在R中有一个cox模型,可以显示患糖尿病的时间。 公式是:
model <- coxph(Surv(days_followed,$outcome_status) ~exposure_status +
+ age + gender + start_follow_date, method="breslow")
我现在想用Ggplot2绘制调整后的生存曲线。两条曲线,一条用于治疗组,一条用于对照组。
所以它就像Kaplan-Meier一样,除了曲线根据年龄,性别,开始跟随日期等进行调整。
有人可以帮我这么做吗?