Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 66 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
66
Dung lượng
0,98 MB
File đính kèm
128. Xu ly va phan tich du kien.rar
(1 MB)
Nội dung
Xử lý phân tích kiện Thái Thanh Trúc Nội dung học • Kiểm định so sánh với kết định tính – Kiểm định so sánh tỉ lệ – Kiểm định Chi bình phương – Kiểm định xác Fisher • Phân tích NC dịch tễ học Quá trình kiểm định ý nghĩa Xây dựng giả thuyết Ho Chọn lựa kiểm định thích hợp Tính giá trị thống kê số liệu thu thập Tính xác suất thống kê - kí hiệu P() gọi giá trị p Nếu P() đủ nhỏ kết luận P(Ho) nhỏ bác bỏ giả thuyết Ho Chọn lựa kiểm định phù hợp Biến phụ thuộc Biến độc lập Nhị giá Danh định Thứ tự -Định lượng Đa biến T-test ANOVA Hồi quy tuyến tính Thứ tự MannWhitney Kruskal-Wallis TQ Spearman Nhị giá Chi bình phương Chi bình phương Hồi quy logistic Hồi quy Poisson Sống Wilcoxon tổng quát Logrank Wilcoxon tổng quát Logrank Hồi quy Cox Định lượng có phân phối bình thường Kiểm định so sánh tỉ lệ • Điều trị bệnh X có phương pháp A B Trong 93 người điều trị theo phương pháp A có 52.7% đáp ứng điều trị Trong 91 người điều trị theo phương pháp B có 19.8% đáp ứng điều trị • Hỏi: phương pháp có tỉ lệ đáp ứng điều trị cao hơn? • Bước 1: Giả thuyết Ho – Tỉ lệ đáp ứng điều trị phương pháp A B – p A = pB Kiểm định so sánh tỉ lệ • Bước 2: lựa chọn kiểm định thích hợp – So sánh tỉ lệ • Bước 3: Tính số thống kê tương ứng p = (p1 x n1 + p2 x n2) / (n1 + n2) p = (0.527 x 93 + 0.198 x 91) / (93 + 91) = 0.364 S.E p(1 - p)( 1 ) n1 n 1 S.E 0.364(1 - 0.364)( ) 0.071 93 91 z = (p1 - p2) / SE z = (0.527 – 0.198) / 0.071 = 4.63 Kiểm định so sánh tỉ lệ • Bước 4: tính xác suất thống kê z: tra bảng với z = 4.63 • Dựa vào tra bảng P(|Z| 1,960) =0,05 p < 0.05 • Bước 5: kết luận p < 0.05 bác bỏ giả thuyết Ho tỉ lệ đáp ứng điều trị nhóm có khác biệt Kiểm định so sánh tỉ lệ prtesti n1 p1 n2 p2 • Ví dụ: prtesti 93 0.527 91 0.198 Ví dụ • Hai nghiên cứu xác định tỉ lệ tăng huyết áp tỉnh Việt Nam Nghiên cứu thứ tiến hành 2000 người cho tỉ lệ 9.6% Nghiên cứu thứ tiến hành 3500 người cho tỉ lệ 8% • Hãy cho biết kết nghiên cứu tương đồng hay có khác biệt? Kiểm định Chi bình phương • Vọng trị – Vọng trị # giá trị kỳ vọng # giá trị mong đợi # kỳ vọng toán học – Là trung bình biến số ngẫu nhiên phép thử lập lại nhiều lần • Để ước lượng thực nghiệm giá trị kỳ vọng biến ngẫu nhiên, người ta liên tục thực quan sát biến tính trung bình cộng kết 10 Các bước phân tích Lượng hóa mối quan hệ (PR) Giả sử mối quan hệ nhóm tuổi tình trạng nhẹ cân có tính chất khuynh hướng glm nhecan nhomtuoi, family(poisson) link(log) robust eform Nhóm tuổi OR (KTC 95%) PR = 0.94 (KTC 95%: 0.72 – 1.22) Ở người mẹ có tuổi cao 0.05) 35 – 39 tuổi 0.88 (0.52 – 1.49) Nhóm tuổi cao tỉ lệ sinh nhẹ cân thấp điều không 40+ tuổi 0.83 (0.37 – 1.82) có ý nghĩa thống kê 52 Các bước phân tích Lượng hóa mối quan hệ (PR) glm nhecan i.nghenghiep, family(poisson) link(log) robust eform mặc định lấy nhóm làm chuẩn glm nhecan ib3.nghenghiep, family(poisson) link(log) robust eform lấy nhóm có giá trị làm chuẩn 53 Các bước phân tích Lượng hóa mối quan hệ (OR) cc biếnphụthuộc biếnđộclập nptrend biếnphụthuộc, by(biếnđộclập) logistic biếnphụthuộc i.biếnđộclập logistic biếnphụthuộc biếnđộclập 54 Các bước phân tích Nhiễu/tương tác • Yếu tố gây nhiễu – Liên quan đến phơi nhiễm (biến độc lập) – Liên quan đến kết (biến phụ thuộc) – Không phải yếu tố ngun nhân kết Số lần khám tiền sản (biến số độc lập) Sanh nhẹ 2500 g (biến số phụ thuộc) Thu nhập gia đình (biến số gây nhiễu) 55 Các bước phân tích Nhiễu/tương tác • Yếu tố tương tác – Mối quan hệ (mức độ hướng) yếu tố phơi nhiễm kết khác theo mức độ yếu tố 56 Các bước phân tích Nhiễu/tương tác PHÂN TÍCH BẢNG 2x2 PHÂN TẦNG KIỂM ĐỊNH TÍNH ĐỒNG NHẤT CỦA OR/PR THEO TẦNG P VALUE < 0.05 CÓ CÓ TƯƠNG TÁC – BÁO CÁO OR/RR TỪNG TẦNG CÓ CÓ GÂY NHIỄU – BÁO CÁO OR/PR KẾT HỢP MH KHÔNG SO SÁNH OR/RR VÀ OR/PR KẾT HỢP MH KHÁC BIỆT >10% KHÔNG KHÔNG NHIỄU VÀ KHÔNG TƯƠNG TÁC – BÁO CÁO OR/PR THƠ 57 Các bước phân tích Nhiễu/tương tác Mối quan tâm chính: Sinh non Nhẹ cân (PRthơ) Yếu tố cần xét: Giới Phân tích bảng 2x2 phân tầng: Trai: PRsinhnonnhecan (PRtrai) Gái: PRsinhnonnhecan (PRgái) Kiểm tra tính đồng PRtrai ≠ PRgái: tương tác PRtrai ≈ PRgái tính PR có giới (PRkếthợp) Xét khác biệt PRthô ≠ PRkếthợp (>10%): gây nhiễu PRthô ≈ PRkếthợp (10%): gây nhiễu báo cáo PR tính Giới (PRkếthợp) PRthơ ≈ PRkếthợp ( 0.05 PR nhóm trai nhóm gái khơng tương tác Sự khác biệt PR thô PR kết hợp = (11.961 – 11.945)/11.945 = 0.0013 ≈ 0,13% yếu tố giới xuất không làm thay đổi đáng kể mối quan hệ sinh non nhẹ cân ( 0.05 PR nhóm huyết áp bình thường nhóm huyết áp tăng khơng tương tác Sự khác biệt PR thô PR kết hợp = (11.961 – 9.918)/9.918 = 0.2060 ≈ 20,60% yếu tố tăng huyết áp xuất làm thay đổi đáng kể mối quan hệ sinh non nhẹ cân (>10%) yếu tố gây nhiễu 61 Các bước phân tích Mơ hình hóa (kiểm sốt nhiều yếu tố tác động) Kiểm soát yếu tố tăng huyết áp nghề nghiệp glm nhecan sinhnon tang_ha i.nghenghiep, family(poisson) link(log) robust eform Sau kiểm soát theo yếu tố tăng huyết áp nghề nghiệp người mẹ sinh non có tỉ lệ sinh nhẹ cân 10.88 lần người không sinh non Sự khác biệt có ý nghĩa thống kê (p