BÀI TẬP
2023-06-04
Đề bài: Mô phỏng ít nhất 5 biến ngẫu nhiên (có phân phối xác suất khác nhau), mô phỏng, vẽ đồ thị, tính toán các đặc trưng đo lường và giải thích ý nghĩa.
Bài làm: Thống kê và theo dõi các bệnh nhân mắc bệnh cao huyết áp tại xã A dựa trên độ tuổi và các chỉ số về cholesterol của 15 bệnh nhân, với các biến đầu vào là: Giới tính, độ tuổi, chỉ số BMI, chỉ số HLD-Cholesterol và chỉ số LDL-Cholesterol.
# Trích dữ liệu từ file và gắn vào biến *ins*
setwd("D:/HN")
ins <- read.csv("bn.csv", header=TRUE)
str(ins)## 'data.frame': 15 obs. of 5 variables:
## $ SEX: chr "Nam" "Nu" "Nu" "Nam" ...
## $ AGE: int 57 64 60 65 47 65 76 61 59 57 ...
## $ BMI: int 17 18 18 18 18 18 19 19 19 19 ...
## $ HDL: int 5000 4380 3360 5920 6250 4150 737 7170 6942 5000 ...
## $ LDL: num 2 3 3 4 2.1 3 3 3 3 2 ...ab <- table(ins$SEX)
ab##
## Nam Nu
## 8 7# Trong 15 bệnh nhân mắc bệnh cao huyết áp, bao gồm 8 bệnh nhân nam và 7 bệnh nhân nữ.
# Thống kê mô tả các biến
summary(ins)## SEX AGE BMI HDL
## Length:15 Min. :45.00 Min. :17.00 Min. : 737
## Class :character 1st Qu.:57.00 1st Qu.:18.00 1st Qu.:4265
## Mode :character Median :61.00 Median :19.00 Median :5920
## Mean :59.87 Mean :20.33 Mean :5495
## 3rd Qu.:64.00 3rd Qu.:24.00 3rd Qu.:7170
## Max. :76.00 Max. :25.00 Max. :7880
## LDL
## Min. :1.000
## 1st Qu.:2.000
## Median :3.000
## Mean :2.713
## 3rd Qu.:3.000
## Max. :4.000# Độ tuổi trung bình của các bệnh nhân là 61 tuổi, nhỏ nhất là bệnh nhân 45 tuổi và cao nhất là 76 tuổi.
# Chỉ số BMI trung bình là 19, nhìn chung hầu hết bệnh nhân đều có mức BMI bình thường.
# Chỉ số HDL trung bình là 5.92, thấp nhất là 0.737 và cao nhất là 7.88
# Chỉ số LDL trung bình là 3, thấp nhất là 1 và cao nhất là 4
barplot(ab, main="GIỚI TÍNH")
hist(ins$AGE, main="TUỔI", xlab="Số tuổi", ylab="Số bệnh nhân")
hist(ins$BMI, main="Chỉ số BMI", ylab="Số bệnh nhân")
hist(ins$HDL, main="Chỉ số HDL", ylab="Số bệnh nhân")
hist(ins$LDL, main="Chỉ số LDL", ylab="Số bệnh nhân")