6.1 Các Phương Pháp Xấp Xỉ Tiệm Cận Không Điều Kiện cho Bảng 2 × 2 Đơn Lẻ
Các số đếm quan sát được cho phân tích không phân tầng (unstratified analysis) được trình bày trong Bảng 4.1. Thực hiện phép thay thế π1 = RR · π2, hàm xác suất đồng thời (4.1) có thể được tái tham số hóa để thu được một hàm hợp lý (likelihood) phụ thuộc vào các tham số RR và π2. Điều này dẫn đến các phương trình hợp lý tối đa không điều kiện (unconditional maximum likelihood equations):
a1 = RR̂ · π̂2 · r1
và
(a1 − RR̂ · π̂2 · r1) / (1 − RR̂ · π̂2) + (a2 − π̂2 · r2) / (1 − π̂2) = 0
trong đó RR̂ ký hiệu ước lượng hợp lý tối đa không điều kiện của RR. Các nghiệm là:
RR̂ = π̂1 / π̂2 = (a1 · r2) / (a2 · r1) (6.1)
và
π̂2 = a2 / r2
trong đó π̂1 = a1 / r1. Ước lượng hợp lý tối đa không điều kiện của phương sai của log(RR) [var(log RR)] là:
vâr(log RR̂) = (1 − π̂1) / (π̂1 · r1) + (1 − π̂2) / (π̂2 · r2) = b1 / (a1 · r1) + b2 / (a2 · r2) (6.2)
và khoảng tin cậy (1 − α) × 100% cho RR thu được bằng cách lấy lũy thừa (exponentiating):
(log RR, log RR) = log(RR̂) ± zα/2 · √(b1 / (a1 · r1) + b2 / (a2 · r2))
Nếu a1 hoặc a2 bằng 0, chúng ta thay thế (6.1) và (6.2) bằng:
RR̂ = (a1 + 0.5) · r2 / [(a2 + 0.5) · r1]
và
vâr(log RR̂) = (b1 + 0.5) / [(a1 + 0.5) · r1] + (b2 + 0.5) / [(a2 + 0.5) · r2]
Vì π1 = π2 tương đương với log(RR) = 0, giả thuyết không có mối liên hệ (no association) có thể được biểu diễn dưới dạng H0: log(RR) = 0. Dưới H0, một ước lượng của var(log RR̂) là:
vâr0(log RR̂) = f̂1 / (ê1 · r1) + f̂2 / (ê2 · r2) = (r · m2) / (r1 · r2 · m1)
Kiểm định Wald về mối liên hệ (Wald test of association) là:
Xw2 = (log RR̂)2 · (r1 · r2 · m1) / (r · m2) (df = 1)
và kiểm định tỉ số hợp lý (likelihood ratio test) về mối liên hệ chính xác là công thức (4.12).
Ví dụ 6.1 (Mức độ thụ thể – Ung thư vú)
Dữ liệu cho ví dụ này được lấy từ Bảng 4.5(a). Ước lượng của tỉ số nguy cơ là RR̂ = 2,23, vâr(log RR̂) = 0,048, và khoảng tin cậy 95% cho RR là [1,45; 3,42]. Kiểm định Wald về mối liên hệ cho kết quả Xw2 = 9,02 (p = 0,003) và, từ Ví dụ 4.2, kiểm định tỉ số hợp lý về mối liên hệ là Xlr2 = 11,68 (p = 0,001). Cách diễn giải các kết quả này tương tự như trong Ví dụ 4.2.
(Nguồn: Newman SC. Biostatistical Methods in Epidemiology. Wiley, 2001. tr.143-144.)