TRABAJO PRACTICO Nº1

Informe: Poblacion por comunas en de la Ciudad de Bs. As.

library (tidyverse)
## -- Attaching packages --------------------------------------- tidyverse 1.3.1 --
## v ggplot2 3.3.5     v purrr   0.3.4
## v tibble  3.1.4     v dplyr   1.0.7
## v tidyr   1.1.3     v stringr 1.4.0
## v readr   2.0.1     v forcats 0.5.1
## -- Conflicts ------------------------------------------ tidyverse_conflicts() --
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library (sf)
## Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1

CARGAR DATASET GEOJSON

test2 <- st_read("https://cdn.buenosaires.gob.ar/datosabiertos/datasets/informacion-censal-por-radio/caba_radios_censales.geojson", stringsAsFactors = TRUE)
## Reading layer `informacion_censal_por_radio_wgs84' from data source 
##   `https://cdn.buenosaires.gob.ar/datosabiertos/datasets/informacion-censal-por-radio/caba_radios_censales.geojson' 
##   using driver `GeoJSON'
## Simple feature collection with 3554 features and 15 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -58.53152 ymin: -34.70529 xmax: -58.33514 ymax: -34.52755
## Geodetic CRS:  WGS 84
  test2 <- test2 %>%
  select(-WKT)

CARGAR DATASET CSV.

censo <- read.csv("informacion-censal-por-radio-2010.csv", stringsAsFactors = TRUE)
  censo <- censo %>%
  select(-WKT)

ENCABEZADOS

head(censo)
##   ID CO_FRAC_RA COMUNA FRACCION RADIO TOTAL_POB T_VARON T_MUJER T_VIVIENDA
## 1  1      1_1_1      1        1     1       336     212     124         82
## 2  2     1_12_1      1       12     1       341     184     157        365
## 3  3    1_12_10      1       12    10       296     162     134        629
## 4  4    1_12_11      1       12    11       528     294     234        375
## 5  5     1_12_2      1       12     2       229     101     128        445
## 6  6     1_12_3      1       12     3       723     374     349        744
##   V_PARTICUL V_COLECTIV T_HOGAR H_CON_NBI H_SIN_NBI
## 1         80          2      65        19        46
## 2        361          4     116        25        91
## 3        627          2     101         1       100
## 4        370          5     136         7       129
## 5        445          0     129        16       113
## 6        739          5     314       104       210
head(test2)
## Simple feature collection with 6 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -58.38652 ymin: -34.60801 xmax: -58.35114 ymax: -34.57821
## Geodetic CRS:  WGS 84
##   ID CO_FRAC_RA COMUNA FRACCION RADIO TOTAL_POB T_VARON T_MUJER T_VIVIENDA
## 1  1      1_1_1      1        1     1       336     212     124         82
## 2  2     1_12_1      1       12     1       341     184     157        365
## 3  3    1_12_10      1       12    10       296     162     134        629
## 4  4    1_12_11      1       12    11       528     294     234        375
## 5  5     1_12_2      1       12     2       229     101     128        445
## 6  6     1_12_3      1       12     3       723     374     349        744
##   V_PARTICUL V_COLECTIV T_HOGAR H_CON_NBI H_SIN_NBI
## 1         80          2      65        19        46
## 2        361          4     116        25        91
## 3        627          2     101         1       100
## 4        370          5     136         7       129
## 5        445          0     129        16       113
## 6        739          5     314       104       210
##                         geometry
## 1 MULTIPOLYGON (((-58.37249 -...
## 2 MULTIPOLYGON (((-58.38652 -...
## 3 MULTIPOLYGON (((-58.37938 -...
## 4 MULTIPOLYGON (((-58.36792 -...
## 5 MULTIPOLYGON (((-58.38513 -...
## 6 MULTIPOLYGON (((-58.38213 -...

RESUMEN ESTADISTICO

