我有一个如下所示的数据框:
ID Time Censor Region
X1 5 4 NorthAmerica
X2 10 12 Europe
...
我使用autoplot通过以下方式绘制K-M曲线:
survivalFit <- survfit(Surv(survivalPlotData$Time, survivalPlotData$Censor)~Region, data = survivalPlotData)
autoplot(survivalFit)
工作正常。
但是,我还想绘制相同情节的整体KM生存曲线,该曲线不按地区区分。最好的方法是什么?
我找到的唯一方法是将数据框基本上加倍,另一组说“ALL”然后绘制,但这看起来非常低效。