Słownik metodologiczny · Statystyka opisowa

Obserwacja cenzurowana (ucięta) w analizie przeżycia

Obserwacja cenzurowana to czas do zdarzenia, o którym wiemy tylko, że jest dłuższy niż okres obserwacji. Pokazujemy rodzaje cenzurowania i przykład w R, jak błędne podejście zaniża medianę.

Strona główna / Słownik / Obserwacja cenzurowana
W skrócie: Obserwacja cenzurowana to czas do zdarzenia znany tylko częściowo: wiemy, że zdarzenie nie nastąpiło do końca obserwacji. Nie usuwaj takich osób i nie licz ich jako zdarzeń, tylko użyj metod analizy przeżycia, np. estymatora Kaplana–Meiera.

Co to jest obserwacja cenzurowana?

Obserwacja cenzurowana (ucięta, ang. censored observation) to pomiar czasu do zdarzenia, o którym wiemy tylko tyle, że zdarzenie nie nastąpiło przed pewnym momentem. Osoba skończyła udział w badaniu bez nawrotu, wyprowadziła się albo wycofała zgodę. Jej czas nie jest więc brakiem danych, lecz informacją niepełną: „co najmniej 20 tygodni”. Leung i in. (1997) wskazują, że właśnie cenzurowanie odróżnia analizę przeżycia od innych działów statystyki.

Rodzaje cenzurowania

  • Prawostronne – najczęstsze. Zdarzenie mogło wystąpić dopiero po końcu obserwacji: koniec badania, utrata kontaktu, rezygnacja.
  • Lewostronne – zdarzenie nastąpiło przed pierwszym pomiarem, ale nie wiadomo kiedy, np. pierwszy epizod objawów przed włączeniem do badania.
  • Przedziałowe – wiadomo tylko, że zdarzenie wystąpiło między dwiema wizytami kontrolnymi.

Kluczowe założenie standardowych metod to cenzurowanie nieinformatywne: osoby ucięte mają takie same dalsze rokowanie jak te, które zostały w obserwacji. Przy tym założeniu powstał estymator Kaplana i Meiera (1958). Leung i in. (1997) przypominają, że w badaniach obserwacyjnych mechanizm cenzurowania jest zwykle nieznany, więc to założenie trzeba uzasadnić, a nie przyjmować milcząco.

Przykład

W symulacji 300 osób obserwowano po terapii uzależnienia przez 52 tygodnie, a prawdziwa mediana czasu do nawrotu wynosiła 30 tygodni. Do końca obserwacji 43.3% osób pozostało bez nawrotu lub odeszło wcześniej. Mediana policzona tylko z zaobserwowanych nawrotów wyniosła 16.8 tygodnia, a mediana, w której ucięcia potraktowano jak nawroty – 19.2. Estymator Kaplana–Meiera dał 26.7 tygodnia, najbliżej prawdy; resztę różnicy tłumaczy zmienność losowa próby.

Jak opisać w pracy

Obserwacje osób, u których do 52. tygodnia nie wystąpił nawrót, oraz osób utraconych z obserwacji potraktowano jako prawostronnie cenzurowane (43.3% próby). Medianę czasu do nawrotu oszacowano estymatorem Kaplana–Meiera (26.7 tygodnia), zakładając, że cenzurowanie nie zależało od rokowania.

Na co uważać

  • Nie usuwaj obserwacji uciętych i nie licz ich jako zdarzeń. Oba skróty zaniżają czas do zdarzenia, co widać w przykładzie. Dane z ucięciem wymagają metod analizy przeżycia (Clark i in., 2003).
  • Raportuj, ile osób ucięto i dlaczego, osobno dla grup. Różnicowe odpadanie (zob. utrata uczestników) może łamać założenie nieinformatywności.
  • Cenzurowanie to nie to samo co braki danych. Ucięty czas wnosi do analizy informację; brak wyniku jej nie wnosi.

Porównanie grup z obserwacjami uciętymi umożliwia test log-rank.

Przykład reprodukowalny

Obserwacje cenzurowane w R

Kod symuluje czasy do nawrotu z końcem badania i wcześniejszymi odejściami, a potem porównuje trzy sposoby liczenia mediany z prawdziwą wartością.

Mediana czasu do nawrotu z cenzurowaniem i bez

Dane wejściowe: 300 osób (dane symulowane, set.seed(212)), prawdziwa mediana 30 tygodni, obserwacja do 52 tygodni.

Co odczytać: odsetek obserwacji cenzurowanych i mediany: tylko ze zdarzeń, z cenzurowaniem liczonym jako zdarzenie, z estymatora Kaplana–Meiera.

