Assunzioni del modello
Hai già dato un'occhiata al modello Cox PH nell'ultimo esercizio di coding. In questo esercizio, scoprirai se il tuo modello è appropriato. Il tuo modello è ancora salvato nell'oggetto fitCPH.
Questo esercizio fa parte del corso
Machine Learning per il Marketing Analytics in R
Istruzioni dell'esercizio
- Verifica l'assunzione di rischi proporzionali del modello
fitCPHusandocox.zph(). Salva il risultato del test in un oggetto chiamatotestCPHe stampalo. - L'assunzione sembra violata per una variabile al livello alfa 0.05. Quale? Traccia il coefficiente beta in funzione del tempo per questa variabile.
- Valida il modello usando la cross validation con la funzione
validate()dal pacchettorms.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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)