r - 如何绘制生存函数的“逆”？

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"or这样的预定义函数"cloglog"，因此您可以轻松地为其赋予您自己的自定义函数。

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