Wie berechnet man die vorhergesagten Gefährdungsraten aus einem Cox PH-Modell?

11

Ich habe das folgende Cox PH-Modell:

(Zeit, Ereignis) ~ X + Y + Z.

Ich mag die vorhergesagte Gefahr erhalten Preise (i Raten über Gefahren sprechen bin nicht bestimmte Werte von bestimmten Hazard Ratios) X, Y, Z. Ich weiß, dass das Muhaz R-Paket die beobachteten Gefährdungsraten berechnen kann, aber ich interessiere mich für das vorhergesagte Modell.

Gibt es eine Möglichkeit, dies in R zu tun?

r survival hazard cox-model

— kein Bestand
quelle

2

Es scheint, dass Ihre "Gleichung" ein binäres Ergebnis modelliert. Nun, Cox PH-Modell befasst sich mit Überlebensdaten und was es auf der linken Seite liegt, sollte eine Gefährdungsrate sein ... Ich denke, Sie sollten klarstellen ...

— ocram

2

Das Cox-Modell schätzt die Grundlinien-Gefährdungsfunktion nicht (normalerweise angegeben als

λ_{0} (t)

$\lambda_0(t)$ ). Um die Inzidenzrate aus einem Modell abzuschätzen, sollten Sie stattdessen ein parametrisches Überlebensmodell verwenden.

— Boscovich

1

Ihre Frage enthält einige Unklarheiten. Die Gefahrenfunktion ist eine Wahrscheinlichkeit, keine Rate. Ist es die Gefahrenfunktion, die Sie wollen, oder etwas anderes?

— Fomite

16

Die Funktion im R-Überlebenspaket zum Abrufen der Basisrisikorate ist basehaz.

Dann müssen Sie es für die verschiedenen multiplizieren $e^{\beta}$ um die spezifische Gefährdungsrate anhand der von Ihnen gefundenen Koeffizienten zu erhalten.

Ein einfaches Beispiel kann helfen:

library(survival) #survival analysis
library(eha) #used for data 
data(oldmort) #create the data

# Create surv data set
mort <- Surv(time=oldmort$enter,time2=oldmort$exit,event=oldmort$event)

reg_fit <- coxph(formula=mort~oldmort$sex)
summary(reg_fit)

# Now get baseline curve
baseline <- basehaz(reg_fit)

# Draw baseline hazard (that's male)
plot(baseline$time, baseline$hazard, type='l',main="Hazard rates") 

# Draw female hazard
lines(baseline$time, exp(-0.1929)*baseline$hazard, col="blue")

das $\exp(-0.1929)=0.8245$ wird aus den zusammenfassenden Ergebnissen der Cox-Regression kopiert und eingefügt summary(reg_fit), was bedeutet, dass die Frauen in den Daten Gefährdungsraten aufweisen $18\%$ niedriger als die Männchen.

— CarrKnight
quelle

13

Die Funktion basehaz(aus der vorherigen Antwort) liefert die kumulative Gefahr , nicht die Gefahrenfunktion (die Rate). Ich glaube, diese Frage betraf die Gefahrenfunktion. Die Schätzung der Gefahrenfunktion würde die Angabe der Art der Glättung erfordern (wie bei der Dichteschätzung). Das MuhazR-Paket kann dies für eine Beispieldaten tun. Mir ist keine Funktion bekannt, die dies für das Grundrisiko im Cox-Modell tun kann. Ich brauche das auch. Ich denke, ich muss die kumulative Gefahr selbst glätten.

— user21413
quelle