1 Variable CENSO

Necesitamos calcular las frecuencias a nivel censal de las respuestas correspondientes a la categoría: “Parquet, piso flotante, cerámico, madera, alfombra, flexit, cubrepiso u otro similar, sobre radier o vigas de madera” del campo P03C del Censo de viviendas. Recordemos que ésta fué la más alta correlación en relación a los ingresos expandidos (ver punto 1.2 aquí).

1.1 Lectura y filtrado de la tabla censal de viviendas

Leemos la tabla Censo 2017 de viviendas que ya tiene integrada la clave zonal:

tabla_con_clave <- readRDS("../censo_viviendas_con_clave_17.rds")
r3_100 <- tabla_con_clave[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
REGION PROVINCIA COMUNA DC AREA ZC_LOC ID_ZONA_LOC NVIV P01 P02 P03A P03B P03C P04 P05 CANT_HOG CANT_PER REGION_15R PROVINCIA_15R COMUNA_15R clave
15 152 15202 1 2 6 13225 1 3 1 5 3 5 1 4 1 1 15 152 15202 15202012006
15 152 15202 1 2 6 13225 2 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 3 1 1 5 3 5 2 3 1 4 15 152 15202 15202012006
15 152 15202 1 2 6 13225 4 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 5 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 6 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 7 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 8 3 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 9 3 1 5 3 5 1 4 1 4 15 152 15202 15202012006
15 152 15202 1 2 6 13225 10 1 1 5 3 4 1 4 1 1 15 152 15202 15202012006
15 152 15202 1 2 6 13225 11 1 2 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 12 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 13 1 1 5 3 4 1 4 1 3 15 152 15202 15202012006
15 152 15202 1 2 6 13225 14 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 15 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 16 1 1 5 3 5 1 3 1 3 15 152 15202 15202012006
15 152 15202 1 2 6 13225 17 1 1 5 3 5 2 4 1 8 15 152 15202 15202012006
15 152 15202 1 2 6 13225 18 3 1 5 3 5 1 1 1 1 15 152 15202 15202012006
15 152 15202 1 2 6 13225 19 3 1 5 3 5 1 3 1 1 15 152 15202 15202012006
15 152 15202 1 2 6 13225 20 1 1 5 3 5 2 1 1 2 15 152 15202 15202012006
15 152 15202 1 2 6 13225 21 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 22 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 23 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 24 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 25 1 1 5 3 5 1 4 1 2 15 152 15202 15202012006
15 152 15202 1 2 6 13225 26 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 27 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 28 3 1 5 3 5 1 4 1 5 15 152 15202 15202012006
15 152 15202 1 2 6 13225 29 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 30 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 31 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 32 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 33 1 1 5 3 5 3 4 1 4 15 152 15202 15202012006
15 152 15202 1 2 6 13225 34 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 35 3 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 36 1 1 5 3 5 3 2 1 9 15 152 15202 15202012006
15 152 15202 1 2 6 13225 37 3 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 38 1 1 5 3 5 99 4 1 1 15 152 15202 15202012006
15 152 15202 1 2 6 13225 39 1 1 5 3 5 1 4 1 2 15 152 15202 15202012006
15 152 15202 1 2 6 13225 40 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 41 1 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 42 3 3 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 43 3 1 5 3 5 2 1 1 5 15 152 15202 15202012006
15 152 15202 1 2 6 13225 44 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 6 13225 45 1 4 98 98 98 98 98 0 0 15 152 15202 15202012006
15 152 15202 1 2 8 13910 1 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 2 1 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 3 1 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 4 3 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 5 3 1 5 3 5 2 3 1 3 15 152 15202 15202012008
15 152 15202 1 2 8 13910 6 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 7 1 1 5 99 5 2 4 1 2 15 152 15202 15202012008
15 152 15202 1 2 8 13910 8 3 1 5 3 5 3 3 1 2 15 152 15202 15202012008
15 152 15202 1 2 8 13910 9 3 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 10 3 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 11 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 12 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 13 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 14 3 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 15 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 16 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 17 3 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 18 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 19 3 1 99 99 99 99 99 1 1 15 152 15202 15202012008
15 152 15202 1 2 8 13910 20 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 21 1 1 5 3 5 1 4 1 2 15 152 15202 15202012008
15 152 15202 1 2 8 13910 22 3 1 5 3 5 2 4 1 1 15 152 15202 15202012008
15 152 15202 1 2 8 13910 23 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 24 3 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 25 3 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 26 1 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 27 1 4 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 28 3 2 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 29 1 3 98 98 98 98 98 0 0 15 152 15202 15202012008
15 152 15202 1 2 8 13910 30 1 1 5 1 4 2 4 1 1 15 152 15202 15202012008
15 152 15202 1 2 12 8394 1 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 2 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 3 3 1 5 3 5 2 3 1 4 15 152 15202 15202012012
15 152 15202 1 2 12 8394 4 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 5 3 3 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 6 3 3 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 7 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 8 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 9 3 1 5 3 5 1 4 1 1 15 152 15202 15202012012
15 152 15202 1 2 12 8394 10 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 11 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 12 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 13 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 14 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 15 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 16 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 17 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 18 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 19 3 1 99 99 99 99 99 1 1 15 152 15202 15202012012
15 152 15202 1 2 12 8394 20 3 1 5 3 5 3 99 1 1 15 152 15202 15202012012
15 152 15202 1 2 12 8394 21 3 1 5 99 5 1 4 1 2 15 152 15202 15202012012
15 152 15202 1 2 12 8394 22 3 2 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 23 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012
15 152 15202 1 2 12 8394 24 3 1 5 3 5 1 2 1 2 15 152 15202 15202012012
15 152 15202 1 2 12 8394 25 3 4 98 98 98 98 98 0 0 15 152 15202 15202012012

Despleguemos los códigos de regiones de nuestra tabla:

regiones <- unique(tabla_con_clave$REGION)

Hagamos un subset con la 1:

tabla_con_clave <- filter(tabla_con_clave, tabla_con_clave$REGION == 1) 
tabla_con_clave <- filter(tabla_con_clave, tabla_con_clave$AREA== 1)

1.2 Cálculo de frecuencias

tabla_con_clave_f <- tabla_con_clave[,-c(1,2,4,5,6,7,8,9,10,11,12,14,15,16,17,18,19,20),drop=FALSE]

aterial de construcción del piso

names(tabla_con_clave_f)[2] <- "Tipo de piso"
# Ahora filtramos por Tipo de techo = 1.
tabla_con_clave_ff <- filter(tabla_con_clave_f, tabla_con_clave_f$`Tipo de piso` == 1)
# Determinamos las frecuencias por zona:
b <- tabla_con_clave_ff$clave
c <- tabla_con_clave_ff$`Tipo de piso`
d <- tabla_con_clave_ff$COMUNA
cross_tab =  xtabs( ~ unlist(b) + unlist(c)+ unlist(d))
tabla <- as.data.frame(cross_tab)
d <-tabla[!(tabla$Freq == 0),]
names(d)[1] <- "zona" 
d$anio <- "2017"

Veamos los primeros 100 registros:

r3_100 <- d[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona unlist.c. unlist.d. Freq anio
1 1101011001 1 1101 807 2017
2 1101011002 1 1101 452 2017
3 1101021001 1 1101 257 2017
4 1101021002 1 1101 11 2017
5 1101021003 1 1101 571 2017
6 1101021004 1 1101 451 2017
7 1101021005 1 1101 1046 2017
8 1101031001 1 1101 653 2017
9 1101031002 1 1101 837 2017
10 1101031003 1 1101 1315 2017
11 1101031004 1 1101 903 2017
12 1101041001 1 1101 417 2017
13 1101041002 1 1101 524 2017
14 1101041003 1 1101 1036 2017
15 1101041004 1 1101 1406 2017
16 1101041005 1 1101 894 2017
17 1101041006 1 1101 588 2017
18 1101051001 1 1101 857 2017
19 1101051002 1 1101 1279 2017
20 1101051003 1 1101 1252 2017
21 1101051004 1 1101 558 2017
22 1101051005 1 1101 1207 2017
23 1101051006 1 1101 858 2017
24 1101061001 1 1101 606 2017
25 1101061002 1 1101 1338 2017
26 1101061003 1 1101 1436 2017
27 1101061004 1 1101 1294 2017
28 1101061005 1 1101 626 2017
29 1101071001 1 1101 691 2017
30 1101071002 1 1101 753 2017
31 1101071003 1 1101 1050 2017
32 1101071004 1 1101 605 2017
33 1101081001 1 1101 1380 2017
34 1101081002 1 1101 852 2017
35 1101081003 1 1101 658 2017
36 1101081004 1 1101 762 2017
37 1101101001 1 1101 594 2017
38 1101101002 1 1101 993 2017
39 1101101003 1 1101 1052 2017
40 1101101004 1 1101 838 2017
41 1101101005 1 1101 1215 2017
42 1101101006 1 1101 940 2017
43 1101111001 1 1101 869 2017
44 1101111002 1 1101 578 2017
45 1101111003 1 1101 1242 2017
46 1101111004 1 1101 1078 2017
47 1101111005 1 1101 1235 2017
48 1101111006 1 1101 1039 2017
49 1101111007 1 1101 827 2017
50 1101111008 1 1101 1352 2017
51 1101111009 1 1101 1344 2017
52 1101111010 1 1101 97 2017
53 1101111011 1 1101 1327 2017
54 1101111012 1 1101 820 2017
55 1101111013 1 1101 1153 2017
56 1101111014 1 1101 842 2017
57 1101991999 1 1101 194 2017
144 1107011001 1 1107 843 2017
145 1107011002 1 1107 884 2017
146 1107011003 1 1107 1149 2017
147 1107021001 1 1107 1132 2017
148 1107021002 1 1107 803 2017
149 1107021003 1 1107 1702 2017
150 1107021004 1 1107 1035 2017
151 1107021005 1 1107 927 2017
152 1107021006 1 1107 861 2017
153 1107021007 1 1107 1178 2017
154 1107021008 1 1107 1037 2017
155 1107031001 1 1107 1256 2017
156 1107031002 1 1107 1329 2017
157 1107031003 1 1107 1185 2017
158 1107041001 1 1107 719 2017
159 1107041002 1 1107 1145 2017
160 1107041003 1 1107 606 2017
161 1107041004 1 1107 1149 2017
162 1107041005 1 1107 762 2017
163 1107041006 1 1107 734 2017
164 1107041007 1 1107 848 2017
165 1107991999 1 1107 162 2017
252 1401011001 1 1401 651 2017
253 1401011002 1 1401 1385 2017
254 1401991999 1 1401 21 2017
341 1404011001 1 1404 224 2017
342 1404991999 1 1404 12 2017
429 1405011001 1 1405 895 2017
430 1405991999 1 1405 14 2017
NA NA NA NA NA NA
NA.1 NA NA NA NA NA
NA.2 NA NA NA NA NA
NA.3 NA NA NA NA NA
NA.4 NA NA NA NA NA
NA.5 NA NA NA NA NA
NA.6 NA NA NA NA NA
NA.7 NA NA NA NA NA
NA.8 NA NA NA NA NA
NA.9 NA NA NA NA NA
NA.10 NA NA NA NA NA
NA.11 NA NA NA NA NA
NA.12 NA NA NA NA NA
NA.13 NA NA NA NA NA

Agregamos un cero a los códigos comunales de cuatro dígitos:

codigos <- d$unlist.d.
rango <- seq(1:nrow(d))
cadena <- paste("0",codigos[rango], sep = "")
cadena <- substr(cadena,(nchar(cadena)[rango])-(4),6)
codigos <- as.data.frame(codigos)
cadena <- as.data.frame(cadena)
comuna_corr <- cbind(d,cadena)
comuna_corr <- comuna_corr[,-c(2,3),drop=FALSE]
names(comuna_corr)[4] <- "código" 
r3_100 <- comuna_corr[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona Freq anio código
1 1101011001 807 2017 01101
2 1101011002 452 2017 01101
3 1101021001 257 2017 01101
4 1101021002 11 2017 01101
5 1101021003 571 2017 01101
6 1101021004 451 2017 01101
7 1101021005 1046 2017 01101
8 1101031001 653 2017 01101
9 1101031002 837 2017 01101
10 1101031003 1315 2017 01101
11 1101031004 903 2017 01101
12 1101041001 417 2017 01101
13 1101041002 524 2017 01101
14 1101041003 1036 2017 01101
15 1101041004 1406 2017 01101
16 1101041005 894 2017 01101
17 1101041006 588 2017 01101
18 1101051001 857 2017 01101
19 1101051002 1279 2017 01101
20 1101051003 1252 2017 01101
21 1101051004 558 2017 01101
22 1101051005 1207 2017 01101
23 1101051006 858 2017 01101
24 1101061001 606 2017 01101
25 1101061002 1338 2017 01101
26 1101061003 1436 2017 01101
27 1101061004 1294 2017 01101
28 1101061005 626 2017 01101
29 1101071001 691 2017 01101
30 1101071002 753 2017 01101
31 1101071003 1050 2017 01101
32 1101071004 605 2017 01101
33 1101081001 1380 2017 01101
34 1101081002 852 2017 01101
35 1101081003 658 2017 01101
36 1101081004 762 2017 01101
37 1101101001 594 2017 01101
38 1101101002 993 2017 01101
39 1101101003 1052 2017 01101
40 1101101004 838 2017 01101
41 1101101005 1215 2017 01101
42 1101101006 940 2017 01101
43 1101111001 869 2017 01101
44 1101111002 578 2017 01101
45 1101111003 1242 2017 01101
46 1101111004 1078 2017 01101
47 1101111005 1235 2017 01101
48 1101111006 1039 2017 01101
49 1101111007 827 2017 01101
50 1101111008 1352 2017 01101
51 1101111009 1344 2017 01101
52 1101111010 97 2017 01101
53 1101111011 1327 2017 01101
54 1101111012 820 2017 01101
55 1101111013 1153 2017 01101
56 1101111014 842 2017 01101
57 1101991999 194 2017 01101
144 1107011001 843 2017 01107
145 1107011002 884 2017 01107
146 1107011003 1149 2017 01107
147 1107021001 1132 2017 01107
148 1107021002 803 2017 01107
149 1107021003 1702 2017 01107
150 1107021004 1035 2017 01107
151 1107021005 927 2017 01107
152 1107021006 861 2017 01107
153 1107021007 1178 2017 01107
154 1107021008 1037 2017 01107
155 1107031001 1256 2017 01107
156 1107031002 1329 2017 01107
157 1107031003 1185 2017 01107
158 1107041001 719 2017 01107
159 1107041002 1145 2017 01107
160 1107041003 606 2017 01107
161 1107041004 1149 2017 01107
162 1107041005 762 2017 01107
163 1107041006 734 2017 01107
164 1107041007 848 2017 01107
165 1107991999 162 2017 01107
252 1401011001 651 2017 01401
253 1401011002 1385 2017 01401
254 1401991999 21 2017 01401
341 1404011001 224 2017 01404
342 1404991999 12 2017 01404
429 1405011001 895 2017 01405
430 1405991999 14 2017 01405
NA NA NA NA NA
NA.1 NA NA NA NA
NA.2 NA NA NA NA
NA.3 NA NA NA NA
NA.4 NA NA NA NA
NA.5 NA NA NA NA
NA.6 NA NA NA NA
NA.7 NA NA NA NA
NA.8 NA NA NA NA
NA.9 NA NA NA NA
NA.10 NA NA NA NA
NA.11 NA NA NA NA
NA.12 NA NA NA NA
NA.13 NA NA NA NA


2 Variable CASEN

2.1 Tabla de ingresos expandidos

Hemos calculado ya éste valor como conclusión del punto 1.1 de aquí

h_y_m_2017_censo <- readRDS("../ingresos_expandidos_urbano_17.rds")
tablamadre <- head(h_y_m_2017_censo,50)
kbl(tablamadre) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
código comuna.x promedio_i año comuna.y personas Ingresos_expandidos
01101 Iquique 375676.9 2017 1101 191468 71930106513
01107 Alto Hospicio 311571.7 2017 1107 108375 33766585496
01401 Pozo Almonte 316138.5 2017 1401 15711 4966851883
01405 Pica 330061.1 2017 1405 9296 3068247619
02101 Antofagasta 368221.4 2017 2101 361873 133249367039
02102 Mejillones 369770.7 2017 2102 13467 4979702302
02104 Taltal 383666.2 2017 2104 13317 5109282942
02201 Calama 434325.1 2017 2201 165731 71981127235
02203 San Pedro de Atacama 442861.0 2017 2203 10996 4869699464
02301 Tocopilla 286187.2 2017 2301 25186 7207910819
02302 María Elena 477748.0 2017 2302 6457 3084818966
03101 Copiapó 343121.0 2017 3101 153937 52819016037
03102 Caldera 318653.2 2017 3102 17662 5628052276
03103 Tierra Amarilla 333194.9 2017 3103 14019 4671058718
03201 Chañaral 286389.3 2017 3201 12219 3499391196
03202 Diego de Almagro 351583.9 2017 3202 13925 4895805596
03301 Vallenar 315981.5 2017 3301 51917 16404810756
03303 Freirina 289049.9 2017 3303 7041 2035200054
03304 Huasco 337414.8 2017 3304 10149 3424422750
04101 La Serena 279340.1 2017 4101 221054 61749247282
04102 Coquimbo 269078.6 2017 4102 227730 61277269093
04103 Andacollo 258539.7 2017 4103 11044 2855312920
04104 La Higuera 214257.0 2017 4104 4241 908664019
04106 Vicuña 254177.0 2017 4106 27771 7058750373
04201 Illapel 282139.3 2017 4201 30848 8703433491
04202 Canela 233397.3 2017 4202 9093 2122281844
04203 Los Vilos 285214.0 2017 4203 21382 6098444926
04204 Salamanca 262056.9 2017 4204 29347 7690585032
04301 Ovalle 280373.5 2017 4301 111272 31197719080
04302 Combarbalá 234537.3 2017 4302 13322 3124505460
04303 Monte Patria 225369.1 2017 4303 30751 6930326684
04304 Punitaqui 212496.1 2017 4304 10956 2328107498
05101 Valparaíso 306572.5 2017 5101 296655 90946261553
05102 Casablanca 348088.6 2017 5102 26867 9352095757
05103 Concón 333932.4 2017 5103 42152 14075920021
05105 Puchuncaví 296035.5 2017 5105 18546 5490274928
05107 Quintero 308224.7 2017 5107 31923 9839456903
05109 Viña del Mar 354715.9 2017 5109 334248 118563074323
05301 Los Andes 355446.2 2017 5301 66708 23711104774
05302 Calle Larga 246387.3 2017 5302 14832 3654416747
05303 Rinconada 279807.9 2017 5303 10207 2855998928
05304 San Esteban 219571.6 2017 5304 18855 4140022481
05401 La Ligua 259482.3 2017 5401 35390 9183080280
05402 Cabildo 262745.9 2017 5402 19388 5094117762
05403 Papudo 302317.1 2017 5403 6356 1921527704
05404 Petorca 237510.8 2017 5404 9826 2333781007
05405 Zapallar 294389.2 2017 5405 7339 2160521991
05501 Quillota 288694.2 2017 5501 90517 26131733924
05502 Calera 282823.6 2017 5502 50554 14297866792
05503 Hijuelas 268449.7 2017 5503 17988 4828872604

3 Unión Censo-Casen

Integramos a la tabla censal de frecuencias la tabla de ingresos expandidos de la Casen.

comunas_con_ing_exp = merge( x = comuna_corr, y = h_y_m_2017_censo, by = "código", all.x = TRUE)
comunas_con_ing_exp <-comunas_con_ing_exp[!(is.na(comunas_con_ing_exp$Ingresos_expandidos)),]
r3_100 <- comunas_con_ing_exp[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
código zona Freq anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos
1 01101 1101021001 257 2017 Iquique 375676.9 2017 1101 191468 71930106513
2 01101 1101021002 11 2017 Iquique 375676.9 2017 1101 191468 71930106513
3 01101 1101011001 807 2017 Iquique 375676.9 2017 1101 191468 71930106513
4 01101 1101011002 452 2017 Iquique 375676.9 2017 1101 191468 71930106513
5 01101 1101021005 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513
6 01101 1101031001 653 2017 Iquique 375676.9 2017 1101 191468 71930106513
7 01101 1101031002 837 2017 Iquique 375676.9 2017 1101 191468 71930106513
8 01101 1101031003 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513
9 01101 1101031004 903 2017 Iquique 375676.9 2017 1101 191468 71930106513
10 01101 1101041001 417 2017 Iquique 375676.9 2017 1101 191468 71930106513
11 01101 1101041002 524 2017 Iquique 375676.9 2017 1101 191468 71930106513
12 01101 1101041003 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513
13 01101 1101041004 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513
14 01101 1101041005 894 2017 Iquique 375676.9 2017 1101 191468 71930106513
15 01101 1101041006 588 2017 Iquique 375676.9 2017 1101 191468 71930106513
16 01101 1101021003 571 2017 Iquique 375676.9 2017 1101 191468 71930106513
17 01101 1101021004 451 2017 Iquique 375676.9 2017 1101 191468 71930106513
18 01101 1101051003 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513
19 01101 1101051004 558 2017 Iquique 375676.9 2017 1101 191468 71930106513
20 01101 1101051005 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513
21 01101 1101051006 858 2017 Iquique 375676.9 2017 1101 191468 71930106513
22 01101 1101061001 606 2017 Iquique 375676.9 2017 1101 191468 71930106513
23 01101 1101061002 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513
24 01101 1101061003 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513
25 01101 1101061004 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513
26 01101 1101061005 626 2017 Iquique 375676.9 2017 1101 191468 71930106513
27 01101 1101071001 691 2017 Iquique 375676.9 2017 1101 191468 71930106513
28 01101 1101071002 753 2017 Iquique 375676.9 2017 1101 191468 71930106513
29 01101 1101051001 857 2017 Iquique 375676.9 2017 1101 191468 71930106513
30 01101 1101051002 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513
31 01101 1101081001 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513
32 01101 1101081002 852 2017 Iquique 375676.9 2017 1101 191468 71930106513
33 01101 1101081003 658 2017 Iquique 375676.9 2017 1101 191468 71930106513
34 01101 1101081004 762 2017 Iquique 375676.9 2017 1101 191468 71930106513
35 01101 1101101001 594 2017 Iquique 375676.9 2017 1101 191468 71930106513
36 01101 1101101002 993 2017 Iquique 375676.9 2017 1101 191468 71930106513
37 01101 1101101003 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513
38 01101 1101101004 838 2017 Iquique 375676.9 2017 1101 191468 71930106513
39 01101 1101101005 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513
40 01101 1101101006 940 2017 Iquique 375676.9 2017 1101 191468 71930106513
41 01101 1101111001 869 2017 Iquique 375676.9 2017 1101 191468 71930106513
42 01101 1101071003 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513
43 01101 1101071004 605 2017 Iquique 375676.9 2017 1101 191468 71930106513
44 01101 1101111004 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513
45 01101 1101111005 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513
46 01101 1101111006 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513
47 01101 1101111007 827 2017 Iquique 375676.9 2017 1101 191468 71930106513
48 01101 1101111008 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513
49 01101 1101111009 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513
50 01101 1101111010 97 2017 Iquique 375676.9 2017 1101 191468 71930106513
51 01101 1101111011 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513
52 01101 1101111012 820 2017 Iquique 375676.9 2017 1101 191468 71930106513
53 01101 1101111013 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513
54 01101 1101111014 842 2017 Iquique 375676.9 2017 1101 191468 71930106513
55 01101 1101111002 578 2017 Iquique 375676.9 2017 1101 191468 71930106513
56 01101 1101111003 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513
57 01101 1101991999 194 2017 Iquique 375676.9 2017 1101 191468 71930106513
58 01107 1107011002 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
59 01107 1107011003 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
60 01107 1107021001 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
61 01107 1107021003 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
62 01107 1107021004 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
63 01107 1107021005 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
64 01107 1107021002 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
65 01107 1107021007 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
66 01107 1107021008 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
67 01107 1107031001 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
68 01107 1107021006 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
69 01107 1107011001 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
70 01107 1107041001 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
71 01107 1107041002 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
72 01107 1107041003 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
73 01107 1107041004 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
74 01107 1107041005 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
75 01107 1107041006 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
76 01107 1107041007 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
77 01107 1107991999 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
78 01107 1107031003 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
79 01107 1107031002 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
80 01401 1401011002 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
81 01401 1401991999 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
82 01401 1401011001 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
85 01405 1405011001 895 2017 Pica 330061.1 2017 1405 9296 3068247619
86 01405 1405991999 14 2017 Pica 330061.1 2017 1405 9296 3068247619
NA NA NA NA NA NA NA NA NA NA NA
NA.1 NA NA NA NA NA NA NA NA NA NA
NA.2 NA NA NA NA NA NA NA NA NA NA
NA.3 NA NA NA NA NA NA NA NA NA NA
NA.4 NA NA NA NA NA NA NA NA NA NA
NA.5 NA NA NA NA NA NA NA NA NA NA
NA.6 NA NA NA NA NA NA NA NA NA NA
NA.7 NA NA NA NA NA NA NA NA NA NA
NA.8 NA NA NA NA NA NA NA NA NA NA
NA.9 NA NA NA NA NA NA NA NA NA NA
NA.10 NA NA NA NA NA NA NA NA NA NA
NA.11 NA NA NA NA NA NA NA NA NA NA
NA.12 NA NA NA NA NA NA NA NA NA NA
NA.13 NA NA NA NA NA NA NA NA NA NA
NA.14 NA NA NA NA NA NA NA NA NA NA
NA.15 NA NA NA NA NA NA NA NA NA NA


4 Proporción poblacional zonal respecto a la comunal

Del censo obtenemos la cantidad de población a nivel de zona y estimamos su proporción a nivel comunal. Ya hemos calculado ésta proporción aquí.

prop_pob <- readRDS("../tabla_de_prop_pob.rds")
names(prop_pob)[1] <- "zona"
names(prop_pob)[3] <- "p_poblacional"

Veamos los 100 primeros registros:

r3_100 <- prop_pob[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona Freq p_poblacional código
1101011001 2491 0.0130100 01101
1101011002 1475 0.0077036 01101
1101021001 1003 0.0052385 01101
1101021002 54 0.0002820 01101
1101021003 2895 0.0151200 01101
1101021004 2398 0.0125243 01101
1101021005 4525 0.0236332 01101
1101031001 2725 0.0142321 01101
1101031002 3554 0.0185618 01101
1101031003 5246 0.0273988 01101
1101031004 3389 0.0177001 01101
1101041001 1800 0.0094010 01101
1101041002 2538 0.0132555 01101
1101041003 3855 0.0201339 01101
1101041004 5663 0.0295767 01101
1101041005 4162 0.0217373 01101
1101041006 2689 0.0140441 01101
1101051001 3296 0.0172144 01101
1101051002 4465 0.0233198 01101
1101051003 4656 0.0243174 01101
1101051004 2097 0.0109522 01101
1101051005 3569 0.0186402 01101
1101051006 2741 0.0143157 01101
1101061001 1625 0.0084871 01101
1101061002 4767 0.0248971 01101
1101061003 4826 0.0252053 01101
1101061004 4077 0.0212934 01101
1101061005 2166 0.0113126 01101
1101071001 2324 0.0121378 01101
1101071002 2801 0.0146291 01101
1101071003 3829 0.0199981 01101
1101071004 1987 0.0103777 01101
1101081001 5133 0.0268087 01101
1101081002 3233 0.0168853 01101
1101081003 2122 0.0110828 01101
1101081004 2392 0.0124929 01101
1101092001 57 0.0002977 01101
1101092004 247 0.0012900 01101
1101092005 76 0.0003969 01101
1101092006 603 0.0031494 01101
1101092007 84 0.0004387 01101
1101092010 398 0.0020787 01101
1101092012 58 0.0003029 01101
1101092014 23 0.0001201 01101
1101092016 20 0.0001045 01101
1101092017 8 0.0000418 01101
1101092018 74 0.0003865 01101
1101092019 25 0.0001306 01101
1101092021 177 0.0009244 01101
1101092022 23 0.0001201 01101
1101092023 288 0.0015042 01101
1101092024 14 0.0000731 01101
1101092901 30 0.0001567 01101
1101101001 2672 0.0139553 01101
1101101002 4398 0.0229699 01101
1101101003 4524 0.0236280 01101
1101101004 3544 0.0185096 01101
1101101005 4911 0.0256492 01101
1101101006 3688 0.0192617 01101
1101111001 3886 0.0202958 01101
1101111002 2312 0.0120751 01101
1101111003 4874 0.0254560 01101
1101111004 4543 0.0237272 01101
1101111005 4331 0.0226200 01101
1101111006 3253 0.0169898 01101
1101111007 4639 0.0242286 01101
1101111008 4881 0.0254925 01101
1101111009 5006 0.0261454 01101
1101111010 366 0.0019115 01101
1101111011 4351 0.0227244 01101
1101111012 2926 0.0152819 01101
1101111013 3390 0.0177053 01101
1101111014 2940 0.0153550 01101
1101112003 33 0.0001724 01101
1101112013 104 0.0005432 01101
1101112019 34 0.0001776 01101
1101112025 21 0.0001097 01101
1101112901 6 0.0000313 01101
1101991999 1062 0.0055466 01101
1107011001 4104 0.0378685 01107
1107011002 4360 0.0402307 01107
1107011003 8549 0.0788835 01107
1107012003 3 0.0000277 01107
1107012901 17 0.0001569 01107
1107021001 6701 0.0618316 01107
1107021002 3971 0.0366413 01107
1107021003 6349 0.0585836 01107
1107021004 5125 0.0472895 01107
1107021005 4451 0.0410704 01107
1107021006 3864 0.0356540 01107
1107021007 5235 0.0483045 01107
1107021008 4566 0.0421315 01107
1107031001 4195 0.0387082 01107
1107031002 7099 0.0655040 01107
1107031003 4720 0.0435525 01107
1107032005 38 0.0003506 01107
1107032006 2399 0.0221361 01107
1107032008 4 0.0000369 01107
1107041001 3630 0.0334948 01107
1107041002 5358 0.0494394 01107


5 Ingreso medio

Deseamos el valor del ingreso promedio a nivel comunal, pero expandido a nivel zonal. Ésta información está contenida en el campo promedio_i de la tabla obtenida en el punto 3.

r3_100 <- comunas_con_ing_exp[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
código zona Freq anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos
1 01101 1101021001 257 2017 Iquique 375676.9 2017 1101 191468 71930106513
2 01101 1101021002 11 2017 Iquique 375676.9 2017 1101 191468 71930106513
3 01101 1101011001 807 2017 Iquique 375676.9 2017 1101 191468 71930106513
4 01101 1101011002 452 2017 Iquique 375676.9 2017 1101 191468 71930106513
5 01101 1101021005 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513
6 01101 1101031001 653 2017 Iquique 375676.9 2017 1101 191468 71930106513
7 01101 1101031002 837 2017 Iquique 375676.9 2017 1101 191468 71930106513
8 01101 1101031003 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513
9 01101 1101031004 903 2017 Iquique 375676.9 2017 1101 191468 71930106513
10 01101 1101041001 417 2017 Iquique 375676.9 2017 1101 191468 71930106513
11 01101 1101041002 524 2017 Iquique 375676.9 2017 1101 191468 71930106513
12 01101 1101041003 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513
13 01101 1101041004 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513
14 01101 1101041005 894 2017 Iquique 375676.9 2017 1101 191468 71930106513
15 01101 1101041006 588 2017 Iquique 375676.9 2017 1101 191468 71930106513
16 01101 1101021003 571 2017 Iquique 375676.9 2017 1101 191468 71930106513
17 01101 1101021004 451 2017 Iquique 375676.9 2017 1101 191468 71930106513
18 01101 1101051003 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513
19 01101 1101051004 558 2017 Iquique 375676.9 2017 1101 191468 71930106513
20 01101 1101051005 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513
21 01101 1101051006 858 2017 Iquique 375676.9 2017 1101 191468 71930106513
22 01101 1101061001 606 2017 Iquique 375676.9 2017 1101 191468 71930106513
23 01101 1101061002 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513
24 01101 1101061003 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513
25 01101 1101061004 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513
26 01101 1101061005 626 2017 Iquique 375676.9 2017 1101 191468 71930106513
27 01101 1101071001 691 2017 Iquique 375676.9 2017 1101 191468 71930106513
28 01101 1101071002 753 2017 Iquique 375676.9 2017 1101 191468 71930106513
29 01101 1101051001 857 2017 Iquique 375676.9 2017 1101 191468 71930106513
30 01101 1101051002 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513
31 01101 1101081001 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513
32 01101 1101081002 852 2017 Iquique 375676.9 2017 1101 191468 71930106513
33 01101 1101081003 658 2017 Iquique 375676.9 2017 1101 191468 71930106513
34 01101 1101081004 762 2017 Iquique 375676.9 2017 1101 191468 71930106513
35 01101 1101101001 594 2017 Iquique 375676.9 2017 1101 191468 71930106513
36 01101 1101101002 993 2017 Iquique 375676.9 2017 1101 191468 71930106513
37 01101 1101101003 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513
38 01101 1101101004 838 2017 Iquique 375676.9 2017 1101 191468 71930106513
39 01101 1101101005 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513
40 01101 1101101006 940 2017 Iquique 375676.9 2017 1101 191468 71930106513
41 01101 1101111001 869 2017 Iquique 375676.9 2017 1101 191468 71930106513
42 01101 1101071003 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513
43 01101 1101071004 605 2017 Iquique 375676.9 2017 1101 191468 71930106513
44 01101 1101111004 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513
45 01101 1101111005 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513
46 01101 1101111006 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513
47 01101 1101111007 827 2017 Iquique 375676.9 2017 1101 191468 71930106513
48 01101 1101111008 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513
49 01101 1101111009 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513
50 01101 1101111010 97 2017 Iquique 375676.9 2017 1101 191468 71930106513
51 01101 1101111011 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513
52 01101 1101111012 820 2017 Iquique 375676.9 2017 1101 191468 71930106513
53 01101 1101111013 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513
54 01101 1101111014 842 2017 Iquique 375676.9 2017 1101 191468 71930106513
55 01101 1101111002 578 2017 Iquique 375676.9 2017 1101 191468 71930106513
56 01101 1101111003 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513
57 01101 1101991999 194 2017 Iquique 375676.9 2017 1101 191468 71930106513
58 01107 1107011002 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
59 01107 1107011003 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
60 01107 1107021001 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
61 01107 1107021003 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
62 01107 1107021004 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
63 01107 1107021005 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
64 01107 1107021002 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
65 01107 1107021007 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
66 01107 1107021008 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
67 01107 1107031001 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
68 01107 1107021006 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
69 01107 1107011001 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
70 01107 1107041001 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
71 01107 1107041002 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
72 01107 1107041003 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
73 01107 1107041004 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
74 01107 1107041005 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
75 01107 1107041006 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
76 01107 1107041007 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
77 01107 1107991999 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
78 01107 1107031003 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
79 01107 1107031002 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496
80 01401 1401011002 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
81 01401 1401991999 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
82 01401 1401011001 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883
85 01405 1405011001 895 2017 Pica 330061.1 2017 1405 9296 3068247619
86 01405 1405991999 14 2017 Pica 330061.1 2017 1405 9296 3068247619
NA NA NA NA NA NA NA NA NA NA NA
NA.1 NA NA NA NA NA NA NA NA NA NA
NA.2 NA NA NA NA NA NA NA NA NA NA
NA.3 NA NA NA NA NA NA NA NA NA NA
NA.4 NA NA NA NA NA NA NA NA NA NA
NA.5 NA NA NA NA NA NA NA NA NA NA
NA.6 NA NA NA NA NA NA NA NA NA NA
NA.7 NA NA NA NA NA NA NA NA NA NA
NA.8 NA NA NA NA NA NA NA NA NA NA
NA.9 NA NA NA NA NA NA NA NA NA NA
NA.10 NA NA NA NA NA NA NA NA NA NA
NA.11 NA NA NA NA NA NA NA NA NA NA
NA.12 NA NA NA NA NA NA NA NA NA NA
NA.13 NA NA NA NA NA NA NA NA NA NA
NA.14 NA NA NA NA NA NA NA NA NA NA
NA.15 NA NA NA NA NA NA NA NA NA NA


6 Ingreso promedio expandido por zona (multi_pob)

En éste momento vamos a construir nuestra variable dependiente de regresión aplicando la siguiente fórmula:

\[ multi\_pob = promedio\_i \cdot personas \cdot p\_poblacional \]

Para ello integramos a la tabla de ingresos expandidos a nivel zonal (punto 3) la tabla de proporciones poblacionales zonales respecto al total comunal (punto 4) :

h_y_m_comuna_corr_01 = merge( x = comunas_con_ing_exp, y = prop_pob, by = "zona", all.x = TRUE)
tablamadre <- head(h_y_m_comuna_corr_01,100)
kbl(tablamadre) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona código.x Freq.x anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos Freq.y p_poblacional código.y
1101011001 01101 807 2017 Iquique 375676.9 2017 1101 191468 71930106513 2491 0.0130100 01101
1101011002 01101 452 2017 Iquique 375676.9 2017 1101 191468 71930106513 1475 0.0077036 01101
1101021001 01101 257 2017 Iquique 375676.9 2017 1101 191468 71930106513 1003 0.0052385 01101
1101021002 01101 11 2017 Iquique 375676.9 2017 1101 191468 71930106513 54 0.0002820 01101
1101021003 01101 571 2017 Iquique 375676.9 2017 1101 191468 71930106513 2895 0.0151200 01101
1101021004 01101 451 2017 Iquique 375676.9 2017 1101 191468 71930106513 2398 0.0125243 01101
1101021005 01101 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513 4525 0.0236332 01101
1101031001 01101 653 2017 Iquique 375676.9 2017 1101 191468 71930106513 2725 0.0142321 01101
1101031002 01101 837 2017 Iquique 375676.9 2017 1101 191468 71930106513 3554 0.0185618 01101
1101031003 01101 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513 5246 0.0273988 01101
1101031004 01101 903 2017 Iquique 375676.9 2017 1101 191468 71930106513 3389 0.0177001 01101
1101041001 01101 417 2017 Iquique 375676.9 2017 1101 191468 71930106513 1800 0.0094010 01101
1101041002 01101 524 2017 Iquique 375676.9 2017 1101 191468 71930106513 2538 0.0132555 01101
1101041003 01101 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513 3855 0.0201339 01101
1101041004 01101 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513 5663 0.0295767 01101
1101041005 01101 894 2017 Iquique 375676.9 2017 1101 191468 71930106513 4162 0.0217373 01101
1101041006 01101 588 2017 Iquique 375676.9 2017 1101 191468 71930106513 2689 0.0140441 01101
1101051001 01101 857 2017 Iquique 375676.9 2017 1101 191468 71930106513 3296 0.0172144 01101
1101051002 01101 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513 4465 0.0233198 01101
1101051003 01101 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513 4656 0.0243174 01101
1101051004 01101 558 2017 Iquique 375676.9 2017 1101 191468 71930106513 2097 0.0109522 01101
1101051005 01101 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513 3569 0.0186402 01101
1101051006 01101 858 2017 Iquique 375676.9 2017 1101 191468 71930106513 2741 0.0143157 01101
1101061001 01101 606 2017 Iquique 375676.9 2017 1101 191468 71930106513 1625 0.0084871 01101
1101061002 01101 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513 4767 0.0248971 01101
1101061003 01101 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513 4826 0.0252053 01101
1101061004 01101 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513 4077 0.0212934 01101
1101061005 01101 626 2017 Iquique 375676.9 2017 1101 191468 71930106513 2166 0.0113126 01101
1101071001 01101 691 2017 Iquique 375676.9 2017 1101 191468 71930106513 2324 0.0121378 01101
1101071002 01101 753 2017 Iquique 375676.9 2017 1101 191468 71930106513 2801 0.0146291 01101
1101071003 01101 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513 3829 0.0199981 01101
1101071004 01101 605 2017 Iquique 375676.9 2017 1101 191468 71930106513 1987 0.0103777 01101
1101081001 01101 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513 5133 0.0268087 01101
1101081002 01101 852 2017 Iquique 375676.9 2017 1101 191468 71930106513 3233 0.0168853 01101
1101081003 01101 658 2017 Iquique 375676.9 2017 1101 191468 71930106513 2122 0.0110828 01101
1101081004 01101 762 2017 Iquique 375676.9 2017 1101 191468 71930106513 2392 0.0124929 01101
1101101001 01101 594 2017 Iquique 375676.9 2017 1101 191468 71930106513 2672 0.0139553 01101
1101101002 01101 993 2017 Iquique 375676.9 2017 1101 191468 71930106513 4398 0.0229699 01101
1101101003 01101 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513 4524 0.0236280 01101
1101101004 01101 838 2017 Iquique 375676.9 2017 1101 191468 71930106513 3544 0.0185096 01101
1101101005 01101 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513 4911 0.0256492 01101
1101101006 01101 940 2017 Iquique 375676.9 2017 1101 191468 71930106513 3688 0.0192617 01101
1101111001 01101 869 2017 Iquique 375676.9 2017 1101 191468 71930106513 3886 0.0202958 01101
1101111002 01101 578 2017 Iquique 375676.9 2017 1101 191468 71930106513 2312 0.0120751 01101
1101111003 01101 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513 4874 0.0254560 01101
1101111004 01101 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513 4543 0.0237272 01101
1101111005 01101 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513 4331 0.0226200 01101
1101111006 01101 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513 3253 0.0169898 01101
1101111007 01101 827 2017 Iquique 375676.9 2017 1101 191468 71930106513 4639 0.0242286 01101
1101111008 01101 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513 4881 0.0254925 01101
1101111009 01101 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513 5006 0.0261454 01101
1101111010 01101 97 2017 Iquique 375676.9 2017 1101 191468 71930106513 366 0.0019115 01101
1101111011 01101 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513 4351 0.0227244 01101
1101111012 01101 820 2017 Iquique 375676.9 2017 1101 191468 71930106513 2926 0.0152819 01101
1101111013 01101 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513 3390 0.0177053 01101
1101111014 01101 842 2017 Iquique 375676.9 2017 1101 191468 71930106513 2940 0.0153550 01101
1101991999 01101 194 2017 Iquique 375676.9 2017 1101 191468 71930106513 1062 0.0055466 01101
1107011001 01107 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4104 0.0378685 01107
1107011002 01107 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4360 0.0402307 01107
1107011003 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 8549 0.0788835 01107
1107021001 01107 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6701 0.0618316 01107
1107021002 01107 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3971 0.0366413 01107
1107021003 01107 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6349 0.0585836 01107
1107021004 01107 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5125 0.0472895 01107
1107021005 01107 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4451 0.0410704 01107
1107021006 01107 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3864 0.0356540 01107
1107021007 01107 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5235 0.0483045 01107
1107021008 01107 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4566 0.0421315 01107
1107031001 01107 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4195 0.0387082 01107
1107031002 01107 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 7099 0.0655040 01107
1107031003 01107 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4720 0.0435525 01107
1107041001 01107 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3630 0.0334948 01107
1107041002 01107 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5358 0.0494394 01107
1107041003 01107 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4520 0.0417070 01107
1107041004 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5317 0.0490611 01107
1107041005 01107 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3725 0.0343714 01107
1107041006 01107 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4009 0.0369919 01107
1107041007 01107 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5217 0.0481384 01107
1107991999 01107 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 819 0.0075571 01107
1401011001 01401 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 2771 0.1763732 01401
1401011002 01401 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 6506 0.4141048 01401
1401991999 01401 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 818 0.0520654 01401
1405011001 01405 895 2017 Pica 330061.1 2017 1405 9296 3068247619 3876 0.4169535 01405
1405991999 01405 14 2017 Pica 330061.1 2017 1405 9296 3068247619 36 0.0038726 01405


Hacemos la multiplicación que queda almacenada en la variable multi_pob:

h_y_m_comuna_corr_01$multi_pob <- h_y_m_comuna_corr_01$promedio_i * h_y_m_comuna_corr_01$personas * h_y_m_comuna_corr_01$p_poblacional
tablamadre <- head(h_y_m_comuna_corr_01,100)
kbl(tablamadre) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona código.x Freq.x anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos Freq.y p_poblacional código.y multi_pob
1101011001 01101 807 2017 Iquique 375676.9 2017 1101 191468 71930106513 2491 0.0130100 01101 935811182
1101011002 01101 452 2017 Iquique 375676.9 2017 1101 191468 71930106513 1475 0.0077036 01101 554123442
1101021001 01101 257 2017 Iquique 375676.9 2017 1101 191468 71930106513 1003 0.0052385 01101 376803940
1101021002 01101 11 2017 Iquique 375676.9 2017 1101 191468 71930106513 54 0.0002820 01101 20286553
1101021003 01101 571 2017 Iquique 375676.9 2017 1101 191468 71930106513 2895 0.0151200 01101 1087584653
1101021004 01101 451 2017 Iquique 375676.9 2017 1101 191468 71930106513 2398 0.0125243 01101 900873229
1101021005 01101 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513 4525 0.0236332 01101 1699938016
1101031001 01101 653 2017 Iquique 375676.9 2017 1101 191468 71930106513 2725 0.0142321 01101 1023719578
1101031002 01101 837 2017 Iquique 375676.9 2017 1101 191468 71930106513 3554 0.0185618 01101 1335155736
1101031003 01101 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513 5246 0.0273988 01101 1970801067
1101031004 01101 903 2017 Iquique 375676.9 2017 1101 191468 71930106513 3389 0.0177001 01101 1273169046
1101041001 01101 417 2017 Iquique 375676.9 2017 1101 191468 71930106513 1800 0.0094010 01101 676218437
1101041002 01101 524 2017 Iquique 375676.9 2017 1101 191468 71930106513 2538 0.0132555 01101 953467996
1101041003 01101 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513 3855 0.0201339 01101 1448234486
1101041004 01101 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513 5663 0.0295767 01101 2127458339
1101041005 01101 894 2017 Iquique 375676.9 2017 1101 191468 71930106513 4162 0.0217373 01101 1563567297
1101041006 01101 588 2017 Iquique 375676.9 2017 1101 191468 71930106513 2689 0.0140441 01101 1010195210
1101051001 01101 857 2017 Iquique 375676.9 2017 1101 191468 71930106513 3296 0.0172144 01101 1238231094
1101051002 01101 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513 4465 0.0233198 01101 1677397401
1101051003 01101 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513 4656 0.0243174 01101 1749151691
1101051004 01101 558 2017 Iquique 375676.9 2017 1101 191468 71930106513 2097 0.0109522 01101 787794479
1101051005 01101 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513 3569 0.0186402 01101 1340790890
1101051006 01101 858 2017 Iquique 375676.9 2017 1101 191468 71930106513 2741 0.0143157 01101 1029730409
1101061001 01101 606 2017 Iquique 375676.9 2017 1101 191468 71930106513 1625 0.0084871 01101 610474978
1101061002 01101 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513 4767 0.0248971 01101 1790851828
1101061003 01101 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513 4826 0.0252053 01101 1813016765
1101061004 01101 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513 4077 0.0212934 01101 1531634760
1101061005 01101 626 2017 Iquique 375676.9 2017 1101 191468 71930106513 2166 0.0113126 01101 813716186
1101071001 01101 691 2017 Iquique 375676.9 2017 1101 191468 71930106513 2324 0.0121378 01101 873073138
1101071002 01101 753 2017 Iquique 375676.9 2017 1101 191468 71930106513 2801 0.0146291 01101 1052271024
1101071003 01101 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513 3829 0.0199981 01101 1438466887
1101071004 01101 605 2017 Iquique 375676.9 2017 1101 191468 71930106513 1987 0.0103777 01101 746470019
1101081001 01101 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513 5133 0.0268087 01101 1928349577
1101081002 01101 852 2017 Iquique 375676.9 2017 1101 191468 71930106513 3233 0.0168853 01101 1214563448
1101081003 01101 658 2017 Iquique 375676.9 2017 1101 191468 71930106513 2122 0.0110828 01101 797186402
1101081004 01101 762 2017 Iquique 375676.9 2017 1101 191468 71930106513 2392 0.0124929 01101 898619168
1101101001 01101 594 2017 Iquique 375676.9 2017 1101 191468 71930106513 2672 0.0139553 01101 1003808702
1101101002 01101 993 2017 Iquique 375676.9 2017 1101 191468 71930106513 4398 0.0229699 01101 1652227048
1101101003 01101 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513 4524 0.0236280 01101 1699562339
1101101004 01101 838 2017 Iquique 375676.9 2017 1101 191468 71930106513 3544 0.0185096 01101 1331398967
1101101005 01101 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513 4911 0.0256492 01101 1844949303
1101101006 01101 940 2017 Iquique 375676.9 2017 1101 191468 71930106513 3688 0.0192617 01101 1385496442
1101111001 01101 869 2017 Iquique 375676.9 2017 1101 191468 71930106513 3886 0.0202958 01101 1459880470
1101111002 01101 578 2017 Iquique 375676.9 2017 1101 191468 71930106513 2312 0.0120751 01101 868565015
1101111003 01101 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513 4874 0.0254560 01101 1831049257
1101111004 01101 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513 4543 0.0237272 01101 1706700200
1101111005 01101 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513 4331 0.0226200 01101 1627056695
1101111006 01101 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513 3253 0.0169898 01101 1222076987
1101111007 01101 827 2017 Iquique 375676.9 2017 1101 191468 71930106513 4639 0.0242286 01101 1742765183
1101111008 01101 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513 4881 0.0254925 01101 1833678995
1101111009 01101 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513 5006 0.0261454 01101 1880638609
1101111010 01101 97 2017 Iquique 375676.9 2017 1101 191468 71930106513 366 0.0019115 01101 137497749
1101111011 01101 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513 4351 0.0227244 01101 1634570233
1101111012 01101 820 2017 Iquique 375676.9 2017 1101 191468 71930106513 2926 0.0152819 01101 1099230637
1101111013 01101 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513 3390 0.0177053 01101 1273544723
1101111014 01101 842 2017 Iquique 375676.9 2017 1101 191468 71930106513 2940 0.0153550 01101 1104490114
1101991999 01101 194 2017 Iquique 375676.9 2017 1101 191468 71930106513 1062 0.0055466 01101 398968878
1107011001 01107 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4104 0.0378685 01107 1278690352
1107011002 01107 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4360 0.0402307 01107 1358452713
1107011003 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 8549 0.0788835 01107 2663626661
1107021001 01107 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6701 0.0618316 01107 2087842117
1107021002 01107 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3971 0.0366413 01107 1237251313
1107021003 01107 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6349 0.0585836 01107 1978168870
1107021004 01107 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5125 0.0472895 01107 1596805081
1107021005 01107 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4451 0.0410704 01107 1386805740
1107021006 01107 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3864 0.0356540 01107 1203913138
1107021007 01107 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5235 0.0483045 01107 1631077971
1107021008 01107 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4566 0.0421315 01107 1422636488
1107031001 01107 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4195 0.0387082 01107 1307043379
1107031002 01107 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 7099 0.0655040 01107 2211847663
1107031003 01107 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4720 0.0435525 01107 1470618533
1107041001 01107 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3630 0.0334948 01107 1131005355
1107041002 01107 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5358 0.0494394 01107 1669401293
1107041003 01107 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4520 0.0417070 01107 1408304189
1107041004 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5317 0.0490611 01107 1656626852
1107041005 01107 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3725 0.0343714 01107 1160604669
1107041006 01107 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4009 0.0369919 01107 1249091038
1107041007 01107 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5217 0.0481384 01107 1625469680
1107991999 01107 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 819 0.0075571 01107 255177241
1401011001 01401 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 2771 0.1763732 01401 876019768
1401011002 01401 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 6506 0.4141048 01401 2056797044
1401991999 01401 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 818 0.0520654 01401 258601288
1405011001 01405 895 2017 Pica 330061.1 2017 1405 9296 3068247619 3876 0.4169535 01405 1279316671
1405991999 01405 14 2017 Pica 330061.1 2017 1405 9296 3068247619 36 0.0038726 01405 11882198

7 Análisis de regresión

Aplicaremos un análisis de regresión donde:

\[ Y(dependiente) = ingreso \ expandido \ por \ zona \ (multi\_pob)\]

\[ X(independiente) = frecuencia \ de \ población \ que \ posee \ la \ variable \ Censal \ respecto \ a \ la \ zona \ (Freq.x) \]

7.1 Diagrama de dispersión loess

scatter.smooth(x=h_y_m_comuna_corr_01$Freq.x, y=h_y_m_comuna_corr_01$multi_pob, main="multi_pob ~ Freq.x",
     xlab = "Freq.x",
     ylab = "multi_pob",
           col = 2)

7.2 Outliers

Hemos demostrado en el punto 5.7.2 de aquí que la exclusión de ouliers no genera ninguna mejora en el modelo de regresión.

7.3 Modelo lineal

Aplicaremos un análisis de regresión lineal del ingreso expandido por zona sobre las frecuencias de respuestas zonales.

linearMod <- lm( multi_pob~(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 
## 
## Call:
## lm(formula = multi_pob ~ (Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -467801362 -131533867  -19308985  105879621 1027825825 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 142703202   64860152    2.20   0.0306 *  
## Freq.x        1299476      67952   19.12   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.24e+08 on 82 degrees of freedom
## Multiple R-squared:  0.8168, Adjusted R-squared:  0.8146 
## F-statistic: 365.7 on 1 and 82 DF,  p-value: < 2.2e-16

7.4 Gráfica de la recta de regresión lineal

ggplot(h_y_m_comuna_corr_01, aes(x = Freq.x , y = multi_pob)) + 
  geom_point() +
  stat_smooth(method = "lm", col = "red")

Si bien obtenemos nuestro modelo lineal da cuenta del 0.8214 de la variabilidad de los datos de respuesta en torno a su media, modelos alternativos pueden ofrecernos una explicación de la variable dependiente aún mayor.

8 Modelos alternativos

### 8.1 Modelo cuadrático

linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "cuadrático"
sintaxis <- "linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = '\'beta_0 + '\'beta_1  X^2  $$"
modelos1 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.2 Modelo cúbico
 
linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "cúbico"
sintaxis <- "linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = '\'beta_0 + '\'beta_1  X^3  $$"
modelos2 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.3 Modelo logarítmico
 
linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "logarítmico"
sintaxis <- "linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = '\'beta_0 + '\'beta_1 ln X  $$"
modelos3 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.5 Modelo con raíz cuadrada 
 
linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "raíz cuadrada"
sintaxis <- "linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = '\'beta_0 + '\'beta_1  '\'sqrt {X}  $$"
modelos5 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.6 Modelo raíz-raíz
 
linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "raíz-raíz"
sintaxis <- "linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = {'\'beta_0}^2 + 2  '\'beta_0  '\'beta_1 '\'sqrt{X}+  '\'beta_1^2 X  $$"
modelos6 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.7 Modelo log-raíz
 
linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "log-raíz"
sintaxis <- "linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = e^{'\'beta_0 + '\'beta_1 '\'sqrt{X}} $$"
modelos7 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.8 Modelo raíz-log
 
linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "raíz-log"
sintaxis <- "linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = {'\'beta_0}^2 + 2  '\'beta_0  '\'beta_1 '\'ln{X}+  '\'beta_1^2 ln^2X  $$"
modelos8 <- cbind(modelo,dato,sintaxis,latex)
 
### 8.9 Modelo log-log
 
linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)
datos <- summary(linearMod)
dato <- datos$adj.r.squared
modelo <- "log-log"
sintaxis <- "linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)"
latex <- "$$ '\'hat Y = e^{'\'beta_0+'\'beta_1  ln{X}} $$"
modelos9 <- cbind(modelo,dato,sintaxis,latex)
 
modelos_bind <- rbind(modelos1,modelos2,modelos3,modelos5,modelos6,modelos7,modelos8,modelos9)
modelos_bind <- as.data.frame(modelos_bind)

modelos_bind <<- modelos_bind[order(modelos_bind$dato ),]
h_y_m_comuna_corr_01 <<- h_y_m_comuna_corr_01

kbl(modelos_bind) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
modelo dato sintaxis latex
3 logarítmico 0.608464579130019 linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = ''beta_0 + ''beta_1 ln X \]
6 log-raíz 0.782054751860565 linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = e^{''beta_0 + ''beta_1 ''sqrt{X}} \]
7 raíz-log 0.792192908719751 linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = {''beta_0}^2 + 2 ''beta_0 ''beta_1 ''ln{X}+ ''beta_1^2 ln^2X \]
4 raíz cuadrada 0.793054820442292 linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = ''beta_0 + ''beta_1 ''sqrt {X} \]
1 cuadrático 0.814612334906238 linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr_01) \[ ''hat Y = ''beta_0 + ''beta_1 X^2 \]
2 cúbico 0.814612334906238 linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01) \[ ''hat Y = ''beta_0 + ''beta_1 X^3 \]
5 raíz-raíz 0.878472394176248 linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = {''beta_0}^2 + 2 ''beta_0 ''beta_1 ''sqrt{X}+ ''beta_1^2 X \]
8 log-log 0.891995389342853 linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = e^{''beta_0+''beta_1 ln{X}} \] 
h_y_m_comuna_corr <- h_y_m_comuna_corr_01
metodo <- 8


