Arbol de decisión - house pricing
Pablo Eduardo Salazar Sánchez A01666311
2026-02-19
Instalar paquetes y llamar librerías
# install.packages("rpart")
library(rpart)
# install.packages("rpart.plot")
library(rpart.plot)Importar la base de datos
#file.choose()
house2<- read.csv("/Users/edu_sssedu/Desktop/Concentración/HousePriceData.csv")Entender la base de datos
summary(house2)## Observation Dist_Taxi Dist_Market Dist_Hospital
## Min. : 1.0 Min. : 146 Min. : 1666 Min. : 3227
## 1st Qu.:237.0 1st Qu.: 6477 1st Qu.: 9367 1st Qu.:11302
## Median :469.0 Median : 8228 Median :11149 Median :13189
## Mean :468.4 Mean : 8235 Mean :11022 Mean :13091
## 3rd Qu.:700.0 3rd Qu.: 9939 3rd Qu.:12675 3rd Qu.:14855
## Max. :932.0 Max. :20662 Max. :20945 Max. :23294
##
## Carpet Builtup Parking City_Category
## Min. : 775 Min. : 932 Length:905 Length:905
## 1st Qu.: 1317 1st Qu.: 1579 Class :character Class :character
## Median : 1478 Median : 1774 Mode :character Mode :character
## Mean : 1511 Mean : 1794
## 3rd Qu.: 1654 3rd Qu.: 1985
## Max. :24300 Max. :12730
## NA's :7
## Rainfall House_Price
## Min. :-110.0 Min. : 1492000
## 1st Qu.: 600.0 1st Qu.: 4623000
## Median : 780.0 Median : 5860000
## Mean : 786.9 Mean : 6083992
## 3rd Qu.: 970.0 3rd Qu.: 7200000
## Max. :1560.0 Max. :150000000
## str(house2)## 'data.frame': 905 obs. of 10 variables:
## $ Observation : int 1 2 3 4 5 6 7 8 9 10 ...
## $ Dist_Taxi : int 9796 8294 11001 8301 10510 6665 13153 5882 7495 8233 ...
## $ Dist_Market : int 5250 8186 14399 11188 12629 5142 11869 9948 11589 7067 ...
## $ Dist_Hospital: int 10703 12694 16991 12289 13921 9972 17811 13315 13370 11400 ...
## $ Carpet : int 1659 1461 1340 1451 1770 1442 1542 1261 1090 1030 ...
## $ Builtup : int 1961 1752 1609 1748 2111 1733 1858 1507 1321 1235 ...
## $ Parking : chr "Open" "Not Provided" "Not Provided" "Covered" ...
## $ City_Category: chr "CAT B" "CAT B" "CAT A" "CAT B" ...
## $ Rainfall : int 530 210 720 620 450 760 1030 1020 680 1130 ...
## $ House_Price : int 6649000 3982000 5401000 5373000 4662000 4526000 7224000 3772000 4631000 4415000 ...head(house2)## Observation Dist_Taxi Dist_Market Dist_Hospital Carpet Builtup Parking
## 1 1 9796 5250 10703 1659 1961 Open
## 2 2 8294 8186 12694 1461 1752 Not Provided
## 3 3 11001 14399 16991 1340 1609 Not Provided
## 4 4 8301 11188 12289 1451 1748 Covered
## 5 5 10510 12629 13921 1770 2111 Not Provided
## 6 6 6665 5142 9972 1442 1733 Open
## City_Category Rainfall House_Price
## 1 CAT B 530 6649000
## 2 CAT B 210 3982000
## 3 CAT A 720 5401000
## 4 CAT B 620 5373000
## 5 CAT B 450 4662000
## 6 CAT B 760 4526000Crear árbol de decisión
house2 <- house2[-348,]
house2$Parking <- as.factor(house2$Parking)
house2$City_Category <- as.factor(house2$City_Category)
str(house2)## 'data.frame': 904 obs. of 10 variables:
## $ Observation : int 1 2 3 4 5 6 7 8 9 10 ...
## $ Dist_Taxi : int 9796 8294 11001 8301 10510 6665 13153 5882 7495 8233 ...
## $ Dist_Market : int 5250 8186 14399 11188 12629 5142 11869 9948 11589 7067 ...
## $ Dist_Hospital: int 10703 12694 16991 12289 13921 9972 17811 13315 13370 11400 ...
## $ Carpet : int 1659 1461 1340 1451 1770 1442 1542 1261 1090 1030 ...
## $ Builtup : int 1961 1752 1609 1748 2111 1733 1858 1507 1321 1235 ...
## $ Parking : Factor w/ 4 levels "Covered","No Parking",..: 4 3 3 1 3 4 2 4 3 4 ...
## $ City_Category: Factor w/ 3 levels "CAT A","CAT B",..: 2 2 1 2 2 2 1 3 2 3 ...
## $ Rainfall : int 530 210 720 620 450 760 1030 1020 680 1130 ...
## $ House_Price : int 6649000 3982000 5401000 5373000 4662000 4526000 7224000 3772000 4631000 4415000 ...arbol_house2 <- rpart(House_Price~., data=house2)
options(scipen=999)
rpart.plot(arbol_house2)
Conclusiones
Si la casa es categoría A el precio sube a 7.4 millones de
dólares.
Si la casa es categoría B y tiene una distancia mayor al sitio de taxis
el precio sube a 6.5 millones de dólares y si está cerca del sitio el
precio se reduce a 5.3 millones de dólares.
Si la casa es categoría C el precio es de 4.5 millones de dólares.
LS0tCnRpdGxlOiAiQXJib2wgZGUgZGVjaXNpw7NuIC0gaG91c2UgcHJpY2luZyIKYXV0aG9yOiAiUGFibG8gRWR1YXJkbyBTYWxhemFyIFPDoW5jaGV6IEEwMTY2NjMxMSIKZGF0ZTogIjIwMjYtMDItMTkiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDogCiAgICB0b2M6IFRSVUUKICAgIHRvY19mbG9hdDogVFJVRQogICAgY29kZV9kb3dubG9hZDogVFJVRQogICAgdGhlbWU6IGNvc21vCi0tLQoKIVtdKGh0dHBzOi8vbWVkaWEuaXN0b2NrcGhvdG8uY29tL2lkLzk2MTY0NDg5Ni9lcy9mb3RvL2VtYnJ1amFkYS1ob3VzZS5qcGc/cz02MTJ4NjEyJnc9MCZrPTIwJmM9TFNIZzRIMkFDSWJPUF9hNGhmcnh4Q2pkUmxBZlVfSFA3TEt6MHZOU0tHZz0pCgojIDxzcGFuIHN0eWxlPSJjb2xvcjpibHVlIj4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hcyA8L3NwYW4+CmBgYHtyIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0V9CiMgaW5zdGFsbC5wYWNrYWdlcygicnBhcnQiKQpsaWJyYXJ5KHJwYXJ0KQojIGluc3RhbGwucGFja2FnZXMoInJwYXJ0LnBsb3QiKQpsaWJyYXJ5KHJwYXJ0LnBsb3QpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZSI+IEltcG9ydGFyIGxhIGJhc2UgZGUgZGF0b3MgPC9zcGFuPgpgYGB7cn0KI2ZpbGUuY2hvb3NlKCkKaG91c2UyPC0gcmVhZC5jc3YoIi9Vc2Vycy9lZHVfc3NzZWR1L0Rlc2t0b3AvQ29uY2VudHJhY2lvzIFuL0hvdXNlUHJpY2VEYXRhLmNzdiIpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZSI+IEVudGVuZGVyIGxhIGJhc2UgZGUgZGF0b3MgPC9zcGFuPgpgYGB7cn0Kc3VtbWFyeShob3VzZTIpCnN0cihob3VzZTIpCmhlYWQoaG91c2UyKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmJsdWUiPiBDcmVhciDDoXJib2wgZGUgZGVjaXNpw7NuIDwvc3Bhbj4KYGBge3J9CmhvdXNlMiA8LSBob3VzZTJbLTM0OCxdCmhvdXNlMiRQYXJraW5nIDwtIGFzLmZhY3Rvcihob3VzZTIkUGFya2luZykKaG91c2UyJENpdHlfQ2F0ZWdvcnkgPC0gYXMuZmFjdG9yKGhvdXNlMiRDaXR5X0NhdGVnb3J5KQpzdHIoaG91c2UyKQphcmJvbF9ob3VzZTIgPC0gcnBhcnQoSG91c2VfUHJpY2V+LiwgZGF0YT1ob3VzZTIpCm9wdGlvbnMoc2NpcGVuPTk5OSkKcnBhcnQucGxvdChhcmJvbF9ob3VzZTIpCmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZSI+IENvbmNsdXNpb25lcyA8L3NwYW4+CgpTaSBsYSBjYXNhIGVzIGNhdGVnb3LDrWEgQSBlbCBwcmVjaW8gc3ViZSBhIDcuNCBtaWxsb25lcyBkZSBkw7NsYXJlcy4gIApTaSBsYSBjYXNhIGVzIGNhdGVnb3LDrWEgQiB5IHRpZW5lIHVuYSBkaXN0YW5jaWEgbWF5b3IgYWwgc2l0aW8gZGUgdGF4aXMgZWwgcHJlY2lvIHN1YmUgYSA2LjUgbWlsbG9uZXMgZGUgZMOzbGFyZXMgeSBzaSBlc3TDoSBjZXJjYSBkZWwgc2l0aW8gZWwgcHJlY2lvIHNlIHJlZHVjZSBhIDUuMyBtaWxsb25lcyBkZSBkw7NsYXJlcy4gIApTaSBsYSBjYXNhIGVzIGNhdGVnb3LDrWEgQyBlbCBwcmVjaW8gZXMgZGUgNC41IG1pbGxvbmVzIGRlIGTDs2xhcmVzLiAg