R Überlebenspaket; Log (log (Überleben)) gegen Log (Zeit)

Ich versuche, das R-Survival-Paket zu verwenden, um ein Diagramm von log(-log(survival)) gegen Log (Zeit)

_{(Dies wird manchmal als eine Möglichkeit zur visuellen Prüfung auf beschleunigte Lebensdauer oder proportionale Gefahreneigenschaften empfohlen).}

Die Option " fun=cloglog " in plot.survfit produziert nicht das, was ich von ihr erwarte.

Verwendung der "Gehan" -Daten aus der "MASS" -Bibliothek:

Erstens, hier ist ein einfaches Diagramm des Überlebens gegen die Zeit, für die Behandlung und Kontrollgruppen:

OK bis jetzt. Wenn ich nun die Option fun=cloglog verwende, dann lässt die Dokumentation von plot.survfit mich denken, dass ich eine Zeichnung von log(-log(survival)) gegen log(time) bekommen werde:

  

Spaß

     

eine beliebige Funktion, die eine Transformation der Überlebenskurve definiert. Zum Beispiel fun = log   ist eine alternative Möglichkeit, eine logarithmische Überlebenskurve zu zeichnen (aber die Achse ist mit log (S) gekennzeichnet)   values), und fun = sqrt würde eine Kurve auf der Quadratwurzelskala erzeugen. Vier oft benutzt   Transformationen können stattdessen mit einem Zeichenargument angegeben werden: "log" ist dasselbe wie   Mit der Option log = T zeichnet "Ereignis" kumulative Ereignisse (f (y) = 1-y), "cumhaz" zeichnet die   kumulative Hazard-Funktion (f (y) = -log (y)) und "cloglog" erstellt ein komplementäres log-log   Überlebensdiagramm (f (y) = log (-log (y)) zusammen mit der logarithmischen Skala für die x-Achse).

Wenn ich dies versuche, scheint es jedoch nicht die log(-log(y)) -Funktion zu verwenden, da die angezeigte Kurve immer noch abnimmt (da die ursprüngliche Überlebenskurve abnimmt und die angewendete f(y)=log(-log(y)) -Funktion eine abnehmende Funktion ist , die resultierende log(-log(survival)) -Kurve sollte zunehmen).

Außerdem ist die x-Achse nicht logarithmisch skaliert:

Ich kann bekommen, was ich will, indem ich meine eigene log(-log()) Funktion definiere und die Option log="x" verwende:

Also: Was mache ich oben falsch (oder wie interpretiere ich die plot.survfit -Dokumentation falsch)?

Ergänzende Frage: Wie würde eine Option "fun=" die Skalierung auf der horizontalen Achse ändern, wie die Dokumentation behauptet, dass "fun=cloglog" , wenn das Argument auf "fun" eine Funktion ist, die auf die vertikale Variable angewendet wird ?