About

Qualitative Descriptive Analytics aims to gather an in-depth understanding of the underlying reasons and motivations for an event or observation. It is typically represented with visuals or charts.

Quantitative Descriptive Analytics focuses on investigating a phenomenon via statistical, mathematical, and computational techniques. It aims to quantify an event with metrics and numbers.

In this project, we will explore quantitative analytics using the babies data set provided. This data is collected to understand different factors effecting baby weights at birth. Below are the explanations of the variables:

Bwt: Baby weights in ounces Gestation: Length of pregnancy in days Height: Mother’s height in inches Smoke: =1 if mother is smoker, = 0 Nonsmoker Age: Mother’s age in years Weight: Mother’s pregnancy weight Parity: = 0 if the baby is first born, =1 otherwise babies: number of babies before this birth

Read this worksheet carefully and follow the instructions to complete the tasks and answer any questions. Submit your work as an HTML or PDF or Word document to Sakai. If you are submitting as a team project only one submission from one of the team members will be sufficient but in this case you need to input the names of the team members in the above Title section.

Task 1: Testing for missing data (1 point)

babies_old data set will be used to check if it contains any missing data. First read the babies_old file and make sure you correctly read the file with the above mentioned variables.

#Read babies_old file and check the variables
setwd("C:/Users/kkobayashi/Downloads/LAB 02") # Change working directory
mydata_old <- read.csv("babies_old.csv") 
dim(mydata_old)
## [1] 1174    8

Find NA’s in babies_old file

#Check your file for missing data
missing_data <- is.na(mydata_old)
missing_count <- which(is.na(mydata_old$Count))
missing_bwt <- which(is.na(mydata_old$bwt))
missing_gest <- which(is.na(mydata_old$gestation))
missing_parity <- which(is.na(mydata_old$parity))
missing_height <- which(is.na(mydata_old$height))
missing_weight <- which(is.na(mydata_old$weight))
missing_smoke <- which(is.na(mydata_old$smoke))
missing_age <- which(is.na(mydata_old$age))

Implement omission strategy to remove observations with missing values. (0.5 point)

#implement omission strategy and check if all the observations with NA's were omitted.
CompletedData <-na.omit(mydata_old)
dim(CompletedData)
## [1] 1163    8

Implement mean imputation strategy (0.5 point)

