Expansión de la CASEN sobre el CENSO (Región 15)

De ingresos sobre una categoría de respuesta

DataIntelligence

Abstract

Expandiremos los ingresos promedios comunales obtenidos de la CASEN sobre la categoría de respuesta: “tejas o tejuelas de arcilla, metálicas, de cemento, de madera, asfálticas o plásticas” del campo P03B del CENSO de viviendas, que fue la categoría de respuesta que más alto correlaciona con los ingresos expandidos (obtenidos de la multiplicación del ingreso promedio y los habitantes), ambos a nivel comunal.

Haremos el análisis sobre la región 15.

Ensayaremos diferentes modelos dentro del análisis de regresión cuya variable independiente será: “frecuencia de población que posee la variable Censal respecto a la zona” y la dependiente: “ingreso expandido por zona”

Lo anterior para elegir el que posea el mayor coeficiente de determinación y así contruir una tabla de valores predichos.

1 Variable CENSO

Necesitamos calcular las frecuencias a nivel censal de las respuestas correspondientes a la categoría: “tejas o tejuelas de arcilla, metálicas, de cemento, de madera, asfálticas o plásticas” del campo P03B 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 Casen 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)
regiones

##  [1] 15 14 13 12 11 10  9 16  8  7  6  5  4  3  2  1

Hagamos un subset con la 1:

tabla_con_clave <- filter(tabla_con_clave, tabla_con_clave$REGION == 15) 

1.2 Cálculo de frecuencias

