Ước tính cỡ mẫu cho nghiên cứu y học

1. Ước tính cỡ mẫu cho nghiên cứu y học Tuan V. Nguyen Garvan Institute of Medical Research Sydney

2. Nội dung • Lí thuyết về ước tính cỡ mẫu • Thông số cần thiết • Ví dụ về ước tính cỡ mẫu

3. Qui trình phản nghiệm (falsificationism) • Qui trình phản nghiệm: • Bước 1: đề ra giả thuyết • Bước 2: thu thập dữ liệu để phản nghiệm giả thuyết (chứng minh rằng giả thuyết sai) • Thống kê: • Bước 1: phát biểu giả thuyết Ho (không có hiệu quả, không có liên hệ, không khác nhau giữa hai nhóm, v.v…) • Bước 2: thu thập dữ liệu D • Bước 3: ước tính xác suất dữ liệu D nếu Ho đúng.

4. Sai sót loại I và II S = “significant” (tức p<0.05); NS = “non-significant” (tức p>0.05). Power = 1 – b = P(S | Ha) b = P(NS | Ha); a = P(S | Ho).

5. Tương quan giữa chẩn đoán y khoa và suy luận trong khoa học

6. Bảng 1. Hằng số liên quan đến sai sót loại I và II

7. Các thông số cần thiết • Mô hình nghiên cứu • Nghiên cứu cắt ngang (cross-sectional study) • Nghiên cứu xuôi thời gian (prospective study) • Nghiên cứu bệnh chứng (case-control study) • Kiểm tra (survey) • Sai lầm I và II • Tiêu chí lâm sàng • Biến phụ thuộc chính (primary dependent variable hay outcome hay endpoint) • Hệ số ảnh hưởng (effect size) • Mức độ ảnh hưởng mà nhà nghiên cứu biết trước hay nằm trong giả thuyết.

8. Tiêu chí lâm sàng • Là biến số mà nghiên cứu đặt trọng tâm để thu thập và phân tích. • Ví dụ: • Nghiên cứu hiệu quả của thuốc chống gãy xương thì tiêu chí lâm sàng có thể là mức độ thay đổi mật độ xương (bone mineral density) hay tỉ lệ gãy xương (incidence of fracture) • Nghiên cứu hiệu quả của rosiglitazone đến bệnh đái tháo đường, tiêu chí có thể là lượng đường trong máu. • Nghiên cứu hiệu quả của rosiglitazone đến bệnh nhồi máu cơ tim, tiêu chí có thể là tỉ lệ nhồi máu cơ tim. • Nghiên cứu ước tính tỉ lệ hút thuốc lá thì tiêu chí cũng chính là tỉ lệ này.

9. Khái niệm “effect size” • Một cách đơn giản, là một chỉ số về độ ảnh hưởng của một thuật can thiệp. • Vấn đề là đơn vị đo lường của tiêu chí lâm sàng khác nhau quá. • Dẫn đến “standardized effect size” – hệ số ảnh hưởng chuẩn hóa. • Chuẩn hóa lấy độ lệch chuẩn (standard deviation) làm đơn vị đo lường chính! X1 là số trung bình của nhóm can thiệp; X0 là số trung bình của nhóm đối chứng; và S0 là độ lệch chuẩn của nhóm đối chứng.

10. Vài ví dụ về “effect size” • Ví dụ 1: Nghiên cứu về bệnh cao huyết áp: bệnh nhân được điều trị bằng một thuốc trong nhóm beta-blocker. • Trước khi điều trị, SBP: 140 mmHg, độ lệch chuẩn (SD): 22 • Sau khi điều trị: 125 mmHg. • Ví dụ 2: Nghiên cứu loãng xương. • MĐX (mật độ xương) lúc ban đầu: 0.68 g/cm2, SD = 0.12 g/cm2. • 6 tháng sau điều trị: MĐX = 0.72 g/cm2, SD=0.13 g/cm2. ES cho nghiên cứu 1: ES = (140 – 125) / 22 = 0.68 ES cho nghiên cứu 2: ES = (0.72 – 0.68) / 0.12 = 0.33

11. Diễn dịch Effect size ES = 0.2 tương đương với độ khác biệt về chiều cao của một em bé 15 tuổi và một em bé 16 tuổi. ES = 0.5 tương đương với độ khác biệt về chiều cao của một em bé 14 tuổi và một em bé 18 tuổi. ES = 0.8 tương đương với độ khác biệt về chỉ số thông minh (IQ) của một sinh viên năm thứ nhất và một tiến sĩ. Đây là ý kiến của các vị tâm lí học – không phải của người soạn bài giảng!!!

12. Công thức cỡ mẫu chung • C là hằng số liên quan đến sai sót I và II trong bảng số 1 trình bày slide trước. • ES là hệ số ảnh hưởng • Nếu nghiên cứu có 2 nhóm: tổng số đối tượng sẽ là 2n. • Nhận xét: • ES càng thấp, số lượng cỡ mẫu càng lớn. • Sai sót I và II càng thấp, số lượng cỡ mẫu càng lớn.

