Giải thích - Dataset (infert):Infertility after Spontaneous and Induced Abortion (Tập dữ liệu cho biết tình trạng của việc phá thai tự nhiên và dùng thuốc)
#Dữ liệu infert
data("infert")
#Gán dữ liệu infert vào i
i <- infert
#Cấu trúc của dataset(infert)
str(infert)## 'data.frame': 248 obs. of 8 variables:
## $ education : Factor w/ 3 levels "0-5yrs","6-11yrs",..: 1 1 1 1 2 2 2 2 2 2 ...
## $ age : num 26 42 39 34 35 36 23 32 21 28 ...
## $ parity : num 6 1 6 4 3 4 1 2 1 2 ...
## $ induced : num 1 1 2 2 1 2 0 0 0 0 ...
## $ case : num 1 1 1 1 1 1 1 1 1 1 ...
## $ spontaneous : num 2 0 0 0 1 1 0 0 1 0 ...
## $ stratum : int 1 2 3 4 5 6 7 8 9 10 ...
## $ pooled.stratum: num 3 1 4 2 32 36 6 22 5 19 ...#Đặt tên cho các biến trong dataset
names(i) <- c('E','A','Pa','I','C','Sp','St','Po')head(i)## E A Pa I C Sp St Po
## 1 0-5yrs 26 6 1 1 2 1 3
## 2 0-5yrs 42 1 1 1 0 2 1
## 3 0-5yrs 39 6 2 1 0 3 4
## 4 0-5yrs 34 4 2 1 0 4 2
## 5 6-11yrs 35 3 1 1 1 5 32
## 6 6-11yrs 36 4 2 1 1 6 36#Gán biến A trong data(i) vào AGE
AGE <- i$A
AGE## [1] 26 42 39 34 35 36 23 32 21 28 29 37 31 29 31 27 30 26 25 44 40 35 28 36 27
## [26] 40 38 34 28 30 32 34 42 32 39 35 36 34 30 28 39 35 41 37 30 37 28 27 26 38
## [51] 24 36 27 28 29 36 28 28 28 27 35 25 34 31 26 32 21 28 37 25 32 25 31 38 26
## [76] 31 31 25 31 34 35 29 23 26 42 39 34 35 36 23 32 21 28 29 37 31 29 31 27 30
## [101] 26 25 44 40 35 28 36 27 40 38 34 28 30 32 34 42 32 39 35 36 34 30 28 39 35
## [126] 41 37 30 37 28 27 26 38 24 36 27 28 29 36 28 28 28 27 35 25 34 31 26 32 21
## [151] 28 37 25 32 25 31 26 31 31 25 31 34 35 29 23 26 42 39 34 35 36 23 32 21 28
## [176] 29 37 31 29 31 27 30 26 25 44 40 35 28 36 27 40 38 34 28 30 32 34 42 32 39
## [201] 35 36 34 30 28 39 35 41 37 30 37 28 27 26 38 24 36 27 28 29 36 28 28 28 27
## [226] 35 25 34 31 26 32 21 28 37 25 32 25 31 38 26 31 31 25 31 34 35 29 23#LỌc dữ liệu có độ tuổi từ 20 đến 40
AGE20_40 <- AGE[AGE>20&AGE<40]
AGE20_40## [1] 26 39 34 35 36 23 32 21 28 29 37 31 29 31 27 30 26 25 35 28 36 27 38 34 28
## [26] 30 32 34 32 39 35 36 34 30 28 39 35 37 30 37 28 27 26 38 24 36 27 28 29 36
## [51] 28 28 28 27 35 25 34 31 26 32 21 28 37 25 32 25 31 38 26 31 31 25 31 34 35
## [76] 29 23 26 39 34 35 36 23 32 21 28 29 37 31 29 31 27 30 26 25 35 28 36 27 38
## [101] 34 28 30 32 34 32 39 35 36 34 30 28 39 35 37 30 37 28 27 26 38 24 36 27 28
## [126] 29 36 28 28 28 27 35 25 34 31 26 32 21 28 37 25 32 25 31 26 31 31 25 31 34
## [151] 35 29 23 26 39 34 35 36 23 32 21 28 29 37 31 29 31 27 30 26 25 35 28 36 27
## [176] 38 34 28 30 32 34 32 39 35 36 34 30 28 39 35 37 30 37 28 27 26 38 24 36 27
## [201] 28 29 36 28 28 28 27 35 25 34 31 26 32 21 28 37 25 32 25 31 38 26 31 31 25
## [226] 31 34 35 29 23#Gán biến Pa trong data(i) vào PAR
PAR <- i$Pa
PAR## [1] 6 1 6 4 3 4 1 2 1 2 2 4 1 3 2 2 5 1 3 1 1 2 2 1 2 2 2 3 4 4 1 2 1 2 1 2 1
## [38] 3 3 1 3 1 1 2 1 1 2 4 2 3 3 5 3 1 2 2 2 2 1 2 2 1 1 2 2 1 1 3 3 1 1 1 1 6
## [75] 2 1 2 1 1 1 2 1 1 6 1 6 4 3 4 1 2 1 2 2 4 1 3 2 2 5 1 3 1 1 2 2 1 2 2 2 3
## [112] 4 4 1 2 1 2 1 2 1 3 3 1 3 1 1 2 1 1 2 4 2 3 3 5 3 1 2 2 2 2 1 2 2 1 1 2 2
## [149] 1 1 3 3 1 1 1 1 2 1 2 1 1 1 2 1 1 6 1 6 4 3 4 1 2 1 2 2 4 1 3 2 2 5 1 3 1
## [186] 1 2 2 1 2 2 2 3 4 4 1 2 1 2 1 2 1 3 3 1 3 1 1 2 1 1 2 4 2 3 3 5 3 1 2 2 2
## [223] 2 1 2 2 1 1 2 2 1 1 3 3 1 1 1 1 6 2 1 2 1 1 1 2 1 1#Lọc dữ liệu có Age lơn hơn 30 và Parity lớn hơn 3
AGE30PAR3 <- i[i$A>30|i$Pa>3,]
AGE30PAR3## E A Pa I C Sp St Po
## 1 0-5yrs 26 6 1 1 2 1 3
## 2 0-5yrs 42 1 1 1 0 2 1
## 3 0-5yrs 39 6 2 1 0 3 4
## 4 0-5yrs 34 4 2 1 0 4 2
## 5 6-11yrs 35 3 1 1 1 5 32
## 6 6-11yrs 36 4 2 1 1 6 36
## 8 6-11yrs 32 2 0 1 0 8 22
## 12 6-11yrs 37 4 2 1 1 12 37
## 13 6-11yrs 31 1 1 1 0 13 9
## 15 6-11yrs 31 2 1 1 1 15 21
## 17 6-11yrs 30 5 2 1 1 17 38
## 20 6-11yrs 44 1 0 1 1 20 17
## 21 6-11yrs 40 1 0 1 1 21 14
## 22 6-11yrs 35 2 2 1 0 22 24
## 24 6-11yrs 36 1 0 1 1 24 12
## 26 6-11yrs 40 2 0 1 2 26 27
## 27 6-11yrs 38 2 0 1 2 27 26
## 28 6-11yrs 34 3 0 1 2 28 31
## 29 6-11yrs 28 4 1 1 2 29 34
## 30 6-11yrs 30 4 2 1 0 30 35
## 31 6-11yrs 32 1 0 1 1 31 10
## 32 6-11yrs 34 2 1 1 0 32 23
## 33 6-11yrs 42 1 1 1 0 33 16
## 34 6-11yrs 32 2 0 1 2 34 22
## 35 6-11yrs 39 1 1 1 0 35 13
## 36 6-11yrs 35 2 0 1 2 36 24
## 37 6-11yrs 36 1 0 1 1 37 12
## 38 6-11yrs 34 3 1 1 2 38 31
## 41 6-11yrs 39 3 0 1 2 41 33
## 42 6-11yrs 35 1 0 1 0 42 11
## 43 6-11yrs 41 1 0 1 0 43 15
## 44 6-11yrs 37 2 1 1 1 44 25
## 46 12+ yrs 37 1 1 1 0 46 48
## 48 12+ yrs 27 4 2 1 0 48 61
## 50 12+ yrs 38 3 0 1 2 50 60
## 52 12+ yrs 36 5 1 1 2 52 62
## 56 12+ yrs 36 2 0 1 2 56 55
## 61 12+ yrs 35 2 0 1 2 61 54
## 63 12+ yrs 34 1 0 1 1 63 47
## 64 12+ yrs 31 2 0 1 2 64 53
## 66 12+ yrs 32 1 0 1 1 66 46
## 69 12+ yrs 37 3 0 1 2 69 59
## 71 12+ yrs 32 1 1 1 0 71 46
## 73 12+ yrs 31 1 0 1 1 73 45
## 74 12+ yrs 38 6 0 1 2 74 63
## 76 12+ yrs 31 1 0 1 1 76 45
## 77 12+ yrs 31 2 0 1 1 77 53
## 79 12+ yrs 31 1 0 1 1 79 45
## 80 12+ yrs 34 1 0 1 1 80 47
## 81 12+ yrs 35 2 2 1 0 81 54
## 84 0-5yrs 26 6 2 0 0 1 3
## 85 0-5yrs 42 1 0 0 0 2 1
## 86 0-5yrs 39 6 2 0 0 3 4
## 87 0-5yrs 34 4 0 0 1 4 2
## 88 6-11yrs 35 3 2 0 0 5 32
## 89 6-11yrs 36 4 1 0 1 6 36
## 91 6-11yrs 32 2 2 0 0 8 22
## 95 6-11yrs 37 4 1 0 1 12 37
## 96 6-11yrs 31 1 0 0 0 13 9
## 98 6-11yrs 31 2 1 0 0 15 21
## 100 6-11yrs 30 5 0 0 2 17 38
## 103 6-11yrs 44 1 0 0 0 20 17
## 104 6-11yrs 40 1 0 0 0 21 14
## 105 6-11yrs 35 2 0 0 0 22 24
## 107 6-11yrs 36 1 0 0 0 24 12
## 109 6-11yrs 40 2 0 0 0 26 27
## 110 6-11yrs 38 2 0 0 0 27 26
## 111 6-11yrs 34 3 0 0 0 28 31
## 112 6-11yrs 28 4 0 0 2 29 34
## 113 6-11yrs 30 4 1 0 1 30 35
## 114 6-11yrs 32 1 0 0 0 31 10
## 115 6-11yrs 34 2 1 0 0 32 23
## 116 6-11yrs 42 1 1 0 0 33 16
## 117 6-11yrs 32 2 0 0 0 34 22
## 118 6-11yrs 39 1 0 0 0 35 13
## 119 6-11yrs 35 2 0 0 0 36 24
## 120 6-11yrs 36 1 0 0 0 37 12
## 121 6-11yrs 34 3 2 0 0 38 31
## 124 6-11yrs 39 3 1 0 0 41 33
## 125 6-11yrs 35 1 0 0 0 42 11
## 126 6-11yrs 41 1 0 0 0 43 15
## 127 6-11yrs 37 2 0 0 0 44 25
## 129 12+ yrs 37 1 0 0 0 46 48
## 131 12+ yrs 27 4 2 0 1 48 61
## 133 12+ yrs 38 3 1 0 0 50 60
## 135 12+ yrs 36 5 1 0 1 52 62
## 139 12+ yrs 36 2 1 0 0 56 55
## 144 12+ yrs 35 2 2 0 0 61 54
## 146 12+ yrs 34 1 0 0 0 63 47
## 147 12+ yrs 31 2 0 0 0 64 53
## 149 12+ yrs 32 1 0 0 0 66 46
## 152 12+ yrs 37 3 1 0 1 69 59
## 154 12+ yrs 32 1 1 0 0 71 46
## 156 12+ yrs 31 1 0 0 1 73 45
## 158 12+ yrs 31 1 0 0 0 76 45
## 159 12+ yrs 31 2 2 0 0 77 53
## 161 12+ yrs 31 1 0 0 0 79 45
## 162 12+ yrs 34 1 0 0 0 80 47
## 163 12+ yrs 35 2 0 0 0 81 54
## 166 0-5yrs 26 6 2 0 0 1 3
## 167 0-5yrs 42 1 0 0 0 2 1
## 168 0-5yrs 39 6 2 0 0 3 4
## 169 0-5yrs 34 4 0 0 2 4 2
## 170 6-11yrs 35 3 0 0 0 5 32
## 171 6-11yrs 36 4 0 0 2 6 36
## 173 6-11yrs 32 2 0 0 1 8 22
## 177 6-11yrs 37 4 0 0 1 12 37
## 178 6-11yrs 31 1 0 0 0 13 9
## 180 6-11yrs 31 2 1 0 0 15 21
## 182 6-11yrs 30 5 1 0 2 17 38
## 185 6-11yrs 44 1 1 0 0 20 17
## 186 6-11yrs 40 1 0 0 0 21 14
## 187 6-11yrs 35 2 0 0 0 22 24
## 189 6-11yrs 36 1 0 0 1 24 12
## 191 6-11yrs 40 2 0 0 0 26 27
## 192 6-11yrs 38 2 0 0 0 27 26
## 193 6-11yrs 34 3 0 0 0 28 31
## 194 6-11yrs 28 4 2 0 1 29 34
## 195 6-11yrs 30 4 1 0 1 30 35
## 196 6-11yrs 32 1 0 0 0 31 10
## 197 6-11yrs 34 2 0 0 0 32 23
## 198 6-11yrs 42 1 0 0 0 33 16
## 199 6-11yrs 32 2 2 0 0 34 22
## 200 6-11yrs 39 1 0 0 0 35 13
## 201 6-11yrs 35 2 0 0 0 36 24
## 202 6-11yrs 36 1 0 0 0 37 12
## 203 6-11yrs 34 3 2 0 0 38 31
## 206 6-11yrs 39 3 0 0 0 41 33
## 207 6-11yrs 35 1 0 0 0 42 11
## 208 6-11yrs 41 1 0 0 0 43 15
## 209 6-11yrs 37 2 0 0 0 44 25
## 211 12+ yrs 37 1 0 0 1 46 48
## 213 12+ yrs 27 4 2 0 0 48 61
## 215 12+ yrs 38 3 1 0 0 50 60
## 217 12+ yrs 36 5 2 0 1 52 62
## 221 12+ yrs 36 2 0 0 1 56 55
## 226 12+ yrs 35 2 1 0 0 61 54
## 228 12+ yrs 34 1 0 0 0 63 47
## 229 12+ yrs 31 2 1 0 0 64 53
## 231 12+ yrs 32 1 1 0 0 66 46
## 234 12+ yrs 37 3 0 0 2 69 59
## 236 12+ yrs 32 1 0 0 0 71 46
## 238 12+ yrs 31 1 0 0 0 73 45
## 239 12+ yrs 38 6 0 0 2 74 63
## 241 12+ yrs 31 1 1 0 0 76 45
## 242 12+ yrs 31 2 0 0 1 77 53
## 244 12+ yrs 31 1 0 0 1 79 45
## 245 12+ yrs 34 1 0 0 0 80 47
## 246 12+ yrs 35 2 2 0 0 81 54#Chia dữ liệu AGE thành 4 tổ và lập bảng
cut(AGE,4)## [1] (21,26.8] (38.2,44] (38.2,44] (32.5,38.2] (32.5,38.2] (32.5,38.2]
## [7] (21,26.8] (26.8,32.5] (21,26.8] (26.8,32.5] (26.8,32.5] (32.5,38.2]
## [13] (26.8,32.5] (26.8,32.5] (26.8,32.5] (26.8,32.5] (26.8,32.5] (21,26.8]
## [19] (21,26.8] (38.2,44] (38.2,44] (32.5,38.2] (26.8,32.5] (32.5,38.2]
## [25] (26.8,32.5] (38.2,44] (32.5,38.2] (32.5,38.2] (26.8,32.5] (26.8,32.5]
## [31] (26.8,32.5] (32.5,38.2] (38.2,44] (26.8,32.5] (38.2,44] (32.5,38.2]
## [37] (32.5,38.2] (32.5,38.2] (26.8,32.5] (26.8,32.5] (38.2,44] (32.5,38.2]
## [43] (38.2,44] (32.5,38.2] (26.8,32.5] (32.5,38.2] (26.8,32.5] (26.8,32.5]
## [49] (21,26.8] (32.5,38.2] (21,26.8] (32.5,38.2] (26.8,32.5] (26.8,32.5]
## [55] (26.8,32.5] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5] (26.8,32.5]
## [61] (32.5,38.2] (21,26.8] (32.5,38.2] (26.8,32.5] (21,26.8] (26.8,32.5]
## [67] (21,26.8] (26.8,32.5] (32.5,38.2] (21,26.8] (26.8,32.5] (21,26.8]
## [73] (26.8,32.5] (32.5,38.2] (21,26.8] (26.8,32.5] (26.8,32.5] (21,26.8]
## [79] (26.8,32.5] (32.5,38.2] (32.5,38.2] (26.8,32.5] (21,26.8] (21,26.8]
## [85] (38.2,44] (38.2,44] (32.5,38.2] (32.5,38.2] (32.5,38.2] (21,26.8]
## [91] (26.8,32.5] (21,26.8] (26.8,32.5] (26.8,32.5] (32.5,38.2] (26.8,32.5]
## [97] (26.8,32.5] (26.8,32.5] (26.8,32.5] (26.8,32.5] (21,26.8] (21,26.8]
## [103] (38.2,44] (38.2,44] (32.5,38.2] (26.8,32.5] (32.5,38.2] (26.8,32.5]
## [109] (38.2,44] (32.5,38.2] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5]
## [115] (32.5,38.2] (38.2,44] (26.8,32.5] (38.2,44] (32.5,38.2] (32.5,38.2]
## [121] (32.5,38.2] (26.8,32.5] (26.8,32.5] (38.2,44] (32.5,38.2] (38.2,44]
## [127] (32.5,38.2] (26.8,32.5] (32.5,38.2] (26.8,32.5] (26.8,32.5] (21,26.8]
## [133] (32.5,38.2] (21,26.8] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5]
## [139] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5] (26.8,32.5] (32.5,38.2]
## [145] (21,26.8] (32.5,38.2] (26.8,32.5] (21,26.8] (26.8,32.5] (21,26.8]
## [151] (26.8,32.5] (32.5,38.2] (21,26.8] (26.8,32.5] (21,26.8] (26.8,32.5]
## [157] (21,26.8] (26.8,32.5] (26.8,32.5] (21,26.8] (26.8,32.5] (32.5,38.2]
## [163] (32.5,38.2] (26.8,32.5] (21,26.8] (21,26.8] (38.2,44] (38.2,44]
## [169] (32.5,38.2] (32.5,38.2] (32.5,38.2] (21,26.8] (26.8,32.5] (21,26.8]
## [175] (26.8,32.5] (26.8,32.5] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5]
## [181] (26.8,32.5] (26.8,32.5] (21,26.8] (21,26.8] (38.2,44] (38.2,44]
## [187] (32.5,38.2] (26.8,32.5] (32.5,38.2] (26.8,32.5] (38.2,44] (32.5,38.2]
## [193] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5] (32.5,38.2] (38.2,44]
## [199] (26.8,32.5] (38.2,44] (32.5,38.2] (32.5,38.2] (32.5,38.2] (26.8,32.5]
## [205] (26.8,32.5] (38.2,44] (32.5,38.2] (38.2,44] (32.5,38.2] (26.8,32.5]
## [211] (32.5,38.2] (26.8,32.5] (26.8,32.5] (21,26.8] (32.5,38.2] (21,26.8]
## [217] (32.5,38.2] (26.8,32.5] (26.8,32.5] (26.8,32.5] (32.5,38.2] (26.8,32.5]
## [223] (26.8,32.5] (26.8,32.5] (26.8,32.5] (32.5,38.2] (21,26.8] (32.5,38.2]
## [229] (26.8,32.5] (21,26.8] (26.8,32.5] (21,26.8] (26.8,32.5] (32.5,38.2]
## [235] (21,26.8] (26.8,32.5] (21,26.8] (26.8,32.5] (32.5,38.2] (21,26.8]
## [241] (26.8,32.5] (26.8,32.5] (21,26.8] (26.8,32.5] (32.5,38.2] (32.5,38.2]
## [247] (26.8,32.5] (21,26.8]
## Levels: (21,26.8] (26.8,32.5] (32.5,38.2] (38.2,44]table(cut(AGE,4))##
## (21,26.8] (26.8,32.5] (32.5,38.2] (38.2,44]
## 45 105 71 27#Lập bảng tần số của AGE được chia thành 4 tổ với lần lượt với các biến Education, Parity, Induced và Case
table(cut(AGE,4),i$E)##
## 0-5yrs 6-11yrs 12+ yrs
## (21,26.8] 3 12 30
## (26.8,32.5] 0 48 57
## (32.5,38.2] 3 39 29
## (38.2,44] 6 21 0table(cut(AGE,4),i$Pa)##
## 1 2 3 4 5 6
## (21,26.8] 27 9 6 0 0 3
## (26.8,32.5] 36 45 12 9 3 0
## (32.5,38.2] 18 24 15 9 3 2
## (38.2,44] 18 3 3 0 0 3table(cut(AGE,4),i$I)##
## 0 1 2
## (21,26.8] 24 15 6
## (26.8,32.5] 56 32 17
## (32.5,38.2] 45 15 11
## (38.2,44] 18 6 3table(cut(AGE,4),i$C)##
## 0 1
## (21,26.8] 30 15
## (26.8,32.5] 70 35
## (32.5,38.2] 47 24
## (38.2,44] 18 9#Chia dữ liệu Parity thành 4 tổ và lập bảng
cut(PAR,4)## [1] (4.75,6] (0.995,2.25] (4.75,6] (3.5,4.75] (2.25,3.5]
## [6] (3.5,4.75] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [11] (0.995,2.25] (3.5,4.75] (0.995,2.25] (2.25,3.5] (0.995,2.25]
## [16] (0.995,2.25] (4.75,6] (0.995,2.25] (2.25,3.5] (0.995,2.25]
## [21] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [26] (0.995,2.25] (0.995,2.25] (2.25,3.5] (3.5,4.75] (3.5,4.75]
## [31] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [36] (0.995,2.25] (0.995,2.25] (2.25,3.5] (2.25,3.5] (0.995,2.25]
## [41] (2.25,3.5] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [46] (0.995,2.25] (0.995,2.25] (3.5,4.75] (0.995,2.25] (2.25,3.5]
## [51] (2.25,3.5] (4.75,6] (2.25,3.5] (0.995,2.25] (0.995,2.25]
## [56] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [61] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [66] (0.995,2.25] (0.995,2.25] (2.25,3.5] (2.25,3.5] (0.995,2.25]
## [71] (0.995,2.25] (0.995,2.25] (0.995,2.25] (4.75,6] (0.995,2.25]
## [76] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [81] (0.995,2.25] (0.995,2.25] (0.995,2.25] (4.75,6] (0.995,2.25]
## [86] (4.75,6] (3.5,4.75] (2.25,3.5] (3.5,4.75] (0.995,2.25]
## [91] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (3.5,4.75]
## [96] (0.995,2.25] (2.25,3.5] (0.995,2.25] (0.995,2.25] (4.75,6]
## [101] (0.995,2.25] (2.25,3.5] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [106] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [111] (2.25,3.5] (3.5,4.75] (3.5,4.75] (0.995,2.25] (0.995,2.25]
## [116] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [121] (2.25,3.5] (2.25,3.5] (0.995,2.25] (2.25,3.5] (0.995,2.25]
## [126] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [131] (3.5,4.75] (0.995,2.25] (2.25,3.5] (2.25,3.5] (4.75,6]
## [136] (2.25,3.5] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [141] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [146] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [151] (2.25,3.5] (2.25,3.5] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [156] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [161] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [166] (4.75,6] (0.995,2.25] (4.75,6] (3.5,4.75] (2.25,3.5]
## [171] (3.5,4.75] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [176] (0.995,2.25] (3.5,4.75] (0.995,2.25] (2.25,3.5] (0.995,2.25]
## [181] (0.995,2.25] (4.75,6] (0.995,2.25] (2.25,3.5] (0.995,2.25]
## [186] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [191] (0.995,2.25] (0.995,2.25] (2.25,3.5] (3.5,4.75] (3.5,4.75]
## [196] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [201] (0.995,2.25] (0.995,2.25] (2.25,3.5] (2.25,3.5] (0.995,2.25]
## [206] (2.25,3.5] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [211] (0.995,2.25] (0.995,2.25] (3.5,4.75] (0.995,2.25] (2.25,3.5]
## [216] (2.25,3.5] (4.75,6] (2.25,3.5] (0.995,2.25] (0.995,2.25]
## [221] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [226] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [231] (0.995,2.25] (0.995,2.25] (2.25,3.5] (2.25,3.5] (0.995,2.25]
## [236] (0.995,2.25] (0.995,2.25] (0.995,2.25] (4.75,6] (0.995,2.25]
## [241] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## [246] (0.995,2.25] (0.995,2.25] (0.995,2.25]
## Levels: (0.995,2.25] (2.25,3.5] (3.5,4.75] (4.75,6]table(cut(PAR,4))##
## (0.995,2.25] (2.25,3.5] (3.5,4.75] (4.75,6]
## 180 36 18 14#Gán biến Po của data(i)vào POO, chia dữ liệu POO thành 4 tổ và lập bảng
POO <- i$Po
cut(POO,4)## [1] (0.938,16.5] (0.938,16.5] (0.938,16.5] (0.938,16.5] (16.5,32]
## [6] (32,47.5] (0.938,16.5] (16.5,32] (0.938,16.5] (16.5,32]
## [11] (16.5,32] (32,47.5] (0.938,16.5] (16.5,32] (16.5,32]
## [16] (16.5,32] (32,47.5] (0.938,16.5] (16.5,32] (16.5,32]
## [21] (0.938,16.5] (16.5,32] (16.5,32] (0.938,16.5] (16.5,32]
## [26] (16.5,32] (16.5,32] (16.5,32] (32,47.5] (32,47.5]
## [31] (0.938,16.5] (16.5,32] (0.938,16.5] (16.5,32] (0.938,16.5]
## [36] (16.5,32] (0.938,16.5] (16.5,32] (16.5,32] (0.938,16.5]
## [41] (32,47.5] (0.938,16.5] (0.938,16.5] (16.5,32] (32,47.5]
## [46] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1]
## [51] (47.5,63.1] (47.5,63.1] (47.5,63.1] (32,47.5] (47.5,63.1]
## [56] (47.5,63.1] (47.5,63.1] (47.5,63.1] (32,47.5] (47.5,63.1]
## [61] (47.5,63.1] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1]
## [66] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1] (32,47.5]
## [71] (32,47.5] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1]
## [76] (32,47.5] (47.5,63.1] (32,47.5] (32,47.5] (32,47.5]
## [81] (47.5,63.1] (32,47.5] (32,47.5] (0.938,16.5] (0.938,16.5]
## [86] (0.938,16.5] (0.938,16.5] (16.5,32] (32,47.5] (0.938,16.5]
## [91] (16.5,32] (0.938,16.5] (16.5,32] (16.5,32] (32,47.5]
## [96] (0.938,16.5] (16.5,32] (16.5,32] (16.5,32] (32,47.5]
## [101] (0.938,16.5] (16.5,32] (16.5,32] (0.938,16.5] (16.5,32]
## [106] (16.5,32] (0.938,16.5] (16.5,32] (16.5,32] (16.5,32]
## [111] (16.5,32] (32,47.5] (32,47.5] (0.938,16.5] (16.5,32]
## [116] (0.938,16.5] (16.5,32] (0.938,16.5] (16.5,32] (0.938,16.5]
## [121] (16.5,32] (16.5,32] (0.938,16.5] (32,47.5] (0.938,16.5]
## [126] (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1] (47.5,63.1]
## [131] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1]
## [136] (47.5,63.1] (32,47.5] (47.5,63.1] (47.5,63.1] (47.5,63.1]
## [141] (47.5,63.1] (32,47.5] (47.5,63.1] (47.5,63.1] (32,47.5]
## [146] (32,47.5] (47.5,63.1] (47.5,63.1] (32,47.5] (32,47.5]
## [151] (47.5,63.1] (47.5,63.1] (32,47.5] (32,47.5] (32,47.5]
## [156] (32,47.5] (47.5,63.1] (32,47.5] (47.5,63.1] (32,47.5]
## [161] (32,47.5] (32,47.5] (47.5,63.1] (32,47.5] (32,47.5]
## [166] (0.938,16.5] (0.938,16.5] (0.938,16.5] (0.938,16.5] (16.5,32]
## [171] (32,47.5] (0.938,16.5] (16.5,32] (0.938,16.5] (16.5,32]
## [176] (16.5,32] (32,47.5] (0.938,16.5] (16.5,32] (16.5,32]
## [181] (16.5,32] (32,47.5] (0.938,16.5] (16.5,32] (16.5,32]
## [186] (0.938,16.5] (16.5,32] (16.5,32] (0.938,16.5] (16.5,32]
## [191] (16.5,32] (16.5,32] (16.5,32] (32,47.5] (32,47.5]
## [196] (0.938,16.5] (16.5,32] (0.938,16.5] (16.5,32] (0.938,16.5]
## [201] (16.5,32] (0.938,16.5] (16.5,32] (16.5,32] (0.938,16.5]
## [206] (32,47.5] (0.938,16.5] (0.938,16.5] (16.5,32] (32,47.5]
## [211] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1] (47.5,63.1]
## [216] (47.5,63.1] (47.5,63.1] (47.5,63.1] (32,47.5] (47.5,63.1]
## [221] (47.5,63.1] (47.5,63.1] (47.5,63.1] (32,47.5] (47.5,63.1]
## [226] (47.5,63.1] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1]
## [231] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1] (32,47.5]
## [236] (32,47.5] (32,47.5] (32,47.5] (47.5,63.1] (47.5,63.1]
## [241] (32,47.5] (47.5,63.1] (32,47.5] (32,47.5] (32,47.5]
## [246] (47.5,63.1] (32,47.5] (32,47.5]
## Levels: (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1]table(cut(POO,4))##
## (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1]
## 51 63 69 65#Lập bảng tần số của dữ liệu Parity với lần lượt các biến Education, Case, Spontaneous và bảng tần số Pooled.stratum
table(cut(PAR,4),i$E)##
## 0-5yrs 6-11yrs 12+ yrs
## (0.995,2.25] 3 84 93
## (2.25,3.5] 0 21 15
## (3.5,4.75] 3 12 3
## (4.75,6] 6 3 5table(cut(PAR,4),i$C)##
## 0 1
## (0.995,2.25] 120 60
## (2.25,3.5] 24 12
## (3.5,4.75] 12 6
## (4.75,6] 9 5table(cut(PAR,4),i$Sp)##
## 0 1 2
## (0.995,2.25] 116 48 16
## (2.25,3.5] 16 10 10
## (3.5,4.75] 4 10 4
## (4.75,6] 5 3 6table(cut(PAR,4),cut(POO,4))##
## (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1]
## (0.995,2.25] 42 45 51 42
## (2.25,3.5] 0 18 3 15
## (3.5,4.75] 3 0 12 3
## (4.75,6] 6 0 3 5#Lập bảng tần số của POO với lần lượt các biến Education, Case, Spontaneous và Parity
table(cut(POO,4),i$E)##
## 0-5yrs 6-11yrs 12+ yrs
## (0.938,16.5] 12 39 0
## (16.5,32] 0 63 0
## (32,47.5] 0 18 51
## (47.5,63.1] 0 0 65table(cut(POO,4),i$C)##
## 0 1
## (0.938,16.5] 34 17
## (16.5,32] 42 21
## (32,47.5] 46 23
## (47.5,63.1] 43 22table(cut(POO,4),i$Sp)##
## 0 1 2
## (0.938,16.5] 39 10 2
## (16.5,32] 38 15 10
## (32,47.5] 32 31 6
## (47.5,63.1] 32 15 18table(cut(POO,4),i$Pa)##
## 1 2 3 4 5 6
## (0.938,16.5] 42 0 0 3 0 6
## (16.5,32] 3 42 18 0 0 0
## (32,47.5] 51 0 3 12 3 0
## (47.5,63.1] 3 39 15 3 3 2#Lập bảng tần số giữa các dữ liệu AGE, POO và PAR với nhau
table(cut(AGE,4),cut(POO,4))##
## (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1]
## (21,26.8] 12 3 18 12
## (26.8,32.5] 9 30 36 30
## (32.5,38.2] 12 24 12 23
## (38.2,44] 18 6 3 0table(cut(AGE,4),cut(PAR,4))##
## (0.995,2.25] (2.25,3.5] (3.5,4.75] (4.75,6]
## (21,26.8] 36 6 0 3
## (26.8,32.5] 81 12 9 3
## (32.5,38.2] 42 15 9 5
## (38.2,44] 21 3 0 3table(cut(PAR,4),cut(POO,4))##
## (0.938,16.5] (16.5,32] (32,47.5] (47.5,63.1]
## (0.995,2.25] 42 45 51 42
## (2.25,3.5] 0 18 3 15
## (3.5,4.75] 3 0 12 3
## (4.75,6] 6 0 3 5#Tính logarit của các biến A, Pa, Po trong data(i)
i$logA <- log(i$A)
i$logA## [1] 3.258097 3.737670 3.663562 3.526361 3.555348 3.583519 3.135494 3.465736
## [9] 3.044522 3.332205 3.367296 3.610918 3.433987 3.367296 3.433987 3.295837
## [17] 3.401197 3.258097 3.218876 3.784190 3.688879 3.555348 3.332205 3.583519
## [25] 3.295837 3.688879 3.637586 3.526361 3.332205 3.401197 3.465736 3.526361
## [33] 3.737670 3.465736 3.663562 3.555348 3.583519 3.526361 3.401197 3.332205
## [41] 3.663562 3.555348 3.713572 3.610918 3.401197 3.610918 3.332205 3.295837
## [49] 3.258097 3.637586 3.178054 3.583519 3.295837 3.332205 3.367296 3.583519
## [57] 3.332205 3.332205 3.332205 3.295837 3.555348 3.218876 3.526361 3.433987
## [65] 3.258097 3.465736 3.044522 3.332205 3.610918 3.218876 3.465736 3.218876
## [73] 3.433987 3.637586 3.258097 3.433987 3.433987 3.218876 3.433987 3.526361
## [81] 3.555348 3.367296 3.135494 3.258097 3.737670 3.663562 3.526361 3.555348
## [89] 3.583519 3.135494 3.465736 3.044522 3.332205 3.367296 3.610918 3.433987
## [97] 3.367296 3.433987 3.295837 3.401197 3.258097 3.218876 3.784190 3.688879
## [105] 3.555348 3.332205 3.583519 3.295837 3.688879 3.637586 3.526361 3.332205
## [113] 3.401197 3.465736 3.526361 3.737670 3.465736 3.663562 3.555348 3.583519
## [121] 3.526361 3.401197 3.332205 3.663562 3.555348 3.713572 3.610918 3.401197
## [129] 3.610918 3.332205 3.295837 3.258097 3.637586 3.178054 3.583519 3.295837
## [137] 3.332205 3.367296 3.583519 3.332205 3.332205 3.332205 3.295837 3.555348
## [145] 3.218876 3.526361 3.433987 3.258097 3.465736 3.044522 3.332205 3.610918
## [153] 3.218876 3.465736 3.218876 3.433987 3.258097 3.433987 3.433987 3.218876
## [161] 3.433987 3.526361 3.555348 3.367296 3.135494 3.258097 3.737670 3.663562
## [169] 3.526361 3.555348 3.583519 3.135494 3.465736 3.044522 3.332205 3.367296
## [177] 3.610918 3.433987 3.367296 3.433987 3.295837 3.401197 3.258097 3.218876
## [185] 3.784190 3.688879 3.555348 3.332205 3.583519 3.295837 3.688879 3.637586
## [193] 3.526361 3.332205 3.401197 3.465736 3.526361 3.737670 3.465736 3.663562
## [201] 3.555348 3.583519 3.526361 3.401197 3.332205 3.663562 3.555348 3.713572
## [209] 3.610918 3.401197 3.610918 3.332205 3.295837 3.258097 3.637586 3.178054
## [217] 3.583519 3.295837 3.332205 3.367296 3.583519 3.332205 3.332205 3.332205
## [225] 3.295837 3.555348 3.218876 3.526361 3.433987 3.258097 3.465736 3.044522
## [233] 3.332205 3.610918 3.218876 3.465736 3.218876 3.433987 3.637586 3.258097
## [241] 3.433987 3.433987 3.218876 3.433987 3.526361 3.555348 3.367296 3.135494i$logPAR <- log(i$Pa)
i$logPAR## [1] 1.7917595 0.0000000 1.7917595 1.3862944 1.0986123 1.3862944 0.0000000
## [8] 0.6931472 0.0000000 0.6931472 0.6931472 1.3862944 0.0000000 1.0986123
## [15] 0.6931472 0.6931472 1.6094379 0.0000000 1.0986123 0.0000000 0.0000000
## [22] 0.6931472 0.6931472 0.0000000 0.6931472 0.6931472 0.6931472 1.0986123
## [29] 1.3862944 1.3862944 0.0000000 0.6931472 0.0000000 0.6931472 0.0000000
## [36] 0.6931472 0.0000000 1.0986123 1.0986123 0.0000000 1.0986123 0.0000000
## [43] 0.0000000 0.6931472 0.0000000 0.0000000 0.6931472 1.3862944 0.6931472
## [50] 1.0986123 1.0986123 1.6094379 1.0986123 0.0000000 0.6931472 0.6931472
## [57] 0.6931472 0.6931472 0.0000000 0.6931472 0.6931472 0.0000000 0.0000000
## [64] 0.6931472 0.6931472 0.0000000 0.0000000 1.0986123 1.0986123 0.0000000
## [71] 0.0000000 0.0000000 0.0000000 1.7917595 0.6931472 0.0000000 0.6931472
## [78] 0.0000000 0.0000000 0.0000000 0.6931472 0.0000000 0.0000000 1.7917595
## [85] 0.0000000 1.7917595 1.3862944 1.0986123 1.3862944 0.0000000 0.6931472
## [92] 0.0000000 0.6931472 0.6931472 1.3862944 0.0000000 1.0986123 0.6931472
## [99] 0.6931472 1.6094379 0.0000000 1.0986123 0.0000000 0.0000000 0.6931472
## [106] 0.6931472 0.0000000 0.6931472 0.6931472 0.6931472 1.0986123 1.3862944
## [113] 1.3862944 0.0000000 0.6931472 0.0000000 0.6931472 0.0000000 0.6931472
## [120] 0.0000000 1.0986123 1.0986123 0.0000000 1.0986123 0.0000000 0.0000000
## [127] 0.6931472 0.0000000 0.0000000 0.6931472 1.3862944 0.6931472 1.0986123
## [134] 1.0986123 1.6094379 1.0986123 0.0000000 0.6931472 0.6931472 0.6931472
## [141] 0.6931472 0.0000000 0.6931472 0.6931472 0.0000000 0.0000000 0.6931472
## [148] 0.6931472 0.0000000 0.0000000 1.0986123 1.0986123 0.0000000 0.0000000
## [155] 0.0000000 0.0000000 0.6931472 0.0000000 0.6931472 0.0000000 0.0000000
## [162] 0.0000000 0.6931472 0.0000000 0.0000000 1.7917595 0.0000000 1.7917595
## [169] 1.3862944 1.0986123 1.3862944 0.0000000 0.6931472 0.0000000 0.6931472
## [176] 0.6931472 1.3862944 0.0000000 1.0986123 0.6931472 0.6931472 1.6094379
## [183] 0.0000000 1.0986123 0.0000000 0.0000000 0.6931472 0.6931472 0.0000000
## [190] 0.6931472 0.6931472 0.6931472 1.0986123 1.3862944 1.3862944 0.0000000
## [197] 0.6931472 0.0000000 0.6931472 0.0000000 0.6931472 0.0000000 1.0986123
## [204] 1.0986123 0.0000000 1.0986123 0.0000000 0.0000000 0.6931472 0.0000000
## [211] 0.0000000 0.6931472 1.3862944 0.6931472 1.0986123 1.0986123 1.6094379
## [218] 1.0986123 0.0000000 0.6931472 0.6931472 0.6931472 0.6931472 0.0000000
## [225] 0.6931472 0.6931472 0.0000000 0.0000000 0.6931472 0.6931472 0.0000000
## [232] 0.0000000 1.0986123 1.0986123 0.0000000 0.0000000 0.0000000 0.0000000
## [239] 1.7917595 0.6931472 0.0000000 0.6931472 0.0000000 0.0000000 0.0000000
## [246] 0.6931472 0.0000000 0.0000000i$logPo <- log(i$Po)
i$logPo## [1] 1.0986123 0.0000000 1.3862944 0.6931472 3.4657359 3.5835189 1.7917595
## [8] 3.0910425 1.6094379 2.9444390 2.9957323 3.6109179 2.1972246 3.3672958
## [15] 3.0445224 2.8903718 3.6375862 1.9459101 3.3322045 2.8332133 2.6390573
## [22] 3.1780538 2.9444390 2.4849066 2.8903718 3.2958369 3.2580965 3.4339872
## [29] 3.5263605 3.5553481 2.3025851 3.1354942 2.7725887 3.0910425 2.5649494
## [36] 3.1780538 2.4849066 3.4339872 3.4011974 2.0794415 3.4965076 2.3978953
## [43] 2.7080502 3.2188758 3.7841896 3.8712010 3.9318256 4.1108739 3.8918203
## [50] 4.0943446 4.0253517 4.1271344 4.0430513 3.7376696 3.9512437 4.0073332
## [57] 3.9318256 3.9318256 3.7376696 3.9120230 3.9889840 3.7135721 3.8501476
## [64] 3.9702919 3.8918203 3.8286414 3.6635616 4.0604430 4.0775374 3.7135721
## [71] 3.8286414 3.7135721 3.8066625 4.1431347 3.8918203 3.8066625 3.9702919
## [78] 3.7135721 3.8066625 3.8501476 3.9889840 3.7612001 3.6888795 1.0986123
## [85] 0.0000000 1.3862944 0.6931472 3.4657359 3.5835189 1.7917595 3.0910425
## [92] 1.6094379 2.9444390 2.9957323 3.6109179 2.1972246 3.3672958 3.0445224
## [99] 2.8903718 3.6375862 1.9459101 3.3322045 2.8332133 2.6390573 3.1780538
## [106] 2.9444390 2.4849066 2.8903718 3.2958369 3.2580965 3.4339872 3.5263605
## [113] 3.5553481 2.3025851 3.1354942 2.7725887 3.0910425 2.5649494 3.1780538
## [120] 2.4849066 3.4339872 3.4011974 2.0794415 3.4965076 2.3978953 2.7080502
## [127] 3.2188758 3.7841896 3.8712010 3.9318256 4.1108739 3.8918203 4.0943446
## [134] 4.0253517 4.1271344 4.0430513 3.7376696 3.9512437 4.0073332 3.9318256
## [141] 3.9318256 3.7376696 3.9120230 3.9889840 3.7135721 3.8501476 3.9702919
## [148] 3.8918203 3.8286414 3.6635616 4.0604430 4.0775374 3.7135721 3.8286414
## [155] 3.7135721 3.8066625 3.8918203 3.8066625 3.9702919 3.7135721 3.8066625
## [162] 3.8501476 3.9889840 3.7612001 3.6888795 1.0986123 0.0000000 1.3862944
## [169] 0.6931472 3.4657359 3.5835189 1.7917595 3.0910425 1.6094379 2.9444390
## [176] 2.9957323 3.6109179 2.1972246 3.3672958 3.0445224 2.8903718 3.6375862
## [183] 1.9459101 3.3322045 2.8332133 2.6390573 3.1780538 2.9444390 2.4849066
## [190] 2.8903718 3.2958369 3.2580965 3.4339872 3.5263605 3.5553481 2.3025851
## [197] 3.1354942 2.7725887 3.0910425 2.5649494 3.1780538 2.4849066 3.4339872
## [204] 3.4011974 2.0794415 3.4965076 2.3978953 2.7080502 3.2188758 3.7841896
## [211] 3.8712010 3.9318256 4.1108739 3.8918203 4.0943446 4.0253517 4.1271344
## [218] 4.0430513 3.7376696 3.9512437 4.0073332 3.9318256 3.9318256 3.7376696
## [225] 3.9120230 3.9889840 3.7135721 3.8501476 3.9702919 3.8918203 3.8286414
## [232] 3.6635616 4.0604430 4.0775374 3.7135721 3.8286414 3.7135721 3.8066625
## [239] 4.1431347 3.8918203 3.8066625 3.9702919 3.7135721 3.8066625 3.8501476
## [246] 3.9889840 3.7612001 3.6888795