14강 / 20강
로그순위 검정 — 곡선의 차이를 숫자 하나로
관찰 − 기대, 그리고 z
두 생존곡선이 다른지 어떻게 판단할까요? "곡선 사이의 면적" 같은 것을 재고 싶어지지만, 중도절단 때문에 곤란합니다.
로그순위 검정의 발상은 다릅니다. 죽음이 일어난 시각마다 이렇게 묻습니다 —
"이 순간 살아 있던 사람이 100명이고 그중 40명이 A군이었다. 그런데 한 명이 죽었다. 두 군이 다르지 않다면 A군에서 죽었을 기대값은 0.4명이다. 실제로는 A군에서 죽었나, B군에서 죽었나?"
이 관찰 − 기대 를 모든 사망 시각에 대해 쌓습니다. 쌓인 값이 0에서 멀면 두 군이 다르다는 뜻이에요.
침윤이 높은 쪽에서는 73명이 죽었습니다 — 차이가 없었다면 102.5명이었을 자리예요. 기대보다 약 30명 적게 죽었습니다.
Python
브라우저에서 실행
12개 암종 중 p < 0.05 가 세 개. 방향은 대부분 음수지만 두 곳은 양수입니다.
Python
브라우저에서 실행
결과를 정직하게 읽어 봅시다.
- SKCM(흑색종) — z = −4.38, p = 1.2×10⁻⁵. 강합니다. 흑색종에서 TIL 이 예후 인자라는 것은 오래전부터 알려진 사실이고, 이 자료도 그렇게 말합니다.
- CESC·BRCA — p 가 0.004, 0.017. 같은 방향입니다.
- 나머지 아홉 — p 가 0.2 ~ 0.87. 아무것도 말하지 못합니다.
- PRAD·PAAD·UCEC — z 가 양수 입니다. 방향이 반대예요.
UVM(포도막흑색종)은 아예 표에 없습니다. 침윤률의 절반 이상이 0이라 중앙값에서 두 무리로 가를 수가 없기 때문이에요. 자를 수 없는 자료도 있습니다.
여기서 멈추면 안 됩니다. 우리는 방금 12번 검정했습니다. "세 개가 유의하다" 는 말이 무슨 뜻인지, 16강에서 따져야 합니다.
🤔 확인 문제
로그순위 검정에서 '기대 사망자 수' 는 어떤 가정 아래 계산됩니까?
해설 귀무가설(두 군이 다르지 않다) 아래에서는 그 순간 죽을 사람이 각 군에 위험집합의 비율대로 나뉘어야 합니다. 그것이 기대값입니다.
🤔 확인 문제
z = −4.38 (침윤 높은 쪽을 1군으로) 을 옳게 읽은 것은?
해설 z = (관찰 − 기대)/√분산 이므로 음수는 관찰이 기대보다 작다는 뜻입니다. 1군(침윤 높음)에서 죽음이 덜 일어났다 — 즉 더 오래 살았다는 말입니다.