Nếu dùng R để phân tích, không dùng các packages mà không thể tải qua lệnh install.packages("tên package"). Đúng sai người khác không đánh giá được!
Ví dụ, 2 biến số là benh và khoilop chứa trong dữ kiện tên là dukien, các lệnh sử dụng ở đây là tích hợp sẵn trong phần mềm R
Kiểm định khuynh hướng của tỷ lệ
Có nhiều dạng kiểm định khuynh hướng khác nhau, kiểm định phi tham số Cuzick, kiểm định Chi bình phương khuynh hướng. do lệnh nptrendr không thấy tài liệu giải thích là kiểm định gì nên không thể trả lời được.
Báo cáo kiểm định từng tầng
Thực chất là bạn đang lập phương trình hồi quy giữa benh và khoilop, thay vì chọn hàm liên kết là logit bạn chọn làm liên kết là log để báo cáo RR, thay vì OR (vì biến số bệnh là nhị giá nên thực chất đây là mô hình hồi quy logistic) Xem Chương 4
#Chuyển biến số thứ tự thành factor
dukien$khoilop = as.factor(dukien$khoilop)
#Lập mô hình
mohinh = glm(benh ~ khoilop, data=dukien, family =binomial(link="log"))
#báo cáo p-value từng nhóm
summary(mohinh)
#báo cáo RR và KTC 95% (lấy mũ exp(hệ số) = RR nếu link="log", là OR nếu link="logit")
exp(cbind(RR = coef(mohinh), confint(mohinh)))