tabla_con_clave_f <- tabla_con_clave[,-c(1,2,4,5,6,7,8,9,10,11,13,14,15,16,17,18,19,20),drop=FALSE] 
names(tabla_con_clave_f)[2] <- "Tipo de techo" 
# Ahora filtramos por Tipo de techo = 1.
tabla_con_clave_ff <- filter(tabla_con_clave_f, tabla_con_clave_f$`Tipo de techo` == 1)
# Determinamos las frecuencias por zona:
b <- tabla_con_clave_ff$clave
c <- tabla_con_clave_ff$`Tipo de techo`
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	150	2017
2	15101011002	1	15101	113	2017
3	15101021001	1	15101	427	2017
4	15101031001	1	15101	308	2017
5	15101031002	1	15101	188	2017
6	15101031003	1	15101	49	2017
7	15101031004	1	15101	490	2017
8	15101031005	1	15101	89	2017
9	15101031006	1	15101	206	2017
10	15101031007	1	15101	146	2017
11	15101032003	1	15101	10	2017
12	15101041001	1	15101	210	2017
13	15101041002	1	15101	178	2017
14	15101041003	1	15101	345	2017
15	15101041004	1	15101	312	2017
16	15101051001	1	15101	156	2017
17	15101051002	1	15101	137	2017
18	15101051003	1	15101	280	2017
19	15101051004	1	15101	189	2017
20	15101061001	1	15101	403	2017
21	15101061002	1	15101	234	2017
22	15101061003	1	15101	87	2017
23	15101061004	1	15101	139	2017
24	15101061005	1	15101	46	2017
25	15101062011	1	15101	38	2017
26	15101062014	1	15101	705	2017
27	15101071001	1	15101	158	2017
28	15101071002	1	15101	115	2017
29	15101071003	1	15101	234	2017
30	15101071004	1	15101	138	2017
31	15101081001	1	15101	284	2017
32	15101081002	1	15101	430	2017
33	15101081003	1	15101	297	2017
34	15101081004	1	15101	198	2017
35	15101091001	1	15101	375	2017
36	15101091002	1	15101	187	2017
37	15101101001	1	15101	235	2017
38	15101101002	1	15101	231	2017
39	15101101003	1	15101	227	2017
40	15101111001	1	15101	198	2017
41	15101111002	1	15101	251	2017
42	15101111003	1	15101	162	2017
43	15101121001	1	15101	288	2017
44	15101121002	1	15101	364	2017
45	15101121003	1	15101	395	2017
46	15101121004	1	15101	214	2017
47	15101121005	1	15101	339	2017
48	15101121006	1	15101	33	2017
49	15101121007	1	15101	292	2017
50	15101121008	1	15101	346	2017
51	15101121009	1	15101	251	2017
52	15101121010	1	15101	303	2017
53	15101132015	1	15101	17	2017
54	15101141001	1	15101	220	2017
55	15101141002	1	15101	230	2017
56	15101152005	1	15101	79	2017
57	15101152007	1	15101	5	2017
58	15101152016	1	15101	128	2017
59	15101152901	1	15101	24	2017
60	15101162016	1	15101	31	2017
61	15101162901	1	15101	5	2017
62	15101171001	1	15101	385	2017
63	15101171002	1	15101	119	2017
64	15101171003	1	15101	135	2017
65	15101171004	1	15101	210	2017
66	15101171005	1	15101	144	2017
67	15101171006	1	15101	457	2017
68	15101171007	1	15101	196	2017
69	15101171008	1	15101	228	2017
70	15101171009	1	15101	362	2017
71	15101171010	1	15101	156	2017
72	15101181001	1	15101	328	2017
73	15101181002	1	15101	384	2017
74	15101181003	1	15101	297	2017
75	15101181004	1	15101	154	2017
76	15101181005	1	15101	380	2017
77	15101181006	1	15101	380	2017
78	15101191001	1	15101	227	2017
79	15101191002	1	15101	386	2017
80	15101191003	1	15101	350	2017
81	15101191004	1	15101	235	2017
82	15101991999	1	15101	36	2017
191	15102012002	1	15102	2	2017
192	15102012003	1	15102	41	2017
193	15102022007	1	15102	2	2017
194	15102032005	1	15102	2	2017
195	15102032901	1	15102	8	2017
196	15102042001	1	15102	21	2017
197	15102052004	1	15102	24	2017
306	15201011001	1	15201	34	2017
307	15201022002	1	15201	2	2017
308	15201032006	1	15201	6	2017
309	15201042021	1	15201	10	2017
310	15201052001	1	15201	1	2017
311	15201052012	1	15201	2	2017
312	15201052013	1	15201	2	2017
313	15201052020	1	15201	1	2017
314	15201052901	1	15201	2	2017
315	15201072901	1	15201	1	2017
424	15202012008	1	15202	1	2017

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	150	2017	15101
2	15101011002	113	2017	15101
3	15101021001	427	2017	15101
4	15101031001	308	2017	15101
5	15101031002	188	2017	15101
6	15101031003	49	2017	15101
7	15101031004	490	2017	15101
8	15101031005	89	2017	15101
9	15101031006	206	2017	15101
10	15101031007	146	2017	15101
11	15101032003	10	2017	15101
12	15101041001	210	2017	15101
13	15101041002	178	2017	15101
14	15101041003	345	2017	15101
15	15101041004	312	2017	15101
16	15101051001	156	2017	15101
17	15101051002	137	2017	15101
18	15101051003	280	2017	15101
19	15101051004	189	2017	15101
20	15101061001	403	2017	15101
21	15101061002	234	2017	15101
22	15101061003	87	2017	15101
23	15101061004	139	2017	15101
24	15101061005	46	2017	15101
25	15101062011	38	2017	15101
26	15101062014	705	2017	15101
27	15101071001	158	2017	15101
28	15101071002	115	2017	15101
29	15101071003	234	2017	15101
30	15101071004	138	2017	15101
31	15101081001	284	2017	15101
32	15101081002	430	2017	15101
33	15101081003	297	2017	15101
34	15101081004	198	2017	15101
35	15101091001	375	2017	15101
36	15101091002	187	2017	15101
37	15101101001	235	2017	15101
38	15101101002	231	2017	15101
39	15101101003	227	2017	15101
40	15101111001	198	2017	15101
41	15101111002	251	2017	15101
42	15101111003	162	2017	15101
43	15101121001	288	2017	15101
44	15101121002	364	2017	15101
45	15101121003	395	2017	15101
46	15101121004	214	2017	15101
47	15101121005	339	2017	15101
48	15101121006	33	2017	15101
49	15101121007	292	2017	15101
50	15101121008	346	2017	15101
51	15101121009	251	2017	15101
52	15101121010	303	2017	15101
53	15101132015	17	2017	15101
54	15101141001	220	2017	15101
55	15101141002	230	2017	15101
56	15101152005	79	2017	15101
57	15101152007	5	2017	15101
58	15101152016	128	2017	15101
59	15101152901	24	2017	15101
60	15101162016	31	2017	15101
61	15101162901	5	2017	15101
62	15101171001	385	2017	15101
63	15101171002	119	2017	15101
64	15101171003	135	2017	15101
65	15101171004	210	2017	15101
66	15101171005	144	2017	15101
67	15101171006	457	2017	15101
68	15101171007	196	2017	15101
69	15101171008	228	2017	15101
70	15101171009	362	2017	15101
71	15101171010	156	2017	15101
72	15101181001	328	2017	15101
73	15101181002	384	2017	15101
74	15101181003	297	2017	15101
75	15101181004	154	2017	15101
76	15101181005	380	2017	15101
77	15101181006	380	2017	15101
78	15101191001	227	2017	15101
79	15101191002	386	2017	15101
80	15101191003	350	2017	15101
81	15101191004	235	2017	15101
82	15101991999	36	2017	15101
191	15102012002	2	2017	15102
192	15102012003	41	2017	15102
193	15102022007	2	2017	15102
194	15102032005	2	2017	15102
195	15102032901	8	2017	15102
196	15102042001	21	2017	15102
197	15102052004	24	2017	15102
306	15201011001	34	2017	15201
307	15201022002	2	2017	15201
308	15201032006	6	2017	15201
309	15201042021	10	2017	15201
310	15201052001	1	2017	15201
311	15201052012	2	2017	15201
312	15201052013	2	2017	15201
313	15201052020	1	2017	15201
314	15201052901	2	2017	15201
315	15201072901	1	2017	15201
424	15202012008	1	2017	15202

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_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	354820.7	2017	1101	191468	67936815240
01107	Alto Hospicio	301933.4	2017	1107	108375	32722034397
01401	Pozo Almonte	285981.8	2017	1401	15711	4493059532
01402	Camiña	262850.3	2017	1402	1250	328562901
01404	Huara	253968.5	2017	1404	2730	693334131
01405	Pica	313007.5	2017	1405	9296	2909717399
02101	Antofagasta	347580.2	2017	2101	361873	125779893517
02102	Mejillones	369770.7	2017	2102	13467	4979702302
02103	Sierra Gorda	403458.5	2017	2103	10186	4109628188
02104	Taltal	364539.1	2017	2104	13317	4854566842
02201	Calama	409671.3	2017	2201	165731	67895226712
02203	San Pedro de Atacama	426592.0	2017	2203	10996	4690805471
02301	Tocopilla	246615.3	2017	2301	25186	6211253937
02302	María Elena	466266.9	2017	2302	6457	3010685220
03101	Copiapó	330075.2	2017	3101	153937	50810778473
03102	Caldera	299314.8	2017	3102	17662	5286498241
03103	Tierra Amarilla	314643.9	2017	3103	14019	4410992711
03201	Chañaral	286389.3	2017	3201	12219	3499391196
03202	Diego de Almagro	336256.8	2017	3202	13925	4682376047
03301	Vallenar	304336.7	2017	3301	51917	15800246795
03302	Alto del Carmen	227130.4	2017	3302	5299	1203563833
03303	Freirina	253086.7	2017	3303	7041	1781983257
03304	Huasco	287406.6	2017	3304	10149	2916889629
04101	La Serena	270221.9	2017	4101	221054	59733627577
04102	Coquimbo	261852.6	2017	4102	227730	59631700074
04103	Andacollo	248209.3	2017	4103	11044	2741223967
04104	La Higuera	228356.8	2017	4104	4241	968461330
04105	Paiguano	205942.1	2017	4105	4497	926121774
04106	Vicuña	211431.9	2017	4106	27771	5871675449
04201	Illapel	238674.4	2017	4201	30848	7362627007
04202	Canela	207933.6	2017	4202	9093	1890740321
04203	Los Vilos	255200.4	2017	4203	21382	5456695139
04204	Salamanca	242879.5	2017	4204	29347	7127783272
04301	Ovalle	266522.9	2017	4301	111272	29656533187
04302	Combarbalá	210409.7	2017	4302	13322	2803077721
04303	Monte Patria	211907.9	2017	4303	30751	6516380780
04304	Punitaqui	194997.8	2017	4304	10956	2136395349
04305	Río Hurtado	182027.2	2017	4305	4278	778712384
05101	Valparaíso	298720.7	2017	5101	296655	88616992249
05102	Casablanca	312802.7	2017	5102	26867	8404070481
05103	Concón	318496.3	2017	5103	42152	13425257057
05105	Puchuncaví	288737.2	2017	5105	18546	5354920887
05107	Quintero	316659.1	2017	5107	31923	10108709691
05109	Viña del Mar	337006.1	2017	5109	334248	112643604611
05301	Los Andes	338182.5	2017	5301	66708	22559476922
05302	Calle Larga	245165.4	2017	5302	14832	3636293159
05303	Rinconada	281633.2	2017	5303	10207	2874630315
05304	San Esteban	220958.4	2017	5304	18855	4166170587
05401	La Ligua	229623.7	2017	5401	35390	8126381563
05402	Cabildo	249717.7	2017	5402	19388	4841527150

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)
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
15101	15101031003	49	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031004	490	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031005	89	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031002	188	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101011001	150	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101011002	113	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101021001	427	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031001	308	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051003	280	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051004	189	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061001	403	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061002	234	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031006	206	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031007	146	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101032003	10	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041001	210	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041002	178	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041003	345	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041004	312	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051001	156	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051002	137	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081001	284	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081002	430	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081003	297	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081004	198	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061003	87	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061004	139	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061005	46	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101062011	38	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101062014	705	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071001	158	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071002	115	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071003	234	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071004	138	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121002	364	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121003	395	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121004	214	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121005	339	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101091001	375	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101091002	187	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101001	235	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101002	231	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101003	227	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111001	198	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111002	251	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111003	162	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121001	288	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152007	5	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152016	128	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152901	24	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101162016	31	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121006	33	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121007	292	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121008	346	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121009	251	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121010	303	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101132015	17	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101141001	220	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101141002	230	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152005	79	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171009	362	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171010	156	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181001	328	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181002	384	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101162901	5	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171001	385	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171002	119	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171003	135	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171004	210	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171005	144	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171006	457	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171007	196	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171008	228	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181003	297	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181004	154	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181005	380	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181006	380	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191001	227	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191002	386	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191003	350	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191004	235	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101991999	36	2017	Arica	292440.1	2017	15101	221364	64735720258
15102	15102012003	41	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102022007	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102032005	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102012002	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102032901	8	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102042001	21	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102052004	24	2017	Camarones	233083.8	2017	15102	1255	292520169
15201	15201072901	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052013	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052020	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052901	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201011001	34	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201022002	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201032006	6	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201042021	10	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052001	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052012	2	2017	Putre	251446.5	2017	15201	2765	695249646
15202	15202022901	3	2017	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
15101	15101031003	49	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031004	490	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031005	89	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031002	188	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101011001	150	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101011002	113	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101021001	427	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031001	308	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051003	280	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051004	189	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061001	403	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061002	234	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031006	206	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101031007	146	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101032003	10	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041001	210	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041002	178	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041003	345	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101041004	312	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051001	156	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101051002	137	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081001	284	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081002	430	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081003	297	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101081004	198	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061003	87	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061004	139	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101061005	46	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101062011	38	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101062014	705	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071001	158	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071002	115	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071003	234	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101071004	138	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121002	364	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121003	395	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121004	214	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121005	339	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101091001	375	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101091002	187	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101001	235	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101002	231	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101101003	227	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111001	198	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111002	251	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101111003	162	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121001	288	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152007	5	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152016	128	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152901	24	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101162016	31	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121006	33	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121007	292	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121008	346	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121009	251	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101121010	303	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101132015	17	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101141001	220	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101141002	230	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101152005	79	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171009	362	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171010	156	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181001	328	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181002	384	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101162901	5	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171001	385	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171002	119	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171003	135	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171004	210	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171005	144	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171006	457	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171007	196	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101171008	228	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181003	297	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181004	154	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181005	380	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101181006	380	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191001	227	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191002	386	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191003	350	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101191004	235	2017	Arica	292440.1	2017	15101	221364	64735720258
15101	15101991999	36	2017	Arica	292440.1	2017	15101	221364	64735720258
15102	15102012003	41	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102022007	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102032005	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102012002	2	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102032901	8	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102042001	21	2017	Camarones	233083.8	2017	15102	1255	292520169
15102	15102052004	24	2017	Camarones	233083.8	2017	15102	1255	292520169
15201	15201072901	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052013	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052020	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052901	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201011001	34	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201022002	2	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201032006	6	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201042021	10	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052001	1	2017	Putre	251446.5	2017	15201	2765	695249646
15201	15201052012	2	2017	Putre	251446.5	2017	15201	2765	695249646
15202	15202022901	3	2017	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	150	2017	Arica	292440.1	2017	15101	221364	64735720258	1545	0.0069795	15101
15101011002	15101	113	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101
15101021001	15101	427	2017	Arica	292440.1	2017	15101	221364	64735720258	4454	0.0201207	15101
15101031001	15101	308	2017	Arica	292440.1	2017	15101	221364	64735720258	6208	0.0280443	15101
15101031002	15101	188	2017	Arica	292440.1	2017	15101	221364	64735720258	3166	0.0143022	15101
15101031003	15101	49	2017	Arica	292440.1	2017	15101	221364	64735720258	623	0.0028144	15101
15101031004	15101	490	2017	Arica	292440.1	2017	15101	221364	64735720258	4041	0.0182550	15101
15101031005	15101	89	2017	Arica	292440.1	2017	15101	221364	64735720258	1010	0.0045626	15101
15101031006	15101	206	2017	Arica	292440.1	2017	15101	221364	64735720258	2848	0.0128657	15101
15101031007	15101	146	2017	Arica	292440.1	2017	15101	221364	64735720258	2762	0.0124772	15101
15101032003	15101	10	2017	Arica	292440.1	2017	15101	221364	64735720258	70	0.0003162	15101
15101041001	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2911	0.0131503	15101
15101041002	15101	178	2017	Arica	292440.1	2017	15101	221364	64735720258	2836	0.0128115	15101
15101041003	15101	345	2017	Arica	292440.1	2017	15101	221364	64735720258	3316	0.0149799	15101
15101041004	15101	312	2017	Arica	292440.1	2017	15101	221364	64735720258	2869	0.0129606	15101
15101051001	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	2041	0.0092201	15101
15101051002	15101	137	2017	Arica	292440.1	2017	15101	221364	64735720258	1655	0.0074764	15101
15101051003	15101	280	2017	Arica	292440.1	2017	15101	221364	64735720258	2345	0.0105934	15101
15101051004	15101	189	2017	Arica	292440.1	2017	15101	221364	64735720258	2196	0.0099203	15101
15101061001	15101	403	2017	Arica	292440.1	2017	15101	221364	64735720258	5901	0.0266575	15101
15101061002	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	2066	0.0093330	15101
15101061003	15101	87	2017	Arica	292440.1	2017	15101	221364	64735720258	615	0.0027782	15101
15101061004	15101	139	2017	Arica	292440.1	2017	15101	221364	64735720258	723	0.0032661	15101
15101061005	15101	46	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101
15101062011	15101	38	2017	Arica	292440.1	2017	15101	221364	64735720258	2587	0.0116866	15101
15101062014	15101	705	2017	Arica	292440.1	2017	15101	221364	64735720258	8708	0.0393379	15101
15101071001	15101	158	2017	Arica	292440.1	2017	15101	221364	64735720258	2416	0.0109142	15101
15101071002	15101	115	2017	Arica	292440.1	2017	15101	221364	64735720258	2443	0.0110361	15101
15101071003	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101
15101071004	15101	138	2017	Arica	292440.1	2017	15101	221364	64735720258	2455	0.0110903	15101
15101081001	15101	284	2017	Arica	292440.1	2017	15101	221364	64735720258	2735	0.0123552	15101
15101081002	15101	430	2017	Arica	292440.1	2017	15101	221364	64735720258	2496	0.0112755	15101
15101081003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	2823	0.0127528	15101
15101081004	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1853	0.0083708	15101
15101091001	15101	375	2017	Arica	292440.1	2017	15101	221364	64735720258	2582	0.0116640	15101
15101091002	15101	187	2017	Arica	292440.1	2017	15101	221364	64735720258	1634	0.0073815	15101
15101101001	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	1944	0.0087819	15101
15101101002	15101	231	2017	Arica	292440.1	2017	15101	221364	64735720258	1732	0.0078242	15101
15101101003	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	1542	0.0069659	15101
15101111001	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1879	0.0084883	15101
15101111002	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	1971	0.0089039	15101
15101111003	15101	162	2017	Arica	292440.1	2017	15101	221364	64735720258	1656	0.0074809	15101
15101121001	15101	288	2017	Arica	292440.1	2017	15101	221364	64735720258	3627	0.0163848	15101
15101121002	15101	364	2017	Arica	292440.1	2017	15101	221364	64735720258	4645	0.0209835	15101
15101121003	15101	395	2017	Arica	292440.1	2017	15101	221364	64735720258	4990	0.0225421	15101
15101121004	15101	214	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101
15101121005	15101	339	2017	Arica	292440.1	2017	15101	221364	64735720258	4822	0.0217831	15101
15101121006	15101	33	2017	Arica	292440.1	2017	15101	221364	64735720258	263	0.0011881	15101
15101121007	15101	292	2017	Arica	292440.1	2017	15101	221364	64735720258	3635	0.0164209	15101
15101121008	15101	346	2017	Arica	292440.1	2017	15101	221364	64735720258	3757	0.0169720	15101
15101121009	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	2707	0.0122287	15101
15101121010	15101	303	2017	Arica	292440.1	2017	15101	221364	64735720258	3363	0.0151922	15101
15101132015	15101	17	2017	Arica	292440.1	2017	15101	221364	64735720258	338	0.0015269	15101
15101141001	15101	220	2017	Arica	292440.1	2017	15101	221364	64735720258	1387	0.0062657	15101
15101141002	15101	230	2017	Arica	292440.1	2017	15101	221364	64735720258	1641	0.0074131	15101
15101152005	15101	79	2017	Arica	292440.1	2017	15101	221364	64735720258	870	0.0039302	15101
15101152007	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	130	0.0005873	15101
15101152016	15101	128	2017	Arica	292440.1	2017	15101	221364	64735720258	1759	0.0079462	15101
15101152901	15101	24	2017	Arica	292440.1	2017	15101	221364	64735720258	207	0.0009351	15101
15101162016	15101	31	2017	Arica	292440.1	2017	15101	221364	64735720258	393	0.0017754	15101
15101162901	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	58	0.0002620	15101
15101171001	15101	385	2017	Arica	292440.1	2017	15101	221364	64735720258	4176	0.0188649	15101
15101171002	15101	119	2017	Arica	292440.1	2017	15101	221364	64735720258	2174	0.0098209	15101
15101171003	15101	135	2017	Arica	292440.1	2017	15101	221364	64735720258	2503	0.0113072	15101
15101171004	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2478	0.0111942	15101
15101171005	15101	144	2017	Arica	292440.1	2017	15101	221364	64735720258	2371	0.0107109	15101
15101171006	15101	457	2017	Arica	292440.1	2017	15101	221364	64735720258	4924	0.0222439	15101
15101171007	15101	196	2017	Arica	292440.1	2017	15101	221364	64735720258	2627	0.0118673	15101
15101171008	15101	228	2017	Arica	292440.1	2017	15101	221364	64735720258	2763	0.0124817	15101
15101171009	15101	362	2017	Arica	292440.1	2017	15101	221364	64735720258	3821	0.0172612	15101
15101171010	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	3487	0.0157523	15101
15101181001	15101	328	2017	Arica	292440.1	2017	15101	221364	64735720258	3763	0.0169992	15101
15101181002	15101	384	2017	Arica	292440.1	2017	15101	221364	64735720258	4627	0.0209022	15101
15101181003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	5211	0.0235404	15101
15101181004	15101	154	2017	Arica	292440.1	2017	15101	221364	64735720258	2346	0.0105979	15101
15101181005	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	3954	0.0178620	15101
15101181006	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	7629	0.0344636	15101
15101191001	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	2482	0.0112123	15101
15101191002	15101	386	2017	Arica	292440.1	2017	15101	221364	64735720258	3364	0.0151967	15101
15101191003	15101	350	2017	Arica	292440.1	2017	15101	221364	64735720258	3067	0.0138550	15101
15101191004	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	2862	0.0129289	15101
15101991999	15101	36	2017	Arica	292440.1	2017	15101	221364	64735720258	1947	0.0087955	15101
15102012002	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	25	0.0199203	15102
15102012003	15102	41	2017	Camarones	233083.8	2017	15102	1255	292520169	565	0.4501992	15102
15102022007	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	39	0.0310757	15102
15102032005	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	112	0.0892430	15102
15102032901	15102	8	2017	Camarones	233083.8	2017	15102	1255	292520169	30	0.0239044	15102
15102042001	15102	21	2017	Camarones	233083.8	2017	15102	1255	292520169	181	0.1442231	15102
15102052004	15102	24	2017	Camarones	233083.8	2017	15102	1255	292520169	259	0.2063745	15102
15201011001	15201	34	2017	Putre	251446.5	2017	15201	2765	695249646	1716	0.6206148	15201
15201022002	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	52	0.0188065	15201
15201032006	15201	6	2017	Putre	251446.5	2017	15201	2765	695249646	31	0.0112116	15201
15201042021	15201	10	2017	Putre	251446.5	2017	15201	2765	695249646	75	0.0271248	15201
15201052001	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	63	0.0227848	15201
15201052012	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	25	0.0090416	15201
15201052013	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	17	0.0061483	15201
15201052020	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	18	0.0065099	15201
15201052901	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	7	0.0025316	15201
15201072901	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	53	0.0191682	15201
15202012008	15202	1	2017	NA	NA	NA	NA	NA	NA	12	0.0175439	15202

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	150	2017	Arica	292440.1	2017	15101	221364	64735720258	1545	0.0069795	15101	451820024
15101011002	15101	113	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585
15101021001	15101	427	2017	Arica	292440.1	2017	15101	221364	64735720258	4454	0.0201207	15101	1302528406
15101031001	15101	308	2017	Arica	292440.1	2017	15101	221364	64735720258	6208	0.0280443	15101	1815468420
15101031002	15101	188	2017	Arica	292440.1	2017	15101	221364	64735720258	3166	0.0143022	15101	925865499
15101031003	15101	49	2017	Arica	292440.1	2017	15101	221364	64735720258	623	0.0028144	15101	182190210
15101031004	15101	490	2017	Arica	292440.1	2017	15101	221364	64735720258	4041	0.0182550	15101	1181750626
15101031005	15101	89	2017	Arica	292440.1	2017	15101	221364	64735720258	1010	0.0045626	15101	295364546
15101031006	15101	206	2017	Arica	292440.1	2017	15101	221364	64735720258	2848	0.0128657	15101	832869533
15101031007	15101	146	2017	Arica	292440.1	2017	15101	221364	64735720258	2762	0.0124772	15101	807719680
15101032003	15101	10	2017	Arica	292440.1	2017	15101	221364	64735720258	70	0.0003162	15101	20470810
15101041001	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2911	0.0131503	15101	851293262
15101041002	15101	178	2017	Arica	292440.1	2017	15101	221364	64735720258	2836	0.0128115	15101	829360251
15101041003	15101	345	2017	Arica	292440.1	2017	15101	221364	64735720258	3316	0.0149799	15101	969731521
15101041004	15101	312	2017	Arica	292440.1	2017	15101	221364	64735720258	2869	0.0129606	15101	839010776
15101051001	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	2041	0.0092201	15101	596870336
15101051002	15101	137	2017	Arica	292440.1	2017	15101	221364	64735720258	1655	0.0074764	15101	483988440
15101051003	15101	280	2017	Arica	292440.1	2017	15101	221364	64735720258	2345	0.0105934	15101	685772140
15101051004	15101	189	2017	Arica	292440.1	2017	15101	221364	64735720258	2196	0.0099203	15101	642198558
15101061001	15101	403	2017	Arica	292440.1	2017	15101	221364	64735720258	5901	0.0266575	15101	1725689296
15101061002	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	2066	0.0093330	15101	604181340
15101061003	15101	87	2017	Arica	292440.1	2017	15101	221364	64735720258	615	0.0027782	15101	179850689
15101061004	15101	139	2017	Arica	292440.1	2017	15101	221364	64735720258	723	0.0032661	15101	211434225
15101061005	15101	46	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585
15101062011	15101	38	2017	Arica	292440.1	2017	15101	221364	64735720258	2587	0.0116866	15101	756542655
15101062014	15101	705	2017	Arica	292440.1	2017	15101	221364	64735720258	8708	0.0393379	15101	2546568783
15101071001	15101	158	2017	Arica	292440.1	2017	15101	221364	64735720258	2416	0.0109142	15101	706535390
15101071002	15101	115	2017	Arica	292440.1	2017	15101	221364	64735720258	2443	0.0110361	15101	714431274
15101071003	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374
15101071004	15101	138	2017	Arica	292440.1	2017	15101	221364	64735720258	2455	0.0110903	15101	717940556
15101081001	15101	284	2017	Arica	292440.1	2017	15101	221364	64735720258	2735	0.0123552	15101	799823797
15101081002	15101	430	2017	Arica	292440.1	2017	15101	221364	64735720258	2496	0.0112755	15101	729930602
15101081003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	2823	0.0127528	15101	825558529
15101081004	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1853	0.0083708	15101	541891589
15101091001	15101	375	2017	Arica	292440.1	2017	15101	221364	64735720258	2582	0.0116640	15101	755080454
15101091002	15101	187	2017	Arica	292440.1	2017	15101	221364	64735720258	1634	0.0073815	15101	477847197
15101101001	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	1944	0.0087819	15101	568503642
15101101002	15101	231	2017	Arica	292440.1	2017	15101	221364	64735720258	1732	0.0078242	15101	506506331
15101101003	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	1542	0.0069659	15101	450942704
15101111001	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1879	0.0084883	15101	549495032
15101111002	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	1971	0.0089039	15101	576399526
15101111003	15101	162	2017	Arica	292440.1	2017	15101	221364	64735720258	1656	0.0074809	15101	484280880
15101121001	15101	288	2017	Arica	292440.1	2017	15101	221364	64735720258	3627	0.0163848	15101	1060680406
15101121002	15101	364	2017	Arica	292440.1	2017	15101	221364	64735720258	4645	0.0209835	15101	1358384474
15101121003	15101	395	2017	Arica	292440.1	2017	15101	221364	64735720258	4990	0.0225421	15101	1459276324
15101121004	15101	214	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374
15101121005	15101	339	2017	Arica	292440.1	2017	15101	221364	64735720258	4822	0.0217831	15101	1410146379
15101121006	15101	33	2017	Arica	292440.1	2017	15101	221364	64735720258	263	0.0011881	15101	76911758
15101121007	15101	292	2017	Arica	292440.1	2017	15101	221364	64735720258	3635	0.0164209	15101	1063019927
15101121008	15101	346	2017	Arica	292440.1	2017	15101	221364	64735720258	3757	0.0169720	15101	1098697625
15101121009	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	2707	0.0122287	15101	791635473
15101121010	15101	303	2017	Arica	292440.1	2017	15101	221364	64735720258	3363	0.0151922	15101	983476208
15101132015	15101	17	2017	Arica	292440.1	2017	15101	221364	64735720258	338	0.0015269	15101	98844769
15101141001	15101	220	2017	Arica	292440.1	2017	15101	221364	64735720258	1387	0.0062657	15101	405614481
15101141002	15101	230	2017	Arica	292440.1	2017	15101	221364	64735720258	1641	0.0074131	15101	479894278
15101152005	15101	79	2017	Arica	292440.1	2017	15101	221364	64735720258	870	0.0039302	15101	254422926
15101152007	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	130	0.0005873	15101	38017219
15101152016	15101	128	2017	Arica	292440.1	2017	15101	221364	64735720258	1759	0.0079462	15101	514402215
15101152901	15101	24	2017	Arica	292440.1	2017	15101	221364	64735720258	207	0.0009351	15101	60535110
15101162016	15101	31	2017	Arica	292440.1	2017	15101	221364	64735720258	393	0.0017754	15101	114928977
15101162901	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	58	0.0002620	15101	16961528
15101171001	15101	385	2017	Arica	292440.1	2017	15101	221364	64735720258	4176	0.0188649	15101	1221230046
15101171002	15101	119	2017	Arica	292440.1	2017	15101	221364	64735720258	2174	0.0098209	15101	635764875
15101171003	15101	135	2017	Arica	292440.1	2017	15101	221364	64735720258	2503	0.0113072	15101	731977683
15101171004	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2478	0.0111942	15101	724666679
15101171005	15101	144	2017	Arica	292440.1	2017	15101	221364	64735720258	2371	0.0107109	15101	693375584
15101171006	15101	457	2017	Arica	292440.1	2017	15101	221364	64735720258	4924	0.0222439	15101	1439975274
15101171007	15101	196	2017	Arica	292440.1	2017	15101	221364	64735720258	2627	0.0118673	15101	768240261
15101171008	15101	228	2017	Arica	292440.1	2017	15101	221364	64735720258	2763	0.0124817	15101	808012121
15101171009	15101	362	2017	Arica	292440.1	2017	15101	221364	64735720258	3821	0.0172612	15101	1117413794
15101171010	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	3487	0.0157523	15101	1019738786
15101181001	15101	328	2017	Arica	292440.1	2017	15101	221364	64735720258	3763	0.0169992	15101	1100452266
15101181002	15101	384	2017	Arica	292440.1	2017	15101	221364	64735720258	4627	0.0209022	15101	1353120551
15101181003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	5211	0.0235404	15101	1523905596
15101181004	15101	154	2017	Arica	292440.1	2017	15101	221364	64735720258	2346	0.0105979	15101	686064580
15101181005	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	3954	0.0178620	15101	1156308333
15101181006	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	7629	0.0344636	15101	2231025866
15101191001	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	2482	0.0112123	15101	725836440
15101191002	15101	386	2017	Arica	292440.1	2017	15101	221364	64735720258	3364	0.0151967	15101	983768648
15101191003	15101	350	2017	Arica	292440.1	2017	15101	221364	64735720258	3067	0.0138550	15101	896913925
15101191004	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	2862	0.0129289	15101	836963695
15101991999	15101	36	2017	Arica	292440.1	2017	15101	221364	64735720258	1947	0.0087955	15101	569380962
15102012002	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	25	0.0199203	15102	5827095
15102012003	15102	41	2017	Camarones	233083.8	2017	15102	1255	292520169	565	0.4501992	15102	131692347
15102022007	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	39	0.0310757	15102	9090268
15102032005	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	112	0.0892430	15102	26105386
15102032901	15102	8	2017	Camarones	233083.8	2017	15102	1255	292520169	30	0.0239044	15102	6992514
15102042001	15102	21	2017	Camarones	233083.8	2017	15102	1255	292520169	181	0.1442231	15102	42188168
15102052004	15102	24	2017	Camarones	233083.8	2017	15102	1255	292520169	259	0.2063745	15102	60368704
15201011001	15201	34	2017	Putre	251446.5	2017	15201	2765	695249646	1716	0.6206148	15201	431482240
15201022002	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	52	0.0188065	15201	13075219
15201032006	15201	6	2017	Putre	251446.5	2017	15201	2765	695249646	31	0.0112116	15201	7794842
15201042021	15201	10	2017	Putre	251446.5	2017	15201	2765	695249646	75	0.0271248	15201	18858490
15201052001	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	63	0.0227848	15201	15841131
15201052012	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	25	0.0090416	15201	6286163
15201052013	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	17	0.0061483	15201	4274591
15201052020	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	18	0.0065099	15201	4526037
15201052901	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	7	0.0025316	15201	1760126
15201072901	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	53	0.0191682	15201	13326666
15202012008	15202	1	2017	NA	NA	NA	NA	NA	NA	12	0.0175439	15202	NA

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

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 
## -682954922  -99093368  -57133592  121698909  974109494 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 71550816   40608532   1.762   0.0812 .  
## Freq.x       3119383     170994  18.243   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 246500000 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.7743, Adjusted R-squared:  0.772 
## F-statistic: 332.8 on 1 and 97 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

