Importar base de datos

base_de_datos <- read.csv("/Users/georginamartinez/Documents/Tec/Séptimo Semestre/Analítica para negocios, de los datos a decisiones/titanic.csv")

Entender base de datos

summary(base_de_datos)
##      pclass         survived         name               sex           
##  Min.   :1.000   Min.   :0.000   Length:1310        Length:1310       
##  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                                        
##  NA's   :1       NA's   :1                                            
##       age              sibsp            parch          ticket         
##  Min.   : 0.1667   Min.   :0.0000   Min.   :0.000   Length:1310       
##  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   :264       NA's   :1        NA's   :1                         
##       fare            cabin             embarked             boat          
##  Min.   :  0.000   Length:1310        Length:1310        Length:1310       
##  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   :2                                                                 
##       body        home.dest        
##  Min.   :  1.0   Length:1310       
##  1st Qu.: 72.0   Class :character  
##  Median :155.0   Mode  :character  
##  Mean   :160.8                     
##  3rd Qu.:256.0                     
##  Max.   :328.0                     
##  NA's   :1189

Filtrar base de datos

Titanic <- base_de_datos[,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':    1310 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/ 3 levels "","female","male": 2 3 2 3 2 3 2 3 2 3 ...
##  $ survived: Factor w/ 2 levels "Murio","Sobrevivió": 2 2 1 1 1 2 2 1 2 1 ...
sum(is.na(Titanic))
## [1] 266
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")
library(rpart.plot)

rpart.plot(arbol)

?as.factor
?str

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

Conclusiones

En este análisis pudimos realizar un árbol de decisión que nos permite conocer quien sobrevivió o quien murió en el Titanic, empezando con el factor de género donde un 18% que era hombre y mayor de 9 años murieron, también murio un 38% de la tripulación de hombres que si eran menores de 9 años, pero de clase baja, solo sobrevivio el 100% de niños menores de 9 años y de clase alta (clase 3). Por el lado de las mujeres, de clase alta sobrevivio un 93% y si tenían más de 1 año sobrevivió un 86%.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0iY29sb3I6R29sZCI+IsOBcmJvbCBkZSBkZWNpc2nDs24gVGl0YW5pYyI8L3NwYW4+IAphdXRob3I6ICJLYXJsYSBHZW9yZ2luYSBNYXJ0w61uZXogR29uesOhbGV6IEEwMDgyNzUwMCIKZGF0ZTogIjIwMjItMDktMDYiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogdHJ1ZQogICAgdG9jX2Zsb2F0OiB0cnVlCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlCi0tLQoKPGltZyBzcmM9ICIvVXNlcnMvZ2VvcmdpbmFtYXJ0aW5lei9Eb2N1bWVudHMvVGVjL1NlzIFwdGltbyBTZW1lc3RyZS9BbmFsacyBdGljYSBwYXJhIG5lZ29jaW9zLCBkZSBsb3MgZGF0b3MgYSBkZWNpc2lvbmVzL1RpdGFuaWMtYmFyY28ucG5nIj4KCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5JbXBvcnRhciBiYXNlIGRlIGRhdG9zPC9zcGFuPiAKYGBge3J9CmJhc2VfZGVfZGF0b3MgPC0gcmVhZC5jc3YoIi9Vc2Vycy9nZW9yZ2luYW1hcnRpbmV6L0RvY3VtZW50cy9UZWMvU2XMgXB0aW1vIFNlbWVzdHJlL0FuYWxpzIF0aWNhIHBhcmEgbmVnb2Npb3MsIGRlIGxvcyBkYXRvcyBhIGRlY2lzaW9uZXMvdGl0YW5pYy5jc3YiKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5FbnRlbmRlciBiYXNlIGRlIGRhdG9zPC9zcGFuPgpgYGB7cn0Kc3VtbWFyeShiYXNlX2RlX2RhdG9zKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5GaWx0cmFyIGJhc2UgZGUgZGF0b3M8L3NwYW4+CmBgYHtyfQpUaXRhbmljIDwtIGJhc2VfZGVfZGF0b3NbLGMoInBjbGFzcyIsICJhZ2UiLCAic2V4IiwgInN1cnZpdmVkIildCgpUaXRhbmljJHN1cnZpdmVkIDwtIGFzLmZhY3RvciAoaWZlbHNlKFRpdGFuaWMkc3Vydml2ZWQ9PTAsICJNdXJpbyIsIlNvYnJldml2acOzIikpClRpdGFuaWMkcGNsYXNzIDwtIGFzLmZhY3RvcihUaXRhbmljJHBjbGFzcykKVGl0YW5pYyRzZXggPC0gYXMuZmFjdG9yKFRpdGFuaWMkc2V4KQpzdHIoVGl0YW5pYykKCnN1bShpcy5uYShUaXRhbmljKSkKVGl0YW5pYyA8LSBuYS5vbWl0KFRpdGFuaWMpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6ZGFya2JsdWUiPkNyZWFyIMOhcmJvbCBkZSBkZWNpc2nDs248L3NwYW4+CmBgYHtyfQpsaWJyYXJ5KHJwYXJ0KQphcmJvbCA8LSBycGFydChmb3JtdWxhID0gc3Vydml2ZWQgfi4sIGRhdGEgPSBUaXRhbmljKQphcmJvbAoKIyBpbnN0YWxsLnBhY2thZ2VzKCJycGFydC5wbG90IikKbGlicmFyeShycGFydC5wbG90KQoKcnBhcnQucGxvdChhcmJvbCkKP2FzLmZhY3Rvcgo/c3RyCgpwcnAoYXJib2wsIGV4dHJhPTcsIHByZWZpeD0iZnJhY2Npw7NuXG4iKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5Db25jbHVzaW9uZXM8L3NwYW4+CkVuIGVzdGUgYW7DoWxpc2lzIHB1ZGltb3MgcmVhbGl6YXIgdW4gPHNwYW4gc3R5bGU9ImNvbG9yOmdyZWVuIj4qKsOhcmJvbCBkZSBkZWNpc2nDs24qKjwvc3Bhbj4gcXVlIG5vcyBwZXJtaXRlIGNvbm9jZXIgcXVpZW4gc29icmV2aXZpw7MgbyBxdWllbiBtdXJpw7MgZW4gZWwgVGl0YW5pYywgZW1wZXphbmRvIGNvbiBlbCBmYWN0b3IgZGUgZ8OpbmVybyBkb25kZSB1biAxOCUgcXVlIGVyYSAqaG9tYnJlKiB5IG1heW9yIGRlIDkgYcOxb3MgbXVyaWVyb24sIHRhbWJpw6luIG11cmlvIHVuIDM4JSBkZSBsYSB0cmlwdWxhY2nDs24gZGUgaG9tYnJlcyBxdWUgc2kgZXJhbiBtZW5vcmVzIGRlIDkgYcOxb3MsIHBlcm8gZGUgY2xhc2UgYmFqYSwgc29sbyBzb2JyZXZpdmlvIGVsIDEwMCUgZGUgbmnDsW9zIG1lbm9yZXMgZGUgOSBhw7FvcyB5IGRlIGNsYXNlIGFsdGEgKGNsYXNlIDMpLiBQb3IgZWwgbGFkbyBkZSBsYXMgKm11amVyZXMqLCBkZSBjbGFzZSBhbHRhIHNvYnJldml2aW8gdW4gOTMlIHkgc2kgdGVuw61hbiBtw6FzIGRlIDEgYcOxbyBzb2JyZXZpdmnDsyB1biA4NiUuCg==