#Implement imputation strategy
mean_weight <- mean(mydata_old$weight, na.rm = TRUE)
mydata_old$weight[is.na(mydata_old$weight)] <- mean_weight
mydata_old
##      Count bwt gestation parity height   weight smoke age
## 1        1 120       284      0     62 100.0000     0  27
## 2        2 113       282      0     64 128.4996     0  33
## 3        3 128       279      0     64 115.0000     1  28
## 4        4 108       282      0     67 125.0000     1  23
## 5        5 136       286      0     62  93.0000     0  25
## 6        6 138       244      0     62 178.0000     0  33
## 7        7 132       245      0     65 140.0000     0  23
## 8        8 120       289      0     62 125.0000     0  25
## 9        9 143       299      0     66 136.0000     1  30
## 10      10 140       351      0     68 120.0000     0  27
## 11      11 144       282      0     64 124.0000     1  32
## 12      12 141       279      0     63 128.0000     1  23
## 13      13 110       281      0     61 128.4996     1  36
## 14      14 114       273      0     63 154.0000     0  30
## 15      15 115       285      0     63 130.0000     0  38
## 16      16  92       255      0     65 125.0000     1  25
## 17      17 115       261      0     60 125.0000     1  33
## 18      18 144       261      0     68 170.0000     0  33
## 19      19 119       288      0     66 142.0000     1  43
## 20      20 105       270      0     56  93.0000     0  22
## 21      21 115       274      0     67 175.0000     1  27
## 22      22 137       287      0     66 145.0000     0  25
## 23      23 122       276      0     68 182.0000     0  30
## 24      24 131       294      0     65 122.0000     0  23
## 25      25 103       261      0     65 112.0000     1  27
## 26      26 146       280      0     58 106.0000     0  26
## 27      27 114       266      0     65 175.0000     1  20
## 28      28 125       292      0     65 125.0000     0  32
## 29      29 114       274      0     66 132.0000     1  28
## 30      30 122       270      0     61 105.0000     0  26
## 31      31  93       278      0     61 146.0000     0  34
## 32      32 130       268      0     66 123.0000     0  30
## 33      33 119       275      0     60 105.0000     0  23
## 34      34 113       281      0     65 120.0000     0  24
## 35      35 134       283      0     67 130.0000     0  22
## 36      36 107       279      0     63 115.0000     0  24
## 37      37 134       288      0     63  92.0000     1  23
## 38      38 122       267      0     65 101.0000     1  27
## 39      39 129       293      0     61 160.0000     0  30
## 40      40 110       278      0     63 177.0000     0  23
## 41      41 111       270      0     61 119.0000     0  27
## 42      42  87       248      0     65 130.0000     1  37
## 43      43 143       274      0     63 110.0000     1  27
## 44      44 155       294      0     66 150.0000     0  32
## 45      45 110       272      0     60  90.0000     0  25
## 46      46 122       275      0     66 147.0000     0  26
## 47      47 145       291      0     63 119.0000     1  26
## 48      48 115       258      0     62 130.0000     0  26
## 49      49 108       283      0     65 148.0000     1  31
## 50      50 102       282      0     61 110.0000     0  28
## 51      51 143       286      0     64 126.0000     0  31
## 52      52 146       267      0     67 132.0000     0  30
## 53      53 124       275      0     60 128.4996     0  22
## 54      54 124       278      0     70 145.0000     1  26
## 55      55 145       257      0     65 140.0000     0  33
## 56      56 106       273      0     60 116.0000     0  28
## 57      57  75       232      0     61 110.0000     0  33
## 58      58 107       273      0     61  96.0000     0  24
## 59      59 124       288      0     67 118.0000     0  22
## 60      60 122       280      0     65 125.0000     1  23
## 61      61 101       245      0     63 130.0000     1  23
## 62      62 128       283      0     63 125.0000     1  28
## 63      63 104       282      0     65 115.0000     1  36
## 64      64  97       246      0     63 150.0000     0  37
## 65      65 137       274      0     69 137.0000     1  26
## 66      66 103       273      0     63 170.0000     1  31
## 67      67 142       276      0     63 170.0000     0  38
## 68      68 130       289      0     66 130.0000     0  27
## 69      69 156       292      0     63 118.0000     0  26
## 70      70 133       284      0     66 125.0000     1  25
## 71      71 120       274      0     62 120.0000     0  24
## 72      72  91       270      0     60 149.0000     1  24
## 73      73 127       274      0     62 110.0000     0  21
## 74      74 153       286      0     63 107.0000     1  26
## 75      75 121       276      0     63 130.0000     0  39
## 76      76 120       277      0     63 126.0000     0  27
## 77      77  99       272      0     62 103.0000     1  27
## 78      78 149       293      0     65 116.0000     0  35
## 79      79 129       280      0     64 104.0000     0  23
## 80      80 139       292      0     68 135.0000     0  25
## 81      81 114       274      0     67 148.0000     1  33
## 82      82 138       287      0     66 145.0000     0  30
## 83      83 138       294      0     65 117.0000     0  32
## 84      84 131       296      0     63 143.0000     0  37
## 85      85 125       305      0     70 196.0000     1  22
## 86      86 128       281      0     59 117.0000     0  33
## 87      87 134       268      0     62 112.0000     0  28
## 88      88 114       271      0     60 104.0000     0  27
## 89      89  85       278      0     61 103.0000     1  23
## 90      90 135       282      0     64 100.0000     0  22
## 91      91  87       255      0     61 128.4996     1  28
## 92      92 125       302      0     62 162.0000     0  37
## 93      93 105       254      0     64 137.0000     0  29
## 94      94 120       279      0     60 121.0000     1  27
## 95      95 119       274      0     64 120.0000     0  33
## 96      96 107       280      0     65 117.0000     1  36
## 97      97 119       273      0     61 108.0000     1  24
## 98      98 133       279      0     66 140.0000     0  37
## 99      99 155       287      0     66 143.0000     0  33
## 100    100 126       273      0     65 150.0000     0  22
## 101    101 129       303      0     64 125.0000     0  27
## 102    102 137       274      0     65 154.0000     0  29
## 103    103 125       302      0     65 125.0000     0  28
## 104    104  91       255      0     67 136.0000     1  19
## 105    105  95       279      0     66 145.0000     1  22
## 106    106 118       276      0     64 114.0000     0  29
## 107    107 141       278      0     66 109.0000     1  33
## 108    108 131       283      0     67 215.0000     0  25
## 109    109 121       264      0     66 145.0000     0  32
## 110    110 100       243      0     65 170.0000     1  39
## 111    111 131       288      0     61 103.0000     0  24
## 112    112 118       284      0     66 133.0000     0  26
## 113    113 152       288      0     67 130.0000     0  35
## 114    114 121       284      0     69 155.0000     0  34
## 115    115 117       276      0     69 150.0000     0  31
## 116    116 115       283      0     61 150.0000     1  25
## 117    117 112       277      0     65 110.0000     0  23
## 118    118  94       267      0     62 120.0000     1  30
## 119    119 109       272      0     66 154.0000     0  35
## 120    120 132       225      0     67 148.0000     0  28
## 121    121 117       278      0     62 103.0000     0  25
## 122    122 101       266      0     67 110.0000     1  20
## 123    123 112       294      0     64 125.0000     1  25
## 124    124 128       283      0     60 100.0000     0  24
## 125    125 128       279      0     66 147.0000     1  25
## 126    126 117       258      0     64 120.0000     0  31
## 127    127 134       278      0     69 135.0000     0  24
## 128    128 127       284      0     65 145.0000     0  28
## 129    129  93       269      0     65 104.0000     1  21
## 130    130 122       275      0     65 165.0000     0  27
## 131    131 100       265      0     62 107.0000     1  39
## 132    132 147       293      0     65 123.0000     0  32
## 133    133 120       299      0     65 110.0000     0  25
## 134    134 144       277      0     63 127.0000     0  30
## 135    135 105       268      0     61 115.0000     1  32
## 136    136 136       276      0     66 155.0000     0  23
## 137    137 102       262      0     63 125.0000     0  24
## 138    138 160       300      0     71 175.0000     1  29
## 139    139 113       275      0     68 140.0000     1  24
## 140    140 126       282      0     66 250.0000     0  38
## 141    141 126       271      0     68 148.0000     0  29
## 142    142 115       278      0     61 128.0000     0  29
## 143    143 119       284      0     66 132.0000     0  20
## 144    144 123       318      0     64 152.0000     0  21
## 145    145 118       282      0     68 135.0000     1  22
## 146    146 133       287      0     60 104.0000     1  24
## 147    147 134       290      0     60 121.0000     0  22
## 148    148 144       288      0     67 111.0000     0  21
## 149    149 111       273      0     62 138.0000     0  43
## 150    150 125       262      0     66 190.0000     0  36
## 151    151 135       296      0     63 123.0000     0  30
## 152    152 134       289      0     63 125.0000     0  22
## 153    153 116       289      0     65 160.0000     1  22
## 154    154 129       291      0     69 123.0000     0  29
## 155    155 113       301      0     67 105.0000     1  26
## 156    156 131       295      0     65 123.0000     1  23
## 157    157 121       272      0     62 109.0000     0  22
## 158    158 121       271      0     68 118.0000     1  25
## 159    159 138       287      0     65 115.0000     0  24
## 160    160 136       278      0     61 105.0000     0  23
## 161    161 120       279      0     66 131.0000     0  30
## 162    162 122       278      0     72 155.0000     1  31
## 163    163 134       267      0     66 170.0000     0  30
## 164    164 101       280      0     65 123.0000     1  25
## 165    165 112       288      0     62 125.0000     0  32
## 166    166 132       290      0     64 120.0000     0  25
## 167    167 136       285      0     62 175.0000     0  23
## 168    168 113       277      0     65 192.0000     1  23
## 169    169  96       271      0     64 116.0000     0  23
## 170    170 124       277      0     63 220.0000     0  29
## 171    171 113       306      0     62 150.0000     0  21
## 172    172 137       258      0     63 117.0000     0  25
## 173    173 133       268      0     61  93.0000     0  24
## 174    174 107       244      0     58  97.0000     0  20
## 175    175  96       265      0     59 135.0000     1  28
## 176    176 142       278      0     66 136.0000     1  35
## 177    177 136       275      0     63 110.0000     0  22
## 178    178  75       239      0     63 124.0000     1  26
## 179    179 104       295      0     65 155.0000     1  26
## 180    180 130       274      0     63 150.0000     0  30
## 181    181  90       290      0     63 168.0000     0  22
## 182    182 118       276      0     66 147.0000     1  22
## 183    183 123       320      0     66 117.0000     0  22
## 184    184 137       291      0     61 110.0000     0  34
## 185    185 101       268      0     63 140.0000     0  19
## 186    186 142       275      0     64 132.0000     0  25
## 187    187  98       282      0     63  97.0000     1  20
## 188    188 124       283      0     63 112.0000     0  23
## 189    189 109       281      0     61 105.0000     0  23
## 190    190 150       285      0     61 110.0000     1  22
## 191    191 119       282      0     68 150.0000     1  26
## 192    192 131       280      0     65 125.0000     0  38
## 193    193 101       272      0     63 150.0000     1  29
## 194    194 113       246      0     62 138.0000     1  19
## 195    195 127       270      0     62 150.0000     0  25
## 196    196  97       260      0     61  99.0000     1  23
## 197    197 117       282      0     64 115.0000     0  28
## 198    198 150       290      0     65 125.0000     0  21
## 199    199  85       234      0     67 130.0000     0  33
## 200    200 128       288      0     70 145.0000     0  27
## 201    201 105       233      0     61 130.0000     0  34
## 202    202 115       274      0     65 130.0000     1  22
## 203    203 107       290      0     62 135.0000     0  28
## 204    204 121       275      0     63 121.0000     1  24
## 205    205 119       286      0     64 180.0000     0  20
## 206    206 117       275      0     64 145.0000     1  20
## 207    207 134       264      0     68 136.0000     0  26
## 208    208 117       288      0     65 142.0000     0  35
## 209    209 115       268      0     66 128.0000     0  28
## 210    210 110       254      0     63 120.0000     1  23
## 211    211 130       282      0     62 106.0000     1  21
## 212    212 140       274      0     63 106.0000     1  23
## 213    213  93       249      0     66 117.0000     0  33
## 214    214 154       292      0     65 116.0000     1  42
## 215    215 125       290      0     64 127.0000     0  19
## 216    216  93       318      0     66 135.0000     0  31
## 217    217 122       277      0     63 135.0000     1  33
## 218    218 129       267      0     63 160.0000     0  22
## 219    219 126       276      0     63 120.0000     0  23
## 220    220  85       274      0     68 155.0000     0  24
## 221    221 173       293      0     63 110.0000     0  30
## 222    222 144       329      0     65 190.0000     1  22
## 223    223 114       278      0     65 140.0000     1  25
## 224    224 154       287      0     65 125.0000     1  27
## 225    225 150       274      0     67 117.0000     1  25
## 226    226 111       278      0     62 125.0000     0  21
## 227    227 126       277      0     66 128.0000     0  32
## 228    228 122       261      0     65 124.0000     0  28
## 229    229 141       282      0     68 169.0000     0  24
## 230    230 142       274      0     63 125.0000     0  24
## 231    231  99       262      0     59 110.0000     1  38
## 232    232 113       286      0     63 105.0000     0  23
## 233    233 149       282      0     61 110.0000     0  21
## 234    234 117       328      0     65 125.0000     1  29
## 235    235 128       290      0     64 118.0000     0  22
## 236    236 125       286      0     64 128.4996     0  21
## 237    237 114       290      0     66 160.0000     0  30
## 238    238 130       285      0     63 128.0000     1  23
## 239    239 116       148      0     66 135.0000     0  28
## 240    240  81       256      0     64 148.0000     1  30
## 241    241 124       287      0     62 105.0000     1  27
## 242    242 125       292      0     65 122.0000     0  22
## 243    243 110       262      0     66 140.0000     0  25
## 244    244 125       279      0     63 104.0000     1  23
## 245    245 138       294      0     64 125.0000     0  40
## 246    246 142       284      0     66 132.0000     0  39
## 247    247 115       278      0     60 102.0000     1  23
## 248    248 102       280      0     67 140.0000     0  38
## 249    249 140       294      0     61 103.0000     0  25
## 250    250 133       276      1     63 119.0000     0  22
## 251    251 127       290      0     66 165.0000     0  35
## 252    252 104       274      1     62 115.0000     1  20
## 253    253 119       275      0     67 156.0000     1  42
## 254    254 152       301      0     65 150.0000     0  29
## 255    255 123       284      1     65 120.0000     1  20
## 256    256 143       273      0     66 135.0000     0  19
## 257    257 131       308      0     65 160.0000     0  40
## 258    258 141       319      1     67 140.0000     1  20
## 259    259 129       277      0     66 142.0000     1  30
## 260    260 113       282      1     59 140.0000     0  36
## 261    261 119       292      0     62 118.0000     1  33
## 262    262 109       295      1     63 103.0000     1  23
## 263    263 104       280      1     68 146.0000     1  27
## 264    264 131       282      1     66 126.0000     0  21
## 265    265 110       293      1     64 135.0000     1  28
## 266    266 148       279      0     71 189.0000     0  27
## 267    267 137       283      1     65 157.0000     0  20
## 268    268 117       283      0     63 108.0000     0  27
## 269    269 115       302      1     67 135.0000     0  22
## 270    270  98       280      0     64 122.0000     1  35
## 271    271 136       303      1     68 148.0000     1  20
## 272    272 121       276      1     71 152.0000     1  23
## 273    273 132       285      1     63 140.0000     0  25
## 274    274  91       264      0     60 100.0000     1  36
## 275    275 119       294      0     59 105.0000     0  34
## 276    276  85       273      0     60 105.0000     1  26
## 277    277 106       271      1     61 110.0000     1  26
## 278    278 132       284      0     64 122.0000     0  29
## 279    279  80       266      1     62 125.0000     0  25
## 280    280 109       286      0     64 125.0000     1  24
## 281    281 111       306      0     61 102.0000     0  27
## 282    282 143       292      1     65 125.0000     0  21
## 283    283 136       290      0     66 135.0000     0  26
## 284    284 110       285      1     64 130.0000     0  19
## 285    285  98       257      0     66 130.0000     1  29
## 286    286 108       305      1     65 112.0000     0  24
## 287    287 101       295      0     62 145.0000     1  18
## 288    288  71       281      0     60 117.0000     1  32
## 289    289 124       292      0     68 176.0000     1  29
## 290    290 106       276      0     66 130.0000     0  30
## 291    291 101       278      0     62 112.0000     1  25
## 292    292 100       277      0     62 100.0000     1  31
## 293    293 104       269      0     63 110.0000     1  35
## 294    294 117       270      0     67 135.0000     1  24
## 295    295 117       267      0     65 120.0000     1  29
## 296    296 149       279      0     67 135.0000     0  25
## 297    297 135       284      0     66 123.0000     0  25
## 298    298 110       283      1     66 129.0000     0  21
## 299    299 121       276      0     67 130.0000     0  31
## 300    300 142       285      1     66 136.0000     0  24
## 301    301 104       260      0     64 145.0000     0  33
## 302    302 138       296      0     66 120.0000     0  34
## 303    303 112       278      1     63 120.0000     0  21
## 304    304 117       293      0     60 120.0000     1  39
## 305    305 109       282      0     62 106.0000     1  25
## 306    306 131       266      1     67 135.0000     0  28
## 307    307 120       273      0     64 130.0000     1  29
## 308    308 116       270      0     63 132.0000     0  29
## 309    309 140       290      0     65 110.0000     0  23
## 310    310 103       273      1     64 110.0000     1  22
## 311    311 120       279      1     67 135.0000     0  23
## 312    312 139       260      1     64 127.0000     0  32
## 313    313 123       254      0     62 130.0000     1  26
## 314    314 104       280      1     64 107.0000     1  23
## 315    315 111       270      0     59 103.0000     0  22
## 316    316 122       277      0     63 157.0000     1  32
## 317    317 116       271      1     67 144.0000     1  30
## 318    318 129       277      0     68 130.0000     1  27
## 319    319 133       292      0     65 112.0000     1  30
## 320    320 110       277      0     61 130.0000     0  25
## 321    321 105       276      0     67 130.0000     0  22
## 322    322  93       246      0     65 130.0000     0  37
## 323    323 122       281      0     63 103.0000     1  42
## 324    324 133       293      0     64 110.0000     1  23
## 325    325 130       296      1     66 117.0000     1  22
## 326    326 104       307      0     59 122.0000     0  24
## 327    327 106       278      0     65 110.0000     1  31
## 328    328 120       281      0     63 113.0000     0  33
## 329    329 118       276      1     63 128.0000     0  18
## 330    330 140       290      1     67 132.0000     1  19
## 331    331 114       268      0     64 104.0000     0  22
## 332    332 116       280      0     62 159.0000     0  40
## 333    333 129       284      0     64 115.0000     0  24
## 334    334 120       286      0     62 115.0000     1  22
## 335    335 127       281      0     63 112.0000     1  24
## 336    336  71       234      0     64 110.0000     1  32
## 337    337  88       274      0     66 130.0000     0  30
## 338    338 122       286      0     64 145.0000     0  23
## 339    339 106       302      1     66 147.0000     0  19
## 340    340 135       285      0     66 130.0000     0  30
## 341    341 107       290      0     63 112.0000     0  26
## 342    342 129       294      0     62 170.0000     1  32
## 343    343 126       274      0     62 122.0000     0  39
## 344    344 116       293      1     64 125.0000     0  26
## 345    345 124       294      0     62 122.0000     0  26
## 346    346 123       281      0     68 136.0000     0  23
## 347    347 145       315      0     67 143.0000     1  39
## 348    348 102       278      0     67 135.0000     1  27
## 349    349 129       293      0     65 130.0000     1  30
## 350    350  98       276      1     61 121.0000     0  22
## 351    351 110       272      0     60 108.0000     0  28
## 352    352 135       282      0     67 128.0000     1  24
## 353    353 101       278      1     62 105.0000     0  20
## 354    354  96       266      0     65 125.0000     0  26
## 355    355 104       276      1     60 109.0000     1  18
## 356    356 100       249      0     67 100.0000     0  24
## 357    357 154       292      0     66 145.0000     0  40
## 358    358 127       293      0     67 137.0000     0  31
## 359    359 126       288      0     62 150.0000     0  31
## 360    360 126       282      1     66 115.0000     1  23
## 361    361 127       279      0     67 155.0000     1  26
## 362    362  98       275      0     65 112.0000     1  25
## 363    363 127       288      1     66 130.0000     0  21
## 364    364 129       299      0     68 145.0000     0  22
## 365    365 131       292      1     64 124.0000     1  22
## 366    366 132       289      1     66 145.0000     0  19
## 367    367 127       280      0     62 118.0000     0  27
## 368    368  99       313      1     59 100.0000     1  34
## 369    369 115       290      0     64 140.0000     1  30
## 370    370 145       290      1     67 125.0000     0  24
## 371    371 102       249      1     67 134.0000     1  23
## 372    372 136       299      0     64 115.0000     0  29
## 373    373 121       282      0     66 133.0000     0  22
## 374    374 120       286      0     62 105.0000     0  25
## 375    375 118       261      0     60 104.0000     0  26
## 376    376 127       304      1     62 105.0000     0  26
## 377    377 132       281      1     63 117.0000     0  24
## 378    378 102       258      1     65 135.0000     0  22
## 379    379 143       279      0     65 129.0000     1  39
## 380    380 118       277      0     62 120.0000     0  25
## 381    381 102       286      1     64 140.0000     0  22
## 382    382 163       280      0     69 139.0000     0  35
## 383    383 132       294      0     64 116.0000     0  32
## 384    384 116       276      0     61 180.0000     0  33
## 385    385 138       288      1     66 124.0000     0  19
## 386    386 139       279      0     64 143.0000     0  20
## 387    387 132       298      1     61 137.0000     0  23
## 388    388  87       282      0     63 104.0000     1  27
## 389    389 131       297      0     67 132.0000     0  30
## 390    390 130       282      0     67 147.0000     1  26
## 391    391 123       290      0     66 107.0000     1  28
## 392    392 115       276      1     63 110.0000     0  18
## 393    393 116       272      0     64 130.0000     1  27
## 394    394 119       286      1     67 130.0000     0  20
## 395    395 125       279      1     67 135.0000     0  19
## 396    396 144       282      0     66 155.0000     1  33
## 397    397 123       269      0     67 132.0000     0  26
## 398    398 120       276      0     66 114.0000     0  23
## 399    399 140       251      0     63 210.0000     0  28
## 400    400 120       271      1     64 142.0000     1  17
## 401    401 120       289      1     59 102.0000     0  31
## 402    402 146       280      0     61 145.0000     0  23
## 403    403 112       283      1     62 102.0000     1  21
## 404    404 132       278      0     64 150.0000     1  20
## 405    405 146       263      0     53 110.0000     1  39
## 406    406 122       275      0     68 140.0000     0  30
## 407    407 128       292      0     66 130.0000     0  32
## 408    408 119       277      0     63 120.0000     1  24
## 409    409 135       278      0     66 148.0000     0  27
## 410    410 129       235      0     66 135.0000     0  24
## 411    411 116       293      1     62 108.0000     0  28
## 412    412 100       275      0     64 111.0000     1  27
## 413    413 138       257      0     67 138.0000     0  38
## 414    414 123       282      0     65 130.0000     0  22
## 415    415 113       288      1     61 120.0000     0  21
## 416    416 129       280      1     65 140.0000     1  24
## 417    417 122       280      0     67 127.0000     1  24
## 418    418 132       281      1     67 140.0000     0  21
## 419    419 120       269      1     63 130.0000     0  40
## 420    420 114       283      1     65 115.0000     0  20
## 421    421 130       280      0     66 135.0000     0  29
## 422    422 117       286      0     66 127.0000     1  32
## 423    423 142       285      0     63 124.0000     0  33
## 424    424 144       273      0     62 118.0000     1  27
## 425    425 127       262      1     64 125.0000     0  32
## 426    426 115       270      0     67 165.0000     1  25
## 427    427  85       258      0     67 137.0000     0  41
## 428    428  99       274      0     66 118.0000     1  28
## 429    429 123       323      1     64 140.0000     0  17
## 430    430 112       281      1     61 150.0000     0  23
## 431    431  68       223      0     66 149.0000     1  32
## 432    432 102       283      1     65 127.0000     1  19
## 433    433 109       273      0     65 138.0000     1  37
## 434    434 102       267      1     60  93.0000     1  25
## 435    435  99       275      0     61 125.0000     1  23
## 436    436  78       256      1     65 123.0000     0  29
## 437    437 128       284      1     66 111.0000     1  19
## 438    438 107       303      1     67 133.0000     0  25
## 439    439 136       295      0     64 147.0000     0  23
## 440    440 101       278      0     61  99.0000     1  27
## 441    441 100       275      1     64 125.0000     0  25
## 442    442 109       272      0     66 154.0000     1  41
## 443    443 117       281      1     70 141.0000     1  21
## 444    444  88       252      1     60 115.0000     1  21
## 445    445  95       270      0     65 135.0000     1  35
## 446    446 127       291      1     66 135.0000     1  24
## 447    447 107       293      0     65 155.0000     1  20
## 448    448 126       262      0     66 135.0000     1  37
## 449    449  98       278      0     63 110.0000     1  27
## 450    450  96       241      0     64 130.0000     1  23
## 451    451 104       282      0     63 144.0000     0  24
## 452    452 133       273      1     63 135.0000     0  33
## 453    453  93       267      0     63 135.0000     1  25
## 454    454 101       280      1     65 123.0000     1  24
## 455    455 118       277      0     64 155.0000     0  21
## 456    456 130       289      0     61 130.0000     1  21
## 457    457 125       288      0     63 128.0000     1  22
## 458    458 140       291      1     65 122.0000     0  19
## 459    459 115       290      1     65 118.0000     0  19
## 460    460 130       293      0     63 123.0000     0  26
## 461    461 114       277      1     64 125.0000     0  31
## 462    462 105       278      0     64 120.0000     0  21
## 463    463 101       289      1     60 125.0000     0  31
## 464    464 132       286      0     67 122.0000     1  26
## 465    465 112       252      0     64 162.0000     0  37
## 466    466  69       232      0     59 128.4996     1  31
## 467    467 114       264      0     63 110.0000     1  26
## 468    468 123       267      0     63 111.0000     1  29
## 469    469 129       284      1     66 130.0000     1  20
## 470    470 114       283      1     64 117.0000     1  15
## 471    471 115       290      0     62  95.0000     0  31
## 472    472  98       272      1     64 129.0000     0  35
## 473    473 128       283      0     67 126.0000     0  27
## 474    474 119       271      0     64 175.0000     1  28
## 475    475 154       288      0     65 147.0000     0  25
## 476    476 127       247      1     63 140.0000     0  21
## 477    477 131       263      0     64 180.0000     1  29
## 478    478 129       288      0     59 102.0000     0  28
## 479    479 114       286      1     64 116.0000     1  22
## 480    480 110       280      0     62 110.0000     1  29
## 481    481 103       268      0     64 150.0000     1  31
## 482    482 117       287      0     65 115.0000     1  20
## 483    483 138       282      0     64 120.0000     0  25
## 484    484 126       280      0     66 147.0000     1  24
## 485    485 124       271      0     66 145.0000     0  23
## 486    486 111       284      0     62 110.0000     0  34
## 487    487 132       282      0     67 200.0000     1  28
## 488    488 103       240      0     65 140.0000     0  26
## 489    489 158       285      0     62 130.0000     0  28
## 490    490 101       286      1     64 117.0000     1  21
## 491    491 132       290      0     66 125.0000     0  26
## 492    492 114       293      1     66 180.0000     1  20
## 493    493  71       277      0     69 135.0000     0  40
## 494    494 116       282      0     64 120.0000     0  19
## 495    495 108       271      0     60 109.0000     1  19
## 496    496 123       298      0     64 113.0000     1  25
## 497    497 129       289      0     63 132.0000     0  37
## 498    498 134       282      0     62 110.0000     0  24
## 499    499 113       298      0     60 124.0000     1  30
## 500    500 123       277      1     65 160.0000     0  20
## 501    501 147       277      0     68 160.0000     0  30
## 502    502 121       270      1     62 103.0000     0  20
## 503    503 125       284      1     67 130.0000     0  19
## 504    504 115       277      1     66 128.0000     0  25
## 505    505 101       289      0     59  96.0000     0  27
## 506    506  93       271      0     65 127.0000     1  30
## 507    507 109       275      0     66 120.0000     0  33
## 508    508 115       276      1     60 106.0000     0  23
## 509    509 130       293      1     65 122.0000     1  23
## 510    510 123       278      0     61  89.0000     0  21
## 511    511 111       300      0     64 108.0000     1  20
## 512    512  97       279      1     64 138.0000     1  24
## 513    513 122       292      1     65 125.0000     0  25
## 514    514 124       300      0     63  95.0000     0  28
## 515    515 129       276      0     66 145.0000     0  26
## 516    516 124       290      0     59 140.0000     0  26
## 517    517 107       280      0     60 107.0000     1  20
## 518    518 142       273      1     62 125.0000     0  22
## 519    519 129       287      1     66 135.0000     0  29
## 520    520 174       281      0     67 155.0000     0  37
## 521    521 105       264      0     65 105.0000     1  30
## 522    522 103       291      1     63 102.0000     0  26
## 523    523 124       285      1     63 114.0000     0  27
## 524    524 105       265      0     65 124.0000     0  43
## 525    525 133       275      0     65 137.0000     1  36
## 526    526 161       302      1     70 170.0000     1  22
## 527    527 105       260      0     64 197.0000     0  23
## 528    528 108       281      0     66 171.0000     0  41
## 529    529 153       297      0     66 145.0000     0  27
## 530    530 133       280      1     61 130.0000     0  25
## 531    531 115       269      0     63 165.0000     1  41
## 532    532 127       254      0     67 146.0000     1  27
## 533    533 128       271      0     65 135.0000     1  41
## 534    534 117       265      0     68 134.0000     1  40
## 535    535 123       274      0     66 135.0000     0  23
## 536    536 119       288      1     64 132.0000     1  22
## 537    537 141       284      1     64 105.0000     0  17
## 538    538  91       260      0     62 110.0000     1  26
## 539    539 116       291      0     65 133.0000     1  29
## 540    540 116       255      0     65 132.0000     0  24
## 541    541 121       273      0     64 112.0000     0  32
## 542    542 111       274      0     67 159.0000     0  36
## 543    543 102       257      0     66 135.0000     0  25
## 544    544 118       283      0     65 150.0000     0  24
## 545    545 126       294      1     65 125.0000     1  22
## 546    546  98       286      0     62 105.0000     1  31
## 547    547 131       288      1     65 125.0000     0  28
## 548    548 115       278      0     60 113.0000     0  21
## 549    549 103       281      1     59 128.4996     1  22
## 550    550 147       301      0     65 130.0000     0  26
## 551    551 123       308      1     65 135.0000     0  19
## 552    552 125       283      0     65 119.0000     0  22
## 553    553 117       270      0     67 130.0000     1  30
## 554    554  99       268      0     71 150.0000     0  29
## 555    555 115       283      0     66 127.0000     1  31
## 556    556 116       265      0     63 120.0000     0  36
## 557    557 118       297      0     68 140.0000     1  35
## 558    558 170       303      1     64 129.0000     0  21
## 559    559 104       270      0     61 110.0000     0  25
## 560    560 108       269      1     62 114.0000     0  20
## 561    561 144       289      1     69 130.0000     1  17
## 562    562  99       250      1     66 115.0000     0  26
## 563    563  97       263      1     63 107.0000     0  25
## 564    564  85       270      1     63 118.0000     1  19
## 565    565 130       285      1     66 126.0000     1  24
## 566    566 117       275      0     62 115.0000     1  22
## 567    567 109       302      0     64 110.0000     0  24
## 568    568 147       285      0     64 137.0000     0  24
## 569    569 105       281      1     64 115.0000     0  23
## 570    570 135       278      1     68 139.0000     1  27
## 571    571 115       273      1     67 215.0000     1  23
## 572    572 123       280      0     65 140.0000     1  23
## 573    573 105       274      1     61 100.0000     0  26
## 574    574 154       271      0     69 160.0000     1  36
## 575    575 110       276      0     63 107.0000     1  25
## 576    576 119       285      1     62 108.0000     0  26
## 577    577 103       292      1     62 132.0000     0  28
## 578    578 117       272      0     64 116.0000     0  25
## 579    579 120       289      0     69 165.0000     0  23
## 580    580 145       278      0     62 109.0000     0  24
## 581    581 104       271      0     62  98.0000     1  20
## 582    582 123       268      1     62 110.0000     1  18
## 583    583 124       272      0     62 110.0000     0  27
## 584    584 129       275      0     64 115.0000     1  26
## 585    585  91       248      0     63 202.0000     0  33
## 586    586 109       295      0     61 135.0000     0  32
## 587    587 108       268      0     58 112.0000     1  22
## 588    588  79       268      0     61 108.0000     0  36
## 589    589 133       301      0     62 108.0000     0  23
## 590    590 114       309      1     62 118.0000     0  27
## 591    591 128       273      0     61 125.0000     0  34
## 592    592 129       280      1     65 126.0000     0  24
## 593    593  97       234      1     65 112.0000     0  26
## 594    594 103       276      1     62 130.0000     1  21
## 595    595 176       293      1     68 180.0000     0  19
## 596    596 143       294      0     65 145.0000     0  44
## 597    597 127       292      1     68 130.0000     1  21
## 598    598 107       256      0     59  90.0000     1  28
## 599    599 113       268      0     62 100.0000     0  31
## 600    600 106       279      1     62 118.0000     1  21
## 601    601 152       285      0     61 120.0000     1  24
## 602    602 150       275      0     65 145.0000     0  29
## 603    603 136       278      0     64 118.0000     1  35
## 604    604 124       279      0     66 129.0000     0  35
## 605    605 123       284      1     64 112.0000     1  18
## 606    606 119       288      0     62 128.0000     0  37
## 607    607 122       291      0     64 155.0000     0  40
## 608    608 112       250      0     67 124.0000     0  34
## 609    609  93       270      0     64 125.0000     1  25
## 610    610 136       274      1     63 165.0000     0  20
## 611    611 150       292      0     64 124.0000     0  26
## 612    612  94       264      1     64 135.0000     0  26
## 613    613 146       306      0     63 112.0000     0  38
## 614    614 129       274      0     65 101.0000     1  19
## 615    615 125       292      0     65 117.0000     1  27
## 616    616 124       273      0     63 115.0000     0  21
## 617    617 141       282      0     63 115.0000     0  27
## 618    618  96       266      0     67 135.0000     1  33
## 619    619 138       297      0     66 133.0000     1  30
## 620    620 127       282      0     67 134.0000     0  28
## 621    621 114       251      0     64 119.0000     1  26
## 622    622 103       297      0     64 125.0000     0  31
## 623    623 127       288      1     65 115.0000     1  20
## 624    624 113       274      0     63 108.0000     1  23
## 625    625  99       249      1     57  98.0000     1  31
## 626    626  97       279      0     61 105.0000     1  33
## 627    627 116       275      1     68 145.0000     0  20
## 628    628 126       297      0     66 120.0000     1  26
## 629    629 119       277      0     66 130.0000     1  28
## 630    630 123       283      0     62 110.0000     0  27
## 631    631 129       287      0     60 107.0000     0  24
## 632    632 117       256      0     65 132.0000     1  37
## 633    633 100       275      0     60 115.0000     0  26
## 634    634 131       274      0     64 118.0000     1  28
## 635    635 146       279      0     64 124.0000     0  27
## 636    636  84       267      0     60  95.0000     0  29
## 637    637 115       302      0     64 116.0000     0  28
## 638    638 115       281      0     60  94.0000     0  25
## 639    639 118       284      0     70 145.0000     1  28
## 640    640  91       292      1     61 125.0000     0  19
## 641    641 112       255      0     60 115.0000     0  39
## 642    642 115       316      1     64 110.0000     0  29
## 643    643 110       269      0     61 102.0000     0  38
## 644    644 117       277      0     66 140.0000     0  34
## 645    645 109       268      1     65 120.0000     1  29
## 646    646  99       267      0     62  94.0000     0  22
## 647    647 131       274      0     62 160.0000     1  27
## 648    648 136       291      1     61 105.0000     0  25
## 649    649 130       298      0     62 120.0000     0  20
## 650    650 134       296      0     60 117.0000     1  35
## 651    651 128       271      0     65 126.0000     1  29
## 652    652 150       286      0     67 175.0000     0  38
## 653    653  86       284      0     65 174.0000     1  39
## 654    654 115       278      0     63 112.0000     1  26
## 655    655  78       237      1     63 144.0000     0  23
## 656    656 100       295      1     68 125.0000     1  21
## 657    657 116       270      0     68 169.0000     0  25
## 658    658 110       271      1     66 135.0000     0  26
## 659    659 109       283      0     64 120.0000     0  34
## 660    660 113       259      0     64 128.0000     0  38
## 661    661 136       297      1     66 135.0000     0  23
## 662    662 121       273      1     61 125.0000     0  34
## 663    663 117       288      1     63 140.0000     0  28
## 664    664 166       299      0     68 140.0000     0  26
## 665    665  87       229      0     62 138.0000     0  27
## 666    666 120       294      1     66 128.0000     1  23
## 667    667  95       286      0     66 118.0000     1  26
## 668    668 132       273      0     62 113.0000     0  28
## 669    669  90       286      0     63 105.0000     1  32
## 670    670 131       308      0     58 150.0000     1  30
## 671    671 103       279      1     65 145.0000     1  22
## 672    672 144       287      1     71 153.0000     1  33
## 673    673 137       299      0     62 115.0000     0  24
## 674    674 124       270      0     64 122.0000     0  20
## 675    675 136       281      1     64 127.0000     0  27
## 676    676 117       298      1     64 160.0000     0  22
## 677    677 121       269      0     62 130.0000     0  23
## 678    678 116       280      0     68 198.0000     0  34
## 679    679 139       275      0     62 118.0000     0  33
## 680    680 110       280      0     67 125.0000     0  39
## 681    681  86       242      0     64 110.0000     1  20
## 682    682 133       287      0     65 165.0000     0  20
## 683    683  81       254      0     62 157.0000     0  23
## 684    684 133       281      0     63 120.0000     0  33
## 685    685 132       284      1     66 140.0000     0  20
## 686    686 132       287      0     64 148.0000     0  29
## 687    687 137       274      0     64 126.0000     0  27
## 688    688  84       279      0     63 190.0000     0  34
## 689    689 136       279      0     69 130.0000     1  30
## 690    690  92       270      0     62 100.0000     1  34
## 691    691 114       298      1     67 114.0000     0  28
## 692    692 129       274      0     69 136.0000     1  33
## 693    693 167       288      1     63 117.0000     0  19
## 694    694 124       282      1     65 118.0000     0  22
## 695    695 105       269      0     62 100.0000     1  27
## 696    696 155       283      1     70 137.0000     0  19
## 697    697 125       279      1     66 126.0000     0  21
## 698    698 125       266      0     62 120.0000     1  21
## 699    699 125       283      1     59  96.0000     0  22
## 700    700 115       315      1     62 110.0000     0  22
## 701    701 174       288      0     61 182.0000     0  25
## 702    702 127       290      0     66 122.0000     0  35
## 703    703 113       262      0     60 105.0000     0  24
## 704    704 115       273      0     66 130.0000     1  22
## 705    705 139       277      0     63 140.0000     0  35
## 706    706 127       275      0     62 125.0000     0  26
## 707    707 112       272      0     60  98.0000     0  26
## 708    708 143       285      0     64 135.0000     1  30
## 709    709 116       286      1     58 105.0000     1  22
## 710    710 155       279      0     61 125.0000     0  33
## 711    711 121       290      0     64 127.0000     0  31
## 712    712 110       282      1     66 125.0000     1  21
## 713    713  87       277      0     62 120.0000     1  31
## 714    714 132       330      0     64 130.0000     1  34
## 715    715 129       277      0     68 128.4996     0  24
## 716    716 123       280      0     62 105.0000     1  20
## 717    717  91       279      1     62 118.0000     0  27
## 718    718 147       286      0     68 147.0000     0  30
## 719    719 144       289      0     62 106.0000     0  20
## 720    720 128       292      0     64 127.0000     0  30
## 721    721 137       318      1     64 110.0000     0  19
## 722    722 104       289      0     60 104.0000     1  24
## 723    723 120       271      0     63 130.0000     0  32
## 724    724 112       277      1     64 118.0000     0  23
## 725    725 138       286      0     63 111.0000     0  26
## 726    726  96       280      0     63 105.0000     1  27
## 727    727 134       285      0     62 134.0000     0  35
## 728    728 126       285      0     64 140.0000     0  24
## 729    729 112       300      0     62 121.0000     0  29
## 730    730 138       313      1     65 111.0000     0  27
## 731    731 110       275      0     63 120.0000     0  25
## 732    732  83       253      0     63 110.0000     1  29
## 733    733 112       288      1     62 110.0000     0  20
## 734    734 148       286      0     68 160.0000     0  38
## 735    735 119       300      1     63 124.0000     0  34
## 736    736  86       246      0     64 113.0000     1  25
## 737    737 110       269      0     63 145.0000     1  38
## 738    738 126       282      0     61 120.0000     0  23
## 739    739 125       272      0     60  96.0000     0  30
## 740    740 136       252      0     63 130.0000     0  27
## 741    741 127       283      1     64 119.0000     0  29
## 742    742  84       272      0     64 150.0000     1  25
## 743    743 131       278      0     66 124.0000     0  22
## 744    744 123       286      1     67 130.0000     1  21
## 745    745  96       282      1     68 127.0000     1  30
## 746    746 110       286      0     62 100.0000     0  26
## 747    747 123       282      0     68 164.0000     0  29
## 748    748 152       286      1     67 135.0000     0  19
## 749    749 127       288      0     65 155.0000     0  28
## 750    750 117       269      1     64 149.0000     1  21
## 751    751 125       277      0     66 139.0000     1  29
## 752    752 139       273      0     68 130.0000     0  29
## 753    753 114       280      0     66 134.0000     1  31
## 754    754  96       280      1     62 127.0000     1  34
## 755    755 124       289      0     63 110.0000     0  29
## 756    756 107       272      0     64 140.0000     1  30
## 757    757 113       277      0     64 108.0000     0  38
## 758    758  98       292      1     65 124.0000     1  20
## 759    759 119       285      1     65 127.0000     0  28
## 760    760 107       268      0     58 112.0000     1  37
## 761    761 117       255      0     61 120.0000     0  26
## 762    762 117       305      0     64 155.0000     0  24
## 763    763 144       276      0     67 129.0000     1  23
## 764    764 136       268      0     63 132.0000     1  30
## 765    765 121       278      0     69 132.0000     0  28
## 766    766 165       282      0     66 145.0000     0  29
## 767    767 120       279      0     64 124.0000     0  38
## 768    768 125       280      0     65 130.0000     1  30
## 769    769 137       285      0     65 110.0000     0  29
## 770    770 100       288      1     61 108.0000     1  28
## 771    771 134       284      0     62 112.0000     0  28
## 772    772  88       262      0     65 118.0000     1  20
## 773    773 108       291      0     65 135.0000     0  39
## 774    774 123       271      0     64 162.0000     0  41
## 775    775 141       277      0     66 162.0000     0  38
## 776    776 130       270      1     66 130.0000     0  19
## 777    777 139       299      1     67 112.0000     0  20
## 778    778 130       283      0     65 118.0000     0  32
## 779    779 113       289      1     59  91.0000     0  26
## 780    780  77       238      1     63 103.0000     1  23
## 781    781  62       228      0     61 107.0000     0  24
## 782    782  93       245      0     61 100.0000     1  33
## 783    783 109       275      1     63 112.0000     1  37
## 784    784 145       283      0     65 125.0000     1  27
## 785    785  92       224      0     63 134.0000     1  19
## 786    786 120       281      0     61 115.0000     0  26
## 787    787 135       284      0     67 141.0000     0  39
## 788    788 113       287      0     63 118.0000     0  36
## 789    789 126       251      1     64 123.0000     0  28
## 790    790 143       270      1     70 148.0000     0  27
## 791    791 128       282      1     64 125.0000     0  25
## 792    792  98       262      0     67 120.0000     0  22
## 793    793 110       306      1     61 122.0000     0  32
## 794    794 162       284      0     64 126.0000     0  27
## 795    795 116       292      1     65 118.0000     0  20
## 796    796 128       284      0     62 110.0000     0  23
## 797    797 111       275      1     61 108.0000     1  18
## 798    798 137       280      0     60 107.0000     0  34
## 799    799 100       264      0     64 120.0000     1  29
## 800    800 160       271      0     67 215.0000     0  32
## 801    801 112       267      1     62 138.0000     0  22
## 802    802 134       297      0     67 170.0000     1  27
## 803    803 145       308      0     64 110.0000     1  35
## 804    804 116       295      0     65 120.0000     0  32
## 805    805 126       278      0     64 150.0000     1  26
## 806    806 111       285      0     65 130.0000     0  29
## 807    807 126       282      0     62 117.0000     0  33
## 808    808 109       291      0     64 107.0000     0  39
## 809    809 136       291      0     66 191.0000     0  41
## 810    810 119       286      0     63 185.0000     1  22
## 811    811 103       267      1     66 150.0000     1  21
## 812    812 124       284      1     62 112.0000     0  17
## 813    813 155       286      0     66 127.0000     0  31
## 814    814 122       282      1     66 110.0000     0  21
## 815    815 113       285      0     66 140.0000     0  26
## 816    816 122       273      0     66 210.0000     0  26
## 817    817 126       293      1     62 111.0000     0  27
## 818    818 116       277      0     64 124.0000     1  41
## 819    819 102       294      0     65 130.0000     1  21
## 820    820 110       181      0     64 133.0000     0  27
## 821    821 133       285      1     64 160.0000     0  30
## 822    822 125       283      0     65 125.0000     0  29
## 823    823 164       286      1     66 143.0000     0  32
## 824    824 133       297      0     61 125.0000     0  36
## 825    825 124       293      1     65 150.0000     0  19
## 826    826 122       306      1     62 100.0000     0  22
## 827    827 121       271      1     63 129.0000     1  34
## 828    828 100       272      0     64 150.0000     1  30
## 829    829  90       266      1     67 135.0000     0  26
## 830    830 128       272      1     67 109.0000     0  18
## 831    831  86       276      1     65 125.0000     1  23
## 832    832 123       282      0     63 118.0000     0  30
## 833    833  87       275      0     63 110.0000     1  28
## 834    834 128       291      1     63 132.0000     0  27
## 835    835 120       288      0     63 125.0000     0  28
## 836    836 125       301      1     68 181.0000     0  35
## 837    837 118       265      0     61 123.0000     0  27
## 838    838 116       284      1     66 117.0000     0  24
## 839    839 131       262      0     67 135.0000     0  22
## 840    840 151       286      1     66 130.0000     0  22
## 841    841  88       273      0     66 110.0000     1  20
## 842    842 137       284      0     67 110.0000     0  30
## 843    843 127       289      0     67 140.0000     0  23
## 844    844  96       278      1     60 120.0000     1  18
## 845    845 129       281      0     67 155.0000     0  31
## 846    846 128       288      1     65 114.0000     0  26
## 847    847  85       255      0     68 159.0000     0  24
## 848    848 111       281      1     64 112.0000     0  27
## 849    849 124       275      0     61 116.0000     0  28
## 850    850 112       292      1     62 110.0000     1  28
## 851    851 115       281      0     61 128.0000     1  28
## 852    852  72       271      0     61 136.0000     0  39
## 853    853 122       281      1     65 137.0000     1  24
## 854    854 116       291      0     66 153.0000     0  26
## 855    855 127       272      0     64 130.0000     1  20
## 856    856  90       266      0     61  99.0000     1  23
## 857    857  99       273      1     59 115.0000     0  27
## 858    858 144       307      1     66 125.0000     0  26
## 859    859 138       280      1     65 175.0000     0  30
## 860    860  58       245      0     64 156.0000     1  34
## 861    861 109       265      1     63 107.0000     1  24
## 862    862 110       277      1     62 160.0000     0  19
## 863    863 129       278      0     63 128.0000     0  27
## 864    864 150       284      0     67 130.0000     0  40
## 865    865 128       279      0     66 135.0000     0  27
## 866    866 142       284      1     66 137.0000     1  31
## 867    867 115       268      1     64 125.0000     0  31
## 868    868 108       283      0     62 108.0000     0  35
## 869    869 139       281      0     63 137.0000     0  27
## 870    870 115       275      0     61 155.0000     1  25
## 871    871 136       288      0     62 128.4996     0  23
## 872    872 163       289      1     64 126.0000     1  25
## 873    873 131       285      0     64 130.0000     0  26
## 874    874  77       238      0     67 135.0000     1  38
## 875    875 124       283      1     67 156.0000     1  33
## 876    876 104       270      1     62 115.0000     0  26
## 877    877 102       267      1     61 109.0000     1  24
## 878    878  94       268      0     62 105.0000     1  30
## 879    879 158       295      1     70 137.0000     0  37
## 880    880 112       275      1     68 143.0000     1  21
## 881    881 119       286      0     64 123.0000     1  26
## 882    882  97       279      0     68 178.0000     1  29
## 883    883  99       252      0     64 120.0000     0  21
## 884    884 115       264      1     67 134.0000     1  23
## 885    885 139       284      0     61 121.0000     0  37
## 886    886 144       304      1     58 102.0000     1  27
## 887    887  99       270      0     63 115.0000     1  22
## 888    888 105       280      1     63 116.0000     0  22
## 889    889  89       275      0     66 170.0000     0  34
## 890    890 129       270      0     67 160.0000     0  43
## 891    891 119       270      1     64 109.0000     0  20
## 892    892 114       291      0     60 112.0000     0  35
## 893    893 106       289      0     67 120.0000     1  28
## 894    894 122       292      1     65 133.0000     0  34
## 895    895 136       261      0     65 110.0000     0  24
## 896    896 121       286      1     69 130.0000     1  22
## 897    897 112       282      0     65 122.0000     0  26
## 898    898 112       266      0     64 122.0000     0  26
## 899    899 123       314      0     61 121.0000     1  22
## 900    900 139       286      0     65 125.0000     1  33
## 901    901 125       290      0     59 105.0000     0  36
## 902    902 105       295      1     64 112.0000     1  20
## 903    903 130       276      0     68 130.0000     0  41
## 904    904 146       294      0     66 145.0000     1  22
## 905    905 133       290      0     64 145.0000     0  21
## 906    906 147       296      1     67 124.0000     0  19
## 907    907 109       269      0     63 113.0000     0  23
## 908    908 122       286      0     64 120.0000     1  23
## 909    909 135       260      0     65 135.0000     0  43
## 910    910 117       272      0     66 118.0000     0  32
## 911    911 138       284      0     66 133.0000     1  30
## 912    912 120       283      0     64 122.0000     1  28
## 913    913 119       273      0     65 125.0000     1  35
## 914    914 118       278      1     62 126.0000     0  19
## 915    915 105       330      0     64 112.0000     1  23
## 916    916 113       306      1     65 137.0000     0  21
## 917    917 148       291      1     63 115.0000     0  21
## 918    918 140       281      1     69 135.0000     0  22
## 919    919 134       287      1     67 131.0000     0  33
## 920    920 120       280      0     61 111.0000     0  31
## 921    921 123       296      1     64 110.0000     1  26
## 922    922 102       275      0     64 160.0000     0  43
## 923    923  55       204      0     65 140.0000     0  35
## 924    924 103       276      1     63 149.0000     1  19
## 925    925 123       283      0     65 110.0000     0  21
## 926    926 105       270      1     65 134.0000     1  27
## 927    927 138       289      0     65 155.0000     0  33
## 928    928 128       281      0     63 150.0000     0  28
## 929    929 139       285      0     65 129.0000     1  30
## 930    930 104       288      1     61 122.0000     1  27
## 931    931 159       296      1     64 112.0000     0  27
## 932    932 118       276      0     62 130.0000     1  29
## 933    933  99       285      0     69 128.0000     1  25
## 934    934 144       281      0     63 120.0000     0  20
## 935    935 121       270      0     62 108.0000     1  25
## 936    936 117       265      1     66  98.0000     0  24
## 937    937 119       293      1     65 127.0000     0  23
## 938    938 105       281      1     61 130.0000     0  19
## 939    939 125       283      0     63 145.0000     1  37
## 940    940 119       259      0     62 130.0000     0  37
## 941    941 101       273      0     60 113.0000     0  39
## 942    942 105       277      1     64 156.0000     0  25
## 943    943 110       281      0     60 110.0000     0  27
## 944    944 100       270      1     65 132.0000     1  21
## 945    945  98       284      0     68 140.0000     0  29
## 946    946 127       276      0     64 159.0000     0  37
## 947    947 117       324      0     62 164.0000     1  22
## 948    948 122       278      0     68 114.0000     0  37
## 949    949 122       273      1     64 130.0000     1  23
## 950    950 118       281      1     66 140.0000     1  36
## 951    951 137       303      1     66 127.0000     1  23
## 952    952 120       275      0     63 115.0000     1  32
## 953    953 143       285      0     68 185.0000     0  27
## 954    954 108       270      0     67 124.0000     1  29
## 955    955 131       284      1     61 114.0000     1  19
## 956    956 110       277      0     61 116.0000     0  36
## 957    957 105       276      0     62 112.0000     1  20
## 958    958 125       255      0     63 133.0000     0  23
## 959    959  78       258      1     66 115.0000     1  24
## 960    960 114       289      0     60 115.0000     0  36
## 961    961 111       278      0     65 145.0000     1  29
## 962    962 103       250      0     59 140.0000     0  40
## 963    963 114       276      0     62 127.0000     0  26
## 964    964  75       247      0     64 120.0000     1  36
## 965    965 169       296      0     67 185.0000     0  33
## 966    966  94       271      0     61 130.0000     1  36
## 967    967 150       287      0     62 135.0000     0  36
## 968    968 144       248      0     70 145.0000     0  30
## 969    969 144       291      0     67 130.0000     0  28
## 970    970 143       313      0     68 150.0000     0  20
## 971    971 145       304      1     63 109.0000     1  25
## 972    972 121       285      0     64 110.0000     0  34
## 973    973 105       256      0     66 142.0000     0  31
## 974    974 134       286      0     64 125.0000     0  25
## 975    975 129       294      1     65 132.0000     0  21
## 976    976 114       276      0     63 110.0000     0  24
## 977    977  97       265      0     61 110.0000     0  30
## 978    978 160       292      0     64 120.0000     0  28
## 979    979  65       237      0     67 130.0000     0  31
## 980    980 145       288      0     64 116.0000     0  28
## 981    981  95       273      0     60  90.0000     0  23
## 982    982 139       293      1     69 130.0000     0  21
## 983    983 123       288      0     63 125.0000     0  27
## 984    984 109       283      0     65 112.0000     1  23
## 985    985 110       268      0     64 127.0000     0  34
## 986    986 122       296      1     65 132.0000     0  24
## 987    987 115       307      0     65 128.0000     1  34
## 988    988 108       279      1     64 115.0000     0  19
## 989    989 120       287      0     67 116.0000     1  23
## 990    990 131       269      0     68 145.0000     0  36
## 991    991 136       283      1     63 119.0000     0  24
## 992    992 125       290      0     63 135.0000     0  32
## 993    993  96       285      1     66 117.0000     1  20
## 994    994 102       282      1     65 125.0000     1  29
## 995    995 102       288      1     65 117.0000     0  18
## 996    996 112       277      1     67 120.0000     0  22
## 997    997 135       272      0     65 130.0000     0  30
## 998    998  91       266      0     60 120.0000     1  23
## 999    999 129       276      0     63 125.0000     0  31
## 1000  1000 155       290      0     66 129.0000     1  26
## 1001  1001 109       274      0     69 144.0000     1  33
## 1002  1002  80       262      1     61 100.0000     1  31
## 1003  1003 125       273      0     64 145.0000     0  30
## 1004  1004  94       284      0     63 104.0000     1  24
## 1005  1005 148       281      0     63 110.0000     1  27
## 1006  1006  73       277      0     65 145.0000     0  29
## 1007  1007 123       267      1     66 132.0000     1  19
## 1008  1008  65       232      0     66 125.0000     1  24
## 1009  1009 118       279      1     64 108.0000     0  21
## 1010  1010 102       283      0     60 119.0000     0  39
## 1011  1011 120       280      0     61 118.0000     0  24
## 1012  1012 108       270      1     65 130.0000     1  21
## 1013  1013 122       280      1     62 128.0000     0  45
## 1014  1014 103       268      0     62  97.0000     1  32
## 1015  1015 105       312      0     61 115.0000     1  41
## 1016  1016 126       273      1     68 135.0000     0  25
## 1017  1017 145       316      0     67 142.0000     0  22
## 1018  1018 139       293      0     66 131.0000     0  34
## 1019  1019 124       290      0     65 165.0000     0  26
## 1020  1020 121       282      0     65 122.0000     0  30
## 1021  1021 126       299      1     60 114.0000     0  21
## 1022  1022 119       286      1     67 137.0000     0  33
## 1023  1023 114       277      1     63 107.0000     0  19
## 1024  1024 118       272      0     64 113.0000     0  23
## 1025  1025 127       295      0     65 145.0000     0  36
## 1026  1026 117       290      1     67 110.0000     0  22
## 1027  1027 137       277      0     65 126.0000     0  41
## 1028  1028 133       292      0     65 135.0000     0  29
## 1029  1029 100       264      0     60 111.0000     1  28
## 1030  1030 107       273      1     65 135.0000     0  26
## 1031  1031 115       276      1     62 105.0000     1  20
## 1032  1032  91       292      1     61 113.0000     1  26
## 1033  1033 112       287      0     64 110.0000     1  27
## 1034  1034 125       289      1     61 120.0000     0  31
## 1035  1035 157       291      0     65 121.0000     0  33
## 1036  1036 108       256      1     67 130.0000     0  26
## 1037  1037 130       279      0     62 122.0000     0  31
## 1038  1038 135       289      0     64 127.0000     0  25
## 1039  1039 123       277      0     66 122.0000     0  24
## 1040  1040 100       281      0     61 115.0000     0  24
## 1041  1041 124       277      1     64 104.0000     0  23
## 1042  1042 174       284      0     65 163.0000     0  39
## 1043  1043 129       278      0     67 146.0000     0  26
## 1044  1044 119       275      0     59 113.0000     1  27
## 1045  1045 126       272      1     61 120.0000     1  35
## 1046  1046 128       267      0     61 142.0000     0  37
## 1047  1047 116       282      1     64 124.0000     0  19
## 1048  1048 100       285      0     68 127.0000     1  18
## 1049  1049  96       285      0     66 135.0000     1  37
## 1050  1050 131       279      1     68 122.0000     1  20
## 1051  1051 110       292      0     62 127.0000     0  35
## 1052  1052 108       278      0     63 125.0000     1  28
## 1053  1053 129       275      0     65 135.0000     0  24
## 1054  1054 141       285      0     67 150.0000     0  23
## 1055  1055 110       276      0     70 155.0000     0  31
## 1056  1056 118       273      0     63 120.0000     0  21
## 1057  1057 111       267      1     60 115.0000     0  24
## 1058  1058 160       297      0     68 136.0000     0  20
## 1059  1059 120       280      0     60 115.0000     0  30
## 1060  1060 121       281      0     63 108.0000     0  29
## 1061  1061 113       282      0     64 118.0000     1  30
## 1062  1062 117       270      0     58 115.0000     0  23
## 1063  1063 158       267      0     64 125.0000     0  35
## 1064  1064 128       277      0     61 120.0000     0  39
## 1065  1065 158       289      0     66 140.0000     0  30
## 1066  1066 133       289      0     65 123.0000     1  22
## 1067  1067 163       298      0     61  98.0000     0  37
## 1068  1068 128       282      1     66 118.0000     0  19
## 1069  1069 126       271      1     60 105.0000     0  21
## 1070  1070 127       283      0     62 154.0000     1  42
## 1071  1071 134       287      0     63 118.0000     0  40
## 1072  1072 140       274      0     63 122.0000     0  41
## 1073  1073 102       285      0     63 128.4996     1  29
## 1074  1074 100       252      0     61 150.0000     0  24
## 1075  1075 120       295      0     59 100.0000     1  29
## 1076  1076  98       279      1     65 115.0000     1  18
## 1077  1077 130       246      0     62 118.0000     0  19
## 1078  1078 104       280      0     63 118.0000     1  41
## 1079  1079 122       285      0     62 102.0000     1  31
## 1080  1080 137       276      1     64 127.0000     0  25
## 1081  1081 114       285      1     61 104.0000     0  20
## 1082  1082  63       236      1     58  99.0000     0  24
## 1083  1083  98       318      0     63 107.0000     0  23
## 1084  1084  99       268      0     63 124.0000     1  32
## 1085  1085  89       238      1     64 136.0000     0  26
## 1086  1086 117       283      0     65 142.0000     1  22
## 1087  1087 143       281      0     67 132.0000     0  29
## 1088  1088 106       279      0     63 125.0000     1  29
## 1089  1089  99       246      0     62 106.0000     0  35
## 1090  1090 156       300      0     65 120.0000     1  27
## 1091  1091  72       266      1     66 200.0000     1  25
## 1092  1092  75       266      0     61 113.0000     1  37
## 1093  1093  97       285      0     61 112.0000     1  35
## 1094  1094 106       264      0     64 114.0000     0  41
## 1095  1095  91       225      0     68 117.0000     1  18
## 1096  1096 117       269      1     61  99.0000     0  28
## 1097  1097 117       284      0     66 177.0000     1  25
## 1098  1098 112       291      0     66 145.0000     0  23
## 1099  1099 112       270      0     61 124.0000     0  29
## 1100  1100 141       293      0     61 125.0000     0  28
## 1101  1101 131       259      0     63 134.0000     0  19
## 1102  1102 130       290      0     65 123.0000     1  19
## 1103  1103 132       270      0     67 140.0000     0  26
## 1104  1104 114       265      0     67 130.0000     1  23
## 1105  1105 160       291      0     64 110.0000     1  34
## 1106  1106 106       283      0     63 119.0000     0  24
## 1107  1107  84       260      1     64 104.0000     1  20
## 1108  1108 112       268      1     59 103.0000     0  25
## 1109  1109 139       311      0     66 135.0000     0  37
## 1110  1110 104       267      0     63 180.0000     0  30
## 1111  1111 130       294      0     63 110.0000     1  32
## 1112  1112  71       254      0     61 145.0000     1  19
## 1113  1113  82       270      0     65 150.0000     1  21
## 1114  1114 119       280      1     64 128.0000     0  21
## 1115  1115 123       353      0     63 115.0000     0  26
## 1116  1116 115       278      0     59  95.0000     0  27
## 1117  1117 124       289      1     67 145.0000     1  21
## 1118  1118 138       292      0     65 130.0000     1  25
## 1119  1119  88       276      0     63 103.0000     1  25
## 1120  1120 128       241      1     64 126.0000     0  17
## 1121  1121  82       274      0     64 101.0000     1  31
## 1122  1122 100       274      0     63 113.0000     0  24
## 1123  1123 114       271      0     61 130.0000     0  32
## 1124  1124  97       269      0     65 137.0000     1  20
## 1125  1125 126       298      0     61 112.0000     0  24
## 1126  1126 122       275      1     65 127.0000     0  20
## 1127  1127 152       295      0     62 140.0000     0  39
## 1128  1128 116       274      0     62 110.0000     1  21
## 1129  1129 132       302      0     63 145.0000     1  36
## 1130  1130  84       260      1     66 140.0000     0  37
## 1131  1131 119       277      1     61  89.0000     1  18
## 1132  1132 106       312      0     62 135.0000     1  24
## 1133  1133 139       291      0     65 160.0000     0  24
## 1134  1134 103       273      0     65 158.0000     1  36
## 1135  1135 112       299      0     67 145.0000     1  24
## 1136  1136  96       276      0     64 127.0000     1  33
## 1137  1137 102       281      1     67 135.0000     1  19
## 1138  1138 120       300      0     63 150.0000     1  34
## 1139  1139 102       338      0     64 170.0000     0  19
## 1140  1140  97       255      1     63 107.0000     1  22
## 1141  1141 113       285      0     70 145.0000     0  22
## 1142  1142 130       297      0     58 130.0000     0  32
## 1143  1143  97       260      1     63 115.0000     1  25
## 1144  1144 116       273      0     61 120.0000     0  31
## 1145  1145 114       266      0     64 113.0000     0  29
## 1146  1146 127       242      0     61 135.0000     1  17
## 1147  1147  87       247      1     66 125.0000     1  18
## 1148  1148 141       281      0     54 156.0000     1  29
## 1149  1149 144       283      1     66 140.0000     0  25
## 1150  1150 116       273      0     66 130.0000     1  33
## 1151  1151  75       265      0     65 103.0000     1  21
## 1152  1152 138       286      1     68 120.0000     0  28
## 1153  1153  99       271      0     69 151.0000     0  39
## 1154  1154 118       293      0     63 103.0000     0  21
## 1155  1155  97       266      0     62 128.4996     0  24
## 1156  1156 146       319      0     66 145.0000     0  28
## 1157  1157  81       285      0     63 150.0000     1  19
## 1158  1158 110       321      0     66 180.0000     0  28
## 1159  1159 135       284      1     60  95.0000     0  19
## 1160  1160 114       290      1     65 120.0000     1  21
## 1161  1161 124       288      1     64 116.0000     1  21
## 1162  1162 115       262      1     64 136.0000     1  23
## 1163  1163 143       281      0     65 135.0000     1  28
## 1164  1164 113       287      1     70 145.0000     1  29
## 1165  1165 109       244      1     63 102.0000     1  21
## 1166  1166 103       278      0     60  87.0000     1  30
## 1167  1167 118       276      0     64 116.0000     0  34
## 1168  1168 127       290      0     65 121.0000     0  27
## 1169  1169 132       270      0     65 126.0000     0  27
## 1170  1170 113       275      1     60 100.0000     0  27
## 1171  1171 128       265      0     67 120.0000     0  24
## 1172  1172 130       291      0     65 150.0000     1  30
## 1173  1173 125       281      1     65 110.0000     0  21
## 1174  1174 117       297      0     65 129.0000     0  38

