Założenia modelu
W poprzednim ćwiczeniu programistycznym zapoznałeś się już z modelem Cox PH. Teraz sprawdzisz, czy twój model jest w ogóle odpowiedni. Model jest nadal zapisany w obiekcie fitCPH.
To ćwiczenie jest częścią kursu
Uczenie maszynowe w analizie marketingowej w R
Instrukcje do ćwiczenia
- Sprawdź założenie proporcjonalności hazardów dla modelu
fitCPH, używając funkcjicox.zph(). Zapisz wynik testu w obiekcie o nazwietestCPHi wyświetl go. - Założenie wydaje się być naruszone dla jednej zmiennej na poziomie istotności alfa = 0,05. Dla której? Narysuj wykres współczynnika beta w zależności od czasu dla tej zmiennej.
- Zwaliduj model za pomocą walidacji krzyżowej, korzystając z funkcji
validate()z pakieturms.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)