Modellantaganden

Du har redan bekantat dig med Cox PH-modellen i den förra kodövningen. I den här övningen ska du undersöka om din modell verkligen är lämplig. Modellen är fortfarande lagrad i objektet fitCPH.

Den här övningen är en del av kursen

Maskininlärning för marknadsanalys i R

Visa kurs

Övningsinstruktioner

Kontrollera antagandet om proportionella hazarder för modellen fitCPH med hjälp av cox.zph(). Lagra testresultatet i ett objekt som heter testCPH och skriv ut det.
Antagandet verkar vara brutet för en variabel på signifikansnivån 0,05. Vilken? Plotta koefficienten beta som funktion av tid för den variabeln.
Validera modellen med korsvalidering via funktionen validate() från paketet rms.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# Check proportional hazard assumption and print result
___ <- ___(fitCPH)
___(testCPH)

# Plot time-dependent beta
plot(___, var = "___")

# Load rms package
library(rms)

# Validate model
___(fitCPH, method = "___",
         B = 10, dxy = TRUE, pr = FALSE)

Redigera och kör kod