Task 2: Testing for Outliers(1 point)

Read the babies.csv file and calculate the statistics required to find out possible outliers for baby weights.

First, read the babies.csv file and calculate the mean, standard deviation, maximum, and minimum for the bwt column using R.

#Read babies file
setwd("C:/Users/kkobayashi/Downloads/LAB 02")
babies=read.csv("babies.csv") 
dim(babies)
## [1] 1174    8
#Calculate and display the average, standard deviation, maximum and minimum babies weight. 

average_weight<-mean(babies$weight, na.rm = TRUE)
stnddeviation_weight <-sd(babies$weight, na.rm = TRUE)
max_weight <- max(babies$weight, na.rm = TRUE)
min_weight <- min(babies$weight, na.rm = TRUE)
cat("average baby weight:", average_weight, "\n")
## average baby weight: 128.4787
cat("stn baby weight:", stnddeviation_weight, "\n")
## stn baby weight: 20.73428
cat("max baby weight:", max_weight, "\n")
## max baby weight: 250
cat("min baby weight:", min_weight, "\n")
## min baby weight: 87

Use the standard approach to find out possible outliers (0.5 points)

Use the formula from class to detect any outliers. An outlier is value that “lies outside” most of the other values in a set of data.

#calculate the upper and lower thresholds
Q1 <- quantile(babies$weight, 0.25, na.rm = TRUE)
Q3 <- quantile(babies$weight, 0.75, na.rm = TRUE)
interquartile_range <- Q3 -Q1
upper_threshold <- Q3 + 1.5 * interquartile_range
lower_threshold <- Q1 - 1.5 * interquartile_range

