4

データは実際には時間の経過に伴うイベントの割合の増加であるため、生存関数の逆数をプロットしようとしています。カプラン・マイヤー生存プロットを作成できますが、これらの「反対」を作成したいと考えています。次を使用して、必要なものを取得できます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 軸の上に伸びている理由がよくわかりません。

これについて少し助けていただければ幸いです。

乾杯

4

1 に答える 1

12

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))
于 2013-04-05T11:52:21.563 に答える