Các mô hình tiên lượng tử vong Hà Tấn Đức Bệnh viện đa khoa trung ương Cần Thơ Đại học Tôn Đức Thắng Thực hành R Ngày 4 1 Câu hỏi 13 Khoảng tin cậy đối với những biến cố chưa xuất hiện một thử nghiệm.
Thực hành R Ngày Hà Tấn Đức Bệnh viện đa khoa trung ương Cần Thơ Đại học Tôn Đức Thắng Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 13 • Khoảng tin cậy biến cố chưa xuất hiện: thử nghiệm lâm sàng theo dõi nhóm 90 bệnh nhân vịng năm, khơng thấy trường hợp bị phản tác dụng thuốc Tính khoảng tin cậy cho tỷ lệ phát sinh trường hợp bị phản tác dụng thuốc Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 13 • Khoảng tin cậy biến cố chưa xuất hiện: thử nghiệm lâm sàng theo dõi nhóm 90 bệnh nhân vịng năm, khơng thấy trường hợp bị phản tác dụng thuốc Tính khoảng tin cậy cho tỷ lệ phát sinh trường hợp bị phản tác dụng thuốc • library(binom) • binom.confint(x=0, n=90, conf.level=0.95) • library(epitools) • pois.exact(x=0, pt=90*3) Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 14 • So sánh tỷ lệ nghiên cứu cắt ngang: nghiên cứu 700 phụ nữ Việt Nam, nhà nghiên cứu phát 148 người bị loãng xương Một nghiên cứu khác Úc phụ nữ độ tuổi (n = 1287) có 345 người bị loãng xương Hai tỷ lệ thật khác nhau? • Hàm CrossTable Package: gmodels Các argument: − x: vector matrix Nếu có y x phải vector − y: vector − chisq: chi-square test − fisher: Fisher exact test Câu hỏi 14 • Hàm matrix Package: base Các argument hàm: − data: vector tùy ý − nrow: số hàng (dòng) mong muốn − ncol: số cột mong muốn − dimnames: tên hàng, tên cột Câu hỏi 14 Loãng xương Việt Nam Australia 148 345 Khơng lỗng xương 700 – 148 1287 - 345 library(gmodels) CrossTable(matrix(c(148,345,700-148,1287-345),nrow=2), fisher=T) Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 15 Nhóm A Nhóm B Số tử vong 26 Số sống sót 15 28 So sánh hai tỷ lệ phát sinh (incidence rates) nghiên cứu đồn hệ Tính RR (95%CI RR) Câu hỏi 15 • Hàm epi2x2 Package: epibasix Các argument: − X: matrix Cột bệnh, cột không bệnh Hàng 1: phơi nhiễm, hàng 2: không phơi nhiễm − alpha: Type I Error Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 15 Nhóm A Nhóm B Số tử vong 26 Số sống sót 15 28 So sánh hai tỷ lệ phát sinh (incidence rates) nghiên cứu đồn hệ Tính RR (95%CI RR) library(epibasix) oc=epi2x2(matrix(c(5,26,15,28),nrow=2)) summary(oc) oc$RR oc$RR.CIL oc$RR.CIU Câu hỏi 16 • File PISA DATA.csv • Đọc liệu vào R gọi đối tượng pisa • Trình bày kết theo bảng số liệu soạn 10 Câu hỏi 16 • Hàm tabular Package: tables Các argument: − table: formula dùng để tính (y ~ x) y: hàng bảng, x: cột bảng • +: concatenation • *: nesting − n: thường thiết lập 11 Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 16 • library(tables) • tabular(pisa$REGION~1*(n=1+Percent('col'))) #Bảng • tabular(pisa$REGION~pisa$AREA*(n=1+Percent('col'))) #Bảng • tabular(pisa$REGION*pisa$TYPE~pisa$AREA*(n=1+Pe rcent('col'))) #Bảng • tabular(pisa$TYPE*pisa$REGION~pisa$AREA*(n=1+Pe rcent('col'))) #Bảng • tabular(pisa$REGION~pisa$PV1MATH*(n=1 + mean +sd)) #Bảng 12 Phân tích mơ tả biến phân nhóm (tables, binomial test) Câu hỏi 16 • tabular(pisa$TYPE*pisa$REGION~pisa$PV1MATH*(n=1 + mean +sd)) #Bảng • tabular(pisa$REGION*pisa$TYPE~pisa$PV1MATH*(n=1 + mean +sd)) #Bảng • tabular(pisa$TYPE*pisa$REGION*pisa$AREA~pisa$PV 1MATH*(n=1 + mean +sd)) #Bảng 13 ... cột Câu hỏi 14 Lỗng xương Việt Nam Australia 148 345 Khơng lỗng xương 700 – 148 1287 - 345 library(gmodels) CrossTable(matrix(c( 148 , 345 ,700- 148 ,1287- 345 ),nrow=2), fisher=T) Phân tích mơ tả biến... test) Câu hỏi 14 • So sánh tỷ lệ nghiên cứu cắt ngang: nghiên cứu 700 phụ nữ Việt Nam, nhà nghiên cứu phát 148 người bị loãng xương Một nghiên cứu khác Úc phụ nữ độ tuổi (n = 1287) có 345 người bị... test Câu hỏi 14 • Hàm matrix Package: base Các argument hàm: − data: vector tùy ý − nrow: số hàng (dòng) mong muốn − ncol: số cột mong muốn − dimnames: tên hàng, tên cột Câu hỏi 14 Lỗng xương