¿Sabias qué… La última sobreviviente murió en mayo de 2009, a los 97 años?

Importar base de datos

#file.choose

bd<- read.csv("/Users/joseramonvazquezguzman/Documents/Tecnológico de Monterrey/Septimo semestre/Analítica de datos /M3/titaniccsv.csv")

Entender base de datos

summary(bd)
##      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.1667   Min.   :0.0000   Min.   :0.000   Length:1309       
##  1st Qu.:21.0000   1st Qu.:0.0000   1st Qu.:0.000   Class :character  
##  Median :28.0000   Median :0.0000   Median :0.000   Mode  :character  
##  Mean   :29.8811   Mean   :0.4989   Mean   :0.385                     
##  3rd Qu.:39.0000   3rd Qu.:1.0000   3rd Qu.:0.000                     
##  Max.   :80.0000   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

Filtrar base de datos

 Titanic <- bd[,c("pclass", "age", "sex", "survived")]

 Titanic$survived <- as.factor(ifelse(Titanic$survived==0, "Murio", "Sobrevivió"))
 
 Titanic$pclass<- as.factor(Titanic$pclass)
 Titanic$sex<- as.factor(Titanic$sex)
 str(Titanic)
## 'data.frame':    1309 obs. of  4 variables:
##  $ pclass  : Factor w/ 3 levels "1","2","3": 1 1 1 1 1 1 1 1 1 1 ...
##  $ age     : num  29 0.917 2 30 25 ...
##  $ sex     : Factor w/ 2 levels "female","male": 1 2 1 2 1 2 1 2 1 2 ...
##  $ survived: Factor w/ 2 levels "Murio","Sobrevivió": 2 2 1 1 1 2 2 1 2 1 ...
 sum(is.na(Titanic)) 
## [1] 263
 Titanic <- na.omit(Titanic) 

Omitir NA

Titanic <- na.omit(Titanic)

Crear árbol de decisión

library(rpart)
arbol <-rpart(formula=survived~ ., data = Titanic)
arbol 
## n= 1046 
## 
## node), split, n, loss, yval, (yprob)
##       * denotes terminal node
## 
##  1) root 1046 427 Murio (0.59177820 0.40822180)  
##    2) sex=male 658 135 Murio (0.79483283 0.20516717)  
##      4) age>=9.5 615 110 Murio (0.82113821 0.17886179) *
##      5) age< 9.5 43  18 Sobrevivió (0.41860465 0.58139535)  
##       10) pclass=3 29  11 Murio (0.62068966 0.37931034) *
##       11) pclass=1,2 14   0 Sobrevivió (0.00000000 1.00000000) *
##    3) sex=female 388  96 Sobrevivió (0.24742268 0.75257732)  
##      6) pclass=3 152  72 Murio (0.52631579 0.47368421)  
##       12) age>=1.5 145  66 Murio (0.54482759 0.45517241) *
##       13) age< 1.5 7   1 Sobrevivió (0.14285714 0.85714286) *
##      7) pclass=1,2 236  16 Sobrevivió (0.06779661 0.93220339) *

install.packages (rpart.plot)

#install.packages("rpart.plot")
library(rpart.plot)

rpart.plot(arbol)

prp(arbol,extra=7, prefix="fracción\n") 

Conclusiones

Se realizo un árbol de decisión, este nos indica quienes fueron las personas que sobrevivieron o murieron en este suceso; Tenemos para empezar el factor genero, que separa en dos las posibilidades.

Podemos deducir en base al árbol que del género hombre, 18% de éstas personas sobrevivieron, esto si eran mayores a 9.5 años de edad, un 38% murió siendo mayor a esta edad de 9.5, pero un 100 de los niños de clase económica 3 sobevivieron.

Las mujeres pertenecientes a la clase 3 sobrevivieron con un porcentaje del 93%, esto sí eras mayor a 1.5 años de edad.

LS0tCnRpdGxlOiAgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWUiPlRpdGFuaWM8L3NwYW4+IAphdXRob3I6ICJKb3PDqSBSYW3Ds24gQTAxNjI1MjAzIgpkYXRlOiAiMjAtMDktMjAyMmAiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogdHJ1ZQogICAgdG9jX2Zsb2F0OiB0cnVlCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCi0tLQoKPGltZyBzcmM9ICIvVXNlcnMvam9zZXJhbW9udmF6cXVlemd1em1hbi9EZXNrdG9wL1RpdGFuaWMtYmFyY28ucG5nIj4KCiMgICoqwr9TYWJpYXMgcXXDqS4uLiBMYSDDumx0aW1hIHNvYnJldml2aWVudGUgbXVyacOzIGVuIG1heW8gZGUgMjAwOSwgYSBsb3MgOTcgYcOxb3M/KioKCgojIyBJbXBvcnRhciBiYXNlIGRlIGRhdG9zCgpgYGB7cn0KI2ZpbGUuY2hvb3NlCgpiZDwtIHJlYWQuY3N2KCIvVXNlcnMvam9zZXJhbW9udmF6cXVlemd1em1hbi9Eb2N1bWVudHMvVGVjbm9sb8yBZ2ljbyBkZSBNb250ZXJyZXkvU2VwdGltbyBzZW1lc3RyZS9BbmFsacyBdGljYSBkZSBkYXRvcyAvTTMvdGl0YW5pY2Nzdi5jc3YiKQoKYGBgCgojIyBFbnRlbmRlciBiYXNlIGRlIGRhdG9zIApgYGB7cn0Kc3VtbWFyeShiZCkKYGBgCgojIyBGaWx0cmFyIGJhc2UgZGUgZGF0b3MgCgpgYGB7cn0KIFRpdGFuaWMgPC0gYmRbLGMoInBjbGFzcyIsICJhZ2UiLCAic2V4IiwgInN1cnZpdmVkIildCgogVGl0YW5pYyRzdXJ2aXZlZCA8LSBhcy5mYWN0b3IoaWZlbHNlKFRpdGFuaWMkc3Vydml2ZWQ9PTAsICJNdXJpbyIsICJTb2JyZXZpdmnDsyIpKQogCiBUaXRhbmljJHBjbGFzczwtIGFzLmZhY3RvcihUaXRhbmljJHBjbGFzcykKIFRpdGFuaWMkc2V4PC0gYXMuZmFjdG9yKFRpdGFuaWMkc2V4KQogc3RyKFRpdGFuaWMpCgogc3VtKGlzLm5hKFRpdGFuaWMpKSAKCiBUaXRhbmljIDwtIG5hLm9taXQoVGl0YW5pYykgCmBgYAoKIyMgT21pdGlyIE5BCgpgYGB7cn0KVGl0YW5pYyA8LSBuYS5vbWl0KFRpdGFuaWMpCmBgYAoKIyMgIENyZWFyIMOhcmJvbCBkZSBkZWNpc2nDs24gCgpgYGB7cn0KbGlicmFyeShycGFydCkKYXJib2wgPC1ycGFydChmb3JtdWxhPXN1cnZpdmVkfiAuLCBkYXRhID0gVGl0YW5pYykKYXJib2wgCmBgYAoKIyMgIGluc3RhbGwucGFja2FnZXMgKHJwYXJ0LnBsb3QpCmBgYHtyfQojaW5zdGFsbC5wYWNrYWdlcygicnBhcnQucGxvdCIpCmxpYnJhcnkocnBhcnQucGxvdCkKCnJwYXJ0LnBsb3QoYXJib2wpCgpwcnAoYXJib2wsZXh0cmE9NywgcHJlZml4PSJmcmFjY2nDs25cbiIpIApgYGAKCiAKIyAqQ29uY2x1c2lvbmVzKgoKU2UgcmVhbGl6byB1biDDoXJib2wgZGUgZGVjaXNpw7NuLCBlc3RlIG5vcyBpbmRpY2EgcXVpZW5lcyBmdWVyb24gbGFzIHBlcnNvbmFzIHF1ZSBzb2JyZXZpdmllcm9uIG8gbXVyaWVyb24gZW4gZXN0ZSBzdWNlc287IFRlbmVtb3MgcGFyYSBlbXBlemFyIGVsIGZhY3RvciBnZW5lcm8sIHF1ZSBzZXBhcmEgZW4gZG9zIGxhcyBwb3NpYmlsaWRhZGVzLiAKClBvZGVtb3MgZGVkdWNpciBlbiBiYXNlIGFsIMOhcmJvbCBxdWUgZGVsIGfDqW5lcm8gaG9tYnJlLCAxOCUgZGUgw6lzdGFzIHBlcnNvbmFzIHNvYnJldml2aWVyb24sIGVzdG8gc2kgZXJhbiBtYXlvcmVzIGEgOS41IGHDsW9zIGRlIGVkYWQsIHVuIDM4JSBtdXJpw7Mgc2llbmRvIG1heW9yIGEgZXN0YSBlZGFkIGRlIDkuNSwgcGVybyB1biAxMDAgZGUgbG9zIG5pw7FvcyBkZSBjbGFzZSBlY29uw7NtaWNhIDMgc29iZXZpdmllcm9uLiAKCkxhcyBtdWplcmVzIHBlcnRlbmVjaWVudGVzIGEgbGEgY2xhc2UgMyBzb2JyZXZpdmllcm9uIGNvbiB1biBwb3JjZW50YWplIGRlbCA5MyUsIGVzdG8gc8OtIGVyYXMgbWF5b3IgYSAxLjUgYcOxb3MgZGUgZWRhZC4KCgoKCg==