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 == 15) 
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 15101011001 1 15101 353 2017
2 15101011002 1 15101 261 2017
3 15101021001 1 15101 1234 2017
4 15101031001 1 15101 2220 2017
5 15101031002 1 15101 1007 2017
6 15101031003 1 15101 162 2017
7 15101031004 1 15101 1102 2017
8 15101031005 1 15101 235 2017
9 15101031006 1 15101 866 2017
10 15101031007 1 15101 824 2017
11 15101041001 1 15101 638 2017
12 15101041002 1 15101 702 2017
13 15101041003 1 15101 665 2017
14 15101041004 1 15101 691 2017
15 15101051001 1 15101 310 2017
16 15101051002 1 15101 327 2017
17 15101051003 1 15101 438 2017
18 15101051004 1 15101 415 2017
19 15101061001 1 15101 965 2017
20 15101061002 1 15101 586 2017
21 15101061003 1 15101 126 2017
22 15101061004 1 15101 137 2017
23 15101061005 1 15101 160 2017
24 15101071001 1 15101 507 2017
25 15101071002 1 15101 632 2017
26 15101071003 1 15101 767 2017
27 15101071004 1 15101 455 2017
28 15101081001 1 15101 539 2017
29 15101081002 1 15101 537 2017
30 15101081003 1 15101 790 2017
31 15101081004 1 15101 445 2017
32 15101091001 1 15101 565 2017
33 15101091002 1 15101 360 2017
34 15101101001 1 15101 430 2017
35 15101101002 1 15101 363 2017
36 15101101003 1 15101 339 2017
37 15101111001 1 15101 462 2017
38 15101111002 1 15101 454 2017
39 15101111003 1 15101 467 2017
40 15101121001 1 15101 737 2017
41 15101121002 1 15101 1263 2017
42 15101121003 1 15101 1083 2017
43 15101121004 1 15101 914 2017
44 15101121005 1 15101 1122 2017
45 15101121006 1 15101 64 2017
46 15101121007 1 15101 1177 2017
47 15101121008 1 15101 893 2017
48 15101121009 1 15101 623 2017
49 15101121010 1 15101 943 2017
50 15101141001 1 15101 318 2017
51 15101141002 1 15101 412 2017
52 15101171001 1 15101 980 2017
53 15101171002 1 15101 575 2017
54 15101171003 1 15101 719 2017
55 15101171004 1 15101 609 2017
56 15101171005 1 15101 559 2017
57 15101171006 1 15101 1117 2017
58 15101171007 1 15101 610 2017
59 15101171008 1 15101 541 2017
60 15101171009 1 15101 715 2017
61 15101171010 1 15101 961 2017
62 15101181001 1 15101 794 2017
63 15101181002 1 15101 1042 2017
64 15101181003 1 15101 1207 2017
65 15101181004 1 15101 643 2017
66 15101181005 1 15101 917 2017
67 15101181006 1 15101 1516 2017
68 15101191001 1 15101 517 2017
69 15101191002 1 15101 705 2017
70 15101191003 1 15101 746 2017
71 15101191004 1 15101 585 2017
72 15101991999 1 15101 100 2017
146 15201011001 1 15201 150 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
NA.14 NA NA NA NA NA
NA.15 NA NA NA NA NA
NA.16 NA NA NA NA NA
NA.17 NA NA NA NA NA
NA.18 NA NA NA NA NA
NA.19 NA NA NA NA NA
NA.20 NA NA NA NA NA
NA.21 NA NA NA NA NA
NA.22 NA NA NA NA NA
NA.23 NA NA NA NA NA
NA.24 NA NA NA NA NA
NA.25 NA NA NA NA NA
NA.26 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 15101011001 353 2017 15101
2 15101011002 261 2017 15101
3 15101021001 1234 2017 15101
4 15101031001 2220 2017 15101
5 15101031002 1007 2017 15101
6 15101031003 162 2017 15101
7 15101031004 1102 2017 15101
8 15101031005 235 2017 15101
9 15101031006 866 2017 15101
10 15101031007 824 2017 15101
11 15101041001 638 2017 15101
12 15101041002 702 2017 15101
13 15101041003 665 2017 15101
14 15101041004 691 2017 15101
15 15101051001 310 2017 15101
16 15101051002 327 2017 15101
17 15101051003 438 2017 15101
18 15101051004 415 2017 15101
19 15101061001 965 2017 15101
20 15101061002 586 2017 15101
21 15101061003 126 2017 15101
22 15101061004 137 2017 15101
23 15101061005 160 2017 15101
24 15101071001 507 2017 15101
25 15101071002 632 2017 15101
26 15101071003 767 2017 15101
27 15101071004 455 2017 15101
28 15101081001 539 2017 15101
29 15101081002 537 2017 15101
30 15101081003 790 2017 15101
31 15101081004 445 2017 15101
32 15101091001 565 2017 15101
33 15101091002 360 2017 15101
34 15101101001 430 2017 15101
35 15101101002 363 2017 15101
36 15101101003 339 2017 15101
37 15101111001 462 2017 15101
38 15101111002 454 2017 15101
39 15101111003 467 2017 15101
40 15101121001 737 2017 15101
41 15101121002 1263 2017 15101
42 15101121003 1083 2017 15101
43 15101121004 914 2017 15101
44 15101121005 1122 2017 15101
45 15101121006 64 2017 15101
46 15101121007 1177 2017 15101
47 15101121008 893 2017 15101
48 15101121009 623 2017 15101
49 15101121010 943 2017 15101
50 15101141001 318 2017 15101
51 15101141002 412 2017 15101
52 15101171001 980 2017 15101
53 15101171002 575 2017 15101
54 15101171003 719 2017 15101
55 15101171004 609 2017 15101
56 15101171005 559 2017 15101
57 15101171006 1117 2017 15101
58 15101171007 610 2017 15101
59 15101171008 541 2017 15101
60 15101171009 715 2017 15101
61 15101171010 961 2017 15101
62 15101181001 794 2017 15101
63 15101181002 1042 2017 15101
64 15101181003 1207 2017 15101
65 15101181004 643 2017 15101
66 15101181005 917 2017 15101
67 15101181006 1516 2017 15101
68 15101191001 517 2017 15101
69 15101191002 705 2017 15101
70 15101191003 746 2017 15101
71 15101191004 585 2017 15101
72 15101991999 100 2017 15101
146 15201011001 150 2017 15201
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
NA.14 NA NA NA NA
NA.15 NA NA NA NA
NA.16 NA NA NA NA
NA.17 NA NA NA NA
NA.18 NA NA NA NA
NA.19 NA NA NA NA
NA.20 NA NA NA NA
NA.21 NA NA NA NA
NA.22 NA NA NA NA
NA.23 NA NA NA NA
NA.24 NA NA NA NA
NA.25 NA NA NA NA
NA.26 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 15101 15101011001 353 2017 Arica 310013.3 2017 15101 221364 68625788545
2 15101 15101011002 261 2017 Arica 310013.3 2017 15101 221364 68625788545
3 15101 15101021001 1234 2017 Arica 310013.3 2017 15101 221364 68625788545
4 15101 15101031001 2220 2017 Arica 310013.3 2017 15101 221364 68625788545
5 15101 15101031002 1007 2017 Arica 310013.3 2017 15101 221364 68625788545
6 15101 15101031003 162 2017 Arica 310013.3 2017 15101 221364 68625788545
7 15101 15101031004 1102 2017 Arica 310013.3 2017 15101 221364 68625788545
8 15101 15101031005 235 2017 Arica 310013.3 2017 15101 221364 68625788545
9 15101 15101031006 866 2017 Arica 310013.3 2017 15101 221364 68625788545
10 15101 15101031007 824 2017 Arica 310013.3 2017 15101 221364 68625788545
11 15101 15101041001 638 2017 Arica 310013.3 2017 15101 221364 68625788545
12 15101 15101041002 702 2017 Arica 310013.3 2017 15101 221364 68625788545
13 15101 15101041003 665 2017 Arica 310013.3 2017 15101 221364 68625788545
14 15101 15101041004 691 2017 Arica 310013.3 2017 15101 221364 68625788545
15 15101 15101051001 310 2017 Arica 310013.3 2017 15101 221364 68625788545
16 15101 15101051002 327 2017 Arica 310013.3 2017 15101 221364 68625788545
17 15101 15101051003 438 2017 Arica 310013.3 2017 15101 221364 68625788545
18 15101 15101051004 415 2017 Arica 310013.3 2017 15101 221364 68625788545
19 15101 15101061001 965 2017 Arica 310013.3 2017 15101 221364 68625788545
20 15101 15101061002 586 2017 Arica 310013.3 2017 15101 221364 68625788545
21 15101 15101061003 126 2017 Arica 310013.3 2017 15101 221364 68625788545
22 15101 15101061004 137 2017 Arica 310013.3 2017 15101 221364 68625788545
23 15101 15101061005 160 2017 Arica 310013.3 2017 15101 221364 68625788545
24 15101 15101071001 507 2017 Arica 310013.3 2017 15101 221364 68625788545
25 15101 15101071002 632 2017 Arica 310013.3 2017 15101 221364 68625788545
26 15101 15101071003 767 2017 Arica 310013.3 2017 15101 221364 68625788545
27 15101 15101071004 455 2017 Arica 310013.3 2017 15101 221364 68625788545
28 15101 15101081001 539 2017 Arica 310013.3 2017 15101 221364 68625788545
29 15101 15101081002 537 2017 Arica 310013.3 2017 15101 221364 68625788545
30 15101 15101081003 790 2017 Arica 310013.3 2017 15101 221364 68625788545
31 15101 15101081004 445 2017 Arica 310013.3 2017 15101 221364 68625788545
32 15101 15101091001 565 2017 Arica 310013.3 2017 15101 221364 68625788545
33 15101 15101091002 360 2017 Arica 310013.3 2017 15101 221364 68625788545
34 15101 15101101001 430 2017 Arica 310013.3 2017 15101 221364 68625788545
35 15101 15101101002 363 2017 Arica 310013.3 2017 15101 221364 68625788545
36 15101 15101101003 339 2017 Arica 310013.3 2017 15101 221364 68625788545
37 15101 15101111001 462 2017 Arica 310013.3 2017 15101 221364 68625788545
38 15101 15101111002 454 2017 Arica 310013.3 2017 15101 221364 68625788545
39 15101 15101111003 467 2017 Arica 310013.3 2017 15101 221364 68625788545
40 15101 15101121001 737 2017 Arica 310013.3 2017 15101 221364 68625788545
41 15101 15101121002 1263 2017 Arica 310013.3 2017 15101 221364 68625788545
42 15101 15101121003 1083 2017 Arica 310013.3 2017 15101 221364 68625788545
43 15101 15101121004 914 2017 Arica 310013.3 2017 15101 221364 68625788545
44 15101 15101121005 1122 2017 Arica 310013.3 2017 15101 221364 68625788545
45 15101 15101121006 64 2017 Arica 310013.3 2017 15101 221364 68625788545
46 15101 15101121007 1177 2017 Arica 310013.3 2017 15101 221364 68625788545
47 15101 15101121008 893 2017 Arica 310013.3 2017 15101 221364 68625788545
48 15101 15101121009 623 2017 Arica 310013.3 2017 15101 221364 68625788545
49 15101 15101121010 943 2017 Arica 310013.3 2017 15101 221364 68625788545
50 15101 15101141001 318 2017 Arica 310013.3 2017 15101 221364 68625788545
51 15101 15101141002 412 2017 Arica 310013.3 2017 15101 221364 68625788545
52 15101 15101171001 980 2017 Arica 310013.3 2017 15101 221364 68625788545
53 15101 15101171002 575 2017 Arica 310013.3 2017 15101 221364 68625788545
54 15101 15101171003 719 2017 Arica 310013.3 2017 15101 221364 68625788545
55 15101 15101171004 609 2017 Arica 310013.3 2017 15101 221364 68625788545
56 15101 15101171005 559 2017 Arica 310013.3 2017 15101 221364 68625788545
57 15101 15101171006 1117 2017 Arica 310013.3 2017 15101 221364 68625788545
58 15101 15101171007 610 2017 Arica 310013.3 2017 15101 221364 68625788545
59 15101 15101171008 541 2017 Arica 310013.3 2017 15101 221364 68625788545
60 15101 15101171009 715 2017 Arica 310013.3 2017 15101 221364 68625788545
61 15101 15101171010 961 2017 Arica 310013.3 2017 15101 221364 68625788545
62 15101 15101181001 794 2017 Arica 310013.3 2017 15101 221364 68625788545
63 15101 15101181002 1042 2017 Arica 310013.3 2017 15101 221364 68625788545
64 15101 15101181003 1207 2017 Arica 310013.3 2017 15101 221364 68625788545
65 15101 15101181004 643 2017 Arica 310013.3 2017 15101 221364 68625788545
66 15101 15101181005 917 2017 Arica 310013.3 2017 15101 221364 68625788545
67 15101 15101181006 1516 2017 Arica 310013.3 2017 15101 221364 68625788545
68 15101 15101191001 517 2017 Arica 310013.3 2017 15101 221364 68625788545
69 15101 15101191002 705 2017 Arica 310013.3 2017 15101 221364 68625788545
70 15101 15101191003 746 2017 Arica 310013.3 2017 15101 221364 68625788545
71 15101 15101191004 585 2017 Arica 310013.3 2017 15101 221364 68625788545
72 15101 15101991999 100 2017 Arica 310013.3 2017 15101 221364 68625788545
73 15201 15201011001 150 2017 Putre 283661.5 2017 15201 2765 784324030
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
NA.16 NA NA NA NA NA NA NA NA NA NA
NA.17 NA NA NA NA NA NA NA NA NA NA
NA.18 NA NA NA NA NA NA NA NA NA NA
NA.19 NA NA NA NA NA NA NA NA NA NA
NA.20 NA NA NA NA NA NA NA NA NA NA
NA.21 NA NA NA NA NA NA NA NA NA NA
NA.22 NA NA NA NA NA NA NA NA NA NA
NA.23 NA NA NA NA NA NA NA NA NA NA
NA.24 NA NA NA NA NA NA NA NA NA NA
NA.25 NA NA NA NA NA NA NA NA NA NA
NA.26 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 15101 15101011001 353 2017 Arica 310013.3 2017 15101 221364 68625788545
2 15101 15101011002 261 2017 Arica 310013.3 2017 15101 221364 68625788545
3 15101 15101021001 1234 2017 Arica 310013.3 2017 15101 221364 68625788545
4 15101 15101031001 2220 2017 Arica 310013.3 2017 15101 221364 68625788545
5 15101 15101031002 1007 2017 Arica 310013.3 2017 15101 221364 68625788545
6 15101 15101031003 162 2017 Arica 310013.3 2017 15101 221364 68625788545
7 15101 15101031004 1102 2017 Arica 310013.3 2017 15101 221364 68625788545
8 15101 15101031005 235 2017 Arica 310013.3 2017 15101 221364 68625788545
9 15101 15101031006 866 2017 Arica 310013.3 2017 15101 221364 68625788545
10 15101 15101031007 824 2017 Arica 310013.3 2017 15101 221364 68625788545
11 15101 15101041001 638 2017 Arica 310013.3 2017 15101 221364 68625788545
12 15101 15101041002 702 2017 Arica 310013.3 2017 15101 221364 68625788545
13 15101 15101041003 665 2017 Arica 310013.3 2017 15101 221364 68625788545
14 15101 15101041004 691 2017 Arica 310013.3 2017 15101 221364 68625788545
15 15101 15101051001 310 2017 Arica 310013.3 2017 15101 221364 68625788545
16 15101 15101051002 327 2017 Arica 310013.3 2017 15101 221364 68625788545
17 15101 15101051003 438 2017 Arica 310013.3 2017 15101 221364 68625788545
18 15101 15101051004 415 2017 Arica 310013.3 2017 15101 221364 68625788545
19 15101 15101061001 965 2017 Arica 310013.3 2017 15101 221364 68625788545
20 15101 15101061002 586 2017 Arica 310013.3 2017 15101 221364 68625788545
21 15101 15101061003 126 2017 Arica 310013.3 2017 15101 221364 68625788545
22 15101 15101061004 137 2017 Arica 310013.3 2017 15101 221364 68625788545
23 15101 15101061005 160 2017 Arica 310013.3 2017 15101 221364 68625788545
24 15101 15101071001 507 2017 Arica 310013.3 2017 15101 221364 68625788545
25 15101 15101071002 632 2017 Arica 310013.3 2017 15101 221364 68625788545
26 15101 15101071003 767 2017 Arica 310013.3 2017 15101 221364 68625788545
27 15101 15101071004 455 2017 Arica 310013.3 2017 15101 221364 68625788545
28 15101 15101081001 539 2017 Arica 310013.3 2017 15101 221364 68625788545
29 15101 15101081002 537 2017 Arica 310013.3 2017 15101 221364 68625788545
30 15101 15101081003 790 2017 Arica 310013.3 2017 15101 221364 68625788545
31 15101 15101081004 445 2017 Arica 310013.3 2017 15101 221364 68625788545
32 15101 15101091001 565 2017 Arica 310013.3 2017 15101 221364 68625788545
33 15101 15101091002 360 2017 Arica 310013.3 2017 15101 221364 68625788545
34 15101 15101101001 430 2017 Arica 310013.3 2017 15101 221364 68625788545
35 15101 15101101002 363 2017 Arica 310013.3 2017 15101 221364 68625788545
36 15101 15101101003 339 2017 Arica 310013.3 2017 15101 221364 68625788545
37 15101 15101111001 462 2017 Arica 310013.3 2017 15101 221364 68625788545
38 15101 15101111002 454 2017 Arica 310013.3 2017 15101 221364 68625788545
39 15101 15101111003 467 2017 Arica 310013.3 2017 15101 221364 68625788545
40 15101 15101121001 737 2017 Arica 310013.3 2017 15101 221364 68625788545
41 15101 15101121002 1263 2017 Arica 310013.3 2017 15101 221364 68625788545
42 15101 15101121003 1083 2017 Arica 310013.3 2017 15101 221364 68625788545
43 15101 15101121004 914 2017 Arica 310013.3 2017 15101 221364 68625788545
44 15101 15101121005 1122 2017 Arica 310013.3 2017 15101 221364 68625788545
45 15101 15101121006 64 2017 Arica 310013.3 2017 15101 221364 68625788545
46 15101 15101121007 1177 2017 Arica 310013.3 2017 15101 221364 68625788545
47 15101 15101121008 893 2017 Arica 310013.3 2017 15101 221364 68625788545
48 15101 15101121009 623 2017 Arica 310013.3 2017 15101 221364 68625788545
49 15101 15101121010 943 2017 Arica 310013.3 2017 15101 221364 68625788545
50 15101 15101141001 318 2017 Arica 310013.3 2017 15101 221364 68625788545
51 15101 15101141002 412 2017 Arica 310013.3 2017 15101 221364 68625788545
52 15101 15101171001 980 2017 Arica 310013.3 2017 15101 221364 68625788545
53 15101 15101171002 575 2017 Arica 310013.3 2017 15101 221364 68625788545
54 15101 15101171003 719 2017 Arica 310013.3 2017 15101 221364 68625788545
55 15101 15101171004 609 2017 Arica 310013.3 2017 15101 221364 68625788545
56 15101 15101171005 559 2017 Arica 310013.3 2017 15101 221364 68625788545
57 15101 15101171006 1117 2017 Arica 310013.3 2017 15101 221364 68625788545
58 15101 15101171007 610 2017 Arica 310013.3 2017 15101 221364 68625788545
59 15101 15101171008 541 2017 Arica 310013.3 2017 15101 221364 68625788545
60 15101 15101171009 715 2017 Arica 310013.3 2017 15101 221364 68625788545
61 15101 15101171010 961 2017 Arica 310013.3 2017 15101 221364 68625788545
62 15101 15101181001 794 2017 Arica 310013.3 2017 15101 221364 68625788545
63 15101 15101181002 1042 2017 Arica 310013.3 2017 15101 221364 68625788545
64 15101 15101181003 1207 2017 Arica 310013.3 2017 15101 221364 68625788545
65 15101 15101181004 643 2017 Arica 310013.3 2017 15101 221364 68625788545
66 15101 15101181005 917 2017 Arica 310013.3 2017 15101 221364 68625788545
67 15101 15101181006 1516 2017 Arica 310013.3 2017 15101 221364 68625788545
68 15101 15101191001 517 2017 Arica 310013.3 2017 15101 221364 68625788545
69 15101 15101191002 705 2017 Arica 310013.3 2017 15101 221364 68625788545
70 15101 15101191003 746 2017 Arica 310013.3 2017 15101 221364 68625788545
71 15101 15101191004 585 2017 Arica 310013.3 2017 15101 221364 68625788545
72 15101 15101991999 100 2017 Arica 310013.3 2017 15101 221364 68625788545
73 15201 15201011001 150 2017 Putre 283661.5 2017 15201 2765 784324030
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
NA.16 NA NA NA NA NA NA NA NA NA NA
NA.17 NA NA NA NA NA NA NA NA NA NA
NA.18 NA NA NA NA NA NA NA NA NA NA
NA.19 NA NA NA NA NA NA NA NA NA NA
NA.20 NA NA NA NA NA NA NA NA NA NA
NA.21 NA NA NA NA NA NA NA NA NA NA
NA.22 NA NA NA NA NA NA NA NA NA NA
NA.23 NA NA NA NA NA NA NA NA NA NA
NA.24 NA NA NA NA NA NA NA NA NA NA
NA.25 NA NA NA NA NA NA NA NA NA NA
NA.26 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
15101011001 15101 353 2017 Arica 310013.3 2017 15101 221364 68625788545 1545 0.0069795 15101
15101011002 15101 261 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101
15101021001 15101 1234 2017 Arica 310013.3 2017 15101 221364 68625788545 4454 0.0201207 15101
15101031001 15101 2220 2017 Arica 310013.3 2017 15101 221364 68625788545 6208 0.0280443 15101
15101031002 15101 1007 2017 Arica 310013.3 2017 15101 221364 68625788545 3166 0.0143022 15101
15101031003 15101 162 2017 Arica 310013.3 2017 15101 221364 68625788545 623 0.0028144 15101
15101031004 15101 1102 2017 Arica 310013.3 2017 15101 221364 68625788545 4041 0.0182550 15101
15101031005 15101 235 2017 Arica 310013.3 2017 15101 221364 68625788545 1010 0.0045626 15101
15101031006 15101 866 2017 Arica 310013.3 2017 15101 221364 68625788545 2848 0.0128657 15101
15101031007 15101 824 2017 Arica 310013.3 2017 15101 221364 68625788545 2762 0.0124772 15101
15101041001 15101 638 2017 Arica 310013.3 2017 15101 221364 68625788545 2911 0.0131503 15101
15101041002 15101 702 2017 Arica 310013.3 2017 15101 221364 68625788545 2836 0.0128115 15101
15101041003 15101 665 2017 Arica 310013.3 2017 15101 221364 68625788545 3316 0.0149799 15101
15101041004 15101 691 2017 Arica 310013.3 2017 15101 221364 68625788545 2869 0.0129606 15101
15101051001 15101 310 2017 Arica 310013.3 2017 15101 221364 68625788545 2041 0.0092201 15101
15101051002 15101 327 2017 Arica 310013.3 2017 15101 221364 68625788545 1655 0.0074764 15101
15101051003 15101 438 2017 Arica 310013.3 2017 15101 221364 68625788545 2345 0.0105934 15101
15101051004 15101 415 2017 Arica 310013.3 2017 15101 221364 68625788545 2196 0.0099203 15101
15101061001 15101 965 2017 Arica 310013.3 2017 15101 221364 68625788545 5901 0.0266575 15101
15101061002 15101 586 2017 Arica 310013.3 2017 15101 221364 68625788545 2066 0.0093330 15101
15101061003 15101 126 2017 Arica 310013.3 2017 15101 221364 68625788545 615 0.0027782 15101
15101061004 15101 137 2017 Arica 310013.3 2017 15101 221364 68625788545 723 0.0032661 15101
15101061005 15101 160 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101
15101071001 15101 507 2017 Arica 310013.3 2017 15101 221364 68625788545 2416 0.0109142 15101
15101071002 15101 632 2017 Arica 310013.3 2017 15101 221364 68625788545 2443 0.0110361 15101
15101071003 15101 767 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101
15101071004 15101 455 2017 Arica 310013.3 2017 15101 221364 68625788545 2455 0.0110903 15101
15101081001 15101 539 2017 Arica 310013.3 2017 15101 221364 68625788545 2735 0.0123552 15101
15101081002 15101 537 2017 Arica 310013.3 2017 15101 221364 68625788545 2496 0.0112755 15101
15101081003 15101 790 2017 Arica 310013.3 2017 15101 221364 68625788545 2823 0.0127528 15101
15101081004 15101 445 2017 Arica 310013.3 2017 15101 221364 68625788545 1853 0.0083708 15101
15101091001 15101 565 2017 Arica 310013.3 2017 15101 221364 68625788545 2582 0.0116640 15101
15101091002 15101 360 2017 Arica 310013.3 2017 15101 221364 68625788545 1634 0.0073815 15101
15101101001 15101 430 2017 Arica 310013.3 2017 15101 221364 68625788545 1944 0.0087819 15101
15101101002 15101 363 2017 Arica 310013.3 2017 15101 221364 68625788545 1732 0.0078242 15101
15101101003 15101 339 2017 Arica 310013.3 2017 15101 221364 68625788545 1542 0.0069659 15101
15101111001 15101 462 2017 Arica 310013.3 2017 15101 221364 68625788545 1879 0.0084883 15101
15101111002 15101 454 2017 Arica 310013.3 2017 15101 221364 68625788545 1971 0.0089039 15101
15101111003 15101 467 2017 Arica 310013.3 2017 15101 221364 68625788545 1656 0.0074809 15101
15101121001 15101 737 2017 Arica 310013.3 2017 15101 221364 68625788545 3627 0.0163848 15101
15101121002 15101 1263 2017 Arica 310013.3 2017 15101 221364 68625788545 4645 0.0209835 15101
15101121003 15101 1083 2017 Arica 310013.3 2017 15101 221364 68625788545 4990 0.0225421 15101
15101121004 15101 914 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101
15101121005 15101 1122 2017 Arica 310013.3 2017 15101 221364 68625788545 4822 0.0217831 15101
15101121006 15101 64 2017 Arica 310013.3 2017 15101 221364 68625788545 263 0.0011881 15101
15101121007 15101 1177 2017 Arica 310013.3 2017 15101 221364 68625788545 3635 0.0164209 15101
15101121008 15101 893 2017 Arica 310013.3 2017 15101 221364 68625788545 3757 0.0169720 15101
15101121009 15101 623 2017 Arica 310013.3 2017 15101 221364 68625788545 2707 0.0122287 15101
15101121010 15101 943 2017 Arica 310013.3 2017 15101 221364 68625788545 3363 0.0151922 15101
15101141001 15101 318 2017 Arica 310013.3 2017 15101 221364 68625788545 1387 0.0062657 15101
15101141002 15101 412 2017 Arica 310013.3 2017 15101 221364 68625788545 1641 0.0074131 15101
15101171001 15101 980 2017 Arica 310013.3 2017 15101 221364 68625788545 4176 0.0188649 15101
15101171002 15101 575 2017 Arica 310013.3 2017 15101 221364 68625788545 2174 0.0098209 15101
15101171003 15101 719 2017 Arica 310013.3 2017 15101 221364 68625788545 2503 0.0113072 15101
15101171004 15101 609 2017 Arica 310013.3 2017 15101 221364 68625788545 2478 0.0111942 15101
15101171005 15101 559 2017 Arica 310013.3 2017 15101 221364 68625788545 2371 0.0107109 15101
15101171006 15101 1117 2017 Arica 310013.3 2017 15101 221364 68625788545 4924 0.0222439 15101
15101171007 15101 610 2017 Arica 310013.3 2017 15101 221364 68625788545 2627 0.0118673 15101
15101171008 15101 541 2017 Arica 310013.3 2017 15101 221364 68625788545 2763 0.0124817 15101
15101171009 15101 715 2017 Arica 310013.3 2017 15101 221364 68625788545 3821 0.0172612 15101
15101171010 15101 961 2017 Arica 310013.3 2017 15101 221364 68625788545 3487 0.0157523 15101
15101181001 15101 794 2017 Arica 310013.3 2017 15101 221364 68625788545 3763 0.0169992 15101
15101181002 15101 1042 2017 Arica 310013.3 2017 15101 221364 68625788545 4627 0.0209022 15101
15101181003 15101 1207 2017 Arica 310013.3 2017 15101 221364 68625788545 5211 0.0235404 15101
15101181004 15101 643 2017 Arica 310013.3 2017 15101 221364 68625788545 2346 0.0105979 15101
15101181005 15101 917 2017 Arica 310013.3 2017 15101 221364 68625788545 3954 0.0178620 15101
15101181006 15101 1516 2017 Arica 310013.3 2017 15101 221364 68625788545 7629 0.0344636 15101
15101191001 15101 517 2017 Arica 310013.3 2017 15101 221364 68625788545 2482 0.0112123 15101
15101191002 15101 705 2017 Arica 310013.3 2017 15101 221364 68625788545 3364 0.0151967 15101
15101191003 15101 746 2017 Arica 310013.3 2017 15101 221364 68625788545 3067 0.0138550 15101
15101191004 15101 585 2017 Arica 310013.3 2017 15101 221364 68625788545 2862 0.0129289 15101
15101991999 15101 100 2017 Arica 310013.3 2017 15101 221364 68625788545 1947 0.0087955 15101
15201011001 15201 150 2017 Putre 283661.5 2017 15201 2765 784324030 1716 0.6206148 15201


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
15101011001 15101 353 2017 Arica 310013.3 2017 15101 221364 68625788545 1545 0.0069795 15101 478970579
15101011002 15101 261 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333
15101021001 15101 1234 2017 Arica 310013.3 2017 15101 221364 68625788545 4454 0.0201207 15101 1380799327
15101031001 15101 2220 2017 Arica 310013.3 2017 15101 221364 68625788545 6208 0.0280443 15101 1924562690
15101031002 15101 1007 2017 Arica 310013.3 2017 15101 221364 68625788545 3166 0.0143022 15101 981502171
15101031003 15101 162 2017 Arica 310013.3 2017 15101 221364 68625788545 623 0.0028144 15101 193138298
15101031004 15101 1102 2017 Arica 310013.3 2017 15101 221364 68625788545 4041 0.0182550 15101 1252763826
15101031005 15101 235 2017 Arica 310013.3 2017 15101 221364 68625788545 1010 0.0045626 15101 313113453
15101031006 15101 866 2017 Arica 310013.3 2017 15101 221364 68625788545 2848 0.0128657 15101 882917935
15101031007 15101 824 2017 Arica 310013.3 2017 15101 221364 68625788545 2762 0.0124772 15101 856256790
15101041001 15101 638 2017 Arica 310013.3 2017 15101 221364 68625788545 2911 0.0131503 15101 902448774
15101041002 15101 702 2017 Arica 310013.3 2017 15101 221364 68625788545 2836 0.0128115 15101 879197775
15101041003 15101 665 2017 Arica 310013.3 2017 15101 221364 68625788545 3316 0.0149799 15101 1028004169
15101041004 15101 691 2017 Arica 310013.3 2017 15101 221364 68625788545 2869 0.0129606 15101 889428215
15101051001 15101 310 2017 Arica 310013.3 2017 15101 221364 68625788545 2041 0.0092201 15101 632737186
15101051002 15101 327 2017 Arica 310013.3 2017 15101 221364 68625788545 1655 0.0074764 15101 513072044
15101051003 15101 438 2017 Arica 310013.3 2017 15101 221364 68625788545 2345 0.0105934 15101 726981235
15101051004 15101 415 2017 Arica 310013.3 2017 15101 221364 68625788545 2196 0.0099203 15101 680789250
15101061001 15101 965 2017 Arica 310013.3 2017 15101 221364 68625788545 5901 0.0266575 15101 1829388601
15101061002 15101 586 2017 Arica 310013.3 2017 15101 221364 68625788545 2066 0.0093330 15101 640487519
15101061003 15101 126 2017 Arica 310013.3 2017 15101 221364 68625788545 615 0.0027782 15101 190658192
15101061004 15101 137 2017 Arica 310013.3 2017 15101 221364 68625788545 723 0.0032661 15101 224139630
15101061005 15101 160 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333
15101071001 15101 507 2017 Arica 310013.3 2017 15101 221364 68625788545 2416 0.0109142 15101 748992181
15101071002 15101 632 2017 Arica 310013.3 2017 15101 221364 68625788545 2443 0.0110361 15101 757362541
15101071003 15101 767 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097
15101071004 15101 455 2017 Arica 310013.3 2017 15101 221364 68625788545 2455 0.0110903 15101 761082700
15101081001 15101 539 2017 Arica 310013.3 2017 15101 221364 68625788545 2735 0.0123552 15101 847886430
15101081002 15101 537 2017 Arica 310013.3 2017 15101 221364 68625788545 2496 0.0112755 15101 773793246
15101081003 15101 790 2017 Arica 310013.3 2017 15101 221364 68625788545 2823 0.0127528 15101 875167602
15101081004 15101 445 2017 Arica 310013.3 2017 15101 221364 68625788545 1853 0.0083708 15101 574454682
15101091001 15101 565 2017 Arica 310013.3 2017 15101 221364 68625788545 2582 0.0116640 15101 800454392
15101091002 15101 360 2017 Arica 310013.3 2017 15101 221364 68625788545 1634 0.0073815 15101 506561765
15101101001 15101 430 2017 Arica 310013.3 2017 15101 221364 68625788545 1944 0.0087819 15101 602665894
15101101002 15101 363 2017 Arica 310013.3 2017 15101 221364 68625788545 1732 0.0078242 15101 536943070
15101101003 15101 339 2017 Arica 310013.3 2017 15101 221364 68625788545 1542 0.0069659 15101 478040539
15101111001 15101 462 2017 Arica 310013.3 2017 15101 221364 68625788545 1879 0.0084883 15101 582515028
15101111002 15101 454 2017 Arica 310013.3 2017 15101 221364 68625788545 1971 0.0089039 15101 611036254
15101111003 15101 467 2017 Arica 310013.3 2017 15101 221364 68625788545 1656 0.0074809 15101 513382058
15101121001 15101 737 2017 Arica 310013.3 2017 15101 221364 68625788545 3627 0.0163848 15101 1124418311
15101121002 15101 1263 2017 Arica 310013.3 2017 15101 221364 68625788545 4645 0.0209835 15101 1440011871
15101121003 15101 1083 2017 Arica 310013.3 2017 15101 221364 68625788545 4990 0.0225421 15101 1546966466
15101121004 15101 914 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097
15101121005 15101 1122 2017 Arica 310013.3 2017 15101 221364 68625788545 4822 0.0217831 15101 1494884229
15101121006 15101 64 2017 Arica 310013.3 2017 15101 221364 68625788545 263 0.0011881 15101 81533503
15101121007 15101 1177 2017 Arica 310013.3 2017 15101 221364 68625788545 3635 0.0164209 15101 1126898418
15101121008 15101 893 2017 Arica 310013.3 2017 15101 221364 68625788545 3757 0.0169720 15101 1164720043
15101121009 15101 623 2017 Arica 310013.3 2017 15101 221364 68625788545 2707 0.0122287 15101 839206057
15101121010 15101 943 2017 Arica 310013.3 2017 15101 221364 68625788545 3363 0.0151922 15101 1042574795
15101141001 15101 318 2017 Arica 310013.3 2017 15101 221364 68625788545 1387 0.0062657 15101 429988475
15101141002 15101 412 2017 Arica 310013.3 2017 15101 221364 68625788545 1641 0.0074131 15101 508731858
15101171001 15101 980 2017 Arica 310013.3 2017 15101 221364 68625788545 4176 0.0188649 15101 1294615624
15101171002 15101 575 2017 Arica 310013.3 2017 15101 221364 68625788545 2174 0.0098209 15101 673968957
15101171003 15101 719 2017 Arica 310013.3 2017 15101 221364 68625788545 2503 0.0113072 15101 775963340
15101171004 15101 609 2017 Arica 310013.3 2017 15101 221364 68625788545 2478 0.0111942 15101 768213007
15101171005 15101 559 2017 Arica 310013.3 2017 15101 221364 68625788545 2371 0.0107109 15101 735041581
15101171006 15101 1117 2017 Arica 310013.3 2017 15101 221364 68625788545 4924 0.0222439 15101 1526505587
15101171007 15101 610 2017 Arica 310013.3 2017 15101 221364 68625788545 2627 0.0118673 15101 814404991
15101171008 15101 541 2017 Arica 310013.3 2017 15101 221364 68625788545 2763 0.0124817 15101 856566803
15101171009 15101 715 2017 Arica 310013.3 2017 15101 221364 68625788545 3821 0.0172612 15101 1184560895
15101171010 15101 961 2017 Arica 310013.3 2017 15101 221364 68625788545 3487 0.0157523 15101 1081016446
15101181001 15101 794 2017 Arica 310013.3 2017 15101 221364 68625788545 3763 0.0169992 15101 1166580123
15101181002 15101 1042 2017 Arica 310013.3 2017 15101 221364 68625788545 4627 0.0209022 15101 1434431631
15101181003 15101 1207 2017 Arica 310013.3 2017 15101 221364 68625788545 5211 0.0235404 15101 1615479410
15101181004 15101 643 2017 Arica 310013.3 2017 15101 221364 68625788545 2346 0.0105979 15101 727291248
15101181005 15101 917 2017 Arica 310013.3 2017 15101 221364 68625788545 3954 0.0178620 15101 1225792667
15101181006 15101 1516 2017 Arica 310013.3 2017 15101 221364 68625788545 7629 0.0344636 15101 2365091617
15101191001 15101 517 2017 Arica 310013.3 2017 15101 221364 68625788545 2482 0.0112123 15101 769453060
15101191002 15101 705 2017 Arica 310013.3 2017 15101 221364 68625788545 3364 0.0151967 15101 1042884808
15101191003 15101 746 2017 Arica 310013.3 2017 15101 221364 68625788545 3067 0.0138550 15101 950810852
15101191004 15101 585 2017 Arica 310013.3 2017 15101 221364 68625788545 2862 0.0129289 15101 887258122
15101991999 15101 100 2017 Arica 310013.3 2017 15101 221364 68625788545 1947 0.0087955 15101 603595934
15201011001 15201 150 2017 Putre 283661.5 2017 15201 2765 784324030 1716 0.6206148 15201 486763123

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 
## -589664137 -105119054  -25404291  100546723  634283330 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 180800649   42775909   4.227 6.95e-05 ***
## Freq.x        1051093      56478  18.611  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 176800000 on 71 degrees of freedom
## Multiple R-squared:  0.8299, Adjusted R-squared:  0.8275 
## F-statistic: 346.4 on 1 and 71 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.720561980229055 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.807082554558394 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.815246548318143 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.822545914371923 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.827484831376365 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.827484831376365 linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01) \[ ''hat Y = ''beta_0 + ''beta_1 X^3 \]
8 log-log 0.845859747949471 linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01) \[ ''hat Y = e^{''beta_0+''beta_1 ln{X}} \]
5 raíz-raíz 0.857606968067538 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 \] 
h_y_m_comuna_corr <- h_y_m_comuna_corr_01
metodo <- 5


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 = sqrt(multi_pob) ~ sqrt(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5841.2 -1723.5  -310.7  1801.9  9652.3 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   5541.38    1157.74   4.786 8.97e-06 ***
## sqrt(Freq.x)   937.45      44.97  20.848  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2753 on 71 degrees of freedom
## Multiple R-squared:  0.8596, Adjusted R-squared:  0.8576 
## F-statistic: 434.6 on 1 and 71 DF,  p-value: < 2.2e-16
aa <- linearMod$coefficients[1]
aa
## (Intercept) 
##    5541.383
bb <- linearMod$coefficients[2]
bb
## sqrt(Freq.x) 
##     937.4488

9 Modelo raíz-raíz (raíz-raíz)

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

9.1 Diagrama de dispersión sobre raíz-raíz

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

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

9.2 Modelo raíz-raíz

Observemos nuevamente el resultado sobre raíz-raíz.

linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr)
summary(linearMod) 
## 
## Call:
## lm(formula = sqrt(multi_pob) ~ sqrt(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5841.2 -1723.5  -310.7  1801.9  9652.3 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   5541.38    1157.74   4.786 8.97e-06 ***
## sqrt(Freq.x)   937.45      44.97  20.848  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2753 on 71 degrees of freedom
## Multiple R-squared:  0.8596, Adjusted R-squared:  0.8576 
## F-statistic: 434.6 on 1 and 71 DF,  p-value: < 2.2e-16
ggplot(h_y_m_comuna_corr, aes(x = sqrt(Freq.x) , y = sqrt(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 = {-309.2197}^2 + 2 -309.2197 942.8527 \sqrt{X}+ 942.8527^2 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 <- ((aa^2)+(2*(aa*bb)) * sqrt(h_y_m_comuna_corr$Freq.x)+((bb^2)*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 15101011001 15101 353 2017 Arica 310013.3 2017 15101 221364 68625788545 1545 0.0069795 15101 478970579 536128381
2 15101011002 15101 261 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333 427924301
3 15101021001 15101 1234 2017 Arica 310013.3 2017 15101 221364 68625788545 4454 0.0201207 15101 1380799327 1480125454
4 15101031001 15101 2220 2017 Arica 310013.3 2017 15101 221364 68625788545 6208 0.0280443 15101 1924562690 2471187560
5 15101031002 15101 1007 2017 Arica 310013.3 2017 15101 221364 68625788545 3166 0.0143022 15101 981502171 1245362341
6 15101031003 15101 162 2017 Arica 310013.3 2017 15101 221364 68625788545 623 0.0028144 15101 193138298 305311241
7 15101031004 15101 1102 2017 Arica 310013.3 2017 15101 221364 68625788545 4041 0.0182550 15101 1252763826 1344050454
8 15101031005 15101 235 2017 Arica 310013.3 2017 15101 221364 68625788545 1010 0.0045626 15101 313113453 396495725
9 15101031006 15101 866 2017 Arica 310013.3 2017 15101 221364 68625788545 2848 0.0128657 15101 882917935 1097498242
10 15101031007 15101 824 2017 Arica 310013.3 2017 15101 221364 68625788545 2762 0.0124772 15101 856256790 1053082012
11 15101041001 15101 638 2017 Arica 310013.3 2017 15101 221364 68625788545 2911 0.0131503 15101 902448774 853813338
12 15101041002 15101 702 2017 Arica 310013.3 2017 15101 221364 68625788545 2836 0.0128115 15101 879197775 922905094
13 15101041003 15101 665 2017 Arica 310013.3 2017 15101 221364 68625788545 3316 0.0149799 15101 1028004169 883036567
14 15101041004 15101 691 2017 Arica 310013.3 2017 15101 221364 68625788545 2869 0.0129606 15101 889428215 911072966
15 15101051001 15101 310 2017 Arica 310013.3 2017 15101 221364 68625788545 2041 0.0092201 15101 632737186 486064553
16 15101051002 15101 327 2017 Arica 310013.3 2017 15101 221364 68625788545 1655 0.0074764 15101 513072044 505953111
17 15101051003 15101 438 2017 Arica 310013.3 2017 15101 221364 68625788545 2345 0.0105934 15101 726981235 633062446
18 15101051004 15101 415 2017 Arica 310013.3 2017 15101 221364 68625788545 2196 0.0099203 15101 680789250 607063877
19 15101061001 15101 965 2017 Arica 310013.3 2017 15101 221364 68625788545 5901 0.0266575 15101 1829388601 1201503650
20 15101061002 15101 586 2017 Arica 310013.3 2017 15101 221364 68625788545 2066 0.0093330 15101 640487519 797193479
21 15101061003 15101 126 2017 Arica 310013.3 2017 15101 221364 68625788545 615 0.0027782 15101 190658192 258059138
22 15101061004 15101 137 2017 Arica 310013.3 2017 15101 221364 68625788545 723 0.0032661 15101 224139630 272710194
23 15101061005 15101 160 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333 302734808
24 15101071001 15101 507 2017 Arica 310013.3 2017 15101 221364 68625788545 2416 0.0109142 15101 748992181 710201108
25 15101071002 15101 632 2017 Arica 310013.3 2017 15101 221364 68625788545 2443 0.0110361 15101 757362541 847303586
26 15101071003 15101 767 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097 992489809
27 15101071004 15101 455 2017 Arica 310013.3 2017 15101 221364 68625788545 2455 0.0110903 15101 761082700 652181712
28 15101081001 15101 539 2017 Arica 310013.3 2017 15101 221364 68625788545 2735 0.0123552 15101 847886430 745592721
29 15101081002 15101 537 2017 Arica 310013.3 2017 15101 221364 68625788545 2496 0.0112755 15101 773793246 743387176
30 15101081003 15101 790 2017 Arica 310013.3 2017 15101 221364 68625788545 2823 0.0127528 15101 875167602 1016984734
31 15101081004 15101 445 2017 Arica 310013.3 2017 15101 221364 68625788545 1853 0.0083708 15101 574454682 640944738
32 15101091001 15101 565 2017 Arica 310013.3 2017 15101 221364 68625788545 2582 0.0116640 15101 800454392 774190883
33 15101091002 15101 360 2017 Arica 310013.3 2017 15101 221364 68625788545 1634 0.0073815 15101 506561765 544205976
34 15101101001 15101 430 2017 Arica 310013.3 2017 15101 221364 68625788545 1944 0.0087819 15101 602665894 624037090
35 15101101002 15101 363 2017 Arica 310013.3 2017 15101 221364 68625788545 1732 0.0078242 15101 536943070 547662067
36 15101101003 15101 339 2017 Arica 310013.3 2017 15101 221364 68625788545 1542 0.0069659 15101 478040539 519915027
37 15101111001 15101 462 2017 Arica 310013.3 2017 15101 221364 68625788545 1879 0.0084883 15101 582515028 660031616
38 15101111002 15101 454 2017 Arica 310013.3 2017 15101 221364 68625788545 1971 0.0089039 15101 611036254 651059234
39 15101111003 15101 467 2017 Arica 310013.3 2017 15101 221364 68625788545 1656 0.0074809 15101 513382058 665630827
40 15101121001 15101 737 2017 Arica 310013.3 2017 15101 221364 68625788545 3627 0.0163848 15101 1124418311 960442215
41 15101121002 15101 1263 2017 Arica 310013.3 2017 15101 221364 68625788545 4645 0.0209835 15101 1440011871 1509874553
42 15101121003 15101 1083 2017 Arica 310013.3 2017 15101 221364 68625788545 4990 0.0225421 15101 1546966466 1324366903
43 15101121004 15101 914 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097 1148040076
44 15101121005 15101 1122 2017 Arica 310013.3 2017 15101 221364 68625788545 4822 0.0217831 15101 1494884229 1364742302
45 15101121006 15101 64 2017 Arica 310013.3 2017 15101 221364 68625788545 263 0.0011881 15101 81533503 170066976
46 15101121007 15101 1177 2017 Arica 310013.3 2017 15101 221364 68625788545 3635 0.0164209 15101 1126898418 1421504484
47 15101121008 15101 893 2017 Arica 310013.3 2017 15101 221364 68625788545 3757 0.0169720 15101 1164720043 1125955717
48 15101121009 15101 623 2017 Arica 310013.3 2017 15101 221364 68625788545 2707 0.0122287 15101 839206057 837527897
49 15101121010 15101 943 2017 Arica 310013.3 2017 15101 221364 68625788545 3363 0.0151922 15101 1042574795 1178469658
50 15101141001 15101 318 2017 Arica 310013.3 2017 15101 221364 68625788545 1387 0.0062657 15101 429988475 495440341
51 15101141002 15101 412 2017 Arica 310013.3 2017 15101 221364 68625788545 1641 0.0074131 15101 508731858 603661056
52 15101171001 15101 980 2017 Arica 310013.3 2017 15101 221364 68625788545 4176 0.0188649 15101 1294615624 1217184511
53 15101171002 15101 575 2017 Arica 310013.3 2017 15101 221364 68625788545 2174 0.0098209 15101 673968957 785154853
54 15101171003 15101 719 2017 Arica 310013.3 2017 15101 221364 68625788545 2503 0.0113072 15101 775963340 941158012
55 15101171004 15101 609 2017 Arica 310013.3 2017 15101 221364 68625788545 2478 0.0111942 15101 768213007 822294265
56 15101171005 15101 559 2017 Arica 310013.3 2017 15101 221364 68625788545 2371 0.0107109 15101 735041581 767603251
57 15101171006 15101 1117 2017 Arica 310013.3 2017 15101 221364 68625788545 4924 0.0222439 15101 1526505587 1359571961
58 15101171007 15101 610 2017 Arica 310013.3 2017 15101 221364 68625788545 2627 0.0118673 15101 814404991 823383492
59 15101171008 15101 541 2017 Arica 310013.3 2017 15101 221364 68625788545 2763 0.0124817 15101 856566803 747797436
60 15101171009 15101 715 2017 Arica 310013.3 2017 15101 221364 68625788545 3821 0.0172612 15101 1184560895 936866762
61 15101171010 15101 961 2017 Arica 310013.3 2017 15101 221364 68625788545 3487 0.0157523 15101 1081016446 1197318814
62 15101181001 15101 794 2017 Arica 310013.3 2017 15101 221364 68625788545 3763 0.0169992 15101 1166580123 1021238326
63 15101181002 15101 1042 2017 Arica 310013.3 2017 15101 221364 68625788545 4627 0.0209022 15101 1434431631 1281801290
64 15101181003 15101 1207 2017 Arica 310013.3 2017 15101 221364 68625788545 5211 0.0235404 15101 1615479410 1452382747
65 15101181004 15101 643 2017 Arica 310013.3 2017 15101 221364 68625788545 2346 0.0105979 15101 727291248 859233696
66 15101181005 15101 917 2017 Arica 310013.3 2017 15101 221364 68625788545 3954 0.0178620 15101 1225792667 1151191567
67 15101181006 15101 1516 2017 Arica 310013.3 2017 15101 221364 68625788545 7629 0.0344636 15101 2365091617 1767508139
68 15101191001 15101 517 2017 Arica 310013.3 2017 15101 221364 68625788545 2482 0.0112123 15101 769453060 721285020
69 15101191002 15101 705 2017 Arica 310013.3 2017 15101 221364 68625788545 3364 0.0151967 15101 1042884808 926129088
70 15101191003 15101 746 2017 Arica 310013.3 2017 15101 221364 68625788545 3067 0.0138550 15101 950810852 970068445
71 15101191004 15101 585 2017 Arica 310013.3 2017 15101 221364 68625788545 2862 0.0129289 15101 887258122 796099983
72 15101991999 15101 100 2017 Arica 310013.3 2017 15101 221364 68625788545 1947 0.0087955 15101 603595934 222483194
73 15201011001 15201 150 2017 Putre 283661.5 2017 15201 2765 784324030 1716 0.6206148 15201 486763123 289773629
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
NA.16 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.17 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.18 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.19 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.20 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.21 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.22 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.23 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.24 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.25 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.26 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 15101011001 15101 353 2017 Arica 310013.3 2017 15101 221364 68625788545 1545 0.0069795 15101 478970579 536128381 347008.7
2 15101011002 15101 261 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333 427924301 427496.8
3 15101021001 15101 1234 2017 Arica 310013.3 2017 15101 221364 68625788545 4454 0.0201207 15101 1380799327 1480125454 332313.8
4 15101031001 15101 2220 2017 Arica 310013.3 2017 15101 221364 68625788545 6208 0.0280443 15101 1924562690 2471187560 398065.0
5 15101031002 15101 1007 2017 Arica 310013.3 2017 15101 221364 68625788545 3166 0.0143022 15101 981502171 1245362341 393355.1
6 15101031003 15101 162 2017 Arica 310013.3 2017 15101 221364 68625788545 623 0.0028144 15101 193138298 305311241 490066.2
7 15101031004 15101 1102 2017 Arica 310013.3 2017 15101 221364 68625788545 4041 0.0182550 15101 1252763826 1344050454 332603.4
8 15101031005 15101 235 2017 Arica 310013.3 2017 15101 221364 68625788545 1010 0.0045626 15101 313113453 396495725 392570.0
9 15101031006 15101 866 2017 Arica 310013.3 2017 15101 221364 68625788545 2848 0.0128657 15101 882917935 1097498242 385357.5
10 15101031007 15101 824 2017 Arica 310013.3 2017 15101 221364 68625788545 2762 0.0124772 15101 856256790 1053082012 381275.2
11 15101041001 15101 638 2017 Arica 310013.3 2017 15101 221364 68625788545 2911 0.0131503 15101 902448774 853813338 293305.9
12 15101041002 15101 702 2017 Arica 310013.3 2017 15101 221364 68625788545 2836 0.0128115 15101 879197775 922905094 325424.9
13 15101041003 15101 665 2017 Arica 310013.3 2017 15101 221364 68625788545 3316 0.0149799 15101 1028004169 883036567 266295.7
14 15101041004 15101 691 2017 Arica 310013.3 2017 15101 221364 68625788545 2869 0.0129606 15101 889428215 911072966 317557.7
15 15101051001 15101 310 2017 Arica 310013.3 2017 15101 221364 68625788545 2041 0.0092201 15101 632737186 486064553 238150.2
16 15101051002 15101 327 2017 Arica 310013.3 2017 15101 221364 68625788545 1655 0.0074764 15101 513072044 505953111 305711.8
17 15101051003 15101 438 2017 Arica 310013.3 2017 15101 221364 68625788545 2345 0.0105934 15101 726981235 633062446 269962.7
18 15101051004 15101 415 2017 Arica 310013.3 2017 15101 221364 68625788545 2196 0.0099203 15101 680789250 607063877 276440.7
19 15101061001 15101 965 2017 Arica 310013.3 2017 15101 221364 68625788545 5901 0.0266575 15101 1829388601 1201503650 203610.2
20 15101061002 15101 586 2017 Arica 310013.3 2017 15101 221364 68625788545 2066 0.0093330 15101 640487519 797193479 385863.3
21 15101061003 15101 126 2017 Arica 310013.3 2017 15101 221364 68625788545 615 0.0027782 15101 190658192 258059138 419608.4
22 15101061004 15101 137 2017 Arica 310013.3 2017 15101 221364 68625788545 723 0.0032661 15101 224139630 272710194 377192.5
23 15101061005 15101 160 2017 Arica 310013.3 2017 15101 221364 68625788545 1001 0.0045220 15101 310323333 302734808 302432.4
24 15101071001 15101 507 2017 Arica 310013.3 2017 15101 221364 68625788545 2416 0.0109142 15101 748992181 710201108 293957.4
25 15101071002 15101 632 2017 Arica 310013.3 2017 15101 221364 68625788545 2443 0.0110361 15101 757362541 847303586 346829.1
26 15101071003 15101 767 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097 992489809 269259.3
27 15101071004 15101 455 2017 Arica 310013.3 2017 15101 221364 68625788545 2455 0.0110903 15101 761082700 652181712 265654.5
28 15101081001 15101 539 2017 Arica 310013.3 2017 15101 221364 68625788545 2735 0.0123552 15101 847886430 745592721 272611.6
29 15101081002 15101 537 2017 Arica 310013.3 2017 15101 221364 68625788545 2496 0.0112755 15101 773793246 743387176 297831.4
30 15101081003 15101 790 2017 Arica 310013.3 2017 15101 221364 68625788545 2823 0.0127528 15101 875167602 1016984734 360249.6
31 15101081004 15101 445 2017 Arica 310013.3 2017 15101 221364 68625788545 1853 0.0083708 15101 574454682 640944738 345895.7
32 15101091001 15101 565 2017 Arica 310013.3 2017 15101 221364 68625788545 2582 0.0116640 15101 800454392 774190883 299841.6
33 15101091002 15101 360 2017 Arica 310013.3 2017 15101 221364 68625788545 1634 0.0073815 15101 506561765 544205976 333051.4
34 15101101001 15101 430 2017 Arica 310013.3 2017 15101 221364 68625788545 1944 0.0087819 15101 602665894 624037090 321006.7
35 15101101002 15101 363 2017 Arica 310013.3 2017 15101 221364 68625788545 1732 0.0078242 15101 536943070 547662067 316202.1
36 15101101003 15101 339 2017 Arica 310013.3 2017 15101 221364 68625788545 1542 0.0069659 15101 478040539 519915027 337169.3
37 15101111001 15101 462 2017 Arica 310013.3 2017 15101 221364 68625788545 1879 0.0084883 15101 582515028 660031616 351267.5
38 15101111002 15101 454 2017 Arica 310013.3 2017 15101 221364 68625788545 1971 0.0089039 15101 611036254 651059234 330319.2
39 15101111003 15101 467 2017 Arica 310013.3 2017 15101 221364 68625788545 1656 0.0074809 15101 513382058 665630827 401951.0
40 15101121001 15101 737 2017 Arica 310013.3 2017 15101 221364 68625788545 3627 0.0163848 15101 1124418311 960442215 264803.5
41 15101121002 15101 1263 2017 Arica 310013.3 2017 15101 221364 68625788545 4645 0.0209835 15101 1440011871 1509874553 325053.7
42 15101121003 15101 1083 2017 Arica 310013.3 2017 15101 221364 68625788545 4990 0.0225421 15101 1546966466 1324366903 265404.2
43 15101121004 15101 914 2017 Arica 310013.3 2017 15101 221364 68625788545 3686 0.0166513 15101 1142709097 1148040076 311459.6
44 15101121005 15101 1122 2017 Arica 310013.3 2017 15101 221364 68625788545 4822 0.0217831 15101 1494884229 1364742302 283024.1
45 15101121006 15101 64 2017 Arica 310013.3 2017 15101 221364 68625788545 263 0.0011881 15101 81533503 170066976 646642.5
46 15101121007 15101 1177 2017 Arica 310013.3 2017 15101 221364 68625788545 3635 0.0164209 15101 1126898418 1421504484 391060.4
47 15101121008 15101 893 2017 Arica 310013.3 2017 15101 221364 68625788545 3757 0.0169720 15101 1164720043 1125955717 299695.4
48 15101121009 15101 623 2017 Arica 310013.3 2017 15101 221364 68625788545 2707 0.0122287 15101 839206057 837527897 309393.4
49 15101121010 15101 943 2017 Arica 310013.3 2017 15101 221364 68625788545 3363 0.0151922 15101 1042574795 1178469658 350422.1
50 15101141001 15101 318 2017 Arica 310013.3 2017 15101 221364 68625788545 1387 0.0062657 15101 429988475 495440341 357202.8
51 15101141002 15101 412 2017 Arica 310013.3 2017 15101 221364 68625788545 1641 0.0074131 15101 508731858 603661056 367861.7
52 15101171001 15101 980 2017 Arica 310013.3 2017 15101 221364 68625788545 4176 0.0188649 15101 1294615624 1217184511 291471.4
53 15101171002 15101 575 2017 Arica 310013.3 2017 15101 221364 68625788545 2174 0.0098209 15101 673968957 785154853 361156.8
54 15101171003 15101 719 2017 Arica 310013.3 2017 15101 221364 68625788545 2503 0.0113072 15101 775963340 941158012 376012.0
55 15101171004 15101 609 2017 Arica 310013.3 2017 15101 221364 68625788545 2478 0.0111942 15101 768213007 822294265 331837.9
56 15101171005 15101 559 2017 Arica 310013.3 2017 15101 221364 68625788545 2371 0.0107109 15101 735041581 767603251 323746.6
57 15101171006 15101 1117 2017 Arica 310013.3 2017 15101 221364 68625788545 4924 0.0222439 15101 1526505587 1359571961 276111.3
58 15101171007 15101 610 2017 Arica 310013.3 2017 15101 221364 68625788545 2627 0.0118673 15101 814404991 823383492 313431.1
59 15101171008 15101 541 2017 Arica 310013.3 2017 15101 221364 68625788545 2763 0.0124817 15101 856566803 747797436 270646.9
60 15101171009 15101 715 2017 Arica 310013.3 2017 15101 221364 68625788545 3821 0.0172612 15101 1184560895 936866762 245188.9
61 15101171010 15101 961 2017 Arica 310013.3 2017 15101 221364 68625788545 3487 0.0157523 15101 1081016446 1197318814 343366.5
62 15101181001 15101 794 2017 Arica 310013.3 2017 15101 221364 68625788545 3763 0.0169992 15101 1166580123 1021238326 271389.4
63 15101181002 15101 1042 2017 Arica 310013.3 2017 15101 221364 68625788545 4627 0.0209022 15101 1434431631 1281801290 277026.4
64 15101181003 15101 1207 2017 Arica 310013.3 2017 15101 221364 68625788545 5211 0.0235404 15101 1615479410 1452382747 278714.8
65 15101181004 15101 643 2017 Arica 310013.3 2017 15101 221364 68625788545 2346 0.0105979 15101 727291248 859233696 366254.8
66 15101181005 15101 917 2017 Arica 310013.3 2017 15101 221364 68625788545 3954 0.0178620 15101 1225792667 1151191567 291146.1
67 15101181006 15101 1516 2017 Arica 310013.3 2017 15101 221364 68625788545 7629 0.0344636 15101 2365091617 1767508139 231682.8
68 15101191001 15101 517 2017 Arica 310013.3 2017 15101 221364 68625788545 2482 0.0112123 15101 769453060 721285020 290606.4
69 15101191002 15101 705 2017 Arica 310013.3 2017 15101 221364 68625788545 3364 0.0151967 15101 1042884808 926129088 275305.9
70 15101191003 15101 746 2017 Arica 310013.3 2017 15101 221364 68625788545 3067 0.0138550 15101 950810852 970068445 316292.3
71 15101191004 15101 585 2017 Arica 310013.3 2017 15101 221364 68625788545 2862 0.0129289 15101 887258122 796099983 278162.1
72 15101991999 15101 100 2017 Arica 310013.3 2017 15101 221364 68625788545 1947 0.0087955 15101 603595934 222483194 114269.7
73 15201011001 15201 150 2017 Putre 283661.5 2017 15201 2765 784324030 1716 0.6206148 15201 486763123 289773629 168865.8
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
NA.16 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.17 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.18 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.19 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.20 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.21 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.22 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.23 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.24 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.25 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
NA.26 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_15_P03C_u.rds")