Första Kaplan-Meier-skattningen
I den här övningen använder vi samma data som visades i videon. Vi tittar närmare på funktionen survfit() och det objekt den genererar. Övningen hjälper dig att utforska survfit-objektet.
Paketet survival är redan inläst för dig i den här övningen.
Den här övningen är en del av kursen
Överlevnadsanalys i R
Övningsinstruktioner
- Utforska funktionen
survfit()genom att ange?survfiti konsolen. - Beräkna Kaplan-Meier-skattningen med hjälp av
survfit(). - Undersök strukturen hos
survfit-objektet medstr(). - Skapa en
data.framemed de fyra tidpunkterna, motsvarande antal i riskzonen (n.risk), antal observationer med en händelse (n.event), antal censurerade observationer (n.censor) och värdet på överlevnadskurvan (surv). Hämta all information frånsurvfit-objektet.
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Create time and event data
time <- c(5, 6, 2, 4, 4)
event <- c(1, 0, 0, 1, 1)
# Compute Kaplan-Meier estimate
km <- survfit(___(___, ___) ~ ___)
km
# Take a look at the structure
str(___)
# Create data.frame
data.frame(time = km$time, n.risk = ____, n.event = ____,
n.censor = ____, surv = ___)