PHÂN TÍCH DỮ LIỆU ĐỊNH TÍNH: BÀI TẬP THỰC HÀNH
Trần Khánh Vy
2024-06-02
Giới thiệu về bệnh AIDS
Bệnh AIDS (HIV/AIDS) là một trong những vấn đề sức khỏe quan trọng ở Úc. AIDS là tình trạng suy giảm hệ miễn dịch do virus HIV gây ra, làm suy yếu cơ thể và tăng nguy cơ mắc các bệnh nặng. Ở Úc, chính sách và chương trình phòng chống HIV/AIDS được triển khai mạnh mẽ, bao gồm giáo dục, xét nghiệm miễn phí và truy cứu người nhiễm HIV. Tuy nhiên, vẫn còn những thách thức về việc giảm thiểu nguy cơ lây nhiễm và cung cấp dịch vụ hỗ trợ cho những người sống với HIV/AIDS.
Dữ liệu sống sót của bệnh AIDS ở Úc
Dữ liệu về bệnh nhân được chẩn đoán mắc bệnh AIDS ở Úc trước ngày 1 tháng 7 năm 1991. Dữ liệu gồm 2843 quan sát và 7 biến. Các biến bao gồm:
State: Nhóm nước xuất xứ
Sex: Giới tính
Diag: Ngày chuẩn đoán
Death: Ngày chết hoặc nagyf kết thúc quan sát
Status: Tình trạng (còn sống hay đã chết)
T.categ: Thể loại truyền tải bái cáo
Age: Tuổi lúc chuẩn đoán
## # A tibble: 2,843 × 8
## rownames state sex diag death status T.categ age
## <dbl> <chr> <chr> <dbl> <dbl> <chr> <chr> <dbl>
## 1 1 NSW M 10905 11081 D hs 35
## 2 2 NSW M 11029 11096 D hs 53
## 3 3 NSW M 9551 9983 D hs 42
## 4 4 NSW M 9577 9654 D haem 44
## 5 5 NSW M 10015 10290 D hs 39
## 6 6 NSW M 9971 10344 D hs 36
## 7 7 NSW M 10746 11135 D other 36
## 8 8 NSW M 10042 11069 D hs 31
## 9 9 NSW M 10464 10956 D hs 26
## 10 10 NSW M 10439 10873 D hsid 27
## # ℹ 2,833 more rowsBảng tần số
table <- table(aids$sex, aids$status)
addmargins(table) ##
## A D Sum
## F 36 53 89
## M 1046 1708 2754
## Sum 1082 1761 2843- Số bệnh nhân mắc bệnh AIDS ở nữ giới còn sống sót là 36 người
- Số bệnh nhân mắc bệnh AIDS ở nữ giới đã chết là 53 người
- Số bệnh nhân mắc bệnh AIDS ở nam giới còn sống sót là 1046 người
- Số bệnh nhân mắc bệnh AIDS ở nam giới đã chết là 1708 người
Bảng tần suất
table2 <- ftable(aids$sex, aids$status)
table2 <- prop.table(table2)
round(addmargins(table2),3)## Sum
## 0.013 0.019 0.031
## 0.368 0.601 0.969
## Sum 0.381 0.619 1.000Đồ thị
library("ggplot2")
ggplot( aids, aes(sex, fill = status)) + geom_bar(position = 'dodge')
Ước lượng tỷ lệ
addmargins(prop.table(table(aids$sex, aids$status)))##
## A D Sum
## F 0.01266268 0.01864228 0.03130496
## M 0.36792121 0.60077383 0.96869504
## Sum 0.38058389 0.61941611 1.00000000NHẬN XÉT: - Tổng tỷ lệ bệnh nhân đã chết cao hơn so với tổng tỷ lệ sống sót. - Tỷ lệ nam giới mắc aids cao hơn rõ rệt so với giới tính nữ.
Rủi ro tương đối (Relative risk)
RelRisk(table(aids$sex,aids$status))## [1] 1.064988m <- matrix(c(36,53,1046,1708),nrow=2)
RelRisk(m, conf.level=.95)## rel. risk lwr.ci upr.ci
## 1.1054999 0.7306622 1.6706692Vậy rủi ro tương đối về khả năng sống sót của giới tính nữ so với giới tính nam là 1 (rủi ro tương đối trong khoảng 0,) 1,1 (rủi ro tương đối trong khoảng 0,73 đến 1,67 với alpha=5%).
Tỷ lệ chênh (odd Ratio)
OddsRatio(table(aids$sex, aids$status), conf.level=.95)## odds ratio lwr.ci upr.ci
## 1.1091309 0.7213243 1.7054345Vì OR là 1.109 nên # Hồi quy
chisq.test(table(aids$sex,aids$status))##
## Pearson's Chi-squared test with Yates' continuity correction
##
## data: table(aids$sex, aids$status)
## X-squared = 0.13041, df = 1, p-value = 0.718