Tài liệu hướng dẫn thống kê y học sử dụng SPSS SPSS Phân tích hiệp phương sai (ANCOVA) Lê Đông Nhật Nam Lời nói đầu Tôi nhớ kỷ niệm lần làm nghiên cứu mình, báo cáo thực tập cuối năm thứ Master Pháp Thầy giao cho đề tài khảo sát biomarker khí thở bệnh nhân tiểu đường Mọi việc suôn sẻ, trừ … kết quả; giả thuyết ban đầu muốn tìm thấy thay đổi so với nhóm chứng, kết lại cho thấy khác biệt ý nghĩa Quả thật báo cáo kết âm tính điều không dễ chịu, điều làm tiếc nuối vài năm sau khám phá loại biomarker ngày trước khảo sát đại lượng bất định, có nhiều yếu tố gây nhiễu tác động lên điều kiện sinh lý bình thường Tôi tha thứ cho ngây thơ ngày không xét toàn diện tất yếu tố gây nhiễu Đó lý 10 năm sau, muốn gửi tặng tất bạn trẻ tập làm nghiên cứu đầu đời, tài liệu phương pháp thống kê cho phép bạn chủ động sửa chữa sai lầm tiềm ẩn nghiên cứu để tiếp cận với thật; phương pháp ANCOVA hay phân tích hiệp phương sai Thông điệp quan trọng tài liệu hồ nghi tò mò, nhiều lúc giá trị p, chí khoảng tin cậy với Bootstrap mà bạn tìm với ANOVA, test t hay mô hình hồi quy dối lừa chúng ta, chưa phải thật Có lẽ bạn bắt đầu … tò mò muốn biết ? Mời bạn đọc thí dụ tài liệu để biết nguyên nhân Chúc bạn thành công tìm thật cho câu hỏi nghiên cứu Quy ước trình bày Trong tài liệu làm quen với nhân vật Bác sĩ Nguyễn Văn Thái Bác sĩ Thái cao thủ thống kê y học giang hồ Anh ta xuất trước vấn đề khó khăn, nhằm đưa câu trả lời xác ôn lại cho bạn kiến thức chuyên sâu lý thuyết thống kê Tuy nhiên BS Thái lại bận rộn nên không đủ thời gian vào cách thực chi tiết Anh ta hay trình bày lý thuyết túy sử dụng nhiều công thức toán học nên không gần gũi với sinh viên gây không trở ngại cho bạn vốn dị ứng với thống kê … Bác sĩ Lê Ngọc Khả Nhi Khả Nhi nữ bác sĩ trẻ dễ thương sử dụng thành thạo SPSS Như tên gọi mình, BS Nhi có tính cách hồn nhiên ngây thơ trẻ con, cô có khuynh hướng đơn giản hóa tối đa vấn đề Khả Nhi hướng dẫn bạn sử dụng SPSS qua bước cụ thể, chia sẻ mẹo vặt, thủ thuật để giúp bạn đến kết nhanh dễ dàng Sinh viên Trần Quốc Bảo Bảo sinh viên y khoa năm thứ sáu bắt đầu làm quen với nghiên cứu khoa học Đây cậu sinh viên tò mò đặt nhiều câu hỏi liên quan đến thống kê Mặc dù đế tài Bảo thực đơn giản, đồng hành với cậu ta, bạn có hội tích lũy cho nhiều kinh nghiệm công việc phân tích số liệu thiết kế nghiên cứu 1.1 Tình thí dụ Dòng máu lưu thông Xơ vữa ĐM Cảnh Bề dày nội trung mô ĐM Cảnh (CIMT) đo Siêu âm Doppler Siêu âm Doppler động mạch cảnh phương pháp đánh giá không xâm lấn đoạn sọ động mạch cảnh Chiều dày nội trung mô động mạch cảnh (CIMT) số chứng vữa xơ đo siêu âm Doppler Một chị bác sĩ nội trú khoa hình ảnh học thực đề tài khảo sát tượng dày hóa thành động mạch cảnh tác động xạ trị ung thư vùng hầu họng; với CIMT biến số Thiết kế nghiên cứu bao gồm: + Nhóm chứng âm: 15 người bình thường tiền sử cao huyết áp xơ vữa động mạch + Nhóm chứng dương: 15 bệnh nhân tiểu đường type 2, yếu tố nguy xơ vữa động mạch + Nhóm nghiên cứu: 15 bệnh nhân ung thư carcinoma hầu họng xạ trị trước từ 3-4 năm Khác biệt ý nghĩa với nhóm chứng dương ? Tuy nhiên ý đồ bị thất bại, sau kết ANOVA cho thấy: Dù CIMT cao nhóm xạ trị ung thư so với nhóm chứng âm, nhiên khác biệt ý nghĩa so với nhóm chứng dương, chí CIMT bệnh nhân tiểu đường type cao ? Chúng ta phải làm để giúp chị ta ? Bảo thân mến, cô bác sĩ phạm sai lầm từ đầu dùng ANOVA đơn biến để khảo sát CIMT… Theo biết, trình lão hóa tự nhiên gây tăng CIMT, kết quan sát nhóm bị gây nhiễu Tuổi bệnh nhân Chúng ta phải hiệu chỉnh kết CIMT theo tuổi trước so sánh Tôi đề nghị dùng ANCOVA thay ANOVA… 1.2 ANCOVA phương pháp kẻ tò mò hoài nghi Nghiên cứu sinh « ngây thơ » Điều kiện liệu kết ANOVA post-hoc test cho kết có ý nghĩa thống kê… ANOVA post-hoc test cho kết âm tính: Không có khác biệt ý nghĩa Vui mừng, thỏa mãn với giá trị p[...]... ,280 Quy trình phân tích phương sai bắt đầu bằng việc kiểm tra giả định có sự đồng nhất phương sai của biến số cần khảo sát bằng test Levene Ở đây ta làm ANOVA song song cho cả Y (Bề dày nội trung mạc) và C (Tuổi) nên bảng Levene test có 2 hàng tương ứng cho C và Y Thông tin cần quan tâm là cột Sig hay giá trị p Kết quả trông đợi là p>0,05 cho thấy có sự đồng nhất về phương sai ở các phân nhóm Nếu p0,05 cho thấy có sự đồng nhất về phương sai, giả định được thỏa Tuy nhiên điều quan trọng ở đây mà ta muốn kiểm tra không phải là đồng nhất phương sai, mà là đồng nhất về RSD (giá trị thặng dư chuẩn hóa), nên test Levene là điều kiện cần nhưng chưa... để kiểm tra giả định về đồng nhất phương sai của C và Y giữa các phân nhóm + Thực hiện F test cho Y,C theo X (ANOVA) Chú ý: + So sánh bắt cặp tuần tự giữa các phân nhóm bằng Posthoc test Bonferroni và Games-Howell Để khái quát hóa bộ Syntax, tác giả cho quy trình ANOVA làm cùng lúc 2 loại post-hoc test là : Ý nghĩa của quy trình là để : Bonferroni: nếu giả định phương sai đồng nhất thỏa mãn (Levene test... cách chỉnh sửa số liệu, chuyển dạng biến số… rồi kiểm tra lại trước khi đi tiếp 2 2.3 Phân tích nội dung Syntax cho ANCOVA Khối lệnh thứ 6: * Bước 6: Kiểm tra giả định Homoscedasticity và giả định đồng nhất phương sai của giá trị thặng dư Kiểm tra giả định cho mô hình hồi quy tuyến tính ANCOVA UNIANOVA Y BY X WITH C /METHOD=SSTYPE(3) /INTERCEPT=INCLUDE /SAVE=PRED (PRE) ZRESID (RSD) /PLOT=PROFILE(X) /PLOT=RESIDUALS... label("Giá trị thặng dư chuẩn hóa (RSD)")) ELEMENT: point(position(PRE*RSD)) END GPL Thực hiện ANCOVA lần thứ 2 với mô hình chứa 2 biến số C và X, không kèm theo bootstrap hay phân tích sâu Mục đích duy nhất là để sao lưu kết quả dự báo và sai số thặng dư, từ đó kiểm tra các giả định về homoscedasticity và đồng nhất phương sai của giá trị thặng dư chuẩn hóa (RSD) Các giả định này rất quan trọng để kiểm tra mức... mình hiệp biến số C (ở đây là tuổi) X: Ảnh hưởng riêng phần của 1 mình biến số X hay yếu tố chính (ở đây là phân nhóm bệnh nhân) Error: Phần biến thiên không thể giải thích được trong mô hình nguyên thủy (không bao gồm biến số C) Total: Tổng bình phương (SST) của mô hình nguyên thủy (không chứa C) Corrected Total: Tổng bình phương (SST) của mô hình sau khi hiệu chỉnh (bao gồm biến số C) a Tổng bình phương. .. của ANCOVA không tốt 3 Diễn giải kết quả 3.8 Đọc kết quả chính của ANCOVA 3 Levene's Test of Equality of Error Variancesa Dependent Variable: Bề dày nội trung mạc F df1 df2 Sig 1,987 2 42 ,150 Tests the null hypothesis that the error variance of the dependent variable is equal across groups a Design: Intercept + C + X Cũng như ANOVA, quy trình ANCOVA bắt đầu bằng test Levene để kiểm tra giả định phương. .. chọn 1) Chuyển dạng biến số (ví dụ thang đo Logarit), với hy vọng sẽ có phân phối chuẩn 2) Sử dụng phương pháp khác ví dụ mô hình hồi quy đa cấp có Bootstrap để thay thế cho GLM -ANCOVA 3) Bỏ qua vi phạm và vẫn làm phân tích ANCOVA (Thực ra điều kiện phân phối không phải là tuyệt đối bắt buộc) 4) Làm ANCOVA song song cho 2 mẫu số liệu: nguyên thủy và đã chuyển dạng (logarit hóa), so sánh kết quả của ... thích ANCOVA ANCOVAcũng dùng để kiểm tra giả định quy luật mà ta tìm với ANOVA KHÔNG cho MỌI trường hợp Lúc ANCOVA cho phép đặt giả thuyết mới, quy luật mới, hướng nghiên cứu mới… Tóm lại, ANCOVA. .. hài lòng với giá trị p, Effectsize… mô hình gốc, muốn dùng ANCOVA với hy vọng làm cho kết đẹp chút 1 1.6 Những giả định ANCOVA Bản chất ANCOVA mô hình hồi quy tuyến tính , nên biến số cần phân... đơn biến cho X Do X C gây ANCOVA lý tưởng (X C độc lập) Không giải thích Do X C gây Không giải thích ANCOVA không lý tưởng (X C có tương tác) 2 2.1 Chuẩn bị quy trình ANCOVA Kiểm tra giả định: