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?

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)). 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

Antworten:

16

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

Dann müssen Sie es für die verschiedenen multiplizieren eβ 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,8245wird aus den zusammenfassenden Ergebnissen der Cox-Regression kopiert und eingefügt summary(reg_fit), was bedeutet, dass die Frauen in den Daten Gefährdungsraten aufweisen18%. 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