Dữ liệu bị khuyết góc (censored data / biến cố kiểm duyệt) trong phân tích sống sót là gì?

Là hiện tượng đối tượng nghiên cứu chưa gặp biến cố quan tâm tại thời điểm kết thúc theo dõi, hoặc bị mất dấu theo dõi (lost to follow-up), hoặc tử vong do nguyên nhân cạnh tranh khác; người phân tích chỉ biết đối tượng đã sống sót ít nhất đến thời điểm kiểm duyệt đó.

Đường cong sống sót Kaplan-Meier và phép kiểm Log-rank được áp dụng như thế nào?

Ước lượng Kaplan-Meier là phương pháp phi tham số từng bước tính xác suất tích lũy sống sót theo thời gian; phép kiểm Log-rank so sánh thống kê xem có sự khác biệt có ý nghĩa giữa hai hay nhiều đường cong sống sót của các nhóm điều trị hay không.

Mô hình hồi quy rủi ro tỷ lệ Cox (Cox proportional hazards model) có ý nghĩa gì trong phân tích đa biến?

Mô hình Cox đánh giá tác động đồng thời của nhiều biến số độc lập (tuổi, giai đoạn bệnh, phương pháp điều trị) lên hàm rủi ro (hazard function) mà không cần giả định dạng phân phối cụ thể của thời gian sống sót nền, cho ra tỷ số rủi ro (Hazard Ratio - HR).

Thời gian sống sót là gì? Ý nghĩa và ứng dụng chuyên sâu

Thời gian sống sót (Survival time / Survival analysis) là biến số thời gian liên tục đo khoảng cách từ một thời điểm mốc xác định (như thời điểm chẩn đoán, bắt đầu can thiệp phẫu thuật hoặc dùng thuốc) cho đến khi xuất hiện một biến cố lâm sàng quan tâm (như tử vong, bệnh tái phát hoặc thất bại điều trị) trong nghiên cứu y học và dịch tễ học.

Định nghĩa thời gian sống sót

Thời gian sống sót (survival time) là khoảng thời gian từ thời điểm bắt đầu theo dõi một cá thể cho đến khi xảy ra một sự kiện xác định, thường là tử vong, tái phát bệnh, hoặc một điểm kết thúc nghiên cứu có ý nghĩa lâm sàng. Biến này được sử dụng phổ biến trong nghiên cứu y sinh, đặc biệt trong các thử nghiệm lâm sàng để đánh giá hiệu quả điều trị hoặc dự báo tiên lượng.

Khác với các biến đầu ra truyền thống, thời gian sống sót là một biến ngẫu nhiên liên tục và có đặc điểm quan trọng là có thể bị kiểm duyệt. Thời gian sống sót không chỉ mô tả thời điểm xảy ra sự kiện mà còn phản ánh xác suất sống sót tích lũy theo thời gian, thường được trình bày bằng các hàm toán học và đồ thị để diễn giải một cách trực quan.

Trong các nghiên cứu lâm sàng, thời gian sống sót được sử dụng để tính các chỉ số như thời gian sống trung vị (median survival time), tỷ lệ sống sau 1, 3 hoặc 5 năm, và ước lượng xác suất sống còn tổng thể (overall survival – OS). Những chỉ số này giúp các bác sĩ lâm sàng và nhà nghiên cứu ra quyết định điều trị, theo dõi và đánh giá kết quả can thiệp.

Đặc điểm dữ liệu sống sót

Dữ liệu sống sót có một số đặc điểm thống kê riêng biệt cần được xử lý cẩn trọng trong quá trình phân tích. Đặc điểm nổi bật nhất là sự hiện diện của dữ liệu kiểm duyệt (censoring), thường gặp nhất là kiểm duyệt phải (right censoring), tức là sự kiện chưa xảy ra tại thời điểm kết thúc quan sát.

Ví dụ, nếu một bệnh nhân vẫn còn sống khi kết thúc nghiên cứu, ta không biết chính xác thời gian tử vong của họ mà chỉ biết họ đã sống ít nhất đến thời điểm đó. Đây là thông tin không đầy đủ nhưng vẫn có giá trị thống kê nếu được xử lý đúng cách.

Một số dạng kiểm duyệt phổ biến:

Kiểm duyệt phải: sự kiện chưa xảy ra tại thời điểm kết thúc nghiên cứu.
Kiểm duyệt trái: sự kiện xảy ra trước khi bắt đầu quan sát.
Kiểm duyệt khoảng: sự kiện xảy ra trong một khoảng thời gian nhưng không rõ thời điểm chính xác.