Do you think there are possible outliers? Why?

If you decided that you have possible outliers can you identify the observations with possible outliers?

#Observations with possible outliers at the upper end
possible_outliers_up<-babies$weight[babies$weight>upper_threshold]
#observations with possible outliers at the lower end
possible_outliers_low<-babies$weight[babies$weight>lower_threshold]

What is the number of outliers at the lower end and upper end of the data?

Use the IQR approach to identify outliers (0.5 point)

Another similar method to find the upper and lower thresholds discussed in class involves finding the interquartile range. Find the interquantile range using the following chunk.

#Find the inter quantile range 
Q1 <- quantile(babies$weight, 0.25, na.rm = TRUE)
Q3 <- quantile(babies$weight, 0.75, na.rm = TRUE)
interquartile_range <- Q3 - Q1

The threshold is the boundaries that determine if a value is an outlier. If the value falls above the upper threshold or below the lower threshold, it is a possible outlier.

#calculate the upper threshold
upper_threshold <- Q3+1.5*interquartile_range
#calculate the lower threshold
lower_threshold <- Q1-1.5*interquartile_range

Do you think there are possible outliers and why?

If you think there are possible outliers can you identify them? If so how many?

possible_outliers_upper <- babies$weight[babies$weight > upper_threshold]
possible_outliers_low <- babies$weight[babies$weight > lower_threshold]
num_possible_outliers <- length(possible_outliers_upper)
length(possible_outliers_low)
## [1] 1174

It can also be useful to visualize the data using a box plot.

boxplot(babies$weight,main="Baby Weight Bowplot")

Can you identify the outliers from the boxplot? If so how many?

Task 3: Calculating Z-Value (0.25 points)

Given a birth weight of 182, calculate the corresponding z-value or z-score.

#  Calculate the z-value and display it
birth_weight <- 182
mean_weight <- mean(babies$weight, na.rm = TRUE)
sd_weight <- sd(babies$weight, na.rm = TRUE)
z_value <- (birth_weight - mean_weight) / sd_weight
cat("Z-value for birth weight:", z_value, "\n")
## Z-value for birth weight: 2.581295

Based on the z-value, how would you evaluate the weight of this baby. Explain your logic.

Task 4:Relationships with the variables (0.75 points)

How can you describe the relationship between baby weights and other variables?

#correlation between baby weight and other variables
correlation <- cor(babies[, -1], babies$weight, use = "complete.obs")

Your comments on these relationships: A correlation value close to 1 means a strong positive connection, while a value near -1 suggest a strong negative link. If thew value is close to 0, then it means there is no clear relationship between baby weight and the other factors in your dataset.