\[ \hat Y = \beta_0 + \beta_1 X^2 \]

linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = multi_pob ~ (Freq.x^2), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -682954922  -99093368  -57133592  121698909  974109494 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 71550816   40608532   1.762   0.0812 .  
## Freq.x       3119383     170994  18.243   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 246500000 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.7743, Adjusted R-squared:  0.772 
## F-statistic: 332.8 on 1 and 97 DF,  p-value: < 2.2e-16

8.2 Modelo cúbico

\[ \hat Y = \beta_0 + \beta_1 X^3 \]

linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = multi_pob ~ (Freq.x^3), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -682954922  -99093368  -57133592  121698909  974109494 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 71550816   40608532   1.762   0.0812 .  
## Freq.x       3119383     170994  18.243   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 246500000 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.7743, Adjusted R-squared:  0.772 
## F-statistic: 332.8 on 1 and 97 DF,  p-value: < 2.2e-16

8.3 Modelo logarítmico

\[ \hat Y = \beta_0 + \beta_1 ln X \]

linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = multi_pob ~ log(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -547163837 -271658885  -41447337  201132079 1421472876 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -355183718   94927848  -3.742  0.00031 ***
## log(Freq.x)  225714391   19732271  11.439  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 338500000 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.5743, Adjusted R-squared:  0.5699 
## F-statistic: 130.8 on 1 and 97 DF,  p-value: < 2.2e-16

8.4 Modelo exponencial

\[ \hat Y = \beta_0 + \beta_1 e^X \]

No es aplicable sin una transformación pues los valores elevados a \(e\) de Freq.x tienden a infinito.

8.5 Modelo con raíz cuadrada

\[ \hat Y = \beta_0 + \beta_1 \sqrt {X} \]

linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = multi_pob ~ sqrt(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -528066559 -172052962   -6157175  106936390 1060111681 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  -194943368   57739148  -3.376  0.00106 ** 
## sqrt(Freq.x)   70067014    4209138  16.646  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 264200000 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.7407, Adjusted R-squared:  0.738 
## F-statistic: 277.1 on 1 and 97 DF,  p-value: < 2.2e-16

8.6 Modelo raíz-raíz

\[ \hat Y = {\beta_0}^2 + 2 \beta_0 \beta_1 \sqrt{X}+ \beta_1^2 X \]

linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = sqrt(multi_pob) ~ sqrt(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -10509.0  -2819.7   -741.2   2233.7  14993.8 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   1929.52     977.71   1.974   0.0513 .  
## sqrt(Freq.x)  1716.62      71.27  24.085   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4474 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.8567, Adjusted R-squared:  0.8553 
## F-statistic: 580.1 on 1 and 97 DF,  p-value: < 2.2e-16

8.7 Modelo log-raíz

\[ \hat Y = e^{\beta_0 + \beta_1 \sqrt{X}} \]

linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = log(multi_pob) ~ sqrt(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -2.5905 -0.3921  0.0514  0.4962  2.3253 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  16.62974    0.17348   95.86   <2e-16 ***
## sqrt(Freq.x)  0.24159    0.01265   19.10   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.7938 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:   0.79,  Adjusted R-squared:  0.7878 
## F-statistic: 364.9 on 1 and 97 DF,  p-value: < 2.2e-16

8.8 Modelo raíz-log

\[ \hat Y = {\beta_0}^2 + 2 \beta_0 \beta_1 \ln{X}+ \beta_1^2 ln^2X \]

linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = sqrt(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -10981.3  -4413.4    117.5   2864.4  15626.0 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  -4337.2     1499.4  -2.893  0.00472 ** 
## log(Freq.x)   6051.2      311.7  19.416  < 2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5347 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.7953, Adjusted R-squared:  0.7932 
## F-statistic:   377 on 1 and 97 DF,  p-value: < 2.2e-16

8.9 Modelo log-log

\[ \hat Y = e^{\beta_0+\beta_1 ln{X}} \]

linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)
summary(linearMod) 

## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.59840 -0.29069 -0.02254  0.26634  1.67931 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  15.3235     0.1496  102.42   <2e-16 ***
## log(Freq.x)   0.9461     0.0311   30.42   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5336 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.9051, Adjusted R-squared:  0.9042 
## F-statistic: 925.5 on 1 and 97 DF,  p-value: < 2.2e-16

9 Modelo log-log (log-log)

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

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_01$Freq.x), y=log(h_y_m_comuna_corr_01$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_01)
summary(linearMod) 

## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -1.59840 -0.29069 -0.02254  0.26634  1.67931 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)  15.3235     0.1496  102.42   <2e-16 ***
## log(Freq.x)   0.9461     0.0311   30.42   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.5336 on 97 degrees of freedom
##   (9 observations deleted due to missingness)
## Multiple R-squared:  0.9051, Adjusted R-squared:  0.9042 
## F-statistic: 925.5 on 1 and 97 DF,  p-value: < 2.2e-16

ggplot(h_y_m_comuna_corr_01, 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.3235+0.9461} \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_01$est_ing <- exp(15.3235+0.9461 * log(h_y_m_comuna_corr_01$Freq.x))

r3_100 <- h_y_m_comuna_corr_01[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
15101011001	15101	150	2017	Arica	292440.1	2017	15101	221364	64735720258	1545	0.0069795	15101	451820024	517263023
15101011002	15101	113	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585	395666252
15101021001	15101	427	2017	Arica	292440.1	2017	15101	221364	64735720258	4454	0.0201207	15101	1302528406	1391743834
15101031001	15101	308	2017	Arica	292440.1	2017	15101	221364	64735720258	6208	0.0280443	15101	1815468420	1021713958
15101031002	15101	188	2017	Arica	292440.1	2017	15101	221364	64735720258	3166	0.0143022	15101	925865499	640460330
15101031003	15101	49	2017	Arica	292440.1	2017	15101	221364	64735720258	623	0.0028144	15101	182190210	179475844
15101031004	15101	490	2017	Arica	292440.1	2017	15101	221364	64735720258	4041	0.0182550	15101	1181750626	1585280088
15101031005	15101	89	2017	Arica	292440.1	2017	15101	221364	64735720258	1010	0.0045626	15101	295364546	315667142
15101031006	15101	206	2017	Arica	292440.1	2017	15101	221364	64735720258	2848	0.0128657	15101	832869533	698330918
15101031007	15101	146	2017	Arica	292440.1	2017	15101	221364	64735720258	2762	0.0124772	15101	807719680	504203354
15101032003	15101	10	2017	Arica	292440.1	2017	15101	221364	64735720258	70	0.0003162	15101	20470810	39903546
15101041001	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2911	0.0131503	15101	851293262	711153199
15101041002	15101	178	2017	Arica	292440.1	2017	15101	221364	64735720258	2836	0.0128115	15101	829360251	608182413
15101041003	15101	345	2017	Arica	292440.1	2017	15101	221364	64735720258	3316	0.0149799	15101	969731521	1137475731
15101041004	15101	312	2017	Arica	292440.1	2017	15101	221364	64735720258	2869	0.0129606	15101	839010776	1034263396
15101051001	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	2041	0.0092201	15101	596870336	536817513
15101051002	15101	137	2017	Arica	292440.1	2017	15101	221364	64735720258	1655	0.0074764	15101	483988440	474747648
15101051003	15101	280	2017	Arica	292440.1	2017	15101	221364	64735720258	2345	0.0105934	15101	685772140	933614755
15101051004	15101	189	2017	Arica	292440.1	2017	15101	221364	64735720258	2196	0.0099203	15101	642198558	643682951
15101061001	15101	403	2017	Arica	292440.1	2017	15101	221364	64735720258	5901	0.0266575	15101	1725689296	1317621272
15101061002	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	2066	0.0093330	15101	604181340	787819300
15101061003	15101	87	2017	Arica	292440.1	2017	15101	221364	64735720258	615	0.0027782	15101	179850689	308951749
15101061004	15101	139	2017	Arica	292440.1	2017	15101	221364	64735720258	723	0.0032661	15101	211434225	481302144
15101061005	15101	46	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585	169062263
15101062011	15101	38	2017	Arica	292440.1	2017	15101	221364	64735720258	2587	0.0116866	15101	756542655	141105764
15101062014	15101	705	2017	Arica	292440.1	2017	15101	221364	64735720258	8708	0.0393379	15101	2546568783	2236573717
15101071001	15101	158	2017	Arica	292440.1	2017	15101	221364	64735720258	2416	0.0109142	15101	706535390	543326595
15101071002	15101	115	2017	Arica	292440.1	2017	15101	221364	64735720258	2443	0.0110361	15101	714431274	402288595
15101071003	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374	787819300
15101071004	15101	138	2017	Arica	292440.1	2017	15101	221364	64735720258	2455	0.0110903	15101	717940556	478025536
15101081001	15101	284	2017	Arica	292440.1	2017	15101	221364	64735720258	2735	0.0123552	15101	799823797	946228392
15101081002	15101	430	2017	Arica	292440.1	2017	15101	221364	64735720258	2496	0.0112755	15101	729930602	1400993109
15101081003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	2823	0.0127528	15101	825558529	987157320
15101081004	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1853	0.0083708	15101	541891589	672645792
15101091001	15101	375	2017	Arica	292440.1	2017	15101	221364	64735720258	2582	0.0116640	15101	755080454	1230842478
15101091002	15101	187	2017	Arica	292440.1	2017	15101	221364	64735720258	1634	0.0073815	15101	477847197	637236784
15101101001	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	1944	0.0087819	15101	568503642	791004215
15101101002	15101	231	2017	Arica	292440.1	2017	15101	221364	64735720258	1732	0.0078242	15101	506506331	778260140
15101101003	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	1542	0.0069659	15101	450942704	765504163
15101111001	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1879	0.0084883	15101	549495032	672645792
15101111002	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	1971	0.0089039	15101	576399526	841865673
15101111003	15101	162	2017	Arica	292440.1	2017	15101	221364	64735720258	1656	0.0074809	15101	484280880	556331497
15101121001	15101	288	2017	Arica	292440.1	2017	15101	221364	64735720258	3627	0.0163848	15101	1060680406	958832456
15101121002	15101	364	2017	Arica	292440.1	2017	15101	221364	64735720258	4645	0.0209835	15101	1358384474	1196656523
15101121003	15101	395	2017	Arica	292440.1	2017	15101	221364	64735720258	4990	0.0225421	15101	1459276324	1292861508
15101121004	15101	214	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374	723962322
15101121005	15101	339	2017	Arica	292440.1	2017	15101	221364	64735720258	4822	0.0217831	15101	1410146379	1118750978
15101121006	15101	33	2017	Arica	292440.1	2017	15101	221364	64735720258	263	0.0011881	15101	76911758	123474573
15101121007	15101	292	2017	Arica	292440.1	2017	15101	221364	64735720258	3635	0.0164209	15101	1063019927	971427087
15101121008	15101	346	2017	Arica	292440.1	2017	15101	221364	64735720258	3757	0.0169720	15101	1098697625	1140594808
15101121009	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	2707	0.0122287	15101	791635473	841865673
15101121010	15101	303	2017	Arica	292440.1	2017	15101	221364	64735720258	3363	0.0151922	15101	983476208	1006014787
15101132015	15101	17	2017	Arica	292440.1	2017	15101	221364	64735720258	338	0.0015269	15101	98844769	65923342
15101141001	15101	220	2017	Arica	292440.1	2017	15101	221364	64735720258	1387	0.0062657	15101	405614481	743151899
15101141002	15101	230	2017	Arica	292440.1	2017	15101	221364	64735720258	1641	0.0074131	15101	479894278	775072270
15101152005	15101	79	2017	Arica	292440.1	2017	15101	221364	64735720258	870	0.0039302	15101	254422926	282004790
15101152007	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	130	0.0005873	15101	38017219	20711284
15101152016	15101	128	2017	Arica	292440.1	2017	15101	221364	64735720258	1759	0.0079462	15101	514402215	445187381
15101152901	15101	24	2017	Arica	292440.1	2017	15101	221364	64735720258	207	0.0009351	15101	60535110	91354375
15101162016	15101	31	2017	Arica	292440.1	2017	15101	221364	64735720258	393	0.0017754	15101	114928977	116382797
15101162901	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	58	0.0002620	15101	16961528	20711284
15101171001	15101	385	2017	Arica	292440.1	2017	15101	221364	64735720258	4176	0.0188649	15101	1221230046	1261873702
15101171002	15101	119	2017	Arica	292440.1	2017	15101	221364	64735720258	2174	0.0098209	15101	635764875	415514777
15101171003	15101	135	2017	Arica	292440.1	2017	15101	221364	64735720258	2503	0.0113072	15101	731977683	468187993
15101171004	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2478	0.0111942	15101	724666679	711153199
15101171005	15101	144	2017	Arica	292440.1	2017	15101	221364	64735720258	2371	0.0107109	15101	693375584	497666316
15101171006	15101	457	2017	Arica	292440.1	2017	15101	221364	64735720258	4924	0.0222439	15101	1439975274	1484083068
15101171007	15101	196	2017	Arica	292440.1	2017	15101	221364	64735720258	2627	0.0118673	15101	768240261	666215853
15101171008	15101	228	2017	Arica	292440.1	2017	15101	221364	64735720258	2763	0.0124817	15101	808012121	768694285
15101171009	15101	362	2017	Arica	292440.1	2017	15101	221364	64735720258	3821	0.0172612	15101	1117413794	1190434959
15101171010	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	3487	0.0157523	15101	1019738786	536817513
15101181001	15101	328	2017	Arica	292440.1	2017	15101	221364	64735720258	3763	0.0169992	15101	1100452266	1084375601
15101181002	15101	384	2017	Arica	292440.1	2017	15101	221364	64735720258	4627	0.0209022	15101	1353120551	1258772553
15101181003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	5211	0.0235404	15101	1523905596	987157320
15101181004	15101	154	2017	Arica	292440.1	2017	15101	221364	64735720258	2346	0.0105979	15101	686064580	530303932
15101181005	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	3954	0.0178620	15101	1156308333	1246363592
15101181006	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	7629	0.0344636	15101	2231025866	1246363592
15101191001	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	2482	0.0112123	15101	725836440	765504163
15101191002	15101	386	2017	Arica	292440.1	2017	15101	221364	64735720258	3364	0.0151967	15101	983768648	1264974417
15101191003	15101	350	2017	Arica	292440.1	2017	15101	221364	64735720258	3067	0.0138550	15101	896913925	1153066275
15101191004	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	2862	0.0129289	15101	836963695	791004215
15101991999	15101	36	2017	Arica	292440.1	2017	15101	221364	64735720258	1947	0.0087955	15101	569380962	134069284
15102012002	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	25	0.0199203	15102	5827095	8703942
15102012003	15102	41	2017	Camarones	233083.8	2017	15102	1255	292520169	565	0.4501992	15102	131692347	151623424
15102022007	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	39	0.0310757	15102	9090268	8703942
15102032005	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	112	0.0892430	15102	26105386	8703942
15102032901	15102	8	2017	Camarones	233083.8	2017	15102	1255	292520169	30	0.0239044	15102	6992514	32309105
15102042001	15102	21	2017	Camarones	233083.8	2017	15102	1255	292520169	181	0.1442231	15102	42188168	80512473
15102052004	15102	24	2017	Camarones	233083.8	2017	15102	1255	292520169	259	0.2063745	15102	60368704	91354375
15201011001	15201	34	2017	Putre	251446.5	2017	15201	2765	695249646	1716	0.6206148	15201	431482240	127011691
15201022002	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	52	0.0188065	15201	13075219	8703942
15201032006	15201	6	2017	Putre	251446.5	2017	15201	2765	695249646	31	0.0112116	15201	7794842	24610497
15201042021	15201	10	2017	Putre	251446.5	2017	15201	2765	695249646	75	0.0271248	15201	18858490	39903546
15201052001	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	63	0.0227848	15201	15841131	4517639
15201052012	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	25	0.0090416	15201	6286163	8703942
15201052013	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	17	0.0061483	15201	4274591	8703942
15201052020	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	18	0.0065099	15201	4526037	4517639
15201052901	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	7	0.0025316	15201	1760126	8703942
15201072901	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	53	0.0191682	15201	13326666	4517639
15202012008	15202	1	2017	NA	NA	NA	NA	NA	NA	12	0.0175439	15202	NA	4517639

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_01$ing_medio_zona <- h_y_m_comuna_corr_01$est_ing  /( h_y_m_comuna_corr_01$personas  * h_y_m_comuna_corr_01$p_poblacional)

r3_100 <- h_y_m_comuna_corr_01[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
15101011001	15101	150	2017	Arica	292440.1	2017	15101	221364	64735720258	1545	0.0069795	15101	451820024	517263023	334798.07
15101011002	15101	113	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585	395666252	395270.98
15101021001	15101	427	2017	Arica	292440.1	2017	15101	221364	64735720258	4454	0.0201207	15101	1302528406	1391743834	312470.55
15101031001	15101	308	2017	Arica	292440.1	2017	15101	221364	64735720258	6208	0.0280443	15101	1815468420	1021713958	164580.21
15101031002	15101	188	2017	Arica	292440.1	2017	15101	221364	64735720258	3166	0.0143022	15101	925865499	640460330	202293.22
15101031003	15101	49	2017	Arica	292440.1	2017	15101	221364	64735720258	623	0.0028144	15101	182190210	179475844	288083.22
15101031004	15101	490	2017	Arica	292440.1	2017	15101	221364	64735720258	4041	0.0182550	15101	1181750626	1585280088	392298.96
15101031005	15101	89	2017	Arica	292440.1	2017	15101	221364	64735720258	1010	0.0045626	15101	295364546	315667142	312541.73
15101031006	15101	206	2017	Arica	292440.1	2017	15101	221364	64735720258	2848	0.0128657	15101	832869533	698330918	245200.46
15101031007	15101	146	2017	Arica	292440.1	2017	15101	221364	64735720258	2762	0.0124772	15101	807719680	504203354	182550.09
15101032003	15101	10	2017	Arica	292440.1	2017	15101	221364	64735720258	70	0.0003162	15101	20470810	39903546	570050.66
15101041001	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2911	0.0131503	15101	851293262	711153199	244298.59
15101041002	15101	178	2017	Arica	292440.1	2017	15101	221364	64735720258	2836	0.0128115	15101	829360251	608182413	214450.78
15101041003	15101	345	2017	Arica	292440.1	2017	15101	221364	64735720258	3316	0.0149799	15101	969731521	1137475731	343026.46
15101041004	15101	312	2017	Arica	292440.1	2017	15101	221364	64735720258	2869	0.0129606	15101	839010776	1034263396	360496.13
15101051001	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	2041	0.0092201	15101	596870336	536817513	263016.91
15101051002	15101	137	2017	Arica	292440.1	2017	15101	221364	64735720258	1655	0.0074764	15101	483988440	474747648	286856.59
15101051003	15101	280	2017	Arica	292440.1	2017	15101	221364	64735720258	2345	0.0105934	15101	685772140	933614755	398129.96
15101051004	15101	189	2017	Arica	292440.1	2017	15101	221364	64735720258	2196	0.0099203	15101	642198558	643682951	293116.10
15101061001	15101	403	2017	Arica	292440.1	2017	15101	221364	64735720258	5901	0.0266575	15101	1725689296	1317621272	223287.79
15101061002	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	2066	0.0093330	15101	604181340	787819300	381325.90
15101061003	15101	87	2017	Arica	292440.1	2017	15101	221364	64735720258	615	0.0027782	15101	179850689	308951749	502360.57
15101061004	15101	139	2017	Arica	292440.1	2017	15101	221364	64735720258	723	0.0032661	15101	211434225	481302144	665701.44
15101061005	15101	46	2017	Arica	292440.1	2017	15101	221364	64735720258	1001	0.0045220	15101	292732585	169062263	168893.37
15101062011	15101	38	2017	Arica	292440.1	2017	15101	221364	64735720258	2587	0.0116866	15101	756542655	141105764	54544.17
15101062014	15101	705	2017	Arica	292440.1	2017	15101	221364	64735720258	8708	0.0393379	15101	2546568783	2236573717	256841.26
15101071001	15101	158	2017	Arica	292440.1	2017	15101	221364	64735720258	2416	0.0109142	15101	706535390	543326595	224886.84
15101071002	15101	115	2017	Arica	292440.1	2017	15101	221364	64735720258	2443	0.0110361	15101	714431274	402288595	164669.91
15101071003	15101	234	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374	787819300	213732.85
15101071004	15101	138	2017	Arica	292440.1	2017	15101	221364	64735720258	2455	0.0110903	15101	717940556	478025536	194715.09
15101081001	15101	284	2017	Arica	292440.1	2017	15101	221364	64735720258	2735	0.0123552	15101	799823797	946228392	345970.16
15101081002	15101	430	2017	Arica	292440.1	2017	15101	221364	64735720258	2496	0.0112755	15101	729930602	1400993109	561295.32
15101081003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	2823	0.0127528	15101	825558529	987157320	349683.78
15101081004	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1853	0.0083708	15101	541891589	672645792	363003.67
15101091001	15101	375	2017	Arica	292440.1	2017	15101	221364	64735720258	2582	0.0116640	15101	755080454	1230842478	476701.19
15101091002	15101	187	2017	Arica	292440.1	2017	15101	221364	64735720258	1634	0.0073815	15101	477847197	637236784	389985.79
15101101001	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	1944	0.0087819	15101	568503642	791004215	406895.17
15101101002	15101	231	2017	Arica	292440.1	2017	15101	221364	64735720258	1732	0.0078242	15101	506506331	778260140	449341.88
15101101003	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	1542	0.0069659	15101	450942704	765504163	496435.90
15101111001	15101	198	2017	Arica	292440.1	2017	15101	221364	64735720258	1879	0.0084883	15101	549495032	672645792	357980.73
15101111002	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	1971	0.0089039	15101	576399526	841865673	427126.17
15101111003	15101	162	2017	Arica	292440.1	2017	15101	221364	64735720258	1656	0.0074809	15101	484280880	556331497	335948.97
15101121001	15101	288	2017	Arica	292440.1	2017	15101	221364	64735720258	3627	0.0163848	15101	1060680406	958832456	264359.65
15101121002	15101	364	2017	Arica	292440.1	2017	15101	221364	64735720258	4645	0.0209835	15101	1358384474	1196656523	257622.50
15101121003	15101	395	2017	Arica	292440.1	2017	15101	221364	64735720258	4990	0.0225421	15101	1459276324	1292861508	259090.48
15101121004	15101	214	2017	Arica	292440.1	2017	15101	221364	64735720258	3686	0.0166513	15101	1077934374	723962322	196408.66
15101121005	15101	339	2017	Arica	292440.1	2017	15101	221364	64735720258	4822	0.0217831	15101	1410146379	1118750978	232009.74
15101121006	15101	33	2017	Arica	292440.1	2017	15101	221364	64735720258	263	0.0011881	15101	76911758	123474573	469485.07
15101121007	15101	292	2017	Arica	292440.1	2017	15101	221364	64735720258	3635	0.0164209	15101	1063019927	971427087	267242.67
15101121008	15101	346	2017	Arica	292440.1	2017	15101	221364	64735720258	3757	0.0169720	15101	1098697625	1140594808	303591.91
15101121009	15101	251	2017	Arica	292440.1	2017	15101	221364	64735720258	2707	0.0122287	15101	791635473	841865673	310995.82
15101121010	15101	303	2017	Arica	292440.1	2017	15101	221364	64735720258	3363	0.0151922	15101	983476208	1006014787	299142.07
15101132015	15101	17	2017	Arica	292440.1	2017	15101	221364	64735720258	338	0.0015269	15101	98844769	65923342	195039.47
15101141001	15101	220	2017	Arica	292440.1	2017	15101	221364	64735720258	1387	0.0062657	15101	405614481	743151899	535798.05
15101141002	15101	230	2017	Arica	292440.1	2017	15101	221364	64735720258	1641	0.0074131	15101	479894278	775072270	472317.04
15101152005	15101	79	2017	Arica	292440.1	2017	15101	221364	64735720258	870	0.0039302	15101	254422926	282004790	324143.44
15101152007	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	130	0.0005873	15101	38017219	20711284	159317.57
15101152016	15101	128	2017	Arica	292440.1	2017	15101	221364	64735720258	1759	0.0079462	15101	514402215	445187381	253091.18
15101152901	15101	24	2017	Arica	292440.1	2017	15101	221364	64735720258	207	0.0009351	15101	60535110	91354375	441325.48
15101162016	15101	31	2017	Arica	292440.1	2017	15101	221364	64735720258	393	0.0017754	15101	114928977	116382797	296139.43
15101162901	15101	5	2017	Arica	292440.1	2017	15101	221364	64735720258	58	0.0002620	15101	16961528	20711284	357091.10
15101171001	15101	385	2017	Arica	292440.1	2017	15101	221364	64735720258	4176	0.0188649	15101	1221230046	1261873702	302172.82
15101171002	15101	119	2017	Arica	292440.1	2017	15101	221364	64735720258	2174	0.0098209	15101	635764875	415514777	191129.15
15101171003	15101	135	2017	Arica	292440.1	2017	15101	221364	64735720258	2503	0.0113072	15101	731977683	468187993	187050.74
15101171004	15101	210	2017	Arica	292440.1	2017	15101	221364	64735720258	2478	0.0111942	15101	724666679	711153199	286986.76
15101171005	15101	144	2017	Arica	292440.1	2017	15101	221364	64735720258	2371	0.0107109	15101	693375584	497666316	209897.22
15101171006	15101	457	2017	Arica	292440.1	2017	15101	221364	64735720258	4924	0.0222439	15101	1439975274	1484083068	301397.86
15101171007	15101	196	2017	Arica	292440.1	2017	15101	221364	64735720258	2627	0.0118673	15101	768240261	666215853	253603.29
15101171008	15101	228	2017	Arica	292440.1	2017	15101	221364	64735720258	2763	0.0124817	15101	808012121	768694285	278210.02
15101171009	15101	362	2017	Arica	292440.1	2017	15101	221364	64735720258	3821	0.0172612	15101	1117413794	1190434959	311550.63
15101171010	15101	156	2017	Arica	292440.1	2017	15101	221364	64735720258	3487	0.0157523	15101	1019738786	536817513	153948.24
15101181001	15101	328	2017	Arica	292440.1	2017	15101	221364	64735720258	3763	0.0169992	15101	1100452266	1084375601	288167.84
15101181002	15101	384	2017	Arica	292440.1	2017	15101	221364	64735720258	4627	0.0209022	15101	1353120551	1258772553	272049.40
15101181003	15101	297	2017	Arica	292440.1	2017	15101	221364	64735720258	5211	0.0235404	15101	1523905596	987157320	189437.21
15101181004	15101	154	2017	Arica	292440.1	2017	15101	221364	64735720258	2346	0.0105979	15101	686064580	530303932	226046.01
15101181005	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	3954	0.0178620	15101	1156308333	1246363592	315215.88
15101181006	15101	380	2017	Arica	292440.1	2017	15101	221364	64735720258	7629	0.0344636	15101	2231025866	1246363592	163371.82
15101191001	15101	227	2017	Arica	292440.1	2017	15101	221364	64735720258	2482	0.0112123	15101	725836440	765504163	308422.31
15101191002	15101	386	2017	Arica	292440.1	2017	15101	221364	64735720258	3364	0.0151967	15101	983768648	1264974417	376032.82
15101191003	15101	350	2017	Arica	292440.1	2017	15101	221364	64735720258	3067	0.0138550	15101	896913925	1153066275	375959.01
15101191004	15101	235	2017	Arica	292440.1	2017	15101	221364	64735720258	2862	0.0129289	15101	836963695	791004215	276381.63
15101991999	15101	36	2017	Arica	292440.1	2017	15101	221364	64735720258	1947	0.0087955	15101	569380962	134069284	68859.42
15102012002	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	25	0.0199203	15102	5827095	8703942	348157.67
15102012003	15102	41	2017	Camarones	233083.8	2017	15102	1255	292520169	565	0.4501992	15102	131692347	151623424	268360.04
15102022007	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	39	0.0310757	15102	9090268	8703942	223177.99
15102032005	15102	2	2017	Camarones	233083.8	2017	15102	1255	292520169	112	0.0892430	15102	26105386	8703942	77713.76
15102032901	15102	8	2017	Camarones	233083.8	2017	15102	1255	292520169	30	0.0239044	15102	6992514	32309105	1076970.17
15102042001	15102	21	2017	Camarones	233083.8	2017	15102	1255	292520169	181	0.1442231	15102	42188168	80512473	444820.29
15102052004	15102	24	2017	Camarones	233083.8	2017	15102	1255	292520169	259	0.2063745	15102	60368704	91354375	352719.59
15201011001	15201	34	2017	Putre	251446.5	2017	15201	2765	695249646	1716	0.6206148	15201	431482240	127011691	74016.14
15201022002	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	52	0.0188065	15201	13075219	8703942	167383.49
15201032006	15201	6	2017	Putre	251446.5	2017	15201	2765	695249646	31	0.0112116	15201	7794842	24610497	793887.01
15201042021	15201	10	2017	Putre	251446.5	2017	15201	2765	695249646	75	0.0271248	15201	18858490	39903546	532047.28
15201052001	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	63	0.0227848	15201	15841131	4517639	71708.55
15201052012	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	25	0.0090416	15201	6286163	8703942	348157.67
15201052013	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	17	0.0061483	15201	4274591	8703942	511996.57
15201052020	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	18	0.0065099	15201	4526037	4517639	250979.93
15201052901	15201	2	2017	Putre	251446.5	2017	15201	2765	695249646	7	0.0025316	15201	1760126	8703942	1243420.24
15201072901	15201	1	2017	Putre	251446.5	2017	15201	2765	695249646	53	0.0191682	15201	13326666	4517639	85238.47
15202012008	15202	1	2017	NA	NA	NA	NA	NA	NA	12	0.0175439	15202	NA	4517639	NA

Guardamos:

saveRDS(h_y_m_comuna_corr_01, "casen_censo_region_15.rds")

12 Referencias

https://rpubs.com/osoramirez/316691

https://dataintelligencechile.shinyapps.io/casenfinal

Manual_de_usuario_Censo_2017_16R.pdf

http://www.censo2017.cl/microdatos/

Censo de Población y Vivienda

https://www.ine.cl/estadisticas/sociales/censos-de-poblacion-y-vivienda/poblacion-y-vivienda

http://r-statistics.co/Linear-Regression.html