Ipotezele modelului

Ai explorat deja modelul Cox PH în ultimul exercițiu de programare. În acest exercițiu, vei verifica dacă modelul tău este potrivit. Modelul este în continuare stocat în obiectul fitCPH.

Acest exercițiu face parte din cursul

Machine Learning pentru Analiză de Marketing în R

Vezi cursul

Instrucțiuni pentru exercițiu

Verifică ipoteza hazardurilor proporționale pentru modelul fitCPH folosind cox.zph(). Stochează rezultatul testului într-un obiect numit testCPH și afișează-l.
Ipoteza pare să fie încălcată pentru o variabilă la nivelul alpha de 0,05. Care anume? Reprezintă grafic coeficientul beta în funcție de timp pentru această variabilă.
Validează modelul prin validare încrucișată folosind funcția validate() din pachetul rms.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

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

Editează și rulează codul