První odhad Kaplan-Meier

V tomto cvičení použijeme stejná data jako ve videu. Podíváme se na funkci survfit() a na objekt, který vytváří. Cvičení ti pomůže lépe porozumět struktuře objektu survfit.

Balíček survival je v tomto cvičení již načtený.

Toto cvičení je součástí kurzu

Survival Analysis in R

Zobrazit kurz

Pokyny k cvičení

Prozkoumej funkci survfit() zadáním ?survfit do konzole.
Pomocí survfit() vypočítej odhad Kaplan-Meier.
Prohlédni si strukturu objektu survfit pomocí str().
Vytvoř data.frame se čtyřmi časovými body a odpovídajícími hodnotami: počet ohrožených (n.risk), počet pozorování s událostí (n.event), počet cenzorovaných pozorování (n.censor) a hodnota křivky přežití (surv). Všechny informace čerpej z objektu survfit.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

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

Upravit a spustit kód