13. Nghiên cứu một nhóm: ước tính chỉ số trung bình • Mục tiêu: ước tính chiều cao ở đàn ông người Việt • Thông số: • Khoảng tin cậy 0.95 (tức a=0.05) và power = 0.8 (hay b = 0.2). Do đó, hằng số C = 7.85 (Bảng 1). • Chấp nhận sai số trong vòng 1 cm (d = 1) • Độ lệch chuẩn chiều cao ở người Việt khoảng 4.6 cm. • Hệ số ảnh hưởng là: ES = 1/4.6 = 0.217 • Số cỡ mẫu cần thiết:

14. Nghiên cứu một nhóm: ước tính tỉ lệ lưu hành (prevalence rate) • Mục tiêu: ước tính tỉ lệ loãng xương ở phụ nữ Việt tuổi 60+. • Thông số: • Tỉ lệ (trong y văn): 20% (tức p = 0.20) • Sai số ước tính khoảng 2% (m = 0.02) • Sai sót I = 5% (tức khoảng tin cậy 95%) • Số cỡ mẫu cần thiết:

15. Nghiên cứu so sánh 2 số trung bình • Mục tiêu: Thử nghiệm thuốc alendronate trong việc điều trị loãng xương ở phụ nữ sau thời kì mãn kinh. Hai nhóm bệnh nhân được tuyền: nhóm can thiệp và nhóm đối chứng. • Thông số: • Tiêu chí để đánh giá hiệu quả của thuốc là mật độ xương (BMD). • BMD sau thời kì mãn kinh: 0.80 g/cm2, SD 0.12 g/cm2. • Giả thuyết: sau 12 tháng MĐX tăng 5% (tức 0.80 x 1.05 = 0.84) • ES = (0.84 – 0.80) / 0.12 = 0.3333 • Power = 0.90, alpha = 0.05, C = 10.51 • Số cỡ mẫu cho mỗi nhóm:

16. Nomogram để ước tính cỡ mẫu cho 2 nhóm POWER ES Cách sử dụng: Lấy ví dụ trước, chúng ta có ES = 0.33, power = 0.9. Đánh dấu 0.33 ở cột trái, 0.9 cột phải; Kẻ nối hai điểm đã đánh dấu bằng một thước thẳng. Điểm giao chéo giữa đường kẻ thẳng và cột giữa chính là số cỡ mẫu cần thiết cho a = 0.05 hay 0.01.

17. Nghiên cứu so sánh 2 tỉ lệ • Mục tiêu: Nghiên cứu lâm sàng đánh giá hiệu quả của một loại thuốc chống gãy xương sống. Hai nhóm bệnh nhân sẽ được tuyển: can thiệp và đối chứng. • Thông số: • Tiêu chí để đánh giá hiệu quả của thuốc là tỉ lệ gãy xương đốt sống. • Tỉ lệ gãy xương trong nhóm chứng là khoảng 10%, và thuốc có thể làm giảm tỉ lệ này xuống còn khoảng 6%. • power = 0.90, alpha = 0.01. • Số cỡ mẫu ?

18. Lí thuyết • Gọi tỉ lệ 2 nhóm là p1 và p2. • D = | p1 + p2 | • p = (p1 + p2) / 2 • Za/2 là hằng số cho sai sót loại I • Zb là hằng số cho sai sót loại II (power) • Công thức ước tính cỡ mẫu là:

19. Nghiên cứu so sánh 2 tỉ lệ • Thông số: • p1 = 0.10, p2 = 0.06, D = 0.04 • p = (0.10 + 0.06) / 2 = 0.08 • power = 0.90, alpha = 0.01. So, Z0.01/2 = 2.57, Z0.9 = 1.28 • Số cỡ mẫu: Công trình nghiên cứu này cần phải tuyển ít nhất là 2722 (1361 x 2) bệnh nhân để kiểm định giả thiết trên.

20. Nghiên cứu bệnh chứng • Mục tiêu: Tìm hiểu mối liên hệ giữa hút thuốc lá và nguy cơ gãy xương cột sống (vertebral fracture). • Thiết kế: Nhóm 1 là những bệnh nhân mới gãy xương cột sống, và nhóm 2 là những đối tượng không gãy xương. • Thông số: • Khoảng 25% đàn ông hút thuốc lá. p = 0.25 • Odds ratio giữa hút thuốc lá gãy xương: OR = 2 • Sai sót I là a = 0.05 và power = 0.80

21. Nghiên cứu bệnh chứng • Thông số: • p = 0.25 • OR = 2 • C = 7.85 (sai sót I là a = 0.05 và power = 0.80) Công trình nghiên cứu cần phải tuyển ít nhất là 350 đối tượng (175 bệnh nhân và 175 đối chứng) để kiểm định giả thiết trên.

22. Các phương pháp ước tính cỡ mẫu khác • Để “chứng minh” hai nhóm tương đương nhau • Để ước tính hệ số tương quan • Nghiên cứu chẩn đoán • Nghiên cứu về sống sót (survival study)

23. Những vấn đề cần lưu ý • Điều chỉnh cho hiện tượng “bỏ cuộc”. • Điều chỉnh cho trường hợp không cân đối giữa hai nhóm. • Ước tính cỡ mẫu là bước đầu trong nghiên cứu.