switch (metodo,
        case = linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr)
)
summary(linearMod)
## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.03806 -0.12241 -0.00788  0.08472  1.68524 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  15.0053     0.2233    67.2   <2e-16 ***
## log(Freq.x)   0.8803     0.0336    26.2   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.2738 on 82 degrees of freedom
## Multiple R-squared:  0.8933, Adjusted R-squared:  0.892 
## F-statistic: 686.5 on 1 and 82 DF,  p-value: < 2.2e-16
aa <- linearMod$coefficients[1]
aa
## (Intercept) 
##    15.00532
bb <- linearMod$coefficients[2]
bb
## log(Freq.x) 
##   0.8803488

9 Modelo log-log (log-log)

Es éste el modelo que nos entrega el mayor coeficiente de determinación de todos (0.892).

9.1 Diagrama de dispersión sobre log-log

Desplegamos una curva suavizada por loess en el diagrama de dispersión.

scatter.smooth(x=log(h_y_m_comuna_corr$Freq.x), y=log(h_y_m_comuna_corr$multi_pob), lpars = list(col = "red", lwd = 2, lty = 1), main="multi_pob ~ Freq.x")

9.2 Modelo log-log

Observemos nuevamente el resultado sobre log-log.

linearMod <- lm(log( multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr)
summary(linearMod) 
## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.03806 -0.12241 -0.00788  0.08472  1.68524 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  15.0053     0.2233    67.2   <2e-16 ***
## log(Freq.x)   0.8803     0.0336    26.2   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.2738 on 82 degrees of freedom
## Multiple R-squared:  0.8933, Adjusted R-squared:  0.892 
## F-statistic: 686.5 on 1 and 82 DF,  p-value: < 2.2e-16
ggplot(h_y_m_comuna_corr, aes(x = log(Freq.x) , y = log(multi_pob))) + geom_point() + stat_smooth(method = "lm", col = "red")

9.3 Análisis de residuos

par(mfrow = c (2,2))
plot(linearMod)

9.4 Ecuación del modelo


\[ \hat Y = e^{15.00532 +0.8803488 \cdot ln{X}} \]


10 Aplicación la regresión a los valores de la variable a nivel de zona

Esta nueva variable se llamará: est_ing

h_y_m_comuna_corr$est_ing <- exp(aa+bb * log(h_y_m_comuna_corr$Freq.x))

r3_100 <- h_y_m_comuna_corr[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona código.x Freq.x anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos Freq.y p_poblacional código.y multi_pob est_ing
1 1101011001 01101 807 2017 Iquique 375676.9 2017 1101 191468 71930106513 2491 0.0130100 01101 935811182 1190667824
2 1101011002 01101 452 2017 Iquique 375676.9 2017 1101 191468 71930106513 1475 0.0077036 01101 554123442 714785823
3 1101021001 01101 257 2017 Iquique 375676.9 2017 1101 191468 71930106513 1003 0.0052385 01101 376803940 434820223
4 1101021002 01101 11 2017 Iquique 375676.9 2017 1101 191468 71930106513 54 0.0002820 01101 20286553 27134188
5 1101021003 01101 571 2017 Iquique 375676.9 2017 1101 191468 71930106513 2895 0.0151200 01101 1087584653 878070290
6 1101021004 01101 451 2017 Iquique 375676.9 2017 1101 191468 71930106513 2398 0.0125243 01101 900873229 713393469
7 1101021005 01101 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513 4525 0.0236332 01101 1699938016 1496129140
8 1101031001 01101 653 2017 Iquique 375676.9 2017 1101 191468 71930106513 2725 0.0142321 01101 1023719578 988174023
9 1101031002 01101 837 2017 Iquique 375676.9 2017 1101 191468 71930106513 3554 0.0185618 01101 1335155736 1229548994
10 1101031003 01101 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513 5246 0.0273988 01101 1970801067 1830081910
11 1101031004 01101 903 2017 Iquique 375676.9 2017 1101 191468 71930106513 3389 0.0177001 01101 1273169046 1314510781
12 1101041001 01101 417 2017 Iquique 375676.9 2017 1101 191468 71930106513 1800 0.0094010 01101 676218437 665827346
13 1101041002 01101 524 2017 Iquique 375676.9 2017 1101 191468 71930106513 2538 0.0132555 01101 953467996 814119260
14 1101041003 01101 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513 3855 0.0201339 01101 1448234486 1483529987
15 1101041004 01101 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513 5663 0.0295767 01101 2127458339 1941123124
16 1101041005 01101 894 2017 Iquique 375676.9 2017 1101 191468 71930106513 4162 0.0217373 01101 1563567297 1302970045
17 1101041006 01101 588 2017 Iquique 375676.9 2017 1101 191468 71930106513 2689 0.0140441 01101 1010195210 901043999
18 1101051001 01101 857 2017 Iquique 375676.9 2017 1101 191468 71930106513 3296 0.0172144 01101 1238231094 1255376908
19 1101051002 01101 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513 4465 0.0233198 01101 1677397401 1785902480
20 1101051003 01101 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513 4656 0.0243174 01101 1749151691 1752670346
21 1101051004 01101 558 2017 Iquique 375676.9 2017 1101 191468 71930106513 2097 0.0109522 01101 787794479 860446979
22 1101051005 01101 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513 3569 0.0186402 01101 1340790890 1697091595
23 1101051006 01101 858 2017 Iquique 375676.9 2017 1101 191468 71930106513 2741 0.0143157 01101 1029730409 1256666398
24 1101061001 01101 606 2017 Iquique 375676.9 2017 1101 191468 71930106513 1625 0.0084871 01101 610474978 925282670
25 1101061002 01101 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513 4767 0.0248971 01101 1790851828 1858231732
26 1101061003 01101 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513 4826 0.0252053 01101 1813016765 1977539223
27 1101061004 01101 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513 4077 0.0212934 01101 1531634760 1804328423
28 1101061005 01101 626 2017 Iquique 375676.9 2017 1101 191468 71930106513 2166 0.0113126 01101 813716186 952113779
29 1101071001 01101 691 2017 Iquique 375676.9 2017 1101 191468 71930106513 2324 0.0121378 01101 873073138 1038625743
30 1101071002 01101 753 2017 Iquique 375676.9 2017 1101 191468 71930106513 2801 0.0146291 01101 1052271024 1120239792
31 1101071003 01101 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513 3829 0.0199981 01101 1438466887 1501164760
32 1101071004 01101 605 2017 Iquique 375676.9 2017 1101 191468 71930106513 1987 0.0103777 01101 746470019 923938360
33 1101081001 01101 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513 5133 0.0268087 01101 1928349577 1909487281
34 1101081002 01101 852 2017 Iquique 375676.9 2017 1101 191468 71930106513 3233 0.0168853 01101 1214563448 1248926756
35 1101081003 01101 658 2017 Iquique 375676.9 2017 1101 191468 71930106513 2122 0.0110828 01101 797186402 994832066
36 1101081004 01101 762 2017 Iquique 375676.9 2017 1101 191468 71930106513 2392 0.0124929 01101 898619168 1132018673
37 1101101001 01101 594 2017 Iquique 375676.9 2017 1101 191468 71930106513 2672 0.0139553 01101 1003808702 909133290
38 1101101002 01101 993 2017 Iquique 375676.9 2017 1101 191468 71930106513 4398 0.0229699 01101 1652227048 1429185699
39 1101101003 01101 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513 4524 0.0236280 01101 1699562339 1503681709
40 1101101004 01101 838 2017 Iquique 375676.9 2017 1101 191468 71930106513 3544 0.0185096 01101 1331398967 1230842130
41 1101101005 01101 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513 4911 0.0256492 01101 1844949303 1706990130
42 1101101006 01101 940 2017 Iquique 375676.9 2017 1101 191468 71930106513 3688 0.0192617 01101 1385496442 1361813158
43 1101111001 01101 869 2017 Iquique 375676.9 2017 1101 191468 71930106513 3886 0.0202958 01101 1459880470 1270838965
44 1101111002 01101 578 2017 Iquique 375676.9 2017 1101 191468 71930106513 2312 0.0120751 01101 868565015 887539828
45 1101111003 01101 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513 4874 0.0254560 01101 1831049257 1740340468
46 1101111004 01101 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513 4543 0.0237272 01101 1706700200 1536350385
47 1101111005 01101 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513 4331 0.0226200 01101 1627056695 1731702489
48 1101111006 01101 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513 3253 0.0169898 01101 1222076987 1487311254
49 1101111007 01101 827 2017 Iquique 375676.9 2017 1101 191468 71930106513 4639 0.0242286 01101 1742765183 1216607427
50 1101111008 01101 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513 4881 0.0254925 01101 1833678995 1875338016
51 1101111009 01101 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513 5006 0.0261454 01101 1880638609 1865565606
52 1101111010 01101 97 2017 Iquique 375676.9 2017 1101 191468 71930106513 366 0.0019115 01101 137497749 184408155
53 1101111011 01101 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513 4351 0.0227244 01101 1634570233 1844776059
54 1101111012 01101 820 2017 Iquique 375676.9 2017 1101 191468 71930106513 2926 0.0152819 01101 1099230637 1207537197
55 1101111013 01101 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513 3390 0.0177053 01101 1273544723 1630068077
56 1101111014 01101 842 2017 Iquique 375676.9 2017 1101 191468 71930106513 2940 0.0153550 01101 1104490114 1236012829
57 1101991999 01101 194 2017 Iquique 375676.9 2017 1101 191468 71930106513 1062 0.0055466 01101 398968878 339462270
58 1107011001 01107 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4104 0.0378685 01107 1278690352 1237305044
59 1107011002 01107 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4360 0.0402307 01107 1358452713 1290130683
60 1107011003 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 8549 0.0788835 01107 2663626661 1625088626
61 1107021001 01107 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6701 0.0618316 01107 2087842117 1603902718
62 1107021002 01107 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3971 0.0366413 01107 1237251313 1185470726
63 1107021003 01107 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6349 0.0585836 01107 1978168870 2296673977
64 1107021004 01107 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5125 0.0472895 01107 1596805081 1482269273
65 1107021005 01107 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4451 0.0410704 01107 1386805740 1345219244
66 1107021006 01107 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3864 0.0356540 01107 1203913138 1260533788
67 1107021007 01107 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5235 0.0483045 01107 1631077971 1661143138
68 1107021008 01107 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4566 0.0421315 01107 1422636488 1484790555
69 1107031001 01107 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4195 0.0387082 01107 1307043379 1757598993
70 1107031002 01107 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 7099 0.0655040 01107 2211847663 1847223535
71 1107031003 01107 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4720 0.0435525 01107 1470618533 1669829952
72 1107041001 01107 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3630 0.0334948 01107 1131005355 1075587735
73 1107041002 01107 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5358 0.0494394 01107 1669401293 1620107101
74 1107041003 01107 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4520 0.0417070 01107 1408304189 925282670
75 1107041004 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5317 0.0490611 01107 1656626852 1625088626
76 1107041005 01107 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3725 0.0343714 01107 1160604669 1132018673
77 1107041006 01107 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4009 0.0369919 01107 1249091038 1095317630
78 1107041007 01107 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5217 0.0481384 01107 1625469680 1243763373
79 1107991999 01107 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 819 0.0075571 01107 255177241 289648909
80 1401011001 01401 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 2771 0.1763732 01401 876019768 985509100
81 1401011002 01401 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 6506 0.4141048 01401 2056797044 1915576595
82 1401991999 01401 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 818 0.0520654 01401 258601288 47944876
83 1405011001 01405 895 2017 Pica 330061.1 2017 1405 9296 3068247619 3876 0.4169535 01405 1279316671 1304253033
84 1405991999 01405 14 2017 Pica 330061.1 2017 1405 9296 3068247619 36 0.0038726 01405 11882198 33552159
NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.1 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.3 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.4 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.5 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.6 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.7 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.8 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.9 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.10 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.11 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.12 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.13 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.14 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.15 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA


11 División del valor estimado entre la población total de la zona para obtener el ingreso medio por zona


\[ Ingreso \_ Medio\_zona = est\_ing / (personas * p\_poblacional) \]


h_y_m_comuna_corr$ing_medio_zona <- h_y_m_comuna_corr$est_ing  /( h_y_m_comuna_corr$personas  * h_y_m_comuna_corr$p_poblacional)

r3_100 <- h_y_m_comuna_corr[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")
zona código.x Freq.x anio comuna.x promedio_i año comuna.y personas Ingresos_expandidos Freq.y p_poblacional código.y multi_pob est_ing ing_medio_zona
1 1101011001 01101 807 2017 Iquique 375676.9 2017 1101 191468 71930106513 2491 0.0130100 01101 935811182 1190667824 477987.89
2 1101011002 01101 452 2017 Iquique 375676.9 2017 1101 191468 71930106513 1475 0.0077036 01101 554123442 714785823 484600.56
3 1101021001 01101 257 2017 Iquique 375676.9 2017 1101 191468 71930106513 1003 0.0052385 01101 376803940 434820223 433519.66
4 1101021002 01101 11 2017 Iquique 375676.9 2017 1101 191468 71930106513 54 0.0002820 01101 20286553 27134188 502484.97
5 1101021003 01101 571 2017 Iquique 375676.9 2017 1101 191468 71930106513 2895 0.0151200 01101 1087584653 878070290 303305.80
6 1101021004 01101 451 2017 Iquique 375676.9 2017 1101 191468 71930106513 2398 0.0125243 01101 900873229 713393469 297495.19
7 1101021005 01101 1046 2017 Iquique 375676.9 2017 1101 191468 71930106513 4525 0.0236332 01101 1699938016 1496129140 330636.27
8 1101031001 01101 653 2017 Iquique 375676.9 2017 1101 191468 71930106513 2725 0.0142321 01101 1023719578 988174023 362632.67
9 1101031002 01101 837 2017 Iquique 375676.9 2017 1101 191468 71930106513 3554 0.0185618 01101 1335155736 1229548994 345962.01
10 1101031003 01101 1315 2017 Iquique 375676.9 2017 1101 191468 71930106513 5246 0.0273988 01101 1970801067 1830081910 348852.82
11 1101031004 01101 903 2017 Iquique 375676.9 2017 1101 191468 71930106513 3389 0.0177001 01101 1273169046 1314510781 387875.71
12 1101041001 01101 417 2017 Iquique 375676.9 2017 1101 191468 71930106513 1800 0.0094010 01101 676218437 665827346 369904.08
13 1101041002 01101 524 2017 Iquique 375676.9 2017 1101 191468 71930106513 2538 0.0132555 01101 953467996 814119260 320771.97
14 1101041003 01101 1036 2017 Iquique 375676.9 2017 1101 191468 71930106513 3855 0.0201339 01101 1448234486 1483529987 384832.68
15 1101041004 01101 1406 2017 Iquique 375676.9 2017 1101 191468 71930106513 5663 0.0295767 01101 2127458339 1941123124 342772.93
16 1101041005 01101 894 2017 Iquique 375676.9 2017 1101 191468 71930106513 4162 0.0217373 01101 1563567297 1302970045 313063.44
17 1101041006 01101 588 2017 Iquique 375676.9 2017 1101 191468 71930106513 2689 0.0140441 01101 1010195210 901043999 335085.16
18 1101051001 01101 857 2017 Iquique 375676.9 2017 1101 191468 71930106513 3296 0.0172144 01101 1238231094 1255376908 380878.92
19 1101051002 01101 1279 2017 Iquique 375676.9 2017 1101 191468 71930106513 4465 0.0233198 01101 1677397401 1785902480 399978.16
20 1101051003 01101 1252 2017 Iquique 375676.9 2017 1101 191468 71930106513 4656 0.0243174 01101 1749151691 1752670346 376432.63
21 1101051004 01101 558 2017 Iquique 375676.9 2017 1101 191468 71930106513 2097 0.0109522 01101 787794479 860446979 410322.83
22 1101051005 01101 1207 2017 Iquique 375676.9 2017 1101 191468 71930106513 3569 0.0186402 01101 1340790890 1697091595 475508.99
23 1101051006 01101 858 2017 Iquique 375676.9 2017 1101 191468 71930106513 2741 0.0143157 01101 1029730409 1256666398 458470.05
24 1101061001 01101 606 2017 Iquique 375676.9 2017 1101 191468 71930106513 1625 0.0084871 01101 610474978 925282670 569404.72
25 1101061002 01101 1338 2017 Iquique 375676.9 2017 1101 191468 71930106513 4767 0.0248971 01101 1790851828 1858231732 389811.57
26 1101061003 01101 1436 2017 Iquique 375676.9 2017 1101 191468 71930106513 4826 0.0252053 01101 1813016765 1977539223 409767.76
27 1101061004 01101 1294 2017 Iquique 375676.9 2017 1101 191468 71930106513 4077 0.0212934 01101 1531634760 1804328423 442562.77
28 1101061005 01101 626 2017 Iquique 375676.9 2017 1101 191468 71930106513 2166 0.0113126 01101 813716186 952113779 439572.38
29 1101071001 01101 691 2017 Iquique 375676.9 2017 1101 191468 71930106513 2324 0.0121378 01101 873073138 1038625743 446912.97
30 1101071002 01101 753 2017 Iquique 375676.9 2017 1101 191468 71930106513 2801 0.0146291 01101 1052271024 1120239792 399942.80
31 1101071003 01101 1050 2017 Iquique 375676.9 2017 1101 191468 71930106513 3829 0.0199981 01101 1438466887 1501164760 392051.39
32 1101071004 01101 605 2017 Iquique 375676.9 2017 1101 191468 71930106513 1987 0.0103777 01101 746470019 923938360 464991.63
33 1101081001 01101 1380 2017 Iquique 375676.9 2017 1101 191468 71930106513 5133 0.0268087 01101 1928349577 1909487281 372002.20
34 1101081002 01101 852 2017 Iquique 375676.9 2017 1101 191468 71930106513 3233 0.0168853 01101 1214563448 1248926756 386305.83
35 1101081003 01101 658 2017 Iquique 375676.9 2017 1101 191468 71930106513 2122 0.0110828 01101 797186402 994832066 468818.13
36 1101081004 01101 762 2017 Iquique 375676.9 2017 1101 191468 71930106513 2392 0.0124929 01101 898619168 1132018673 473251.95
37 1101101001 01101 594 2017 Iquique 375676.9 2017 1101 191468 71930106513 2672 0.0139553 01101 1003808702 909133290 340244.49
38 1101101002 01101 993 2017 Iquique 375676.9 2017 1101 191468 71930106513 4398 0.0229699 01101 1652227048 1429185699 324962.64
39 1101101003 01101 1052 2017 Iquique 375676.9 2017 1101 191468 71930106513 4524 0.0236280 01101 1699562339 1503681709 332378.80
40 1101101004 01101 838 2017 Iquique 375676.9 2017 1101 191468 71930106513 3544 0.0185096 01101 1331398967 1230842130 347303.08
41 1101101005 01101 1215 2017 Iquique 375676.9 2017 1101 191468 71930106513 4911 0.0256492 01101 1844949303 1706990130 347585.04
42 1101101006 01101 940 2017 Iquique 375676.9 2017 1101 191468 71930106513 3688 0.0192617 01101 1385496442 1361813158 369255.19
43 1101111001 01101 869 2017 Iquique 375676.9 2017 1101 191468 71930106513 3886 0.0202958 01101 1459880470 1270838965 327030.10
44 1101111002 01101 578 2017 Iquique 375676.9 2017 1101 191468 71930106513 2312 0.0120751 01101 868565015 887539828 383884.01
45 1101111003 01101 1242 2017 Iquique 375676.9 2017 1101 191468 71930106513 4874 0.0254560 01101 1831049257 1740340468 357066.16
46 1101111004 01101 1078 2017 Iquique 375676.9 2017 1101 191468 71930106513 4543 0.0237272 01101 1706700200 1536350385 338179.70
47 1101111005 01101 1235 2017 Iquique 375676.9 2017 1101 191468 71930106513 4331 0.0226200 01101 1627056695 1731702489 399838.95
48 1101111006 01101 1039 2017 Iquique 375676.9 2017 1101 191468 71930106513 3253 0.0169898 01101 1222076987 1487311254 457212.19
49 1101111007 01101 827 2017 Iquique 375676.9 2017 1101 191468 71930106513 4639 0.0242286 01101 1742765183 1216607427 262256.40
50 1101111008 01101 1352 2017 Iquique 375676.9 2017 1101 191468 71930106513 4881 0.0254925 01101 1833678995 1875338016 384211.85
51 1101111009 01101 1344 2017 Iquique 375676.9 2017 1101 191468 71930106513 5006 0.0261454 01101 1880638609 1865565606 372665.92
52 1101111010 01101 97 2017 Iquique 375676.9 2017 1101 191468 71930106513 366 0.0019115 01101 137497749 184408155 503847.42
53 1101111011 01101 1327 2017 Iquique 375676.9 2017 1101 191468 71930106513 4351 0.0227244 01101 1634570233 1844776059 423988.98
54 1101111012 01101 820 2017 Iquique 375676.9 2017 1101 191468 71930106513 2926 0.0152819 01101 1099230637 1207537197 412692.14
55 1101111013 01101 1153 2017 Iquique 375676.9 2017 1101 191468 71930106513 3390 0.0177053 01101 1273544723 1630068077 480846.04
56 1101111014 01101 842 2017 Iquique 375676.9 2017 1101 191468 71930106513 2940 0.0153550 01101 1104490114 1236012829 420412.53
57 1101991999 01101 194 2017 Iquique 375676.9 2017 1101 191468 71930106513 1062 0.0055466 01101 398968878 339462270 319644.32
58 1107011001 01107 843 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4104 0.0378685 01107 1278690352 1237305044 301487.58
59 1107011002 01107 884 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4360 0.0402307 01107 1358452713 1290130683 295901.53
60 1107011003 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 8549 0.0788835 01107 2663626661 1625088626 190091.08
61 1107021001 01107 1132 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6701 0.0618316 01107 2087842117 1603902718 239352.74
62 1107021002 01107 803 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3971 0.0366413 01107 1237251313 1185470726 298532.04
63 1107021003 01107 1702 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 6349 0.0585836 01107 1978168870 2296673977 361737.91
64 1107021004 01107 1035 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5125 0.0472895 01107 1596805081 1482269273 289223.27
65 1107021005 01107 927 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4451 0.0410704 01107 1386805740 1345219244 302228.54
66 1107021006 01107 861 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3864 0.0356540 01107 1203913138 1260533788 326225.10
67 1107021007 01107 1178 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5235 0.0483045 01107 1631077971 1661143138 317314.83
68 1107021008 01107 1037 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4566 0.0421315 01107 1422636488 1484790555 325184.09
69 1107031001 01107 1256 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4195 0.0387082 01107 1307043379 1757598993 418974.73
70 1107031002 01107 1329 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 7099 0.0655040 01107 2211847663 1847223535 260208.98
71 1107031003 01107 1185 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4720 0.0435525 01107 1470618533 1669829952 353777.53
72 1107041001 01107 719 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3630 0.0334948 01107 1131005355 1075587735 296305.16
73 1107041002 01107 1145 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5358 0.0494394 01107 1669401293 1620107101 302371.61
74 1107041003 01107 606 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4520 0.0417070 01107 1408304189 925282670 204708.56
75 1107041004 01107 1149 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5317 0.0490611 01107 1656626852 1625088626 305640.14
76 1107041005 01107 762 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 3725 0.0343714 01107 1160604669 1132018673 303897.63
77 1107041006 01107 734 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 4009 0.0369919 01107 1249091038 1095317630 273214.67
78 1107041007 01107 848 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 5217 0.0481384 01107 1625469680 1243763373 238405.86
79 1107991999 01107 162 2017 Alto Hospicio 311571.7 2017 1107 108375 33766585496 819 0.0075571 01107 255177241 289648909 353661.67
80 1401011001 01401 651 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 2771 0.1763732 01401 876019768 985509100 355651.06
81 1401011002 01401 1385 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 6506 0.4141048 01401 2056797044 1915576595 294432.31
82 1401991999 01401 21 2017 Pozo Almonte 316138.5 2017 1401 15711 4966851883 818 0.0520654 01401 258601288 47944876 58612.32
83 1405011001 01405 895 2017 Pica 330061.1 2017 1405 9296 3068247619 3876 0.4169535 01405 1279316671 1304253033 336494.59
84 1405991999 01405 14 2017 Pica 330061.1 2017 1405 9296 3068247619 36 0.0038726 01405 11882198 33552159 932004.42
NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.1 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.2 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.3 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.4 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.5 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.6 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.7 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.8 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.9 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.10 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.11 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.12 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.13 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.14 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.15 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA


Guardamos:

saveRDS(h_y_m_comuna_corr, "P03C/region_01_P03C_u.rds")