summary(censo)
##        ID           CO_FRAC_RA       COMUNA          FRACCION    
##  Min.   :   1.0   1_1_1  :   1   Min.   : 1.000   Min.   : 1.00  
##  1st Qu.: 889.2   1_10_1 :   1   1st Qu.: 4.000   1st Qu.: 6.00  
##  Median :1777.5   1_10_10:   1   Median : 8.000   Median :12.00  
##  Mean   :1777.5   1_10_11:   1   Mean   : 8.008   Mean   :12.38  
##  3rd Qu.:2665.8   1_10_12:   1   3rd Qu.:12.000   3rd Qu.:18.00  
##  Max.   :3554.0   1_10_13:   1   Max.   :15.000   Max.   :31.00  
##                   (Other):3548                                   
##      RADIO          TOTAL_POB         T_VARON          T_MUJER      
##  Min.   : 1.000   Min.   :   0.0   Min.   :   0.0   Min.   :   0.0  
##  1st Qu.: 3.000   1st Qu.: 646.2   1st Qu.: 294.0   1st Qu.: 350.2  
##  Median : 6.000   Median : 786.0   Median : 360.0   Median : 425.0  
##  Mean   : 5.839   Mean   : 813.2   Mean   : 374.1   Mean   : 439.1  
##  3rd Qu.: 8.000   3rd Qu.: 928.0   3rd Qu.: 428.0   3rd Qu.: 504.0  
##  Max.   :23.000   Max.   :3945.0   Max.   :2012.0   Max.   :1933.0  
##                                                                     
##    T_VIVIENDA       V_PARTICUL       V_COLECTIV         T_HOGAR      
##  Min.   :   0.0   Min.   :   0.0   Min.   : 0.0000   Min.   :   0.0  
##  1st Qu.: 311.2   1st Qu.: 311.0   1st Qu.: 0.0000   1st Qu.: 259.0  
##  Median : 377.0   Median : 377.0   Median : 0.0000   Median : 310.0  
##  Mean   : 401.4   Mean   : 400.8   Mean   : 0.5253   Mean   : 323.6  
##  3rd Qu.: 462.0   3rd Qu.: 461.0   3rd Qu.: 1.0000   3rd Qu.: 371.0  
##  Max.   :1405.0   Max.   :1405.0   Max.   :14.0000   Max.   :1093.0  
##                                                                      
##    H_CON_NBI        H_SIN_NBI     
##  Min.   :  0.00   Min.   :   0.0  
##  1st Qu.:  2.00   1st Qu.: 242.0  
##  Median :  6.00   Median : 295.0  
##  Mean   : 19.35   Mean   : 304.3  
##  3rd Qu.: 23.00   3rd Qu.: 352.0  
##  Max.   :403.00   Max.   :1052.0  
## 
summary(test2)
##        ID         CO_FRAC_RA       COMUNA        FRACCION        RADIO     
##  1      :   1   1_1_1  :   1   1      : 330   3      : 175   1      : 351  
##  10     :   1   1_10_1 :   1   13     : 306   7      : 165   2      : 349  
##  100    :   1   1_10_10:   1   14     : 294   2      : 160   3      : 347  
##  1000   :   1   1_10_11:   1   3      : 254   4      : 159   4      : 346  
##  1001   :   1   1_10_12:   1   4      : 252   9      : 159   5      : 345  
##  1002   :   1   1_10_13:   1   7      : 250   13     : 157   6      : 342  
##  (Other):3548   (Other):3548   (Other):1868   (Other):2579   (Other):1474  
##    TOTAL_POB       T_VARON        T_MUJER       T_VIVIENDA     V_PARTICUL  
##  726    :  16   371    :  27   419    :  24   324    :  22   324    :  22  
##  720    :  13   377    :  24   442    :  23   316    :  19   333    :  22  
##  735    :  13   324    :  22   426    :  21   341    :  19   361    :  22  
##  753    :  12   356    :  21   402    :  19   347    :  19   329    :  20  
##  801    :  12   331    :  20   406    :  19   309    :  18   341    :  20  
##  822    :  12   349    :  20   418    :  19   312    :  18   359    :  20  
##  (Other):3476   (Other):3420   (Other):3429   (Other):3439   (Other):3428  
##    V_COLECTIV      T_HOGAR       H_CON_NBI      H_SIN_NBI   
##  0      :2435   297    :  29   1      : 418   290    :  27  
##  1      : 734   267    :  25   2      : 408   281    :  26  
##  2      : 228   290    :  25   0      : 295   309    :  26  
##  3      :  80   320    :  25   3      : 269   244    :  24  
##  4      :  33   249    :  24   4      : 199   269    :  24  
##  5      :  13   293    :  24   5      : 161   279    :  23  
##  (Other):  31   (Other):3402   (Other):1804   (Other):3404  
##           geometry   
##  MULTIPOLYGON :3554  
##  epsg:4326    :   0  
##  +proj=long...:   0  
##                      
##                      
##                      
## 

FILTRAR Y AGRUPAR

censo.agrupado <- censo %>%
  group_by(COMUNA) %>%
summarise(HOG.TOTAL = sum(T_HOGAR), POB.TOTAL = sum(TOTAL_POB), H_NBI = sum(H_CON_NBI))

GRAFICOS

ggplot()+
  geom_bar(data= censo.agrupado, aes (x= COMUNA, weight=POB.TOTAL))+
  labs(title = "Analisis censo 2010")+
  scale_x_continuous(breaks = seq(1,15,1))

ggplot(censo.agrupado)+
  geom_bar(aes(x=reorder(COMUNA, POB.TOTAL), weight=POB.TOTAL, fill=H_NBI))+
  labs(title="Censo 2010:Poblacion por comuna, con indice de NBI ",
       fill="NBI",
       x="Comunas",
       y="Poblacion",
      caption="Fuente:https://data.buenosaires.gob.ar/dataset/informacion-censal-por-radio")

MAPEO

ggplot () +
  geom_sf (data = test2, aes (fill = COMUNA))