R
library(survival)
set.seed(212)
n <- 300
t_zd  <- rexp(n, rate = log(2) / 30)                 # prawdziwa mediana: 30 tygodni do nawrotu
t_cen <- pmin(52, rexp(n, rate = 1 / 80))            # koniec badania (52 tyg.) albo wcześniejsze odejście
czas  <- pmin(t_zd, t_cen)
zdarzenie <- as.integer(t_zd <= t_cen)               # 0 = obserwacja cenzurowana
round(mean(zdarzenie == 0), 3)                       # odsetek obserwacji cenzurowanych
round(c(
  tylko_zdarzenia   = median(czas[zdarzenie == 1]),  # błąd 1: usunięcie cenzurowanych
  cenzura_jak_zdarz = median(czas),                  # błąd 2: cenzurowanie = zdarzenie
  kaplan_meier      = unname(summary(survfit(Surv(czas, zdarzenie) ~ 1))$table["median"]),
  prawda            = 30), 1)
Wynik w R
[1] 0.433
  tylko_zdarzenia cenzura_jak_zdarz      kaplan_meier            prawda
             16.8              19.2              26.7              30.0

Bibliografia oprogramowania

Wyniki pochodzą z uruchomienia pokazanego kodu w R 4.5.3 (pakiet survival).

  1. R: R Core Team. (2026). R: A language and environment for statistical computing (Version 4.5.3) [Computer software]. R Foundation for Statistical Computing. https://www.R-project.org/
  2. survival (wersja 3.8.6): Therneau, T. M. (2026). A package for survival analysis in R (R package version 3.8-6) [Computer software]. https://CRAN.R-project.org/package=survival

Najczęstsze pytania

Czy osobę, która wycofała się z badania, traktuje się jako cenzurowaną?

Tak, jej czas obserwacji kończy się w dniu wycofania. Opisz jednak powód odejścia, bo rezygnacja związana ze stanem zdrowia może łamać założenie cenzurowania nieinformatywnego.

Co oznacza cenzurowanie informatywne?

To sytuacja, w której powód ucięcia wiąże się z rokowaniem, np. odchodzą osoby z nasilonymi objawami. Wtedy krzywa przeżycia i test log-rank są obciążone.

Jak zakodować obserwacje cenzurowane w bazie danych?

Dwie kolumny: czas obserwacji (do zdarzenia albo do ucięcia) i wskaźnik zdarzenia, zwykle 1 = zdarzenie, 0 = cenzurowanie. W R łączy je funkcja Surv(czas, zdarzenie).

Powiązane hasła

Poradniki, w których użyjesz tego pojęcia

Potrzebujesz wsparcia w analizie danych?

Wybierz ofertę dopasowaną do pracy dyplomowej albo profesjonalnego projektu naukowego.

Wyceń analizęOferta dla studentów

Bibliografia i dalsza literatura naukowa (APA 7)

  1. Clark, T. G., Bradburn, M. J., Love, S. B., & Altman, D. G. (2003). Survival analysis part I: Basic concepts and first analyses. British Journal of Cancer, 89(2), 232–238. https://doi.org/10.1038/sj.bjc.6601118 ● Otwarty dostęp
  2. Kaplan, E. L., & Meier, P. (1958). Nonparametric estimation from incomplete observations. Journal of the American Statistical Association, 53(282), 457–481. https://doi.org/10.1080/01621459.1958.10501452
  3. Leung, K.-M., Elashoff, R. M., & Afifi, A. A. (1997). Censoring issues in survival analysis. Annual Review of Public Health, 18(1), 83–104. https://doi.org/10.1146/annurev.publhealth.18.1.83

Cytowania BibTeX

@article{clark2003survival, author={Clark, T. G. and Bradburn, M. J. and Love, S. B. and Altman, D. G.}, title={{Survival analysis part I: Basic concepts and first analyses}}, journal={British Journal of Cancer}, year={2003}, volume={89}, number={2}, pages={232--238}, doi={10.1038/sj.bjc.6601118}}
@article{kaplan1958nonparametric, author={Kaplan, E. L. and Meier, Paul}, title={{Nonparametric estimation from incomplete observations}}, journal={Journal of the American Statistical Association}, year={1958}, volume={53}, number={282}, pages={457--481}, doi={10.1080/01621459.1958.10501452}}
@article{leung1997censoring, author={Leung, Kwan-Moon and Elashoff, Robert M. and Afifi, Abdelmonem A.}, title={{Censoring issues in survival analysis}}, journal={Annual Review of Public Health}, year={1997}, volume={18}, number={1}, pages={83--104}, doi={10.1146/annurev.publhealth.18.1.83}}
R