Importar base de datos

base_de_datos <- read.csv("/Users/isaacdiazruizdechavez/Downloads/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,"Murió","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 "Murió","Sobrevivió": 2 2 1 1 1 2 2 1 2 1 ...
sum(is.na(Titanic))
## [1] 266
Titanic<- na.omit(Titanic)

Crear arbol 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 Murió (0.59177820 0.40822180)  
##    2) sex=male 658 135 Murió (0.79483283 0.20516717)  
##      4) age>=9.5 615 110 Murió (0.82113821 0.17886179) *
##      5) age< 9.5 43  18 Sobrevivió (0.41860465 0.58139535)  
##       10) pclass=3 29  11 Murió (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 Murió (0.52631579 0.47368421)  
##       12) age>=1.5 145  66 Murió (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)

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

Conclusiones

En el árbol de decisión podemos ver los siguientes escenarios respecto al análisis de la base de datos de Titanic.

  1. Ser hombre y morir tenía una probabilidad de ocurrencia del 41% y un 100% de los casos que presentaron estas características murieron.

  2. Ser hombre y tener más o igual de 9.5 años representaba un 18% de probabilidad de morir. El 59% de los casos con estas 2 características, murieron.

  3. Ser hombre, tener menos de 9.5 años y ser de tercera clase adquiría el 38% de probabilidad de morir. El 3% de pasajeros con estos 3 atributos murieron.

  4. Ser hombre, tener menos de 9.5 años y NO ser de tercera clase te otorgaba el 100 de probabilidad de sobrevivir. Únicamente el 1% de los pasajeros cumplieron con dicho parámetro.

  5. No ser hombre, es de cir, ser mujer; Sí ser de tercera clase y tener más o 1.5 años te otorgaba el 46% de probabilidad de morir. 14% de las mujeres con estos atributos murieron en el Titanic.

  6. Ser mujer, ser de tercera clase pero NO tener menos de 1.5 años te otorgaba el 86% de probabilidad de sobrevivir. Únicamente el 1% de las pasajeras estuvieron dentro de dicho resultado.

  7. Finalmente, el ser mujer y NO SER de tercera clase te otorgaba el 93% de probabilidades de sobrevivir. Este resultado representa la mayor población de sobrevivientes con un 23% del total de los pasajeros.

El árbol de decisión en este escenario del Titanic nos permitió ver de manera global los diferentes escenarios respecto a los pasajeros y cómo a pesar de las probabilidades que había respecto a la incidencia de cada uno de los escenarios, los porcentajes reales de incidencia muchas veces contradicen los pronósticos.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0gIkNvbG9yOlB1cnBsZSIgPiJUaXRhbmljIiA8L3NwYW4+CmF1dGhvcjogIklzYWFjIETDrWF6IEEwMTU0MDU0MyIKZGF0ZTogIjIwMjItMDktMDYiCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiB5ZXMKICAgIHRvY19mbG9hdDogeWVzCiAgICBjb2RlX2Rvd25sb2FkOiB5ZXMKLS0tCjxpbWcgc3JjPSIvVXNlcnMvaXNhYWNkaWF6cnVpemRlY2hhdmV6L0Rvd25sb2Fkcy84NDBfNTYwLmpwZWciPgoKIyMjIDxzcGFuIHN0eWxlPSAiQ29sb3I6UHVycGxlIiA+IEltcG9ydGFyIGJhc2UgZGUgZGF0b3MgPC9zcGFuPgoKYGBge3J9CgpiYXNlX2RlX2RhdG9zIDwtIHJlYWQuY3N2KCIvVXNlcnMvaXNhYWNkaWF6cnVpemRlY2hhdmV6L0Rvd25sb2Fkcy90aXRhbmljLmNzdiIpCgpgYGAKCiMjIyA8c3BhbiBzdHlsZT0gIkNvbG9yOlB1cnBsZSIgPiBFbnRlbmRlciBiYXNlIGRlIGRhdG9zIDwvc3Bhbj4KCmBgYHtyfQpzdW1tYXJ5KGJhc2VfZGVfZGF0b3MpCmBgYAoKIyMjIDxzcGFuIHN0eWxlPSAiQ29sb3I6UHVycGxlIiA+IEZpbHRyYXIgYmFzZSBkZSBkYXRvcyA8L3NwYW4+CmBgYHtyfQoKVGl0YW5pYzwtIGJhc2VfZGVfZGF0b3NbLGMoInBjbGFzcyIsImFnZSIsInNleCIsInN1cnZpdmVkIildCgpUaXRhbmljJHN1cnZpdmVkPC0gYXMuZmFjdG9yKGlmZWxzZShUaXRhbmljJHN1cnZpdmVkPT0wLCJNdXJpw7MiLCJTb2JyZXZpdmnDsyIpKQpUaXRhbmljJHBjbGFzczwtIGFzLmZhY3RvcihUaXRhbmljJHBjbGFzcykKVGl0YW5pYyRzZXg8LSBhcy5mYWN0b3IoVGl0YW5pYyRzZXgpCnN0cihUaXRhbmljKQoKc3VtKGlzLm5hKFRpdGFuaWMpKQoKVGl0YW5pYzwtIG5hLm9taXQoVGl0YW5pYykKCmBgYAoKIyMjIDxzcGFuIHN0eWxlPSAiQ29sb3I6UHVycGxlIiA+IENyZWFyIGFyYm9sIGRlIGRlY2lzacOzbiA8L3NwYW4+CmBgYHtyfQoKbGlicmFyeShycGFydCkKCmFyYm9sPC0gcnBhcnQoZm9ybXVsYSA9IHN1cnZpdmVkIH4gLiAsIGRhdGE9IFRpdGFuaWMpCmFyYm9sCgojaW5zdGFsbC5wYWNrYWdlcygicnBhcnQucGxvdCIpCmxpYnJhcnkocnBhcnQucGxvdCkKCnJwYXJ0LnBsb3QoYXJib2wpCnBycChhcmJvbCxleHRyYT03LHByZWZpeCA9ICJmcmFjY2nDs25cbiIpCgpgYGAKCiMjIyA8c3BhbiBzdHlsZT0gIkNvbG9yOlB1cnBsZSIgPiBDb25jbHVzaW9uZXMgPC9zcGFuPgoKRW4gZWwgw6FyYm9sIGRlIGRlY2lzacOzbiBwb2RlbW9zIHZlciBsb3Mgc2lndWllbnRlcyBlc2NlbmFyaW9zIHJlc3BlY3RvIGFsIGFuw6FsaXNpcyBkZSBsYSBiYXNlIGRlIGRhdG9zIGRlIFRpdGFuaWMuCgphKSBTZXIgaG9tYnJlIHkgbW9yaXIgdGVuw61hIHVuYSBwcm9iYWJpbGlkYWQgZGUgb2N1cnJlbmNpYSBkZWwgNDElIHkgdW4gMTAwJSBkZSBsb3MgY2Fzb3MgcXVlIHByZXNlbnRhcm9uIGVzdGFzIGNhcmFjdGVyw61zdGljYXMgbXVyaWVyb24uCgpiKSBTZXIgaG9tYnJlIHkgdGVuZXIgbcOhcyBvIGlndWFsIGRlIDkuNSBhw7FvcyByZXByZXNlbnRhYmEgdW4gMTglIGRlIHByb2JhYmlsaWRhZCBkZSBtb3Jpci4gRWwgNTklIGRlIGxvcyBjYXNvcyBjb24gZXN0YXMgMiBjYXJhY3RlcsOtc3RpY2FzLCBtdXJpZXJvbi4gCgpjKSBTZXIgaG9tYnJlLCB0ZW5lciBtZW5vcyBkZSA5LjUgYcOxb3MgeSBzZXIgZGUgdGVyY2VyYSBjbGFzZSBhZHF1aXLDrWEgZWwgMzglIGRlIHByb2JhYmlsaWRhZCBkZSBtb3Jpci4gRWwgMyUgZGUgcGFzYWplcm9zIGNvbiBlc3RvcyAzIGF0cmlidXRvcyBtdXJpZXJvbi4gCgpkKSBTZXIgaG9tYnJlLCB0ZW5lciBtZW5vcyBkZSA5LjUgYcOxb3MgeSBOTyBzZXIgZGUgdGVyY2VyYSBjbGFzZSB0ZSBvdG9yZ2FiYSBlbCAxMDAgZGUgcHJvYmFiaWxpZGFkIGRlIHNvYnJldml2aXIuIMOabmljYW1lbnRlIGVsIDElIGRlIGxvcyBwYXNhamVyb3MgY3VtcGxpZXJvbiBjb24gZGljaG8gcGFyw6FtZXRyby4KCmUpIE5vIHNlciBob21icmUsIGVzIGRlIGNpciwgc2VyIG11amVyOyBTw60gc2VyIGRlIHRlcmNlcmEgY2xhc2UgeSB0ZW5lciBtw6FzIG8gMS41IGHDsW9zIHRlIG90b3JnYWJhIGVsIDQ2JSBkZSBwcm9iYWJpbGlkYWQgZGUgbW9yaXIuIDE0JSBkZSBsYXMgbXVqZXJlcyBjb24gZXN0b3MgYXRyaWJ1dG9zIG11cmllcm9uIGVuIGVsIFRpdGFuaWMuCgpmKSBTZXIgbXVqZXIsIHNlciBkZSB0ZXJjZXJhIGNsYXNlIHBlcm8gTk8gdGVuZXIgbWVub3MgZGUgMS41IGHDsW9zIHRlIG90b3JnYWJhIGVsIDg2JSBkZSBwcm9iYWJpbGlkYWQgZGUgc29icmV2aXZpci4gw5puaWNhbWVudGUgZWwgMSUgZGUgbGFzIHBhc2FqZXJhcyBlc3R1dmllcm9uIGRlbnRybyBkZSBkaWNobyByZXN1bHRhZG8uCgpnKSBGaW5hbG1lbnRlLCBlbCBzZXIgbXVqZXIgeSBOTyBTRVIgZGUgdGVyY2VyYSBjbGFzZSB0ZSBvdG9yZ2FiYSBlbCA5MyUgZGUgcHJvYmFiaWxpZGFkZXMgZGUgc29icmV2aXZpci4gRXN0ZSByZXN1bHRhZG8gcmVwcmVzZW50YSBsYSBtYXlvciBwb2JsYWNpw7NuIGRlIHNvYnJldml2aWVudGVzIGNvbiB1biAyMyUgZGVsIHRvdGFsIGRlIGxvcyBwYXNhamVyb3MuIAoKRWwgw6FyYm9sIGRlIGRlY2lzacOzbiBlbiBlc3RlIGVzY2VuYXJpbyBkZWwgVGl0YW5pYyBub3MgcGVybWl0acOzIHZlciBkZSBtYW5lcmEgZ2xvYmFsIGxvcyBkaWZlcmVudGVzIGVzY2VuYXJpb3MgcmVzcGVjdG8gYSBsb3MgcGFzYWplcm9zIHkgY8OzbW8gYSBwZXNhciBkZSBsYXMgcHJvYmFiaWxpZGFkZXMgcXVlIGhhYsOtYSByZXNwZWN0byBhIGxhIGluY2lkZW5jaWEgZGUgY2FkYSB1bm8gZGUgbG9zIGVzY2VuYXJpb3MsIGxvcyBwb3JjZW50YWplcyByZWFsZXMgZGUgaW5jaWRlbmNpYSBtdWNoYXMgdmVjZXMgY29udHJhZGljZW4gbG9zIHByb27Ds3N0aWNvcy4K