Assignment 2 (10%)

[shichao hou]]

[Insert course section & student number]


Instructions

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. Review this website for more details on using R Markdown http://rmarkdown.rstudio.com.

Use RStudio for this assignment. Complete the assignment by inserting your R code wherever you see the string “#INSERT YOUR ANSWER HERE”.

When you click the Knit button, a document (PDF, Word, or HTML format) will be generated that includes both the assignment content as well as the output of any embedded R code chunks.

Submit both the rmd and generated output files. Failing to submit both files will be subject to mark deduction.

Sample Question and Solution

Use seq() to create the vector [3,5, … ,29).

seq(3, 30, 2)
##  [1]  3  5  7  9 11 13 15 17 19 21 23 25 27 29
seq(3, 29, 2)
##  [1]  3  5  7  9 11 13 15 17 19 21 23 25 27 29

Question 1 (40 points)

The Titanic Passenger Survival Data Set provides information on the fate of passengers on the fatal maiden voyage of the ocean liner “Titanic.” The dataset is available from the Department of Biostatistics at the Vanderbilt University School of Medicine (https://biostat.app.vumc.org/wiki/pub/Main/DataSets/titanic3.csv)in several formats. store the Titanic Data Set titanic_train using the following commands.

#install.packages("titanic")
library(titanic)
titanicDataset <- read.csv(file = "https://biostat.app.vumc.org/wiki/pub/Main/DataSets/titanic3.csv", stringsAsFactors = F)
summary(titanicDataset)
##      pclass         survived         name               sex           
##  Min.   :1.000   Min.   :0.000   Length:1309        Length:1309       
##  1st Qu.:2.000   1st Qu.:0.000   Class :character   Class :character  
##  Median :3.000   Median :0.000   Mode  :character   Mode  :character  
##  Mean   :2.295   Mean   :0.382                                        
##  3rd Qu.:3.000   3rd Qu.:1.000                                        
##  Max.   :3.000   Max.   :1.000                                        
##                                                                       
##       age            sibsp            parch          ticket         
##  Min.   : 0.17   Min.   :0.0000   Min.   :0.000   Length:1309       
##  1st Qu.:21.00   1st Qu.:0.0000   1st Qu.:0.000   Class :character  
##  Median :28.00   Median :0.0000   Median :0.000   Mode  :character  
##  Mean   :29.88   Mean   :0.4989   Mean   :0.385                     
##  3rd Qu.:39.00   3rd Qu.:1.0000   3rd Qu.:0.000                     
##  Max.   :80.00   Max.   :8.0000   Max.   :9.000                     
##  NA's   :263                                                        
##       fare            cabin             embarked             boat          
##  Min.   :  0.000   Length:1309        Length:1309        Length:1309       
##  1st Qu.:  7.896   Class :character   Class :character   Class :character  
##  Median : 14.454   Mode  :character   Mode  :character   Mode  :character  
##  Mean   : 33.295                                                           
##  3rd Qu.: 31.275                                                           
##  Max.   :512.329                                                           
##  NA's   :1                                                                 
##       body        home.dest        
##  Min.   :  1.0   Length:1309       
##  1st Qu.: 72.0   Class :character  
##  Median :155.0   Mode  :character  
##  Mean   :160.8                     
##  3rd Qu.:256.0                     
##  Max.   :328.0                     
##  NA's   :1188
library(rlang)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
 #.libPaths()
 
#oldpkgs <- installed.packages(lib.loc = "Frameworks/R.framework/Versions/3.6/Resources/library")
#install.packages(oldpkgs[, 1], lib = "Frameworks/R.framework/Versions/4.0/Resources/library")
  1. Extract and show the columns sex, age, fare, cabin and survived into a new data frame of the name ‘titanicSubset’. (5 points)
#INSERT YOUR ANSWER HERE
titanicSubset <- titanicDataset[,c('sex','age','fare' ,'cabin','survived')]
titanicSubset
##         sex   age     fare           cabin survived
## 1    female 29.00 211.3375              B5        1
## 2      male  0.92 151.5500         C22 C26        1
## 3    female  2.00 151.5500         C22 C26        0
## 4      male 30.00 151.5500         C22 C26        0
## 5    female 25.00 151.5500         C22 C26        0
## 6      male 48.00  26.5500             E12        1
## 7    female 63.00  77.9583              D7        1
## 8      male 39.00   0.0000             A36        0
## 9    female 53.00  51.4792            C101        1
## 10     male 71.00  49.5042                        0
## 11     male 47.00 227.5250         C62 C64        0
## 12   female 18.00 227.5250         C62 C64        1
## 13   female 24.00  69.3000             B35        1
## 14   female 26.00  78.8500                        1
## 15     male 80.00  30.0000             A23        1
## 16     male    NA  25.9250                        0
## 17     male 24.00 247.5208         B58 B60        0
## 18   female 50.00 247.5208         B58 B60        1
## 19   female 32.00  76.2917             D15        1
## 20     male 36.00  75.2417              C6        0
## 21     male 37.00  52.5542             D35        1
## 22   female 47.00  52.5542             D35        1
## 23     male 26.00  30.0000            C148        1
## 24   female 42.00 227.5250                        1
## 25   female 29.00 221.7792             C97        1
## 26     male 25.00  26.0000                        0
## 27     male 25.00  91.0792             B49        1
## 28   female 19.00  91.0792             B49        1
## 29   female 35.00 135.6333             C99        1
## 30     male 28.00  26.5500             C52        1
## 31     male 45.00  35.5000               T        0
## 32     male 40.00  31.0000             A31        1
## 33   female 30.00 164.8667              C7        1
## 34   female 58.00  26.5500            C103        1
## 35     male 42.00  26.5500             D22        0
## 36   female 45.00 262.3750                        1
## 37   female 22.00  55.0000             E33        1
## 38     male    NA  26.5500                        1
## 39     male 41.00  30.5000             A21        0
## 40     male 48.00  50.4958             B10        0
## 41     male    NA  39.6000                        0
## 42   female 44.00  27.7208              B4        1
## 43   female 59.00  51.4792            C101        1
## 44   female 60.00  76.2917             D15        1
## 45   female 41.00 134.5000             E40        1
## 46     male 45.00  26.5500             B38        0
## 47     male    NA  31.0000                        0
## 48     male 42.00  26.2875             E24        1
## 49   female 53.00  27.4458                        1
## 50     male 36.00 512.3292     B51 B53 B55        1
## 51   female 58.00 512.3292     B51 B53 B55        1
## 52     male 33.00   5.0000     B51 B53 B55        0
## 53     male 28.00  47.1000                        0
## 54     male 17.00  47.1000                        0
## 55     male 11.00 120.0000         B96 B98        1
## 56   female 14.00 120.0000         B96 B98        1
## 57     male 36.00 120.0000         B96 B98        1
## 58   female 36.00 120.0000         B96 B98        1
## 59     male 49.00  26.0000                        0
## 60   female    NA  27.7208                        1
## 61     male 36.00  78.8500             C46        0
## 62   female 76.00  78.8500             C46        1
## 63     male 46.00  61.1750             E31        0
## 64   female 47.00  61.1750             E31        1
## 65     male 27.00  53.1000              E8        1
## 66   female 33.00  53.1000              E8        1
## 67   female 36.00 262.3750             B61        1
## 68   female 30.00  86.5000             B77        1
## 69     male 45.00  29.7000              A9        1
## 70   female    NA  55.0000             E33        1
## 71     male    NA   0.0000                        0
## 72     male 27.00 136.7792             C89        0
## 73   female 26.00 136.7792             C89        1
## 74   female 22.00 151.5500                        1
## 75     male    NA  52.0000             A14        0
## 76     male 47.00  25.5875             E58        0
## 77   female 39.00  83.1583             E49        1
## 78     male 37.00  83.1583             E52        0
## 79   female 64.00  83.1583             E45        1
## 80   female 55.00  25.7000            C101        1
## 81     male    NA  26.5500                        0
## 82     male 70.00  71.0000             B22        0
## 83   female 36.00  71.0000             B22        1
## 84   female 64.00  26.5500             B26        1
## 85     male 39.00  71.2833             C85        0
## 86   female 38.00  71.2833             C85        1
## 87     male 51.00  26.5500             E17        1
## 88     male 27.00  30.5000                        1
## 89   female 33.00 151.5500                        1
## 90     male 31.00  52.0000             B71        0
## 91   female 27.00  52.0000             B71        1
## 92     male 31.00  57.0000             B20        1
## 93   female 17.00  57.0000             B20        1
## 94     male 53.00  81.8583             A34        1
## 95     male  4.00  81.8583             A34        1
## 96   female 54.00  81.8583             A34        1
## 97     male 50.00 106.4250             C86        0
## 98   female 27.00 247.5208         B58 B60        1
## 99   female 48.00 106.4250             C86        1
## 100  female 48.00  39.6000             A16        1
## 101    male 49.00  56.9292             A20        1
## 102    male 39.00  29.7000             A18        0
## 103  female 23.00  83.1583             C54        1
## 104  female 38.00 227.5250             C45        1
## 105  female 54.00  78.2667             D20        1
## 106  female 36.00  31.6792             A29        0
## 107    male    NA 221.7792             C95        0
## 108  female    NA  31.6833                        1
## 109  female    NA 110.8833                        1
## 110    male 36.00  26.3875             E25        1
## 111    male 30.00  27.7500            C111        0
## 112  female 24.00 263.0000     C23 C25 C27        1
## 113  female 28.00 263.0000     C23 C25 C27        1
## 114  female 23.00 263.0000     C23 C25 C27        1
## 115    male 19.00 263.0000     C23 C25 C27        0
## 116    male 64.00 263.0000     C23 C25 C27        0
## 117  female 60.00 263.0000     C23 C25 C27        1
## 118  female 30.00  56.9292             E36        1
## 119    male    NA  26.5500             D34        0
## 120    male 50.00 133.6500                        1
## 121    male 43.00  27.7208             D40        1
## 122  female    NA 133.6500                        1
## 123  female 22.00  49.5000             B39        1
## 124    male 60.00  79.2000             B41        1
## 125  female 48.00  79.2000             B41        1
## 126    male    NA   0.0000            B102        0
## 127    male 37.00  53.1000            C123        0
## 128  female 35.00  53.1000            C123        1
## 129    male 47.00  38.5000             E63        0
## 130  female 35.00 211.5000            C130        1
## 131  female 22.00  59.4000                        1
## 132  female 45.00  59.4000                        1
## 133    male 24.00  79.2000             B86        0
## 134    male 49.00  89.1042             C92        1
## 135  female    NA  89.1042             C92        1
## 136    male 71.00  34.6542              A5        0
## 137    male 53.00  28.5000             C51        1
## 138  female 19.00  30.0000             B42        1
## 139    male 38.00 153.4625             C91        0
## 140  female 58.00 153.4625            C125        1
## 141    male 23.00  63.3583         D10 D12        1
## 142  female 45.00  63.3583         D10 D12        1
## 143    male 46.00  79.2000         B82 B84        0
## 144    male 25.00  55.4417             E50        1
## 145  female 25.00  55.4417             E50        1
## 146    male 48.00  76.7292             D33        1
## 147  female 49.00  76.7292             D33        1
## 148    male    NA  42.4000                        0
## 149    male 45.00  83.4750             C83        0
## 150  female 35.00  83.4750             C83        1
## 151    male 40.00   0.0000             B94        0
## 152    male 27.00  76.7292             D49        1
## 153    male    NA  30.0000             D45        1
## 154  female 24.00  83.1583             C54        1
## 155    male 55.00  93.5000             B69        0
## 156  female 52.00  93.5000             B69        1
## 157    male 42.00  42.5000             B11        0
## 158    male    NA  51.8625             E46        0
## 159    male 55.00  50.0000             C39        0
## 160  female 16.00  57.9792             B18        1
## 161  female 44.00  57.9792             B18        1
## 162  female 51.00  77.9583             D11        1
## 163    male 42.00  52.0000                        0
## 164  female 35.00  52.0000                        1
## 165    male 35.00  26.5500                        1
## 166    male 38.00  90.0000             C93        1
## 167    male    NA  30.6958                        0
## 168  female 35.00  90.0000             C93        1
## 169  female 38.00  80.0000             B28        1
## 170  female 50.00  28.7125             C49        0
## 171    male 49.00   0.0000     B52 B54 B56        1
## 172    male 46.00  26.0000                        0
## 173    male 50.00  26.0000             E60        0
## 174    male 32.50 211.5000            C132        0
## 175    male 58.00  29.7000             B37        0
## 176    male 41.00  51.8625             D21        0
## 177  female    NA  51.8625             D21        1
## 178    male 42.00  52.5542             D19        1
## 179  female 45.00  52.5542             D19        1
## 180    male    NA  26.5500            C124        0
## 181  female 39.00 211.3375                        1
## 182  female 49.00  25.9292             D17        1
## 183  female 30.00 106.4250                        1
## 184    male 35.00 512.3292            B101        1
## 185    male    NA  27.7208                        0
## 186    male 42.00  26.5500                        0
## 187  female 55.00  27.7208                        1
## 188  female 16.00  39.4000             D28        1
## 189  female 51.00  39.4000             D28        1
## 190    male 29.00  30.0000              D6        0
## 191  female 21.00  77.9583              D9        1
## 192    male 30.00  45.5000                        0
## 193  female 58.00 146.5208             B80        1
## 194  female 15.00 211.3375              B5        1
## 195    male 30.00  26.0000            C106        0
## 196  female 16.00  86.5000             B79        1
## 197    male    NA  29.7000             C47        1
## 198    male 19.00  53.1000             D30        0
## 199  female 18.00  53.1000             D30        1
## 200  female 24.00  49.5042             C90        1
## 201    male 46.00  75.2417              C6        0
## 202    male 54.00  51.8625             E46        0
## 203    male 36.00  26.2875             E25        1
## 204    male 28.00  82.1708                        0
## 205  female    NA  82.1708                        1
## 206    male 65.00  26.5500             E38        0
## 207    male 44.00  90.0000             C78        0
## 208  female 33.00  90.0000             C78        1
## 209  female 37.00  90.0000             C78        1
## 210    male 30.00  57.7500             C78        1
## 211    male 55.00  30.5000             C30        0
## 212    male 47.00  42.4000                        0
## 213    male 37.00  29.7000            C118        0
## 214  female 31.00 113.2750             D36        1
## 215  female 23.00 113.2750             D36        1
## 216    male 58.00 113.2750             D48        0
## 217  female 19.00  26.2833             D47        1
## 218    male 64.00  26.0000                        0
## 219  female 39.00 108.9000            C105        1
## 220    male    NA  25.7417                        1
## 221  female 22.00  61.9792             B36        1
## 222    male 65.00  61.9792             B30        0
## 223    male 28.50  27.7208             D43        0
## 224    male    NA   0.0000                        0
## 225    male 45.50  28.5000            C124        0
## 226    male 23.00  93.5000             B24        0
## 227    male 29.00  66.6000              C2        0
## 228  female 22.00  66.6000              C2        1
## 229    male 18.00 108.9000             C65        0
## 230  female 17.00 108.9000             C65        1
## 231  female 30.00  93.5000             B73        1
## 232    male 52.00  30.5000            C104        1
## 233    male 47.00  52.0000            C110        0
## 234  female 56.00  83.1583             C50        1
## 235    male 38.00   0.0000                        0
## 236    male    NA  39.6000                        1
## 237    male 22.00 135.6333                        0
## 238    male    NA 227.5250                        0
## 239  female 43.00 211.3375              B3        1
## 240    male 31.00  50.4958             A24        0
## 241    male 45.00  26.5500                        1
## 242    male    NA  50.0000             A32        0
## 243  female 33.00  27.7208             A11        1
## 244    male 46.00  79.2000                        0
## 245    male 36.00  40.1250             A10        0
## 246  female 33.00  86.5000             B77        1
## 247    male 55.00  59.4000                        0
## 248  female 54.00  59.4000                        1
## 249    male 33.00  26.5500                        0
## 250    male 13.00 262.3750 B57 B59 B63 B66        1
## 251  female 18.00 262.3750 B57 B59 B63 B66        1
## 252  female 21.00 262.3750 B57 B59 B63 B66        1
## 253    male 61.00 262.3750 B57 B59 B63 B66        0
## 254  female 48.00 262.3750 B57 B59 B63 B66        1
## 255    male    NA  30.5000            C106        1
## 256  female 24.00  69.3000             B35        1
## 257    male    NA  26.0000                        1
## 258  female 35.00  57.7500             C28        1
## 259  female 30.00  31.0000                        1
## 260    male 34.00  26.5500                        1
## 261  female 40.00 153.4625            C125        1
## 262    male 35.00  26.2875             E24        1
## 263    male 50.00  55.9000             E44        0
## 264  female 39.00  55.9000             E44        1
## 265    male 56.00  35.5000             A26        1
## 266    male 28.00  35.5000              A6        1
## 267    male 56.00  26.5500                        0
## 268    male 56.00  30.6958              A7        0
## 269    male 24.00  60.0000             C31        0
## 270    male    NA  26.0000             A19        0
## 271  female 18.00  60.0000             C31        1
## 272    male 24.00  82.2667             B45        1
## 273  female 23.00  82.2667             B45        1
## 274    male  6.00 134.5000             E34        1
## 275    male 45.00 134.5000             E34        1
## 276  female 40.00 134.5000             E34        1
## 277    male 57.00 146.5208             B78        0
## 278  female    NA 146.5208             B78        1
## 279    male 32.00  30.5000             B50        1
## 280    male 62.00  26.5500             C87        0
## 281    male 54.00  55.4417            C116        1
## 282  female 43.00  55.4417            C116        1
## 283  female 52.00  78.2667             D20        1
## 284    male    NA  27.7208                        0
## 285  female 62.00  80.0000             B28        1
## 286    male 67.00 221.7792         C55 C57        0
## 287  female 63.00 221.7792         C55 C57        0
## 288    male 61.00  32.3208             D50        0
## 289  female 48.00  25.9292             D17        1
## 290  female 18.00  79.6500             E68        1
## 291    male 52.00  79.6500             E67        0
## 292  female 39.00  79.6500             E67        1
## 293    male 48.00  52.0000            C126        1
## 294  female    NA  52.0000            C126        1
## 295    male 49.00 110.8833             C68        0
## 296    male 17.00 110.8833             C70        1
## 297  female 39.00 110.8833             C68        1
## 298  female    NA  79.2000                        1
## 299    male 31.00  28.5375             C53        1
## 300    male 40.00  27.7208                        0
## 301    male 61.00  33.5000             B19        0
## 302    male 47.00  34.0208             D46        0
## 303  female 35.00 512.3292                        1
## 304    male 64.00  75.2500             D37        0
## 305  female 60.00  75.2500             D37        1
## 306    male 60.00  26.5500                        0
## 307    male 54.00  77.2875             D26        0
## 308    male 21.00  77.2875             D26        0
## 309  female 55.00 135.6333             C32        1
## 310  female 31.00 164.8667              C7        1
## 311    male 57.00 164.8667                        0
## 312  female 45.00 164.8667                        1
## 313    male 50.00 211.5000             C80        0
## 314    male 27.00 211.5000             C82        0
## 315  female 50.00 211.5000             C80        1
## 316  female 21.00  26.5500                        1
## 317    male 51.00  61.3792                        0
## 318    male 21.00  61.3792                        1
## 319    male    NA  35.0000            C128        0
## 320  female 31.00 134.5000         E39 E41        1
## 321    male    NA  35.5000             C52        1
## 322    male 62.00  26.5500                        0
## 323  female 36.00 135.6333             C32        1
## 324    male 30.00  24.0000                        0
## 325  female 28.00  24.0000                        1
## 326    male 30.00  13.0000                        0
## 327    male 18.00  11.5000                        0
## 328    male 25.00  10.5000                        0
## 329    male 34.00  26.0000                        0
## 330  female 36.00  26.0000                        1
## 331    male 57.00  13.0000                        0
## 332    male 18.00  11.5000                        0
## 333    male 23.00  10.5000                        0
## 334  female 36.00  13.0000               D        1
## 335    male 28.00  10.5000                        0
## 336    male 51.00  12.5250                        0
## 337    male 32.00  26.0000                        1
## 338  female 19.00  26.0000                        1
## 339    male 28.00  26.0000                        0
## 340    male  1.00  39.0000              F4        1
## 341  female  4.00  39.0000              F4        1
## 342  female 12.00  39.0000              F4        1
## 343  female 36.00  39.0000              F4        1
## 344    male 34.00  13.0000             D56        1
## 345  female 19.00  13.0000                        1
## 346    male 23.00  13.0000                        0
## 347    male 26.00  13.0000                        0
## 348    male 42.00  13.0000                        0
## 349    male 27.00  13.0000                        0
## 350  female 24.00  13.0000             F33        1
## 351  female 15.00  39.0000                        1
## 352    male 60.00  39.0000                        0
## 353  female 40.00  39.0000                        1
## 354  female 20.00  26.0000                        1
## 355    male 25.00  26.0000                        0
## 356  female 36.00  13.0000                        1
## 357    male 25.00  13.0000                        0
## 358    male 42.00  13.0000                        0
## 359  female 42.00  13.0000                        1
## 360    male  0.83  29.0000                        1
## 361    male 26.00  29.0000                        1
## 362  female 22.00  29.0000                        1
## 363  female 35.00  21.0000                        1
## 364    male    NA   0.0000                        0
## 365    male 19.00  13.0000                        0
## 366  female 44.00  26.0000                        0
## 367    male 54.00  26.0000                        0
## 368    male 52.00  13.5000                        0
## 369    male 37.00  26.0000                        0
## 370  female 29.00  26.0000                        0
## 371  female 25.00  30.0000                        1
## 372  female 45.00  30.0000                        1
## 373    male 29.00  26.0000                        0
## 374  female 28.00  26.0000                        1
## 375    male 29.00  10.5000                        0
## 376    male 28.00  13.0000                        0
## 377    male 24.00  10.5000                        1
## 378  female  8.00  26.2500                        1
## 379    male 31.00  26.2500                        0
## 380  female 31.00  26.2500                        1
## 381  female 22.00  10.5000             F33        1
## 382  female 30.00  13.0000                        0
## 383  female    NA  21.0000                        0
## 384    male 21.00  11.5000                        0
## 385    male    NA   0.0000                        0
## 386    male  8.00  36.7500                        1
## 387    male 18.00  73.5000                        0
## 388  female 48.00  36.7500                        1
## 389  female 28.00  13.0000                        1
## 390    male 32.00  13.0000                        0
## 391    male 17.00  73.5000                        0
## 392    male 29.00  27.7208                        0
## 393  female 24.00  27.7208                        1
## 394    male 25.00  31.5000                        0
## 395    male 18.00  73.5000                        0
## 396  female 18.00  23.0000                        1
## 397  female 34.00  23.0000                        1
## 398    male 54.00  26.0000                        0
## 399    male  8.00  32.5000                        1
## 400    male 42.00  32.5000                        0
## 401  female 34.00  32.5000                        1
## 402  female 27.00  13.8583                        1
## 403  female 30.00  13.8583                        1
## 404    male 23.00  13.0000                        0
## 405    male 21.00  13.0000                        0
## 406    male 18.00  13.0000                        0
## 407    male 40.00  26.0000                        0
## 408  female 29.00  26.0000                        1
## 409    male 18.00  10.5000                        0
## 410    male 36.00  13.0000                        0
## 411    male    NA   0.0000                        0
## 412  female 38.00  13.0000                        0
## 413    male 35.00  26.0000                        0
## 414    male 38.00  21.0000                        0
## 415    male 34.00  21.0000                        0
## 416  female 34.00  13.0000                        1
## 417    male 16.00  26.0000                        0
## 418    male 26.00  10.5000                        0
## 419    male 47.00  10.5000                        0
## 420    male 21.00  11.5000                        0
## 421    male 21.00  11.5000                        0
## 422    male 24.00  13.5000                        0
## 423    male 24.00  13.0000                        0
## 424    male 34.00  13.0000                        0
## 425    male 30.00  13.0000                        0
## 426    male 52.00  13.0000                        0
## 427    male 30.00  13.0000                        0
## 428    male  0.67  14.5000                        1
## 429  female 24.00  14.5000                        1
## 430    male 44.00  13.0000                        0
## 431  female  6.00  33.0000                        1
## 432    male 28.00  33.0000                        0
## 433    male 62.00  10.5000                        1
## 434    male 30.00  10.5000                        0
## 435  female  7.00  26.2500                        1
## 436    male 43.00  26.2500                        0
## 437  female 45.00  26.2500                        1
## 438  female 24.00  65.0000                        1
## 439  female 24.00  65.0000                        1
## 440    male 49.00  65.0000                        0
## 441  female 48.00  65.0000                        1
## 442  female 55.00  16.0000                        1
## 443    male 24.00  73.5000                        0
## 444    male 32.00  73.5000                        0
## 445    male 21.00  73.5000                        0
## 446  female 18.00  13.0000                        0
## 447  female 20.00  23.0000                        1
## 448    male 23.00  11.5000                        0
## 449    male 36.00  13.0000                        0
## 450  female 54.00  23.0000                        1
## 451    male 50.00  13.0000                        0
## 452    male 44.00  26.0000                        0
## 453  female 29.00  26.0000                        1
## 454    male 21.00  73.5000                        0
## 455    male 42.00  13.0000                        1
## 456    male 63.00  26.0000                        0
## 457  female 60.00  26.0000                        0
## 458    male 33.00  12.2750                        0
## 459  female 17.00  10.5000                        1
## 460    male 42.00  27.0000                        0
## 461  female 24.00  27.0000                        1
## 462    male 47.00  15.0000                        0
## 463    male 24.00  31.5000                        0
## 464    male 22.00  31.5000                        0
## 465    male 32.00  10.5000                        0
## 466  female 23.00  13.7917               D        1
## 467    male 34.00  26.0000                        0
## 468  female 24.00  26.0000                        1
## 469  female 22.00  21.0000                        0
## 470  female    NA  12.3500            E101        1
## 471    male 35.00  12.3500                        0
## 472  female 45.00  13.5000                        1
## 473    male 57.00  12.3500                        0
## 474    male    NA   0.0000                        0
## 475    male 31.00  10.5000                        0
## 476  female 26.00  26.0000                        0
## 477    male 30.00  26.0000                        0
## 478    male    NA  10.7083                        0
## 479  female  1.00  41.5792                        1
## 480  female  3.00  41.5792                        1
## 481    male 25.00  41.5792                        0
## 482  female 22.00  41.5792                        1
## 483  female 17.00  12.0000                        1
## 484  female    NA  33.0000                        1
## 485  female 34.00  10.5000             F33        1
## 486    male 36.00  12.8750               D        0
## 487    male 24.00  10.5000                        0
## 488    male 61.00  12.3500                        0
## 489    male 50.00  26.0000                        0
## 490  female 42.00  26.0000                        1
## 491  female 57.00  10.5000             E77        0
## 492    male    NA  15.0458               D        0
## 493    male  1.00  37.0042                        1
## 494    male 31.00  37.0042                        0
## 495  female 24.00  37.0042                        1
## 496    male    NA  15.5792                        0
## 497    male 30.00  13.0000                        0
## 498    male 40.00  16.0000                        0
## 499    male 32.00  13.5000                        0
## 500    male 30.00  13.0000                        0
## 501    male 46.00  26.0000                        0
## 502  female 13.00  19.5000                        1
## 503  female 41.00  19.5000                        1
## 504    male 19.00  10.5000                        1
## 505    male 39.00  13.0000                        0
## 506    male 48.00  13.0000                        0
## 507    male 70.00  10.5000                        0
## 508    male 27.00  13.0000                        0
## 509    male 54.00  14.0000                        0
## 510    male 39.00  26.0000                        0
## 511    male 16.00  10.5000                        0
## 512    male 62.00   9.6875                        0
## 513    male 32.50  30.0708                        0
## 514  female 14.00  30.0708                        1
## 515    male  2.00  26.0000              F2        1
## 516    male  3.00  26.0000              F2        1
## 517    male 36.50  26.0000              F2        0
## 518    male 26.00  13.0000              F2        0
## 519    male 19.00  36.7500                        0
## 520    male 28.00  13.5000                        0
## 521    male 20.00  13.8625             D38        1
## 522  female 29.00  10.5000             F33        1
## 523    male 39.00  13.0000                        0
## 524    male 22.00  10.5000                        1
## 525    male    NA  13.8625                        1
## 526    male 23.00  10.5000                        0
## 527    male 29.00  13.8583                        1
## 528    male 28.00  10.5000                        0
## 529    male    NA   0.0000                        0
## 530  female 50.00  26.0000                        1
## 531    male 19.00  10.5000                        0
## 532    male    NA  15.0500                        0
## 533    male 41.00  13.0000                        0
## 534  female 21.00  21.0000                        1
## 535  female 19.00  26.0000                        1
## 536    male 43.00  21.0000                        0
## 537  female 32.00  13.0000                        1
## 538    male 34.00  13.0000                        0
## 539    male 30.00  12.7375                        1
## 540    male 27.00  15.0333                        0
## 541  female  2.00  26.0000                        1
## 542  female  8.00  26.0000                        1
## 543  female 33.00  26.0000                        1
## 544    male 36.00  10.5000                        0
## 545    male 34.00  21.0000                        0
## 546  female 30.00  21.0000                        1
## 547  female 28.00  13.0000                        1
## 548    male 23.00  15.0458                        0
## 549    male  0.83  18.7500                        1
## 550    male  3.00  18.7500                        1
## 551  female 24.00  18.7500                        1
## 552  female 50.00  10.5000                        1
## 553    male 19.00  10.5000                        0
## 554  female 21.00  10.5000                        1
## 555    male 26.00  13.0000                        0
## 556    male 25.00  13.0000                        0
## 557    male 27.00  26.0000                        0
## 558  female 25.00  26.0000                        1
## 559  female 18.00  13.0000                        1
## 560  female 20.00  36.7500                        1
## 561  female 30.00  13.0000                        1
## 562    male 59.00  13.5000                        0
## 563  female 30.00  12.3500                        1
## 564    male 35.00  10.5000                        0
## 565  female 40.00  13.0000                        1
## 566    male 25.00  13.0000                        0
## 567    male 41.00  15.0458                        0
## 568    male 25.00  10.5000                        0
## 569    male 18.50  13.0000               F        0
## 570    male 14.00  65.0000                        0
## 571  female 50.00  10.5000                        1
## 572    male 23.00  13.0000                        0
## 573  female 28.00  12.6500                        1
## 574  female 27.00  10.5000            E101        1
## 575    male 29.00  21.0000                        0
## 576  female 27.00  21.0000                        0
## 577    male 40.00  13.0000                        0
## 578  female 31.00  21.0000                        1
## 579    male 30.00  21.0000                        0
## 580    male 23.00  10.5000                        0
## 581  female 31.00  21.0000                        1
## 582    male    NA   0.0000                        0
## 583  female 12.00  15.7500                        1
## 584  female 40.00  15.7500                        1
## 585  female 32.50  13.0000            E101        1
## 586    male 27.00  26.0000                        0
## 587  female 29.00  26.0000                        1
## 588    male  2.00  23.0000                        1
## 589  female  4.00  23.0000                        1
## 590  female 29.00  23.0000                        1
## 591  female  0.92  27.7500                        1
## 592  female  5.00  27.7500                        1
## 593    male 36.00  27.7500                        0
## 594  female 33.00  27.7500                        1
## 595    male 66.00  10.5000                        0
## 596    male    NA  12.8750                        0
## 597    male 31.00  13.0000                        1
## 598    male    NA  13.0000                        1
## 599  female 26.00  13.5000                        1
## 600  female 24.00  13.0000                        0
## 601    male 42.00   7.5500                        0
## 602    male 13.00  20.2500                        0
## 603    male 16.00  20.2500                        0
## 604  female 35.00  20.2500                        1
## 605  female 16.00   7.6500                        1
## 606    male 25.00   7.6500           F G63        1
## 607    male 20.00   7.9250                        1
## 608  female 18.00   7.2292                        1
## 609    male 30.00   7.2500                        0
## 610    male 26.00   8.0500                        0
## 611  female 40.00   9.4750                        0
## 612    male  0.83   9.3500                        1
## 613  female 18.00   9.3500                        1
## 614    male 26.00  18.7875                        1
## 615    male 26.00   7.8875                        0
## 616    male 20.00   7.9250                        0
## 617    male 24.00   7.0500                        0
## 618    male 25.00   7.0500                        0
## 619    male 35.00   8.0500                        0
## 620    male 18.00   8.3000                        0
## 621    male 32.00  22.5250                        0
## 622  female 19.00   7.8542                        1
## 623    male  4.00  31.2750                        0
## 624  female  6.00  31.2750                        0
## 625  female  2.00  31.2750                        0
## 626  female 17.00   7.9250                        1
## 627  female 38.00   7.7750                        0
## 628  female  9.00  31.2750                        0
## 629  female 11.00  31.2750                        0
## 630    male 39.00  31.2750                        0
## 631    male 27.00   7.7958                        1
## 632    male 26.00   7.7750                        0
## 633  female 39.00  31.2750                        0
## 634    male 20.00   7.8542                        0
## 635    male 26.00   7.8958                        0
## 636    male 25.00  17.8000                        0
## 637  female 18.00  17.8000                        0
## 638    male 24.00   7.7750                        0
## 639    male 35.00   7.0500                        0
## 640    male  5.00  31.3875                        0
## 641    male  9.00  31.3875                        0
## 642    male  3.00  31.3875                        1
## 643    male 13.00  31.3875                        0
## 644  female  5.00  31.3875                        1
## 645    male 40.00  31.3875                        0
## 646    male 23.00   7.7958                        1
## 647  female 38.00  31.3875                        1
## 648  female 45.00   7.2250                        1
## 649    male 21.00   7.2250                        0
## 650    male 23.00   7.0500                        0
## 651  female 17.00  14.4583                        0
## 652    male 30.00   7.2250                        0
## 653    male 23.00   7.8542                        0
## 654  female 13.00   7.2292                        1
## 655    male 20.00   7.2250                        0
## 656    male 32.00  15.8500                        0
## 657  female 33.00  15.8500                        1
## 658  female  0.75  19.2583                        1
## 659  female  0.75  19.2583                        1
## 660  female  5.00  19.2583                        1
## 661  female 24.00  19.2583                        1
## 662  female 18.00   8.0500                        1
## 663    male 40.00   7.2250                        0
## 664    male 26.00   7.8958                        0
## 665    male 20.00   7.2292                        1
## 666  female 18.00  14.4542                        0
## 667  female 45.00  14.4542                        0
## 668  female 27.00   7.8792                        0
## 669    male 22.00   8.0500                        0
## 670    male 19.00   8.0500                        0
## 671    male 26.00   7.7750                        0
## 672    male 22.00   9.3500                        0
## 673    male    NA   7.2292                        0
## 674    male 20.00   4.0125                        0
## 675    male 32.00  56.4958                        1
## 676    male 21.00   7.7750                        0
## 677    male 18.00   7.7500                        0
## 678    male 26.00   7.8958                        0
## 679    male  6.00  15.2458                        0
## 680  female  9.00  15.2458                        0
## 681    male    NA   7.2250                        0
## 682  female    NA  15.2458                        0
## 683  female    NA   7.7500                        0
## 684    male 40.00  15.5000                        0
## 685  female 32.00  15.5000                        0
## 686    male 21.00  16.1000                        0
## 687  female 22.00   7.7250                        1
## 688  female 20.00   7.8542                        0
## 689    male 29.00   7.0458                        0
## 690    male 22.00   7.2500                        0
## 691    male 22.00   7.7958                        0
## 692    male 35.00   8.0500                        0
## 693  female 18.50   7.2833                        0
## 694    male 21.00   7.8208                        1
## 695    male 19.00   6.7500                        0
## 696  female 18.00   7.8792                        0
## 697  female 21.00   8.6625                        0
## 698  female 30.00   8.6625                        0
## 699    male 18.00   8.6625                        0
## 700    male 38.00   8.6625                        0
## 701    male 17.00   8.6625                        0
## 702    male 17.00   8.6625                        0
## 703  female 21.00   7.7500                        0
## 704    male 21.00   7.7500                        0
## 705    male 21.00   8.0500                        0
## 706    male    NA  14.4583                        0
## 707  female    NA  14.4583                        0
## 708    male 28.00   7.7958                        0
## 709    male 24.00   7.8542                        0
## 710  female 16.00   7.7500                        1
## 711  female 37.00   7.7500                        0
## 712    male 28.00   7.2500                        0
## 713    male 24.00   8.0500                        0
## 714    male 21.00   7.7333                        0
## 715    male 32.00  56.4958                        1
## 716    male 29.00   8.0500                        0
## 717    male 26.00  14.4542                        0
## 718    male 18.00  14.4542                        0
## 719    male 20.00   7.0500                        0
## 720    male 18.00   8.0500                        1
## 721    male 24.00   7.2500                        0
## 722    male 36.00   7.4958                        0
## 723    male 24.00   7.4958                        0
## 724    male 31.00   7.7333                        0
## 725    male 31.00   7.7500                        0
## 726  female 22.00   7.7500                        1
## 727  female 30.00   7.6292                        0
## 728    male 70.50   7.7500                        0
## 729    male 43.00   8.0500                        0
## 730    male 35.00   7.8958                        0
## 731    male 27.00   7.8958                        0
## 732    male 19.00   7.8958                        0
## 733    male 30.00   8.0500                        0
## 734    male  9.00  15.9000                        1
## 735    male  3.00  15.9000                        1
## 736  female 36.00  15.9000                        1
## 737    male 59.00   7.2500                        0
## 738    male 19.00   8.1583                        0
## 739  female 17.00  16.1000                        1
## 740    male 44.00  16.1000                        0
## 741    male 17.00   8.6625                        0
## 742    male 22.50   7.2250                        0
## 743    male 45.00   8.0500                        1
## 744  female 22.00  10.5167                        0
## 745    male 19.00  10.1708                        0
## 746  female 30.00   6.9500                        1
## 747    male 29.00   7.7500                        1
## 748    male  0.33  14.4000                        0
## 749    male 34.00  14.4000                        0
## 750  female 28.00  14.4000                        0
## 751    male 27.00   7.8958                        0
## 752    male 25.00   7.8958                        0
## 753    male 24.00  24.1500                        0
## 754    male 22.00   8.0500                        0
## 755    male 21.00  24.1500                        0
## 756    male 17.00   8.0500                        0
## 757    male    NA  16.1000                        0
## 758  female    NA  16.1000                        1
## 759    male 36.50  17.4000                        1
## 760  female 36.00  17.4000                        1
## 761    male 30.00   9.5000                        1
## 762    male 16.00   9.5000                        0
## 763    male  1.00  20.5750                        1
## 764  female  0.17  20.5750                        1
## 765    male 26.00  20.5750                        0
## 766  female 33.00  20.5750                        1
## 767    male 25.00   7.8958                        0
## 768    male    NA   7.8958                        0
## 769    male    NA   7.8958                        0
## 770    male 22.00   7.2500                        0
## 771    male 36.00   7.2500                        0
## 772  female 19.00   7.8792                        1
## 773    male 17.00   7.8958                        0
## 774    male 42.00   8.6625                        0
## 775    male 43.00   7.8958                        0
## 776    male    NA   7.2292                        0
## 777    male 32.00   7.7500                        0
## 778    male 19.00   8.0500                        1
## 779  female 30.00  12.4750                        1
## 780  female 24.00   7.7500                        0
## 781  female 23.00   8.0500                        1
## 782    male 33.00   7.8958                        0
## 783    male 65.00   7.7500                        0
## 784    male 24.00   7.5500                        1
## 785    male 23.00  13.9000                        0
## 786  female 22.00  13.9000                        1
## 787    male 18.00   7.7750                        0
## 788    male 16.00   7.7750                        0
## 789    male 45.00   6.9750                        0
## 790    male    NA   7.2250                        0
## 791    male 39.00   7.2292                        0
## 792    male 17.00   7.2292                        0
## 793    male 15.00   7.2292                        0
## 794    male 47.00   7.2500                        0
## 795  female  5.00  12.4750                        1
## 796    male    NA   7.2250                        0
## 797    male 40.50  15.1000                        0
## 798    male 40.50   7.7500                        0
## 799    male    NA   7.0500                        1
## 800    male 18.00   7.7958                        0
## 801  female    NA   7.7500                        0
## 802    male    NA   7.7500                        0
## 803    male    NA   6.9500                        0
## 804    male 26.00   7.8792                        0
## 805    male    NA   7.7500                        0
## 806    male    NA  56.4958                        1
## 807  female 21.00  34.3750                        0
## 808  female  9.00  34.3750                        0
## 809    male    NA   8.0500                        0
## 810    male 18.00  34.3750                        0
## 811    male 16.00  34.3750                        0
## 812  female 48.00  34.3750                        0
## 813    male    NA   7.7500                        0
## 814    male    NA   7.2500                        0
## 815    male 25.00   7.7417                        0
## 816    male    NA  14.5000                        0
## 817    male    NA   7.8958                        0
## 818    male 22.00   8.0500                        0
## 819  female 16.00   7.7333                        1
## 820  female    NA   7.7500                        1
## 821    male  9.00  20.5250                        1
## 822    male 33.00  20.5250                        0
## 823    male 41.00   7.8500                        0
## 824  female 31.00  20.5250                        1
## 825    male 38.00   7.0500                        0
## 826    male  9.00  46.9000                        0
## 827    male  1.00  46.9000                        0
## 828    male 11.00  46.9000                        0
## 829  female 10.00  46.9000                        0
## 830  female 16.00  46.9000                        0
## 831    male 14.00  46.9000                        0
## 832    male 40.00  46.9000                        0
## 833  female 43.00  46.9000                        0
## 834    male 51.00   8.0500                        0
## 835    male 32.00   8.3625                        0
## 836    male    NA   8.0500                        0
## 837    male 20.00   9.8458                        0
## 838    male 37.00   7.9250                        0
## 839    male 28.00   7.9250                        0
## 840    male 19.00   7.7750                        0
## 841  female 24.00   8.8500                        0
## 842  female 17.00   7.7333                        0
## 843    male    NA  19.9667                        0
## 844    male    NA  19.9667                        0
## 845    male 28.00  15.8500                        0
## 846  female 24.00  15.8500                        1
## 847    male 20.00   9.5000                        0
## 848    male 23.50   7.2292                        0
## 849    male 41.00  14.1083                        0
## 850    male 26.00   7.8542                        0
## 851    male 21.00   7.8542                        0
## 852  female 45.00  14.1083                        1
## 853  female    NA   7.5500                        0
## 854    male 25.00   7.2500                        0
## 855    male    NA   6.8583                        0
## 856    male 11.00  18.7875                        0
## 857  female    NA   7.7500                        1
## 858    male 27.00   6.9750                        1
## 859    male    NA  56.4958                        1
## 860  female 18.00   6.7500                        0
## 861  female 26.00   7.9250                        1
## 862  female 23.00   7.9250                        0
## 863  female 22.00   8.9625                        1
## 864    male 28.00   7.8958                        0
## 865  female 28.00   7.7750                        0
## 866  female    NA   7.7500                        0
## 867  female  2.00  12.2875                        1
## 868  female 22.00  12.2875                        1
## 869    male 43.00   6.4500                        0
## 870    male 28.00  22.5250                        0
## 871  female 27.00   7.9250                        1
## 872    male    NA   7.7500                        0
## 873  female    NA   8.0500                        1
## 874    male 42.00   7.6500           F G63        0
## 875    male    NA   7.8875                        1
## 876    male 30.00   7.2292                        0
## 877    male    NA   7.8958                        0
## 878  female 27.00   7.9250                        0
## 879  female 25.00   7.9250                        0
## 880    male    NA   7.8958                        0
## 881    male 29.00   7.8958                        1
## 882    male 21.00   7.7958                        1
## 883    male    NA   7.0500                        0
## 884    male 20.00   7.8542                        0
## 885    male 48.00   7.8542                        0
## 886    male 17.00   7.0542                        0
## 887  female    NA   7.7500                        1
## 888    male    NA   8.1125                        1
## 889    male 34.00   6.4958                        0
## 890    male 26.00   7.7750                        1
## 891    male 22.00   7.7958                        0
## 892    male 33.00   8.6542                        0
## 893    male 31.00   7.7750                        0
## 894    male 29.00   7.8542                        0
## 895    male  4.00  11.1333                        1
## 896  female  1.00  11.1333                        1
## 897    male 49.00   0.0000                        0
## 898    male 33.00   7.7750                        0
## 899    male 19.00   0.0000                        0
## 900  female 27.00  11.1333                        1
## 901    male    NA  23.4500                        0
## 902  female    NA  23.4500                        0
## 903    male    NA  23.4500                        0
## 904  female    NA  23.4500                        0
## 905    male 23.00   7.8958                        0
## 906    male 32.00   7.8542                        1
## 907    male 27.00   7.8542                        0
## 908  female 20.00   9.8250                        0
## 909  female 21.00   9.8250                        0
## 910    male 32.00   7.9250                        1
## 911    male 17.00   7.1250                        0
## 912    male 21.00   8.4333                        0
## 913    male 30.00   7.8958                        0
## 914    male 21.00   7.7958                        1
## 915    male 33.00   7.8542                        0
## 916    male 22.00   7.5208                        0
## 917  female  4.00  13.4167                        1
## 918    male 39.00  13.4167                        1
## 919    male    NA   7.2292                        0
## 920    male 18.50   7.2292                        0
## 921    male    NA   7.7500                        0
## 922    male    NA   7.2500                        0
## 923  female    NA   7.7500                        1
## 924  female    NA   7.7500                        1
## 925    male 34.50   7.8292                        0
## 926    male 44.00   8.0500                        0
## 927    male    NA   7.7500                        1
## 928    male    NA  14.4542                        0
## 929  female    NA  14.4542                        0
## 930    male    NA   7.7500                        0
## 931    male    NA   7.7500                        0
## 932    male    NA   7.7375                        0
## 933  female 22.00   8.6625                        0
## 934    male 26.00   8.6625                        0
## 935  female  4.00  22.0250                        1
## 936    male 29.00  22.0250                        1
## 937  female 26.00  22.0250                        1
## 938  female  1.00  12.1833                        0
## 939    male 18.00   7.8542                        0
## 940  female 36.00  12.1833                        0
## 941    male    NA   7.8958                        0
## 942    male 25.00   7.2292           F E57        1
## 943    male    NA   7.2250                        0
## 944  female 37.00   9.5875                        0
## 945    male    NA   7.8958                        0
## 946    male    NA  56.4958                        1
## 947    male    NA  56.4958                        0
## 948  female 22.00   7.2500                        1
## 949    male    NA   7.7500                        0
## 950    male 26.00  56.4958                        1
## 951    male 29.00   9.4833                        0
## 952    male 29.00   7.7750                        0
## 953    male 22.00   7.7750                        0
## 954    male 22.00   7.2250                        1
## 955    male    NA  25.4667                        0
## 956  female    NA  25.4667                        0
## 957  female    NA  25.4667                        0
## 958  female    NA  25.4667                        0
## 959  female    NA  25.4667                        0
## 960    male 32.00   7.9250                        0
## 961    male 34.50   6.4375                        0
## 962  female    NA  15.5000                        0
## 963    male    NA  15.5000                        0
## 964    male 36.00   0.0000                        0
## 965    male 39.00  24.1500                        0
## 966    male 24.00   9.5000                        0
## 967  female 25.00   7.7750                        0
## 968  female 45.00   7.7500                        0
## 969    male 36.00  15.5500                        0
## 970  female 30.00  15.5500                        0
## 971    male 20.00   7.9250                        1
## 972    male    NA   7.8792                        0
## 973    male 28.00  56.4958                        0
## 974    male    NA   7.5500                        0
## 975    male 30.00  16.1000                        0
## 976  female 26.00  16.1000                        0
## 977    male    NA   7.8792                        0
## 978    male 20.50   7.2500                        0
## 979    male 27.00   8.6625                        1
## 980    male 51.00   7.0542                        0
## 981  female 23.00   7.8542                        1
## 982    male 32.00   7.5792                        1
## 983    male    NA   7.8958                        0
## 984    male    NA   7.5500                        0
## 985  female    NA   7.7500                        1
## 986    male 24.00   7.1417                        1
## 987    male 22.00   7.1250                        0
## 988  female    NA   7.8792                        0
## 989    male    NA   7.7500                        0
## 990    male    NA   8.0500                        0
## 991    male 29.00   7.9250                        0
## 992    male    NA   7.2292                        1
## 993  female 30.50   7.7500                        0
## 994  female    NA   7.7375                        1
## 995    male    NA   7.2292           F E46        0
## 996    male 35.00   7.8958                        0
## 997    male 33.00   7.8958                        0
## 998  female    NA   7.2250                        1
## 999    male    NA   7.8958                        0
## 1000 female    NA   7.7500                        1
## 1001   male    NA   7.7500                        1
## 1002 female    NA  23.2500                        1
## 1003 female    NA  23.2500                        1
## 1004   male    NA  23.2500                        1
## 1005 female    NA   7.7875                        1
## 1006   male    NA  15.5000                        0
## 1007 female    NA   7.8792                        1
## 1008 female 15.00   8.0292                        1
## 1009 female 35.00   7.7500                        0
## 1010   male    NA   7.7500                        0
## 1011   male 24.00  16.1000                        0
## 1012 female 19.00  16.1000                        0
## 1013 female    NA   7.7500                        0
## 1014 female    NA   8.0500                        0
## 1015 female    NA   8.0500                        0
## 1016   male 55.50   8.0500                        0
## 1017   male    NA   7.7500                        0
## 1018   male 21.00   7.7750                        1
## 1019   male    NA   8.0500                        0
## 1020   male 24.00   7.8958                        0
## 1021   male 21.00   7.8958                        0
## 1022   male 28.00   7.8958                        0
## 1023   male    NA   7.8958                        0
## 1024 female    NA   7.8792                        1
## 1025   male 25.00   7.6500           F G73        0
## 1026   male  6.00  12.4750            E121        1
## 1027 female 27.00  12.4750            E121        1
## 1028   male    NA   8.0500                        0
## 1029 female    NA  24.1500                        1
## 1030   male    NA  24.1500                        0
## 1031   male    NA   8.4583                        0
## 1032   male 34.00   8.0500                        0
## 1033   male    NA   7.7500                        0
## 1034   male    NA   7.7750                        1
## 1035   male    NA  15.2458                        1
## 1036   male    NA  15.2458                        1
## 1037 female    NA  15.2458                        1
## 1038 female    NA   7.2292                        1
## 1039   male    NA   8.0500                        0
## 1040 female    NA   7.7333                        1
## 1041 female 24.00   7.7500                        1
## 1042   male    NA   8.0500                        0
## 1043 female    NA  15.5000                        1
## 1044 female    NA  15.5000                        1
## 1045 female    NA  15.5000                        1
## 1046   male 18.00   7.7500                        0
## 1047   male 22.00   7.8958                        0
## 1048 female 15.00   7.2250                        1
## 1049 female  1.00  15.7417                        1
## 1050   male 20.00  15.7417                        1
## 1051 female 19.00  15.7417                        1
## 1052   male 33.00   8.0500                        0
## 1053   male    NA   7.8958                        0
## 1054   male    NA   7.2292                        0
## 1055 female    NA   7.7500                        0
## 1056   male    NA   7.8958                        0
## 1057   male 12.00  11.2417                        1
## 1058 female 14.00  11.2417                        1
## 1059 female 29.00   7.9250                        0
## 1060   male 28.00   8.0500                        0
## 1061 female 18.00   7.7750                        1
## 1062 female 26.00   7.8542                        1
## 1063   male 21.00   7.8542                        0
## 1064   male 41.00   7.1250                        0
## 1065   male 39.00   7.9250                        1
## 1066   male 21.00   7.8000                        0
## 1067   male 28.50   7.2292                        0
## 1068 female 22.00   7.7500                        1
## 1069   male 61.00   6.2375                        0
## 1070   male    NA  15.5000                        0
## 1071   male    NA   7.8292                        0
## 1072 female    NA  15.5000                        1
## 1073   male    NA   7.7333                        0
## 1074   male    NA   7.7500                        0
## 1075   male    NA   7.7500                        0
## 1076   male 23.00   9.2250                        0
## 1077 female    NA   7.7500                        0
## 1078 female    NA   7.7500                        1
## 1079 female    NA   7.8792                        1
## 1080 female 22.00   7.7750                        1
## 1081   male    NA   7.7500                        1
## 1082 female    NA   7.8292                        1
## 1083   male  9.00   3.1708                        1
## 1084   male 28.00  22.5250                        0
## 1085   male 42.00   8.4042                        0
## 1086   male    NA   7.3125                        0
## 1087 female 31.00   7.8542                        0
## 1088   male 28.00   7.8542                        0
## 1089   male 32.00   7.7750                        1
## 1090   male 20.00   9.2250                        0
## 1091 female 23.00   8.6625                        0
## 1092 female 20.00   8.6625                        0
## 1093   male 20.00   8.6625                        0
## 1094   male 16.00   9.2167                        0
## 1095 female 31.00   8.6833                        1
## 1096 female    NA   7.6292                        0
## 1097   male  2.00  21.0750                        0
## 1098   male  6.00  21.0750                        0
## 1099 female  3.00  21.0750                        0
## 1100 female  8.00  21.0750                        0
## 1101 female 29.00  21.0750                        0
## 1102   male  1.00  39.6875                        0
## 1103   male  7.00  39.6875                        0
## 1104   male  2.00  39.6875                        0
## 1105   male 16.00  39.6875                        0
## 1106   male 14.00  39.6875                        0
## 1107 female 41.00  39.6875                        0
## 1108   male 21.00   8.6625                        0
## 1109   male 19.00  14.5000                        0
## 1110   male    NA   8.7125                        0
## 1111   male 32.00   7.8958                        0
## 1112   male  0.75  13.7750                        0
## 1113 female  3.00  13.7750                        0
## 1114 female 26.00  13.7750                        0
## 1115   male    NA   7.0000                        0
## 1116   male    NA   7.7750                        0
## 1117   male    NA   8.0500                        0
## 1118   male 21.00   7.9250                        0
## 1119   male 25.00   7.9250                        0
## 1120   male 22.00   7.2500                        0
## 1121   male 25.00   7.7750                        1
## 1122   male    NA  22.3583                        1
## 1123 female    NA  22.3583           F E69        1
## 1124 female    NA  22.3583                        1
## 1125 female    NA   8.1375                        0
## 1126   male 24.00   8.0500                        0
## 1127 female 28.00   7.8958                        0
## 1128   male 19.00   7.8958                        0
## 1129   male    NA   7.8958                        0
## 1130   male 25.00   7.7750                        0
## 1131 female 18.00   7.7750                        0
## 1132   male 32.00   8.0500             E10        1
## 1133   male    NA   7.8958                        0
## 1134   male 17.00   8.6625                        0
## 1135   male 24.00   8.6625                        0
## 1136   male    NA   7.8958                        0
## 1137 female    NA   8.1125                        0
## 1138   male    NA   7.2292                        0
## 1139   male    NA   7.2500                        0
## 1140   male 38.00   7.8958                        0
## 1141   male 21.00   8.0500                        0
## 1142   male 10.00  29.1250                        0
## 1143   male  4.00  29.1250                        0
## 1144   male  7.00  29.1250                        0
## 1145   male  2.00  29.1250                        0
## 1146   male  8.00  29.1250                        0
## 1147 female 39.00  29.1250                        0
## 1148 female 22.00  39.6875                        0
## 1149   male 35.00   7.1250                        0
## 1150 female    NA   7.7208                        1
## 1151   male    NA  14.5000                        0
## 1152 female    NA  14.5000                        0
## 1153   male 50.00  14.5000                        0
## 1154 female 47.00  14.5000                        0
## 1155   male    NA   8.0500                        0
## 1156   male    NA   7.7750                        0
## 1157 female  2.00  20.2125                        0
## 1158   male 18.00  20.2125                        0
## 1159 female 41.00  20.2125                        0
## 1160 female    NA   8.0500                        1
## 1161   male 50.00   8.0500                        0
## 1162   male 16.00   8.0500                        0
## 1163   male    NA   7.7500                        1
## 1164   male    NA  24.1500                        0
## 1165   male    NA   7.2292                        0
## 1166   male 25.00   7.2250                        0
## 1167   male    NA   7.2250                        0
## 1168   male    NA   7.7292                        0
## 1169   male    NA   7.5750                        0
## 1170   male 38.50   7.2500                        0
## 1171   male    NA  69.5500                        0
## 1172   male 14.50  69.5500                        0
## 1173 female    NA  69.5500                        0
## 1174 female    NA  69.5500                        0
## 1175 female    NA  69.5500                        0
## 1176 female    NA  69.5500                        0
## 1177   male    NA  69.5500                        0
## 1178   male    NA  69.5500                        0
## 1179   male    NA  69.5500                        0
## 1180   male    NA  69.5500                        0
## 1181 female    NA  69.5500                        0
## 1182   male 24.00   9.3250                        0
## 1183 female 21.00   7.6500                        1
## 1184   male 39.00   7.9250                        0
## 1185   male    NA  21.6792                        0
## 1186   male    NA  21.6792                        0
## 1187   male    NA  21.6792                        0
## 1188 female  1.00  16.7000              G6        1
## 1189 female 24.00  16.7000              G6        1
## 1190 female  4.00  16.7000              G6        1
## 1191   male 25.00   9.5000                        1
## 1192   male 20.00   8.0500                        0
## 1193   male 24.50   8.0500                        0
## 1194   male    NA   7.7250                        0
## 1195   male    NA   7.8958                        0
## 1196   male    NA   7.7500                        0
## 1197   male 29.00   9.5000                        1
## 1198   male    NA  15.1000                        0
## 1199 female    NA   7.7792                        1
## 1200   male    NA   8.0500                        0
## 1201   male    NA   8.0500                        0
## 1202   male 22.00   7.2292                        0
## 1203   male    NA   8.0500                        0
## 1204   male 40.00   7.8958                        0
## 1205   male 21.00   7.9250                        0
## 1206 female 18.00   7.4958                        1
## 1207   male  4.00  27.9000                        0
## 1208   male 10.00  27.9000                        0
## 1209 female  9.00  27.9000                        0
## 1210 female  2.00  27.9000                        0
## 1211   male 40.00  27.9000                        0
## 1212 female 45.00  27.9000                        0
## 1213   male    NA   7.8958                        0
## 1214   male    NA   8.0500                        0
## 1215   male    NA   8.6625                        0
## 1216   male    NA   7.7500                        0
## 1217 female    NA   7.7333                        1
## 1218   male 19.00   7.6500           F G73        0
## 1219   male 30.00   8.0500                        0
## 1220   male    NA   8.0500                        0
## 1221   male 32.00   8.0500                        0
## 1222   male    NA   7.8958                        0
## 1223   male 33.00   8.6625                        0
## 1224 female 23.00   7.5500                        1
## 1225   male 21.00   8.0500                        0
## 1226   male 60.50       NA                        0
## 1227   male 19.00   7.8958                        0
## 1228 female 22.00   9.8375                        0
## 1229   male 31.00   7.9250                        1
## 1230   male 27.00   8.6625                        0
## 1231 female  2.00  10.4625              G6        0
## 1232 female 29.00  10.4625              G6        0
## 1233   male 16.00   8.0500                        1
## 1234   male 44.00   7.9250                        1
## 1235   male 25.00   7.0500                        0
## 1236   male 74.00   7.7750                        0
## 1237   male 14.00   9.2250                        1
## 1238   male 24.00   7.7958                        0
## 1239   male 25.00   7.7958                        1
## 1240   male 34.00   8.0500                        0
## 1241   male  0.42   8.5167                        1
## 1242   male    NA   6.4375                        0
## 1243   male    NA   6.4375                        0
## 1244   male    NA   7.2250                        0
## 1245 female 16.00   8.5167                        1
## 1246   male    NA   8.0500                        0
## 1247   male    NA  16.1000                        0
## 1248 female    NA  16.1000                        1
## 1249   male 32.00   7.9250                        0
## 1250   male    NA   7.7500             F38        0
## 1251   male    NA   7.8958                        0
## 1252   male 30.50   8.0500                        0
## 1253   male 44.00   8.0500                        0
## 1254   male    NA   7.2292                        0
## 1255   male 25.00   0.0000                        1
## 1256   male    NA   7.2292                        0
## 1257   male  7.00  15.2458                        1
## 1258 female  9.00  15.2458                        1
## 1259 female 29.00  15.2458                        1
## 1260   male 36.00   7.8958                        0
## 1261 female 18.00   9.8417                        1
## 1262 female 63.00   9.5875                        1
## 1263   male    NA  14.5000                        0
## 1264   male 11.50  14.5000                        0
## 1265   male 40.50  14.5000                        0
## 1266 female 10.00  24.1500                        0
## 1267   male 36.00  24.1500                        0
## 1268 female 30.00  24.1500                        0
## 1269   male    NA   9.5000                        0
## 1270   male 33.00   9.5000                        0
## 1271   male 28.00   9.5000                        0
## 1272   male 28.00   9.5000                        0
## 1273   male 47.00   9.0000                        0
## 1274 female 18.00  18.0000                        0
## 1275   male 31.00  18.0000                        0
## 1276   male 16.00  18.0000                        0
## 1277 female 31.00  18.0000                        0
## 1278   male 22.00   7.2250                        1
## 1279   male 20.00   7.8542                        0
## 1280 female 14.00   7.8542                        0
## 1281   male 22.00   7.8958                        0
## 1282   male 22.00   9.0000                        0
## 1283   male    NA   8.0500                        0
## 1284   male    NA   7.5500                        0
## 1285   male    NA   8.0500                        0
## 1286   male 32.50   9.5000                        0
## 1287 female 38.00   7.2292                        1
## 1288   male 51.00   7.7500                        0
## 1289   male 18.00   6.4958                        0
## 1290   male 21.00   6.4958                        0
## 1291 female 47.00   7.0000                        1
## 1292   male    NA   8.7125                        0
## 1293   male    NA   7.5500                        0
## 1294   male    NA   8.0500                        0
## 1295   male 28.50  16.1000                        0
## 1296   male 21.00   7.2500                        0
## 1297   male 27.00   8.6625                        0
## 1298   male    NA   7.2500                        0
## 1299   male 36.00   9.5000                        0
## 1300   male 27.00  14.4542                        0
## 1301 female 15.00  14.4542                        1
## 1302   male 45.50   7.2250                        0
## 1303   male    NA   7.2250                        0
## 1304   male    NA  14.4583                        0
## 1305 female 14.50  14.4542                        0
## 1306 female    NA  14.4542                        0
## 1307   male 26.50   7.2250                        0
## 1308   male 27.00   7.2250                        0
## 1309   male 29.00   7.8750                        0
  1. Numerical data: Use the count() function from the dplyr package to display the total number of passengers within each Ticket Class Pclass. (5 points)
#INSERT YOUR ANSWER HERE
#library(dplyr)
titanicDataset %>% count(pclass)
##   pclass   n
## 1      1 323
## 2      2 277
## 3      3 709
  1. Categorical data: Use count() and group_by() functions from the dplyr package to calculate the number of passengers by sex. (5 points) HINT: use group_by() first then pipe the result to count() to calculate the number of passengers per sex.
#INSERT YOUR ANSWER HERE
#library(dplyr)
titanicDataset <- titanicDataset %>% group_by(sex)
count(titanicDataset,pclass)
## # A tibble: 6 × 3
## # Groups:   sex [2]
##   sex    pclass     n
##   <chr>   <int> <int>
## 1 female      1   144
## 2 female      2   106
## 3 female      3   216
## 4 male        1   179
## 5 male        2   171
## 6 male        3   493
  1. Use the aggregate() function to calculate the ‘survivalCount’ of each sex and calculate the survival rate of each sex. Then draw the conclusion on which sex has the higher survival rate. (5 points)
#INSERT YOUR ANSWER HERE
titanicDataset %>%group_by(sex)%>% count()
## # A tibble: 2 × 2
## # Groups:   sex [2]
##   sex        n
##   <chr>  <int>
## 1 female   466
## 2 male     843
titanicSubset %>%group_by(survived)%>% count()
## # A tibble: 2 × 2
## # Groups:   survived [2]
##   survived     n
##      <int> <int>
## 1        0   809
## 2        1   500
#INSERT YOUR ANSWER HERE
aggregate(titanicSubset$survived, by = list(titanicSubset$sex),FUN = sum)
##   Group.1   x
## 1  female 339
## 2    male 161
# female has higher survival rate
  1. Using boxplot to display the distribution of age for each sex and infer which group is older between male and female. (5 points)
#INSERT YOUR ANSWER HERE
titanicDataset$sex <- as.numeric(titanicDataset$sex) 
## Warning: NAs introduced by coercion
titanicDataset$age <- as.numeric(titanicDataset$age) 
boxplot(titanicDataset$sex, titanicDataset$age)

#titanicDataset.boxplot(column='age', by='sex', grid=False)
#bp = titanicDataset.boxplot(column='age', by='sex', grid=False)

#for i in [1,2,3]:
    #y = titanicDataset.age[titanic.sex==i].dropna()

    #x = np.random.normal(i, 0.04, size=len(y))
    #plot(x, y, 'r.', alpha=0.2)
  1. Calculate the average age for male and female and describe if the calculation agrees with the box plot. (5 points)
#INSERT YOUR ANSWER HERE
#library(titanic)
#titanicDataset <- read.csv(file = "https://biostat.app.vumc.org/wiki/pub/Main/DataSets/titanic3.csv", stringsAsFactors = F)
#summary(titanicDataset)
#titanicDataset

#titanicSubset <- titanicDataset[,c('sex','age')]
#titanicSubset

aggregate(titanicSubset$survived, by = list(titanicSubset$sex),FUN = mean)
##   Group.1         x
## 1  female 0.7274678
## 2    male 0.1909846
#it is matching with boxplot relatively
  1. Use the ‘for’ loop and ‘if’ control statements to list the women’s names, age 34 or more that embarked from S (Southampton), on the Titanic. (10 points)
#INSERT YOUR ANSWER HERE
for (x in 1:(length(titanicDataset$sex))){
  if (titanicDataset$age[x]>=34 & !
      is.na(titanicDataset$age[x])){print(titanicDataset$name[x])}}
## [1] "Anderson, Mr. Harry"
## [1] "Andrews, Miss. Kornelia Theodosia"
## [1] "Andrews, Mr. Thomas Jr"
## [1] "Appleton, Mrs. Edward Dale (Charlotte Lamson)"
## [1] "Artagaveytia, Mr. Ramon"
## [1] "Astor, Col. John Jacob"
## [1] "Barkworth, Mr. Algernon Henry Wilson"
## [1] "Baxter, Mrs. James (Helene DeLaudeniere Chaput)"
## [1] "Beattie, Mr. Thomson"
## [1] "Beckwith, Mr. Richard Leonard"
## [1] "Beckwith, Mrs. Richard Leonard (Sallie Monypeny)"
## [1] "Bidois, Miss. Rosalie"
## [1] "Bissette, Miss. Amelia"
## [1] "Blackwell, Mr. Stephen Weart"
## [1] "Blank, Mr. Henry"
## [1] "Bonnell, Miss. Elizabeth"
## [1] "Borebank, Mr. John James"
## [1] "Bowen, Miss. Grace Scott"
## [1] "Brady, Mr. John Bertram"
## [1] "Brandeis, Mr. Emil"
## [1] "Brown, Mrs. James Joseph (Margaret Tobin)"
## [1] "Brown, Mrs. John Murray (Caroline Lane Lamson)"
## [1] "Bucknell, Mrs. William Robert (Emma Eliza Ward)"
## [1] "Burns, Miss. Elizabeth Margaret"
## [1] "Butt, Major. Archibald Willingham"
## [1] "Calderhead, Mr. Edward Pennington"
## [1] "Candee, Mrs. Edward (Helen Churchill Hungerford)"
## [1] "Cardeza, Mr. Thomas Drake Martinez"
## [1] "Cardeza, Mrs. James Warburton Martinez (Charlotte Wardle Drake)"
## [1] "Carter, Mr. William Ernest"
## [1] "Carter, Mrs. William Ernest (Lucile Polk)"
## [1] "Case, Mr. Howard Brown"
## [1] "Cavendish, Mr. Tyrell William"
## [1] "Cavendish, Mrs. Tyrell William (Julia Florence Siegel)"
## [1] "Chaffee, Mr. Herbert Fuller"
## [1] "Chaffee, Mrs. Herbert Fuller (Carrie Constance Toogood)"
## [1] "Chaudanson, Miss. Victorine"
## [1] "Chevre, Mr. Paul Romaine"
## [1] "Colley, Mr. Edward Pomeroy"
## [1] "Compton, Miss. Sara Rebecca"
## [1] "Compton, Mr. Alexander Taylor Jr"
## [1] "Compton, Mrs. Alexander Taylor (Mary Eliza Ingersoll)"
## [1] "Cornell, Mrs. Robert Clifford (Malvina Helen Lamson)"
## [1] "Crosby, Capt. Edward Gifford"
## [1] "Crosby, Miss. Harriet R"
## [1] "Crosby, Mrs. Edward Gifford (Catherine Elizabeth Halstead)"
## [1] "Cumings, Mr. John Bradley"
## [1] "Cumings, Mrs. John Bradley (Florence Briggs Thayer)"
## [1] "Daly, Mr. Peter Denis "
## [1] "Dodge, Dr. Washington"
## [1] "Dodge, Mrs. Washington (Ruth Vidaver)"
## [1] "Douglas, Mr. Walter Donald"
## [1] "Douglas, Mrs. Walter Donald (Mahala Dutton)"
## [1] "Duff Gordon, Lady. (Lucille Christiana Sutherland) (\"Mrs Morgan\")"
## [1] "Duff Gordon, Sir. Cosmo Edmund (\"Mr Morgan\")"
## [1] "Dulles, Mr. William Crothers"
## [1] "Endres, Miss. Caroline Louise"
## [1] "Eustis, Miss. Elizabeth Mussey"
## [1] "Evans, Miss. Edith Corse"
## [1] "Flynn, Mr. John Irwin (\"Irving\")"
## [1] "Fortune, Mr. Mark"
## [1] "Fortune, Mrs. Mark (Mary McDougald)"
## [1] "Frauenthal, Dr. Henry William"
## [1] "Frauenthal, Mr. Isaac Gerald"
## [1] "Frolicher-Stehli, Mr. Maxmillian"
## [1] "Frolicher-Stehli, Mrs. Maxmillian (Margaretha Emerentia Stehli)"
## [1] "Futrelle, Mr. Jacques Heath"
## [1] "Futrelle, Mrs. Jacques Heath (Lily May Peel)"
## [1] "Gee, Mr. Arthur H"
## [1] "Geiger, Miss. Amalie"
## [1] "Gibson, Mrs. Leonard (Pauline C Boeson)"
## [1] "Goldenberg, Mr. Samuel L"
## [1] "Goldschmidt, Mr. George B"
## [1] "Gracie, Col. Archibald IV"
## [1] "Graham, Mr. George Edward"
## [1] "Graham, Mrs. William Thompson (Edith Junkins)"
## [1] "Greenfield, Mrs. Leo David (Blanche Strouse)"
## [1] "Guggenheim, Mr. Benjamin"
## [1] "Harper, Mr. Henry Sleeper"
## [1] "Harper, Mrs. Henry Sleeper (Myna Haxtun)"
## [1] "Harris, Mr. Henry Birkhardt"
## [1] "Harris, Mrs. Henry Birkhardt (Irene Wallach)"
## [1] "Harrison, Mr. William"
## [1] "Hays, Mr. Charles Melville"
## [1] "Hays, Mrs. Charles Melville (Clara Jennings Gregg)"
## [1] "Head, Mr. Christopher"
## [1] "Hipkins, Mr. William Edward"
## [1] "Hippach, Mrs. Louis Albert (Ida Sophia Fischer)"
## [1] "Hogeboom, Mrs. John C (Anna Andrews)"
## [1] "Holverson, Mr. Alexander Oskar"
## [1] "Holverson, Mrs. Alexander Oskar (Mary Aline Towner)"
## [1] "Homer, Mr. Harry (\"Mr E Haven\")"
## [1] "Hoyt, Mr. Frederick Maxfield"
## [1] "Hoyt, Mrs. Frederick Maxfield (Jane Anne Forby)"
## [1] "Icard, Miss. Amelie"
## [1] "Isham, Miss. Ann Elizabeth"
## [1] "Ismay, Mr. Joseph Bruce"
## [1] "Jones, Mr. Charles Cresson"
## [1] "Julian, Mr. Henry Forbes"
## [1] "Kent, Mr. Edward Austin"
## [1] "Kenyon, Mr. Frederick R"
## [1] "Kimball, Mr. Edwin Nelson Jr"
## [1] "Kimball, Mrs. Edwin Nelson Jr (Gertrude Parsons)"
## [1] "Kreuchen, Miss. Emilie"
## [1] "Leader, Dr. Alice (Farnham)"
## [1] "Lesurer, Mr. Gustave J"
## [1] "Lindeberg-Lind, Mr. Erik Gustaf (\"Mr Edward Lingrey\")"
## [1] "Lindstrom, Mrs. Carl Johan (Sigrid Posse)"
## [1] "Lines, Mrs. Ernest H (Elizabeth Lindsey James)"
## [1] "Lurette, Miss. Elise"
## [1] "McCaffry, Mr. Thomas Francis"
## [1] "McCarthy, Mr. Timothy J"
## [1] "McGough, Mr. James Robert"
## [1] "Millet, Mr. Francis Davis"
## [1] "Minahan, Dr. William Edward"
## [1] "Minahan, Mrs. William Edward (Lillian E Thorpe)"
## [1] "Molson, Mr. Harry Markland"
## [1] "Moore, Mr. Clarence Bloomfield"
## [1] "Natsch, Mr. Charles H"
## [1] "Newell, Mr. Arthur Webster"
## [1] "Nicholson, Mr. Arthur Ernest"
## [1] "Oliva y Ocana, Dona. Fermina"
## [1] "Ostby, Mr. Engelhart Cornelius"
## [1] "Partner, Mr. Austen"
## [1] "Peuchen, Major. Arthur Godfrey"
## [1] "Porter, Mr. Walter Chamberlain"
## [1] "Potter, Mrs. Thomas Jr (Lily Alexenia Wilson)"
## [1] "Reuchlin, Jonkheer. John George"
## [1] "Robert, Mrs. Edward Scott (Elisabeth Walton McMillan)"
## [1] "Romaine, Mr. Charles Hallace (\"Mr C Rolmane\")"
## [1] "Rosenshine, Mr. George (\"Mr George Thorne\")"
## [1] "Ross, Mr. John Hugo"
## [1] "Rothschild, Mr. Martin"
## [1] "Rothschild, Mrs. Martin (Elizabeth L. Barrett)"
## [1] "Ryerson, Mr. Arthur Larned"
## [1] "Ryerson, Mrs. Arthur Larned (Emily Maria Borie)"
## [1] "Schabert, Mrs. Paul (Emma Mock)"
## [1] "Seward, Mr. Frederic Kimber"
## [1] "Shutes, Miss. Elizabeth W"
## [1] "Silverthorne, Mr. Spencer Victor"
## [1] "Silvey, Mr. William Baird"
## [1] "Silvey, Mrs. William Baird (Alice Munger)"
## [1] "Simonius-Blumer, Col. Oberst Alfons"
## [1] "Smart, Mr. John Montgomery"
## [1] "Smith, Mr. James Clinch"
## [1] "Spedden, Mr. Frederic Oakley"
## [1] "Spedden, Mrs. Frederic Oakley (Margaretta Corning Stone)"
## [1] "Spencer, Mr. William Augustus"
## [1] "Stead, Mr. William Thomas"
## [1] "Stengel, Mr. Charles Emil Henry"
## [1] "Stengel, Mrs. Charles Emil Henry (Annie May Morris)"
## [1] "Stephenson, Mrs. Walter Bertram (Martha Eustis)"
## [1] "Stone, Mrs. George Nelson (Martha Evelyn)"
## [1] "Straus, Mr. Isidor"
## [1] "Straus, Mrs. Isidor (Rosalie Ida Blun)"
## [1] "Sutton, Mr. Frederick"
## [1] "Swift, Mrs. Frederick Joel (Margaret Welles Barron)"
## [1] "Taussig, Mr. Emil"
## [1] "Taussig, Mrs. Emil (Tillie Mandelbaum)"
## [1] "Taylor, Mr. Elmer Zebley"
## [1] "Thayer, Mr. John Borland"
## [1] "Thayer, Mrs. John Borland (Marian Longstreth Morris)"
## [1] "Uruchurtu, Don. Manuel E"
## [1] "Van der hoef, Mr. Wyckoff"
## [1] "Walker, Mr. William Anderson"
## [1] "Ward, Miss. Anna"
## [1] "Warren, Mr. Frank Manley"
## [1] "Warren, Mrs. Frank Manley (Anna Sophia Atkinson)"
## [1] "Weir, Col. John"
## [1] "White, Mr. Percival Wayland"
## [1] "White, Mrs. John Stuart (Ella Holmes)"
## [1] "Wick, Mr. George Dennick"
## [1] "Wick, Mrs. George Dennick (Mary Hitchcock)"
## [1] "Widener, Mr. George Dunton"
## [1] "Widener, Mrs. George Dunton (Eleanor Elkins)"
## [1] "Williams, Mr. Charles Duane"
## [1] "Wright, Mr. George"
## [1] "Young, Miss. Marie Grice"
## [1] "Angle, Mr. William A"
## [1] "Angle, Mrs. William A (Florence \"Mary\" Agnes Hughes)"
## [1] "Ashby, Mr. John"
## [1] "Ball, Mrs. (Ada E Hall)"
## [1] "Bateman, Rev. Robert James"
## [1] "Becker, Mrs. Allen Oliver (Nellie E Baumgardner)"
## [1] "Beesley, Mr. Lawrence"
## [1] "Bowenur, Mr. Solomon"
## [1] "Brown, Mr. Thomas William Solomon"
## [1] "Brown, Mrs. Thomas William Solomon (Elizabeth Catherine Ford)"
## [1] "Buss, Miss. Kate"
## [1] "Byles, Rev. Thomas Roussel Davids"
## [1] "Bystrom, Mrs. (Karolina)"
## [1] "Cameron, Miss. Clear Annie"
## [1] "Carter, Mrs. Ernest Courtenay (Lilian Hughes)"
## [1] "Carter, Rev. Ernest Courtenay"
## [1] "Chapman, Mr. Charles Henry"
## [1] "Chapman, Mr. John Henry"
## [1] "Christy, Mrs. (Alice Frances)"
## [1] "Davies, Mrs. John Morgan (Elizabeth Agnes Mary White) "
## [1] "Doling, Mrs. John T (Ada Julia Bone)"
## [1] "Downton, Mr. William James"
## [1] "Drew, Mr. James Vivian"
## [1] "Drew, Mrs. James Vivian (Lulu Thorne Christian)"
## [1] "Faunthorpe, Mr. Harry"
## [1] "Fox, Mr. Stanley Hubert"
## [1] "Funk, Miss. Annie Clemmer"
## [1] "Fynney, Mr. Joseph J"
## [1] "Gale, Mr. Harry"
## [1] "Gale, Mr. Shadrach"
## [1] "Garside, Miss. Ethel"
## [1] "Gilbert, Mr. William"
## [1] "Gillespie, Mr. William Henry"
## [1] "Greenberg, Mr. Samuel"
## [1] "Harbeck, Mr. William H"
## [1] "Harris, Mr. George"
## [1] "Hart, Mr. Benjamin"
## [1] "Hart, Mrs. Benjamin (Esther Ada Bloomfield)"
## [1] "Herman, Mr. Samuel"
## [1] "Herman, Mrs. Samuel (Jane Laver)"
## [1] "Hewlett, Mrs. (Mary D Kingcome) "
## [1] "Hocking, Mr. Samuel James Metcalfe"
## [1] "Hocking, Mrs. Elizabeth (Eliza Needs)"
## [1] "Hodges, Mr. Henry Price"
## [1] "Hold, Mr. Stephen"
## [1] "Hosono, Mr. Masabumi"
## [1] "Howard, Mr. Benjamin"
## [1] "Howard, Mrs. Benjamin (Ellen Truelove Arman)"
## [1] "Jacobsohn, Mr. Sidney Samuel"
## [1] "Jarvis, Mr. John Denzil"
## [1] "Kantor, Mr. Sinai"
## [1] "Keane, Mr. Daniel"
## [1] "Kelly, Mrs. Florence \"Fannie\""
## [1] "Kirkland, Rev. Charles Leonard"
## [1] "Lemore, Mrs. (Amelia Milley)"
## [1] "Levy, Mr. Rene Jacques"
## [1] "Lingane, Mr. John"
## [1] "Louch, Mr. Charles Alexander"
## [1] "Louch, Mrs. Charles Alexander (Alice Adelaide Slow)"
## [1] "Mack, Mrs. (Mary)"
## [1] "Maybery, Mr. Frank Hubert"
## [1] "McKane, Mr. Peter David"
## [1] "Mellinger, Mrs. (Elizabeth Anne Maidment)"
## [1] "Meyer, Mr. August"
## [1] "Milling, Mr. Jacob Christian"
## [1] "Mitchell, Mr. Henry Michael"
## [1] "Moraweck, Dr. Ernest"
## [1] "Morley, Mr. Henry Samuel (\"Mr Henry Marshall\")"
## [1] "Myles, Mr. Thomas Francis"
## [1] "Navratil, Mr. Michel (\"Louis M Hoffman\")"
## [1] "Otter, Mr. Richard"
## [1] "Parrish, Mrs. (Lutie Davis)"
## [1] "Peruschitz, Rev. Joseph Maria"
## [1] "Phillips, Mr. Escott Robert"
## [1] "Ponesell, Mr. Martin"
## [1] "Reeves, Mr. David"
## [1] "Renouf, Mr. Peter Henry"
## [1] "Ridsdale, Miss. Lucy"
## [1] "Sjostedt, Mr. Ernst Adolf"
## [1] "Slemen, Mr. Richard James"
## [1] "Smith, Miss. Marion Elsie"
## [1] "Stanton, Mr. Samuel Ward"
## [1] "Toomey, Miss. Ellen"
## [1] "Veal, Mr. James"
## [1] "Watt, Mrs. James (Elizabeth \"Bessie\" Inglis Milne)"
## [1] "West, Mr. Edwy Arthur"
## [1] "Wheadon, Mr. Edward H"
## [1] "Abbing, Mr. Anthony"
## [1] "Abbott, Mrs. Stanton (Rosa Hunt)"
## [1] "Ahlin, Mrs. Johan (Johanna Persdotter Larsson)"
## [1] "Allen, Mr. William Henry"
## [1] "Andersson, Miss. Ida Augusta Margareta"
## [1] "Andersson, Mr. Anders Johan"
## [1] "Andersson, Mrs. Anders Johan (Alfrida Konstantia Brogren)"
## [1] "Asim, Mr. Adola"
## [1] "Asplund, Mr. Carl Oscar Vilhelm Gustafsson"
## [1] "Asplund, Mrs. Carl Oscar (Selma Augusta Emilia Johansson)"
## [1] "Assaf Khalil, Mrs. Mariana (\"Miriam\")"
## [1] "Badt, Mr. Mohamed"
## [1] "Barbara, Mrs. (Catherine David)"
## [1] "Bourke, Mr. John"
## [1] "Brocklebank, Mr. William Alfred"
## [1] "Cacic, Mr. Luka"
## [1] "Carr, Miss. Jeannie"
## [1] "Coleff, Mr. Peju"
## [1] "Connors, Mr. Patrick"
## [1] "Cook, Mr. Jacob"
## [1] "Cor, Mr. Bartol"
## [1] "Coutts, Mrs. William (Winnie \"Minnie\" Treanor)"
## [1] "Coxon, Mr. Daniel"
## [1] "Cribb, Mr. John Hatfield"
## [1] "Dahl, Mr. Karl Edwart"
## [1] "Danbom, Mr. Ernst Gilbert"
## [1] "de Messemaeker, Mr. Guillaume Joseph"
## [1] "de Messemaeker, Mrs. Guillaume Joseph (Emma)"
## [1] "Dennis, Mr. William"
## [1] "Dimic, Mr. Jovan"
## [1] "Dintcheff, Mr. Valtcho"
## [1] "Duane, Mr. Frank"
## [1] "Ekstrom, Mr. Johan"
## [1] "Elias, Mr. Joseph"
## [1] "Elsbury, Mr. William James"
## [1] "Everett, Mr. Thomas James"
## [1] "Farrell, Mr. James"
## [1] "Ford, Mrs. Edward (Margaret Ann Watson)"
## [1] "Goldsmith, Mr. Nathan"
## [1] "Goncalves, Mr. Manuel Estanslas"
## [1] "Goodwin, Mr. Charles Frederick"
## [1] "Goodwin, Mrs. Frederick (Augusta Tyler)"
## [1] "Green, Mr. George Henry"
## [1] "Gustafsson, Mr. Anders Vilhelm"
## [1] "Hansen, Mr. Claus Peter"
## [1] "Hansen, Mrs. Claus Peter (Jennie L Howard)"
## [1] "Holm, Mr. John Fredrik Alexander"
## [1] "Humblen, Mr. Adolf Mathias Nicolai Olsen"
## [1] "Jensen, Mr. Niels Peder"
## [1] "Johanson, Mr. Jakob Alfred"
## [1] "Johnson, Mr. Alfred"
## [1] "Karun, Mr. Franz"
## [1] "Kelly, Mr. James"
## [1] "Kelly, Mr. James"
## [1] "Klasen, Mrs. (Hulda Kristina Eugenia Lofqvist)"
## [1] "Laitinen, Miss. Kristina Sofia"
## [1] "Lemberopolous, Mr. Peter L"
## [1] "Leonard, Mr. Lionel"
## [1] "Lester, Mr. James"
## [1] "Lindblom, Miss. Augusta Charlotta"
## [1] "Lindell, Mr. Edvard Bengtsson"
## [1] "Lundahl, Mr. Johan Svensson"
## [1] "Markoff, Mr. Marin"
## [1] "McGowan, Miss. Katherine"
## [1] "Meo, Mr. Alfonzo"
## [1] "Morley, Mr. William"
## [1] "Nirva, Mr. Iisakki Antino Aijo"
## [1] "Niskanen, Mr. Juha"
## [1] "Nysveen, Mr. Johan Hansen"
## [1] "Olsen, Mr. Karl Siegwart Andreas"
## [1] "Panula, Mrs. Juha (Maria Emilia Ojala)"
## [1] "Rekic, Mr. Tido"
## [1] "Rice, Mrs. William (Margaret Norton)"
## [1] "Rintamaki, Mr. Matti"
## [1] "Robins, Mr. Alexander A"
## [1] "Robins, Mrs. Alexander A (Grace Charity Laury)"
## [1] "Rosblom, Mrs. Viktor (Helena Wilhelmina)"
## [1] "Rouse, Mr. Richard Henry"
## [1] "Saether, Mr. Simon Sivertsen"
## [1] "Salonen, Mr. Johan Werner"
## [1] "Sivic, Mr. Husein"
## [1] "Skoog, Mr. Wilhelm"
## [1] "Skoog, Mrs. William (Anna Bernhardina Karlsson)"
## [1] "Storey, Mr. Thomas"
## [1] "Sundman, Mr. Johan Julian"
## [1] "Svensson, Mr. Johan"
## [1] "Theobald, Mr. Thomas Leonard"
## [1] "Torber, Mr. Ernst William"
## [1] "Turcin, Mr. Stjepan"
## [1] "Turkula, Mrs. (Hedwig)"
## [1] "van Billiard, Mr. Austin Blyler"
## [1] "Van Impe, Mr. Jean Baptiste"
## [1] "Vander Cruyssen, Mr. Victor"
## [1] "Whabee, Mrs. George Joseph (Shawneene Abi-Saab)"
## [1] "Widegren, Mr. Carl/Charles Peter"
## [1] "Wilkes, Mrs. James (Ellen Needs)"
## [1] "Wittevrongel, Mr. Camille"
## [1] "Youseff, Mr. Gerious"

Question 2 (15 points)

In a shipment of 20 engines, history shows that the probability of any one engine proving unsatisfactory is 0.1

  1. What will be the probability that all 20 engines are satisfactory? (5 points)
#INSERT YOUR ANSWER HERE
dbinom(x=20,size =20, prob= 0.9)
## [1] 0.1215767
  1. Use the Binomial approximation to calculate the probability that more than 10 engines are defective? (5 points)
#INSERT YOUR ANSWER HERE
dbinom(x=10, size =20, prob=0.1)
## [1] 6.442043e-06
  1. Use the binomial approximation to calculate the probability that at most three engines are defective? (5 points)
#INSERT YOUR ANSWER HERE
dbinom(3,20, 0.1)
## [1] 0.1901199

Question 3 (25 points)

On average, John goes to his parents’ place twice a week for visits.

  1. What will be the probabilities if John decides not to visit at all and then does 1, 2, 3, 4 and 5 visits to his parents in the subsequent weeks?(5 points)
#INSERT YOUR ANSWER HERE
dbinom(x=0,size=6,prob = 1/3)*dbinom(x=5, size=6,prob=1/3)
## [1] 0.001445127
  1. Find the probability that John goes to his parents’ place 3 times or less in a week? (5 points)
#INSERT YOUR ANSWER HERE
sum(dbinom(x = 0:4, size =6,prob = 1/3))
## [1] 0.9821674
  1. Compare the similarity between Binomial and Poisson distribution. (15 points @ 5 point each)
  1. Create 100,000 samples for a Binomial random variable using parameters described in Question 2
  2. Create 100,000 samples for a Poisson random variable using parameters described in Question 3
  3. then illustrate on how well the Poisson probability distribution approximates the Binomial probability distribution. Please provide your result by visualizing through a comparative histogram. HINT: use multhist() from the ‘plotrix’ package. Please install if plotrix package is not in your R environment.
#INSERT YOUR ANSWER HERE
binom_sample <- rbinom(100000,20,0.1)
hist(binom_sample)

#INSERT YOUR ANSWER HERE
poisson_sample <- rpois(100000,lambda=2)
hist(poisson_sample)

#install.packages("plotrix")
library(plotrix)
multhist(list(poisson_sample,binom_sample))

Question 4 (20 points)

Given the following conditions of a normal random variable having mean = 16 and variance = 9, please write the R codes to find the probabilities and their values when,

  1. lies between 14.4 and 20.3 (inclusive) (5 points)
#INSERT YOUR ANSWER HERE
pnorm(20.3, 16, sd=3)-pnorm(14.4,16,sd=3)
## [1] 0.6272173
  1. is greater than 21.8 (5 points)
#INSERT YOUR ANSWER HERE
1-pnorm(21.8,16,sd=3)
## [1] 0.02659757
  1. is less than or equal to 10.5 (5 points)
#INSERT YOUR ANSWER HERE
pnorm(10.5,16, sd=3)
## [1] 0.03337651
  1. is less than 13 or greater than 19 (5 points)
#INSERT YOUR ANSWER HERE
pnorm(13,16,sd=3)+ 1-pnorm(19,16,sd=3)
## [1] 0.3173105

END of Assignment #2.