如何绘制生存函数的"逆"？

Question

我试图绘制生存函数的倒数,因为数据实际上是一个事件随时间的比例增加.我可以生产Kaplan-Meier生存地块,但我希望生成这些的"相反".我可以通过以下方式得到我想要的东西fun="cloglog":

plot(survfit(Surv(Days_until_workers,Workers)~Queen_Number+Treatment,data=xdata),
     fun="cloglog", lty=c(1:4), lwd=2, ylab="Colonies with Workers",
     xlab="Days", las=1, font.lab=2, bty="n")

但我完全不明白它对时间的影响(即不是从0开始,距离减少？),以及为什么生存线延伸到y轴以上.

真的很感激一些帮助!

干杯

Answer 1

使用fun="event"以获得所需的输出

fit <- survfit(Surv(time, status) ~ x, data = aml)
par(mfrow=1:2, las=1)
plot(fit, col=2:3)
plot(fit, col=2:3, fun="event")

fun="cloglog"拧紧轴的原因是它根本不绘制一个分数.相反,它是根据以下情况绘制的?plot.survfit:

"cloglog"创建一个免费的对数 - 对数生存图(f(y)= log(-log(y))以及x轴的对数标度)

此外,fun参数不仅限于"event"或等预定义函数"cloglog",因此您可以轻松地为其提供自己的自定义函数.

plot(fit, col=2:3, fun=function(y) 3*sqrt(1-y))