Bên cạnh đó, dữ liệu sống sót còn có thể gặp hiện tượng “cạnh tranh rủi ro” (competing risks), khi có nhiều nguyên nhân dẫn đến kết thúc quan sát và các nguyên nhân này loại trừ lẫn nhau. Điều này yêu cầu mô hình hóa phù hợp để tránh sai lệch ước lượng.

Hàm sống sót và hàm nguy cơ

Hai khái niệm then chốt trong phân tích sống sót là hàm sống sót và hàm nguy cơ, giúp mô tả xác suất còn sống tại một thời điểm cụ thể cũng như nguy cơ xảy ra sự kiện trong thời gian ngắn tiếp theo. Đây là cơ sở để xây dựng các mô hình phân tích định lượng.

Hàm sống sót được định nghĩa là xác suất một cá thể sống sót quá thời điểm $t$ :

$S(t) = P(T > t)$

Nó là hàm giảm dần và luôn nằm trong khoảng từ 0 đến 1. Khi $S(t) = 0.5$ , thời điểm $t$ đó được gọi là thời gian sống trung vị.

Hàm nguy cơ là xác suất điều kiện xảy ra sự kiện tại thời điểm $t$ nếu cá thể còn sống đến thời điểm đó:

$h(t) = \lim_{\Delta t \to 0} \frac{P(t \leq T < t + \Delta t \mid T \geq t)}{\Delta t}$

Hàm nguy cơ có thể diễn giải như “tốc độ chết” hoặc nguy cơ tích lũy. Nó không phải xác suất, mà là mật độ xảy ra sự kiện tại thời điểm cụ thể. Từ hàm nguy cơ, ta có thể khôi phục lại hàm sống sót bằng công thức:

$S(t) = \exp\left(-\int_0^t h(u) \, du\right)$

Phương pháp Kaplan-Meier

Kaplan-Meier là một trong những phương pháp kinh điển và phổ biến nhất để ước lượng hàm sống sót $S(t)$ từ dữ liệu kiểm duyệt. Đây là phương pháp phi tham số, không cần giả định phân phối của thời gian sống và cho phép phân tích từng thời điểm sự kiện xảy ra.

Ước lượng Kaplan-Meier được tính bằng công thức:

$S(t) = \prod_{t_i \leq t} \left(1 - \frac{d_i}{n_i}\right)$

$d_i$ : số sự kiện (ví dụ: tử vong) xảy ra tại thời điểm $t_i$
$n_i$ : số cá thể còn “nguy cơ” (at risk) ngay trước thời điểm $t_i$

Ước lượng được cập nhật tại mỗi thời điểm sự kiện và giữ nguyên giữa các thời điểm. Đường cong Kaplan-Meier là một đường bậc thang giảm dần thể hiện xác suất sống còn theo thời gian. Khoảng tin cậy thường được tính theo phương pháp Greenwood.

Bảng dưới minh họa cách xây dựng Kaplan-Meier cho một mẫu nhỏ:

Thời điểm (t)	Số sống sót trước đó (n)	Số tử vong (d)	S(t)
2	10	1	0.9
5	9	2	0.7
8	7	1	0.6

Phương pháp Kaplan-Meier được sử dụng rộng rãi trong nghiên cứu y học, dược học, và kỹ thuật độ tin cậy, nhờ tính linh hoạt và khả năng trực quan hóa kết quả qua đồ thị sống sót.

Kiểm định so sánh thời gian sống

Khi so sánh thời gian sống giữa hai hoặc nhiều nhóm, ví dụ giữa hai phác đồ điều trị khác nhau, người ta sử dụng các kiểm định thống kê chuyên biệt để đánh giá sự khác biệt giữa các đường cong Kaplan-Meier. Phương pháp phổ biến nhất là kiểm định log-rank, còn gọi là Mantel-Cox test.

Kiểm định log-rank là kiểm định không tham số, so sánh số sự kiện thực tế quan sát được với số sự kiện kỳ vọng tại mỗi thời điểm sự kiện, dựa trên giả định rằng các nhóm có phân phối sống sót giống nhau nếu không có ảnh hưởng điều trị.

Biểu thức kiểm định log-rank:

$Z = \frac{\sum (O_i - E_i)}{\sqrt{\sum V_i}}$

$O_i$ : số sự kiện quan sát ở nhóm i
$E_i$ : số sự kiện kỳ vọng ở nhóm i
$V_i$ : phương sai

Ngoài log-rank, một số kiểm định khác như Breslow (Generalized Wilcoxon) hoặc Tarone-Ware áp dụng trọng số lớn hơn cho các sự kiện sớm, phù hợp khi giả định nguy cơ không hằng định theo thời gian.

