如何繪製生存函數的「逆」？

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」（與x軸數刻度沿F（Y）=日誌（-log（Y）））創建了一個免費的雙對數存活圖

此外，fun參數不限於預先定義的功能，如"event"或"cloglog"，這樣你就可以輕鬆地給它自己的自定義功能。

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