7.3. Ước lượng Mantel-Haenszel (MH) cho Hiệu số Nguy cơ (RD)
Giới thiệu:
Ước lượng Mantel-Haenszel (MH) cho hiệu số nguy cơ (RD) là một phương pháp tổng hợp các ước lượng RD riêng theo từng tầng thành một ước lượng chung, dựa trên nguyên lý trung bình có trọng số.
Công thức ước lượng MH cho RD:
RD̂MH = (Σj Rj − Σj Sj) / (Σj Tj) = (R• − S•) / T• (7.5)
trong đó:
- Rj = a1jr2j / rj
- Sj = a2jr1j / rj
- Tj = r1jr2j / rj
(Theo Greenland và Robins, 1985b).
Có thể dễ dàng chứng minh rằng:
RD̂MH = (1/T•) × Σj Tj × RD̂j
nghĩa là RD̂MH là trung bình có trọng số của các ước lượng RD theo tầng, với trọng số Tj = r1jr2j/rj.
Phương sai của ước lượng MH:
Sato (1989) đưa ra ước lượng phương sai cho RD̂MH có hiệu lực trong cả điều kiện mẫu lớn (large-strata) và mẫu thưa (sparse-strata):
vâr(RD̂MH) = [RD̂MH × U• + V•] / (T•)² (7.6)
trong đó:
- Uj = [r²1ja2j − r²2ja1j + r1jr2j(r2j − r1j)/2] / r²j
- Vj = (a1jb2j + a2jb1j) / (2rj)
- U• = Σj Uj
- V• = Σj Vj
Khoảng tin cậy:
[RDMH, RD̅MH] = RD̂MH ∓ zα/2 × √[vâr(RD̂MH)]
Trường hợp đặc biệt — một tầng duy nhất:
Khi chỉ có một tầng (J = 1), công thức (7.5) và (7.6) đơn giản hóa thành (7.1) và (7.2) — các công thức cho bảng 2×2 đơn.
Ví dụ 7.3 (Mức độ thụ thể – Ung thư vú):
Áp dụng các phương pháp MH cho dữ liệu ung thư vú:
- RD̂MH = (14,79 − 9,14) / 34,05 = 0,166
- vâr(RD̂MH) = [0,166 × 1,12 + 6,49] / (34,05)² = (0,0759)²
- Khoảng tin cậy 95%: [0,0171; 0,315]
Kết quả này rất gần với ước lượng từ phương pháp tiệm cận không điều kiện (AU) ở Mục 7.2 (RD̂ = 0,161; KTC 95%: [0,013; 0,310]), cho thấy tính nhất quán giữa các phương pháp.
(Nguồn: Newman SC. Biostatistical Methods in Epidemiology. Wiley, 2001. Chương 7, tr.155-156.)