Khóa học sử dụng R Ngày 1

Việc 3: Đọc dữ liệu

Yêu cầu đọc dữ liệu birthwt.csv vào R

bw = read.csv("C:\\Users\\anhth\\OneDrive\\Attachments\\Desktop\\R\\DỮ LIỆU THỰC HÀNH (TS Thạch gửi)\\birthwt.csv")

Việc 4: Kiểm tra dữ liệu

4.1. Tìm số biến số và quan sát

dim (bw)
## [1] 189  11

4.2. Liệt kê 6 dòng đầu tiên

head (bw)
##   id low age lwt race smoke ptl ht ui ftv  bwt
## 1 85   0  19 182    2     0   0  0  1   0 2523
## 2 86   0  33 155    3     0   0  0  0   3 2551
## 3 87   0  20 105    1     1   0  0  0   1 2557
## 4 88   0  21 108    1     1   0  0  1   2 2594
## 5 89   0  18 107    1     1   0  0  1   0 2600
## 6 91   0  21 124    3     0   0  0  0   0 2622

4.3. Liệt kê 10 dòng cuối cùng

tail(bw,10)
##     id low age lwt race smoke ptl ht ui ftv  bwt
## 180 71   1  17 120    2     0   0  0  0   2 2438
## 181 75   1  26 154    3     0   1  1  0   1 2442
## 182 76   1  20 105    3     0   0  0  0   3 2450
## 183 77   1  26 190    1     1   0  0  0   0 2466
## 184 78   1  14 101    3     1   1  0  0   0 2466
## 185 79   1  28  95    1     1   0  0  0   2 2466
## 186 81   1  14 100    3     0   0  0  0   2 2495
## 187 82   1  23  94    3     1   0  0  0   0 2495
## 188 83   1  17 142    2     0   0  1  0   0 2495
## 189 84   1  21 130    1     1   0  1  0   3 2495

Việc 5. Biên tập dữ liệu

5.1. Tạo biến mwt (kg)

bw$mwt = bw$lwt*0.45

5.2. Tạo biến ethnicity

bw$ethnicity[bw$race == 1] = "White"
bw$ethnicity[bw$race == 2] = "Black"
bw$ethnicity[bw$race == 3] = "Other"

5.3. Tạo 1 tập dữ liệu bw1 chỉ gồm 3 biến số id, low và bwt. Số biến số và quan sát?

bw1 = bw[, c("id","low","bwt")]
dim(bw1)
## [1] 189   3

5.4. Tạo 1 tập dữ liệu bw3 chỉ gồm những thai phụ có cân nặng thấp (low = 1). Số biến số và quan sát?

bw3 = subset(bw, low == 1)
dim(bw3)
## [1] 59 13

5.5. Tạo 1 tập dữ liệu bw4 chỉ gồm những thai phụ có cân nặng thấp (low = 1) và có hút thuốc trong lúc mang thai (smoke = 1). Số biến số và quan sát?

bw4 = subset(bw, low ==1 & smoke ==1)
dim(bw4)
## [1] 30 13

5.6. Tạo biến ethnicity (bw5)

bw$bw5 = factor(bw$race,
                levels = c(1,2,3),
                labels = c("White", "Black", "Other"))
table(bw$bw5, bw$race)
##        
##          1  2  3
##   White 96  0  0
##   Black  0 26  0
##   Other  0  0 67