Mô hình Cox và hồi quy sống sót

Mô hình Cox là phương pháp hồi quy bán tham số được sử dụng rộng rãi để phân tích ảnh hưởng của các biến giải thích đến nguy cơ xảy ra sự kiện, đồng thời xử lý dữ liệu kiểm duyệt một cách hiệu quả. Đây là mô hình trung tâm trong phân tích sống sót hiện đại.

Biểu thức cơ bản của mô hình:

$h(t|x) = h_0(t) \exp(\beta_1 x_1 + \beta_2 x_2 + \ldots + \beta_p x_p)$

Trong đó:

$h(t|x)$ : hàm nguy cơ có điều kiện theo biến độc lập $x$
$h_0(t)$ : hàm nguy cơ cơ sở (không cần giả định cụ thể)
$\beta_i$ : hệ số hồi quy cần ước lượng

Mô hình Cox không yêu cầu biết trước phân phối của thời gian sống, chỉ giả định tỷ số nguy cơ giữa các cá thể là không đổi theo thời gian (proportional hazards assumption). Nếu giả định này bị vi phạm, cần dùng mô hình mở rộng như Cox theo thời gian hoặc mô hình AFT (Accelerated Failure Time).

Ưu điểm của mô hình Cox:

Ước lượng được ảnh hưởng tương đối của các yếu tố nguy cơ (hazard ratio)
Kiểm soát được nhiều biến đồng thời
Cho phép phân tích dữ liệu kiểm duyệt và đa biến

Ứng dụng trong nghiên cứu lâm sàng

Phân tích sống sót đóng vai trò quan trọng trong nghiên cứu lâm sàng, đặc biệt trong ung thư học, ghép tạng, bệnh truyền nhiễm và điều trị tim mạch. Thời gian sống được sử dụng làm tiêu chí đánh giá chính trong các thử nghiệm thuốc, thiết bị y tế và quy trình điều trị.

Các chỉ số phổ biến trong báo cáo kết quả lâm sàng:

Overall Survival (OS): thời gian từ khi bắt đầu điều trị đến khi tử vong do bất kỳ nguyên nhân nào.
Progression-Free Survival (PFS): thời gian từ khởi trị đến khi bệnh tiến triển hoặc tử vong.
Event-Free Survival (EFS): thời gian sống không có biến cố xác định trước (tái phát, chuyển nặng, tử vong,...)

Bảng dưới minh họa vai trò của thời gian sống trong đánh giá hiệu quả điều trị:

Chỉ số	Ý nghĩa	Ứng dụng
Median OS	Thời gian sống trung vị	So sánh hiệu quả điều trị giữa các nhóm
Hazard Ratio	Tỷ số nguy cơ giữa nhóm can thiệp và chứng	Mô hình Cox
Kaplan-Meier Curve	Đồ thị biểu diễn xác suất sống còn	Phân tích định tính và trực quan

Giải thích kết quả và hạn chế

Việc diễn giải kết quả phân tích sống sót đòi hỏi hiểu rõ bối cảnh lâm sàng và các giả định thống kê. Cần phân biệt tỷ lệ sống thực sự với tỷ lệ sống ước lượng, và kiểm tra các giả định mô hình như proportional hazards trong Cox regression.

Một số hạn chế thường gặp:

Dữ liệu kiểm duyệt không ngẫu nhiên dẫn đến thiên lệch
Thời gian theo dõi ngắn khiến ước lượng không ổn định
Thiếu điều chỉnh nhiễu làm sai lệch kết luận về nhân quả

Trong báo cáo, cần trình bày đầy đủ số lượng mẫu, thời gian theo dõi trung bình, tỷ lệ kiểm duyệt, đồ thị Kaplan-Meier kèm khoảng tin cậy và kết quả kiểm định có p-value.

Phân tích sống sót bằng phần mềm

Nhiều phần mềm thống kê hiện đại hỗ trợ phân tích sống sót toàn diện, từ mô hình Kaplan-Meier đến hồi quy Cox:

R: sử dụng gói survival với hàm survfit(), coxph(). Có thể vẽ đồ thị Kaplan-Meier bằng ggsurvplot() từ gói survminer.
Stata: cung cấp lệnh stset, sts graph và stcox với giao diện trực quan.
SAS: có thủ tục PROC LIFETEST cho Kaplan-Meier và PROC PHREG cho hồi quy Cox.

Các phần mềm này cho phép phân tích sâu dữ liệu lớn, kiểm tra giả định mô hình và xuất báo cáo khoa học đạt chuẩn quốc tế.