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

1 Variable CENSO

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

1.1 Lectura y filtrado de la tabla censal de viviendas

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

tabla_con_clave <- readRDS("../censo_viviendas_con_clave_17.rds")
r3_100 <- tabla_con_clave[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

REGION	PROVINCIA	COMUNA	DC	AREA	ZC_LOC	ID_ZONA_LOC	NVIV	P01	P02	P03A	P03B	P03C	P04	P05	CANT_HOG	CANT_PER	REGION_15R	PROVINCIA_15R	COMUNA_15R	clave
15	152	15202	1	2	6	13225	1	3	1	5	3	5	1	4	1	1	15	152	15202	15202012006
15	152	15202	1	2	6	13225	2	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	3	1	1	5	3	5	2	3	1	4	15	152	15202	15202012006
15	152	15202	1	2	6	13225	4	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	5	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	6	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	7	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	8	3	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	9	3	1	5	3	5	1	4	1	4	15	152	15202	15202012006
15	152	15202	1	2	6	13225	10	1	1	5	3	4	1	4	1	1	15	152	15202	15202012006
15	152	15202	1	2	6	13225	11	1	2	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	12	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	13	1	1	5	3	4	1	4	1	3	15	152	15202	15202012006
15	152	15202	1	2	6	13225	14	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	15	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	16	1	1	5	3	5	1	3	1	3	15	152	15202	15202012006
15	152	15202	1	2	6	13225	17	1	1	5	3	5	2	4	1	8	15	152	15202	15202012006
15	152	15202	1	2	6	13225	18	3	1	5	3	5	1	1	1	1	15	152	15202	15202012006
15	152	15202	1	2	6	13225	19	3	1	5	3	5	1	3	1	1	15	152	15202	15202012006
15	152	15202	1	2	6	13225	20	1	1	5	3	5	2	1	1	2	15	152	15202	15202012006
15	152	15202	1	2	6	13225	21	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	22	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	23	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	24	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	25	1	1	5	3	5	1	4	1	2	15	152	15202	15202012006
15	152	15202	1	2	6	13225	26	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	27	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	28	3	1	5	3	5	1	4	1	5	15	152	15202	15202012006
15	152	15202	1	2	6	13225	29	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	30	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	31	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	32	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	33	1	1	5	3	5	3	4	1	4	15	152	15202	15202012006
15	152	15202	1	2	6	13225	34	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	35	3	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	36	1	1	5	3	5	3	2	1	9	15	152	15202	15202012006
15	152	15202	1	2	6	13225	37	3	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	38	1	1	5	3	5	99	4	1	1	15	152	15202	15202012006
15	152	15202	1	2	6	13225	39	1	1	5	3	5	1	4	1	2	15	152	15202	15202012006
15	152	15202	1	2	6	13225	40	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	41	1	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	42	3	3	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	43	3	1	5	3	5	2	1	1	5	15	152	15202	15202012006
15	152	15202	1	2	6	13225	44	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	6	13225	45	1	4	98	98	98	98	98	0	0	15	152	15202	15202012006
15	152	15202	1	2	8	13910	1	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	2	1	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	3	1	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	4	3	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	5	3	1	5	3	5	2	3	1	3	15	152	15202	15202012008
15	152	15202	1	2	8	13910	6	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	7	1	1	5	99	5	2	4	1	2	15	152	15202	15202012008
15	152	15202	1	2	8	13910	8	3	1	5	3	5	3	3	1	2	15	152	15202	15202012008
15	152	15202	1	2	8	13910	9	3	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	10	3	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	11	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	12	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	13	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	14	3	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	15	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	16	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	17	3	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	18	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	19	3	1	99	99	99	99	99	1	1	15	152	15202	15202012008
15	152	15202	1	2	8	13910	20	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	21	1	1	5	3	5	1	4	1	2	15	152	15202	15202012008
15	152	15202	1	2	8	13910	22	3	1	5	3	5	2	4	1	1	15	152	15202	15202012008
15	152	15202	1	2	8	13910	23	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	24	3	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	25	3	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	26	1	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	27	1	4	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	28	3	2	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	29	1	3	98	98	98	98	98	0	0	15	152	15202	15202012008
15	152	15202	1	2	8	13910	30	1	1	5	1	4	2	4	1	1	15	152	15202	15202012008
15	152	15202	1	2	12	8394	1	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	2	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	3	3	1	5	3	5	2	3	1	4	15	152	15202	15202012012
15	152	15202	1	2	12	8394	4	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	5	3	3	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	6	3	3	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	7	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	8	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	9	3	1	5	3	5	1	4	1	1	15	152	15202	15202012012
15	152	15202	1	2	12	8394	10	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	11	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	12	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	13	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	14	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	15	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	16	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	17	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	18	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	19	3	1	99	99	99	99	99	1	1	15	152	15202	15202012012
15	152	15202	1	2	12	8394	20	3	1	5	3	5	3	99	1	1	15	152	15202	15202012012
15	152	15202	1	2	12	8394	21	3	1	5	99	5	1	4	1	2	15	152	15202	15202012012
15	152	15202	1	2	12	8394	22	3	2	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	23	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012
15	152	15202	1	2	12	8394	24	3	1	5	3	5	1	2	1	2	15	152	15202	15202012012
15	152	15202	1	2	12	8394	25	3	4	98	98	98	98	98	0	0	15	152	15202	15202012012

Despleguemos los códigos de regiones de nuestra tabla:

regiones <- unique(tabla_con_clave$REGION)

Hagamos un subset con la 1:

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

1.2 Cálculo de frecuencias

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

aterial de construcción del piso

names(tabla_con_clave_f)[2] <- "Tipo de piso"

# Ahora filtramos por Tipo de techo = 1.
tabla_con_clave_ff <- filter(tabla_con_clave_f, tabla_con_clave_f$`Tipo de piso` == 1)

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

Veamos los primeros 100 registros:

r3_100 <- d[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

	zona	unlist.c.	unlist.d.	Freq	anio
1	7101011001	1	7101	417	2017
2	7101011002	1	7101	687	2017
3	7101021001	1	7101	408	2017
4	7101021002	1	7101	335	2017
5	7101021003	1	7101	262	2017
6	7101031001	1	7101	610	2017
7	7101031002	1	7101	415	2017
8	7101031003	1	7101	318	2017
9	7101041001	1	7101	261	2017
10	7101041002	1	7101	506	2017
11	7101041003	1	7101	557	2017
12	7101041004	1	7101	506	2017
13	7101041005	1	7101	397	2017
14	7101051001	1	7101	969	2017
15	7101051002	1	7101	336	2017
16	7101061001	1	7101	644	2017
17	7101061002	1	7101	383	2017
18	7101061003	1	7101	826	2017
19	7101061004	1	7101	834	2017
20	7101061005	1	7101	842	2017
21	7101061006	1	7101	994	2017
22	7101071001	1	7101	676	2017
23	7101071002	1	7101	1117	2017
24	7101071003	1	7101	539	2017
25	7101071004	1	7101	559	2017
26	7101071005	1	7101	609	2017
27	7101071006	1	7101	854	2017
28	7101071007	1	7101	1028	2017
29	7101071008	1	7101	894	2017
30	7101071009	1	7101	1185	2017
31	7101071010	1	7101	684	2017
32	7101081001	1	7101	331	2017
33	7101081002	1	7101	587	2017
34	7101081003	1	7101	468	2017
35	7101081004	1	7101	664	2017
36	7101081005	1	7101	893	2017
37	7101081006	1	7101	539	2017
38	7101091001	1	7101	40	2017
39	7101101001	1	7101	378	2017
40	7101111001	1	7101	1458	2017
41	7101111002	1	7101	1044	2017
42	7101111003	1	7101	927	2017
43	7101111004	1	7101	102	2017
44	7101111005	1	7101	1674	2017
45	7101111006	1	7101	2449	2017
46	7101111007	1	7101	509	2017
47	7101111008	1	7101	1501	2017
48	7101111009	1	7101	1179	2017
49	7101111010	1	7101	566	2017
50	7101121001	1	7101	1182	2017
51	7101121002	1	7101	540	2017
52	7101121003	1	7101	914	2017
53	7101121004	1	7101	1224	2017
54	7101121005	1	7101	97	2017
55	7101121006	1	7101	1899	2017
56	7101121007	1	7101	715	2017
57	7101121008	1	7101	634	2017
58	7101121009	1	7101	894	2017
59	7101121010	1	7101	643	2017
60	7101131001	1	7101	514	2017
61	7101131002	1	7101	1408	2017
62	7101131003	1	7101	1030	2017
63	7101131004	1	7101	1087	2017
64	7101131005	1	7101	923	2017
65	7101131006	1	7101	607	2017
66	7101131007	1	7101	859	2017
67	7101131008	1	7101	1075	2017
68	7101141001	1	7101	811	2017
69	7101141002	1	7101	835	2017
70	7101141003	1	7101	753	2017
71	7101141004	1	7101	1284	2017
72	7101141005	1	7101	1119	2017
73	7101141006	1	7101	1514	2017
74	7101141007	1	7101	1413	2017
75	7101141008	1	7101	1505	2017
76	7101141009	1	7101	1110	2017
77	7101141010	1	7101	658	2017
78	7101151001	1	7101	79	2017
79	7101161001	1	7101	341	2017
80	7101991999	1	7101	306	2017
396	7102011001	1	7102	535	2017
397	7102011002	1	7102	658	2017
398	7102021001	1	7102	1496	2017
399	7102021002	1	7102	1123	2017
400	7102021003	1	7102	1651	2017
401	7102021004	1	7102	1588	2017
402	7102021005	1	7102	760	2017
403	7102031001	1	7102	1357	2017
404	7102031002	1	7102	1421	2017
405	7102041001	1	7102	615	2017
406	7102071001	1	7102	367	2017
407	7102991999	1	7102	23	2017
723	7103011001	1	7103	1113	2017
724	7103991999	1	7103	10	2017
1040	7104011001	1	7104	895	2017
1041	7104021001	1	7104	9	2017
1042	7104991999	1	7104	2	2017
1358	7105011001	1	7105	1273	2017
1359	7105051001	1	7105	221	2017
1360	7105051002	1	7105	987	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	7101011001	417	2017	07101
2	7101011002	687	2017	07101
3	7101021001	408	2017	07101
4	7101021002	335	2017	07101
5	7101021003	262	2017	07101
6	7101031001	610	2017	07101
7	7101031002	415	2017	07101
8	7101031003	318	2017	07101
9	7101041001	261	2017	07101
10	7101041002	506	2017	07101
11	7101041003	557	2017	07101
12	7101041004	506	2017	07101
13	7101041005	397	2017	07101
14	7101051001	969	2017	07101
15	7101051002	336	2017	07101
16	7101061001	644	2017	07101
17	7101061002	383	2017	07101
18	7101061003	826	2017	07101
19	7101061004	834	2017	07101
20	7101061005	842	2017	07101
21	7101061006	994	2017	07101
22	7101071001	676	2017	07101
23	7101071002	1117	2017	07101
24	7101071003	539	2017	07101
25	7101071004	559	2017	07101
26	7101071005	609	2017	07101
27	7101071006	854	2017	07101
28	7101071007	1028	2017	07101
29	7101071008	894	2017	07101
30	7101071009	1185	2017	07101
31	7101071010	684	2017	07101
32	7101081001	331	2017	07101
33	7101081002	587	2017	07101
34	7101081003	468	2017	07101
35	7101081004	664	2017	07101
36	7101081005	893	2017	07101
37	7101081006	539	2017	07101
38	7101091001	40	2017	07101
39	7101101001	378	2017	07101
40	7101111001	1458	2017	07101
41	7101111002	1044	2017	07101
42	7101111003	927	2017	07101
43	7101111004	102	2017	07101
44	7101111005	1674	2017	07101
45	7101111006	2449	2017	07101
46	7101111007	509	2017	07101
47	7101111008	1501	2017	07101
48	7101111009	1179	2017	07101
49	7101111010	566	2017	07101
50	7101121001	1182	2017	07101
51	7101121002	540	2017	07101
52	7101121003	914	2017	07101
53	7101121004	1224	2017	07101
54	7101121005	97	2017	07101
55	7101121006	1899	2017	07101
56	7101121007	715	2017	07101
57	7101121008	634	2017	07101
58	7101121009	894	2017	07101
59	7101121010	643	2017	07101
60	7101131001	514	2017	07101
61	7101131002	1408	2017	07101
62	7101131003	1030	2017	07101
63	7101131004	1087	2017	07101
64	7101131005	923	2017	07101
65	7101131006	607	2017	07101
66	7101131007	859	2017	07101
67	7101131008	1075	2017	07101
68	7101141001	811	2017	07101
69	7101141002	835	2017	07101
70	7101141003	753	2017	07101
71	7101141004	1284	2017	07101
72	7101141005	1119	2017	07101
73	7101141006	1514	2017	07101
74	7101141007	1413	2017	07101
75	7101141008	1505	2017	07101
76	7101141009	1110	2017	07101
77	7101141010	658	2017	07101
78	7101151001	79	2017	07101
79	7101161001	341	2017	07101
80	7101991999	306	2017	07101
396	7102011001	535	2017	07102
397	7102011002	658	2017	07102
398	7102021001	1496	2017	07102
399	7102021002	1123	2017	07102
400	7102021003	1651	2017	07102
401	7102021004	1588	2017	07102
402	7102021005	760	2017	07102
403	7102031001	1357	2017	07102
404	7102031002	1421	2017	07102
405	7102041001	615	2017	07102
406	7102071001	367	2017	07102
407	7102991999	23	2017	07102
723	7103011001	1113	2017	07103
724	7103991999	10	2017	07103
1040	7104011001	895	2017	07104
1041	7104021001	9	2017	07104
1042	7104991999	2	2017	07104
1358	7105011001	1273	2017	07105
1359	7105051001	221	2017	07105
1360	7105051002	987	2017	07105

2 Variable CASEN

2.1 Tabla de ingresos expandidos

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

h_y_m_2017_censo <- readRDS("../ingresos_expandidos_urbano_17.rds")
tablamadre <- head(h_y_m_2017_censo,50)
kbl(tablamadre) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

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

3 Unión Censo-Casen

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

comunas_con_ing_exp = merge( x = comuna_corr, y = h_y_m_2017_censo, by = "código", all.x = TRUE)
comunas_con_ing_exp <-comunas_con_ing_exp[!(is.na(comunas_con_ing_exp$Ingresos_expandidos)),]
r3_100 <- comunas_con_ing_exp[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

código	zona	Freq	anio	comuna.x	promedio_i	año	comuna.y	personas	Ingresos_expandidos
07101	7101011001	417	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101011002	687	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021001	408	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021002	335	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021003	262	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031001	610	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031002	415	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031003	318	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041001	261	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041002	506	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041003	557	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041004	506	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041005	397	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101051001	969	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101051002	336	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061001	644	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061002	383	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061003	826	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061004	834	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061005	842	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061006	994	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071001	676	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071002	1117	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071003	539	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071004	559	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071005	609	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071006	854	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071007	1028	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071008	894	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071009	1185	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071010	684	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081001	331	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081002	587	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081003	468	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081004	664	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081005	893	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081006	539	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101091001	40	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101101001	378	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111001	1458	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111002	1044	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111003	927	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111004	102	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111005	1674	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111006	2449	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111007	509	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111008	1501	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111009	1179	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111010	566	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121001	1182	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121002	540	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121003	914	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121004	1224	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121005	97	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121006	1899	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121007	715	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121008	634	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121009	894	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121010	643	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131001	514	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131002	1408	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131003	1030	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131004	1087	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131005	923	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131006	607	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131007	859	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131008	1075	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141001	811	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141002	835	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141003	753	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141004	1284	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141005	1119	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141006	1514	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141007	1413	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141008	1505	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141009	1110	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141010	658	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101151001	79	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101161001	341	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101991999	306	2017	Talca	307377.4	2017	7101	220357	67732753814
07102	7102011001	535	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102011002	658	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021001	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021002	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021003	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021004	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021005	760	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102031001	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102031002	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102041001	615	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102071001	367	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102991999	23	2017	Constitución	280736.9	2017	7102	46068	12932986800
07103	7103011001	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120
07103	7103991999	10	2017	Curepto	281855.5	2017	7103	9448	2662971120
07104	7104011001	895	2017	Empedrado	209235.2	2017	7104	4142	866652110
07104	7104021001	9	2017	Empedrado	209235.2	2017	7104	4142	866652110
07104	7104991999	2	2017	Empedrado	209235.2	2017	7104	4142	866652110
07105	7105011001	1273	2017	Maule	245019.7	2017	7105	49721	12182624190
07105	7105051001	221	2017	Maule	245019.7	2017	7105	49721	12182624190
07105	7105051002	987	2017	Maule	245019.7	2017	7105	49721	12182624190

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
07101	7101011001	417	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101011002	687	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021001	408	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021002	335	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101021003	262	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031001	610	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031002	415	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101031003	318	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041001	261	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041002	506	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041003	557	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041004	506	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101041005	397	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101051001	969	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101051002	336	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061001	644	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061002	383	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061003	826	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061004	834	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061005	842	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101061006	994	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071001	676	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071002	1117	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071003	539	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071004	559	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071005	609	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071006	854	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071007	1028	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071008	894	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071009	1185	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101071010	684	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081001	331	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081002	587	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081003	468	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081004	664	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081005	893	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101081006	539	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101091001	40	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101101001	378	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111001	1458	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111002	1044	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111003	927	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111004	102	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111005	1674	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111006	2449	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111007	509	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111008	1501	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111009	1179	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101111010	566	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121001	1182	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121002	540	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121003	914	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121004	1224	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121005	97	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121006	1899	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121007	715	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121008	634	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121009	894	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101121010	643	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131001	514	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131002	1408	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131003	1030	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131004	1087	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131005	923	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131006	607	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131007	859	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101131008	1075	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141001	811	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141002	835	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141003	753	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141004	1284	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141005	1119	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141006	1514	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141007	1413	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141008	1505	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141009	1110	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101141010	658	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101151001	79	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101161001	341	2017	Talca	307377.4	2017	7101	220357	67732753814
07101	7101991999	306	2017	Talca	307377.4	2017	7101	220357	67732753814
07102	7102011001	535	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102011002	658	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021001	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021002	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021003	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021004	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102021005	760	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102031001	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102031002	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102041001	615	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102071001	367	2017	Constitución	280736.9	2017	7102	46068	12932986800
07102	7102991999	23	2017	Constitución	280736.9	2017	7102	46068	12932986800
07103	7103011001	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120
07103	7103991999	10	2017	Curepto	281855.5	2017	7103	9448	2662971120
07104	7104011001	895	2017	Empedrado	209235.2	2017	7104	4142	866652110
07104	7104021001	9	2017	Empedrado	209235.2	2017	7104	4142	866652110
07104	7104991999	2	2017	Empedrado	209235.2	2017	7104	4142	866652110
07105	7105011001	1273	2017	Maule	245019.7	2017	7105	49721	12182624190
07105	7105051001	221	2017	Maule	245019.7	2017	7105	49721	12182624190
07105	7105051002	987	2017	Maule	245019.7	2017	7105	49721	12182624190

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
7101011001	07101	417	2017	Talca	307377.4	2017	7101	220357	67732753814	1693	0.0076830	07101
7101011002	07101	687	2017	Talca	307377.4	2017	7101	220357	67732753814	1889	0.0085725	07101
7101021001	07101	408	2017	Talca	307377.4	2017	7101	220357	67732753814	951	0.0043157	07101
7101021002	07101	335	2017	Talca	307377.4	2017	7101	220357	67732753814	912	0.0041387	07101
7101021003	07101	262	2017	Talca	307377.4	2017	7101	220357	67732753814	857	0.0038891	07101
7101031001	07101	610	2017	Talca	307377.4	2017	7101	220357	67732753814	1465	0.0066483	07101
7101031002	07101	415	2017	Talca	307377.4	2017	7101	220357	67732753814	1428	0.0064804	07101
7101031003	07101	318	2017	Talca	307377.4	2017	7101	220357	67732753814	1037	0.0047060	07101
7101041001	07101	261	2017	Talca	307377.4	2017	7101	220357	67732753814	1306	0.0059267	07101
7101041002	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1682	0.0076331	07101
7101041003	07101	557	2017	Talca	307377.4	2017	7101	220357	67732753814	1627	0.0073835	07101
7101041004	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1684	0.0076421	07101
7101041005	07101	397	2017	Talca	307377.4	2017	7101	220357	67732753814	1196	0.0054276	07101
7101051001	07101	969	2017	Talca	307377.4	2017	7101	220357	67732753814	2870	0.0130243	07101
7101051002	07101	336	2017	Talca	307377.4	2017	7101	220357	67732753814	1049	0.0047605	07101
7101061001	07101	644	2017	Talca	307377.4	2017	7101	220357	67732753814	1756	0.0079689	07101
7101061002	07101	383	2017	Talca	307377.4	2017	7101	220357	67732753814	1371	0.0062217	07101
7101061003	07101	826	2017	Talca	307377.4	2017	7101	220357	67732753814	3032	0.0137595	07101
7101061004	07101	834	2017	Talca	307377.4	2017	7101	220357	67732753814	3215	0.0145900	07101
7101061005	07101	842	2017	Talca	307377.4	2017	7101	220357	67732753814	2884	0.0130879	07101
7101061006	07101	994	2017	Talca	307377.4	2017	7101	220357	67732753814	3367	0.0152798	07101
7101071001	07101	676	2017	Talca	307377.4	2017	7101	220357	67732753814	2136	0.0096934	07101
7101071002	07101	1117	2017	Talca	307377.4	2017	7101	220357	67732753814	3210	0.0145673	07101
7101071003	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1773	0.0080460	07101
7101071004	07101	559	2017	Talca	307377.4	2017	7101	220357	67732753814	1769	0.0080279	07101
7101071005	07101	609	2017	Talca	307377.4	2017	7101	220357	67732753814	1956	0.0088765	07101
7101071006	07101	854	2017	Talca	307377.4	2017	7101	220357	67732753814	3097	0.0140545	07101
7101071007	07101	1028	2017	Talca	307377.4	2017	7101	220357	67732753814	3759	0.0170587	07101
7101071008	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2586	0.0117355	07101
7101071009	07101	1185	2017	Talca	307377.4	2017	7101	220357	67732753814	3775	0.0171313	07101
7101071010	07101	684	2017	Talca	307377.4	2017	7101	220357	67732753814	2262	0.0102652	07101
7101081001	07101	331	2017	Talca	307377.4	2017	7101	220357	67732753814	1170	0.0053096	07101
7101081002	07101	587	2017	Talca	307377.4	2017	7101	220357	67732753814	1841	0.0083546	07101
7101081003	07101	468	2017	Talca	307377.4	2017	7101	220357	67732753814	1281	0.0058133	07101
7101081004	07101	664	2017	Talca	307377.4	2017	7101	220357	67732753814	1562	0.0070885	07101
7101081005	07101	893	2017	Talca	307377.4	2017	7101	220357	67732753814	3349	0.0151981	07101
7101081006	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1555	0.0070567	07101
7101091001	07101	40	2017	Talca	307377.4	2017	7101	220357	67732753814	179	0.0008123	07101
7101101001	07101	378	2017	Talca	307377.4	2017	7101	220357	67732753814	1347	0.0061128	07101
7101111001	07101	1458	2017	Talca	307377.4	2017	7101	220357	67732753814	5385	0.0244376	07101
7101111002	07101	1044	2017	Talca	307377.4	2017	7101	220357	67732753814	3446	0.0156383	07101
7101111003	07101	927	2017	Talca	307377.4	2017	7101	220357	67732753814	3435	0.0155883	07101
7101111004	07101	102	2017	Talca	307377.4	2017	7101	220357	67732753814	347	0.0015747	07101
7101111005	07101	1674	2017	Talca	307377.4	2017	7101	220357	67732753814	5190	0.0235527	07101
7101111006	07101	2449	2017	Talca	307377.4	2017	7101	220357	67732753814	7789	0.0353472	07101
7101111007	07101	509	2017	Talca	307377.4	2017	7101	220357	67732753814	1912	0.0086768	07101
7101111008	07101	1501	2017	Talca	307377.4	2017	7101	220357	67732753814	5035	0.0228493	07101
7101111009	07101	1179	2017	Talca	307377.4	2017	7101	220357	67732753814	3737	0.0169588	07101
7101111010	07101	566	2017	Talca	307377.4	2017	7101	220357	67732753814	1695	0.0076921	07101
7101121001	07101	1182	2017	Talca	307377.4	2017	7101	220357	67732753814	4067	0.0184564	07101
7101121002	07101	540	2017	Talca	307377.4	2017	7101	220357	67732753814	1855	0.0084182	07101
7101121003	07101	914	2017	Talca	307377.4	2017	7101	220357	67732753814	3773	0.0171222	07101
7101121004	07101	1224	2017	Talca	307377.4	2017	7101	220357	67732753814	4426	0.0200856	07101
7101121005	07101	97	2017	Talca	307377.4	2017	7101	220357	67732753814	344	0.0015611	07101
7101121006	07101	1899	2017	Talca	307377.4	2017	7101	220357	67732753814	6821	0.0309543	07101
7101121007	07101	715	2017	Talca	307377.4	2017	7101	220357	67732753814	2524	0.0114541	07101
7101121008	07101	634	2017	Talca	307377.4	2017	7101	220357	67732753814	2216	0.0100564	07101
7101121009	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2511	0.0113951	07101
7101121010	07101	643	2017	Talca	307377.4	2017	7101	220357	67732753814	2406	0.0109186	07101
7101131001	07101	514	2017	Talca	307377.4	2017	7101	220357	67732753814	1494	0.0067799	07101
7101131002	07101	1408	2017	Talca	307377.4	2017	7101	220357	67732753814	4280	0.0194230	07101
7101131003	07101	1030	2017	Talca	307377.4	2017	7101	220357	67732753814	3547	0.0160966	07101
7101131004	07101	1087	2017	Talca	307377.4	2017	7101	220357	67732753814	3431	0.0155702	07101
7101131005	07101	923	2017	Talca	307377.4	2017	7101	220357	67732753814	2875	0.0130470	07101
7101131006	07101	607	2017	Talca	307377.4	2017	7101	220357	67732753814	2089	0.0094801	07101
7101131007	07101	859	2017	Talca	307377.4	2017	7101	220357	67732753814	2855	0.0129562	07101
7101131008	07101	1075	2017	Talca	307377.4	2017	7101	220357	67732753814	3485	0.0158152	07101
7101141001	07101	811	2017	Talca	307377.4	2017	7101	220357	67732753814	3226	0.0146399	07101
7101141002	07101	835	2017	Talca	307377.4	2017	7101	220357	67732753814	2620	0.0118898	07101
7101141003	07101	753	2017	Talca	307377.4	2017	7101	220357	67732753814	2174	0.0098658	07101
7101141004	07101	1284	2017	Talca	307377.4	2017	7101	220357	67732753814	5692	0.0258308	07101
7101141005	07101	1119	2017	Talca	307377.4	2017	7101	220357	67732753814	4770	0.0216467	07101
7101141006	07101	1514	2017	Talca	307377.4	2017	7101	220357	67732753814	4860	0.0220551	07101
7101141007	07101	1413	2017	Talca	307377.4	2017	7101	220357	67732753814	4908	0.0222729	07101
7101141008	07101	1505	2017	Talca	307377.4	2017	7101	220357	67732753814	4971	0.0225588	07101
7101141009	07101	1110	2017	Talca	307377.4	2017	7101	220357	67732753814	3856	0.0174989	07101
7101141010	07101	658	2017	Talca	307377.4	2017	7101	220357	67732753814	2182	0.0099021	07101
7101151001	07101	79	2017	Talca	307377.4	2017	7101	220357	67732753814	282	0.0012797	07101
7101161001	07101	341	2017	Talca	307377.4	2017	7101	220357	67732753814	1644	0.0074606	07101
7101991999	07101	306	2017	Talca	307377.4	2017	7101	220357	67732753814	883	0.0040071	07101
7102011001	07102	535	2017	Constitución	280736.9	2017	7102	46068	12932986800	1744	0.0378571	07102
7102011002	07102	658	2017	Constitución	280736.9	2017	7102	46068	12932986800	2069	0.0449119	07102
7102021001	07102	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800	4504	0.0977685	07102
7102021002	07102	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800	3767	0.0817704	07102
7102021003	07102	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800	5586	0.1212555	07102
7102021004	07102	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800	5245	0.1138534	07102
7102021005	07102	760	2017	Constitución	280736.9	2017	7102	46068	12932986800	2544	0.0552227	07102
7102031001	07102	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800	4021	0.0872840	07102
7102031002	07102	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800	4542	0.0985934	07102
7102041001	07102	615	2017	Constitución	280736.9	2017	7102	46068	12932986800	2083	0.0452158	07102
7102071001	07102	367	2017	Constitución	280736.9	2017	7102	46068	12932986800	1093	0.0237258	07102
7102991999	07102	23	2017	Constitución	280736.9	2017	7102	46068	12932986800	75	0.0016280	07102
7103011001	07103	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120	3368	0.3564776	07103
7103991999	07103	10	2017	Curepto	281855.5	2017	7103	9448	2662971120	21	0.0022227	07103
7104011001	07104	895	2017	Empedrado	209235.2	2017	7104	4142	866652110	2983	0.7201835	07104
7104021001	07104	9	2017	Empedrado	209235.2	2017	7104	4142	866652110	22	0.0053114	07104
7104991999	07104	2	2017	Empedrado	209235.2	2017	7104	4142	866652110	3	0.0007243	07104
7105011001	07105	1273	2017	Maule	245019.7	2017	7105	49721	12182624190	4782	0.0961767	07105
7105051001	07105	221	2017	Maule	245019.7	2017	7105	49721	12182624190	620	0.0124696	07105
7105051002	07105	987	2017	Maule	245019.7	2017	7105	49721	12182624190	3942	0.0792824	07105

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
7101011001	07101	417	2017	Talca	307377.4	2017	7101	220357	67732753814	1693	0.0076830	07101	520389877.4
7101011002	07101	687	2017	Talca	307377.4	2017	7101	220357	67732753814	1889	0.0085725	07101	580635840.7
7101021001	07101	408	2017	Talca	307377.4	2017	7101	220357	67732753814	951	0.0043157	07101	292315873.2
7101021002	07101	335	2017	Talca	307377.4	2017	7101	220357	67732753814	912	0.0041387	07101	280328156.0
7101021003	07101	262	2017	Talca	307377.4	2017	7101	220357	67732753814	857	0.0038891	07101	263422401.0
7101031001	07101	610	2017	Talca	307377.4	2017	7101	220357	67732753814	1465	0.0066483	07101	450307838.4
7101031002	07101	415	2017	Talca	307377.4	2017	7101	220357	67732753814	1428	0.0064804	07101	438934875.9
7101031003	07101	318	2017	Talca	307377.4	2017	7101	220357	67732753814	1037	0.0047060	07101	318750326.5
7101041001	07101	261	2017	Talca	307377.4	2017	7101	220357	67732753814	1306	0.0059267	07101	401434837.5
7101041002	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1682	0.0076331	07101	517008726.4
7101041003	07101	557	2017	Talca	307377.4	2017	7101	220357	67732753814	1627	0.0073835	07101	500102971.3
7101041004	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1684	0.0076421	07101	517623481.1
7101041005	07101	397	2017	Talca	307377.4	2017	7101	220357	67732753814	1196	0.0054276	07101	367623327.4
7101051001	07101	969	2017	Talca	307377.4	2017	7101	220357	67732753814	2870	0.0130243	07101	882173034.9
7101051002	07101	336	2017	Talca	307377.4	2017	7101	220357	67732753814	1049	0.0047605	07101	322438854.9
7101061001	07101	644	2017	Talca	307377.4	2017	7101	220357	67732753814	1756	0.0079689	07101	539754651.3
7101061002	07101	383	2017	Talca	307377.4	2017	7101	220357	67732753814	1371	0.0062217	07101	421414366.1
7101061003	07101	826	2017	Talca	307377.4	2017	7101	220357	67732753814	3032	0.0137595	07101	931968167.9
7101061004	07101	834	2017	Talca	307377.4	2017	7101	220357	67732753814	3215	0.0145900	07101	988218225.5
7101061005	07101	842	2017	Talca	307377.4	2017	7101	220357	67732753814	2884	0.0130879	07101	886476318.0
7101061006	07101	994	2017	Talca	307377.4	2017	7101	220357	67732753814	3367	0.0152798	07101	1034939584.8
7101071001	07101	676	2017	Talca	307377.4	2017	7101	220357	67732753814	2136	0.0096934	07101	656558049.7
7101071002	07101	1117	2017	Talca	307377.4	2017	7101	220357	67732753814	3210	0.0145673	07101	986681338.7
7101071003	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1773	0.0080460	07101	544980066.5
7101071004	07101	559	2017	Talca	307377.4	2017	7101	220357	67732753814	1769	0.0080279	07101	543750557.0
7101071005	07101	609	2017	Talca	307377.4	2017	7101	220357	67732753814	1956	0.0088765	07101	601230124.1
7101071006	07101	854	2017	Talca	307377.4	2017	7101	220357	67732753814	3097	0.0140545	07101	951947696.5
7101071007	07101	1028	2017	Talca	307377.4	2017	7101	220357	67732753814	3759	0.0170587	07101	1155431511.5
7101071008	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2586	0.0117355	07101	794877863.5
7101071009	07101	1185	2017	Talca	307377.4	2017	7101	220357	67732753814	3775	0.0171313	07101	1160349549.4
7101071010	07101	684	2017	Talca	307377.4	2017	7101	220357	67732753814	2262	0.0102652	07101	695287597.5
7101081001	07101	331	2017	Talca	307377.4	2017	7101	220357	67732753814	1170	0.0053096	07101	359631516.0
7101081002	07101	587	2017	Talca	307377.4	2017	7101	220357	67732753814	1841	0.0083546	07101	565881727.3
7101081003	07101	468	2017	Talca	307377.4	2017	7101	220357	67732753814	1281	0.0058133	07101	393750403.4
7101081004	07101	664	2017	Talca	307377.4	2017	7101	220357	67732753814	1562	0.0070885	07101	480123442.7
7101081005	07101	893	2017	Talca	307377.4	2017	7101	220357	67732753814	3349	0.0151981	07101	1029406792.3
7101081006	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1555	0.0070567	07101	477971801.1
7101091001	07101	40	2017	Talca	307377.4	2017	7101	220357	67732753814	179	0.0008123	07101	55020548.2
7101101001	07101	378	2017	Talca	307377.4	2017	7101	220357	67732753814	1347	0.0061128	07101	414037309.4
7101111001	07101	1458	2017	Talca	307377.4	2017	7101	220357	67732753814	5385	0.0244376	07101	1655227105.5
7101111002	07101	1044	2017	Talca	307377.4	2017	7101	220357	67732753814	3446	0.0156383	07101	1059222396.6
7101111003	07101	927	2017	Talca	307377.4	2017	7101	220357	67732753814	3435	0.0155883	07101	1055841245.6
7101111004	07101	102	2017	Talca	307377.4	2017	7101	220357	67732753814	347	0.0015747	07101	106659945.3
7101111005	07101	1674	2017	Talca	307377.4	2017	7101	220357	67732753814	5190	0.0235527	07101	1595288519.5
7101111006	07101	2449	2017	Talca	307377.4	2017	7101	220357	67732753814	7789	0.0353472	07101	2394162288.7
7101111007	07101	509	2017	Talca	307377.4	2017	7101	220357	67732753814	1912	0.0086768	07101	587705520.1
7101111008	07101	1501	2017	Talca	307377.4	2017	7101	220357	67732753814	5035	0.0228493	07101	1547645028.1
7101111009	07101	1179	2017	Talca	307377.4	2017	7101	220357	67732753814	3737	0.0169588	07101	1148669209.5
7101111010	07101	566	2017	Talca	307377.4	2017	7101	220357	67732753814	1695	0.0076921	07101	521004632.1
7101121001	07101	1182	2017	Talca	307377.4	2017	7101	220357	67732753814	4067	0.0184564	07101	1250103739.7
7101121002	07101	540	2017	Talca	307377.4	2017	7101	220357	67732753814	1855	0.0084182	07101	570185010.3
7101121003	07101	914	2017	Talca	307377.4	2017	7101	220357	67732753814	3773	0.0171222	07101	1159734794.6
7101121004	07101	1224	2017	Talca	307377.4	2017	7101	220357	67732753814	4426	0.0200856	07101	1360452213.4
7101121005	07101	97	2017	Talca	307377.4	2017	7101	220357	67732753814	344	0.0015611	07101	105737813.2
7101121006	07101	1899	2017	Talca	307377.4	2017	7101	220357	67732753814	6821	0.0309543	07101	2096621000.3
7101121007	07101	715	2017	Talca	307377.4	2017	7101	220357	67732753814	2524	0.0114541	07101	775820466.9
7101121008	07101	634	2017	Talca	307377.4	2017	7101	220357	67732753814	2216	0.0100564	07101	681148238.8
7101121009	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2511	0.0113951	07101	771824561.2
7101121010	07101	643	2017	Talca	307377.4	2017	7101	220357	67732753814	2406	0.0109186	07101	739549938.0
7101131001	07101	514	2017	Talca	307377.4	2017	7101	220357	67732753814	1494	0.0067799	07101	459221781.9
7101131002	07101	1408	2017	Talca	307377.4	2017	7101	220357	67732753814	4280	0.0194230	07101	1315575118.2
7101131003	07101	1030	2017	Talca	307377.4	2017	7101	220357	67732753814	3547	0.0160966	07101	1090267510.4
7101131004	07101	1087	2017	Talca	307377.4	2017	7101	220357	67732753814	3431	0.0155702	07101	1054611736.1
7101131005	07101	923	2017	Talca	307377.4	2017	7101	220357	67732753814	2875	0.0130470	07101	883709921.7
7101131006	07101	607	2017	Talca	307377.4	2017	7101	220357	67732753814	2089	0.0094801	07101	642111313.5
7101131007	07101	859	2017	Talca	307377.4	2017	7101	220357	67732753814	2855	0.0129562	07101	877562374.4
7101131008	07101	1075	2017	Talca	307377.4	2017	7101	220357	67732753814	3485	0.0158152	07101	1071210113.8
7101141001	07101	811	2017	Talca	307377.4	2017	7101	220357	67732753814	3226	0.0146399	07101	991599376.5
7101141002	07101	835	2017	Talca	307377.4	2017	7101	220357	67732753814	2620	0.0118898	07101	805328693.9
7101141003	07101	753	2017	Talca	307377.4	2017	7101	220357	67732753814	2174	0.0098658	07101	668238389.5
7101141004	07101	1284	2017	Talca	307377.4	2017	7101	220357	67732753814	5692	0.0258308	07101	1749591956.3
7101141005	07101	1119	2017	Talca	307377.4	2017	7101	220357	67732753814	4770	0.0216467	07101	1466190026.6
7101141006	07101	1514	2017	Talca	307377.4	2017	7101	220357	67732753814	4860	0.0220551	07101	1493853989.4
7101141007	07101	1413	2017	Talca	307377.4	2017	7101	220357	67732753814	4908	0.0222729	07101	1508608102.9
7101141008	07101	1505	2017	Talca	307377.4	2017	7101	220357	67732753814	4971	0.0225588	07101	1527972876.8
7101141009	07101	1110	2017	Talca	307377.4	2017	7101	220357	67732753814	3856	0.0174989	07101	1185247115.9
7101141010	07101	658	2017	Talca	307377.4	2017	7101	220357	67732753814	2182	0.0099021	07101	670697408.4
7101151001	07101	79	2017	Talca	307377.4	2017	7101	220357	67732753814	282	0.0012797	07101	86680416.7
7101161001	07101	341	2017	Talca	307377.4	2017	7101	220357	67732753814	1644	0.0074606	07101	505328386.5
7101991999	07101	306	2017	Talca	307377.4	2017	7101	220357	67732753814	883	0.0040071	07101	271414212.5
7102011001	07102	535	2017	Constitución	280736.9	2017	7102	46068	12932986800	1744	0.0378571	07102	489605126.8
7102011002	07102	658	2017	Constitución	280736.9	2017	7102	46068	12932986800	2069	0.0449119	07102	580844614.3
7102021001	07102	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800	4504	0.0977685	07102	1264438928.3
7102021002	07102	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800	3767	0.0817704	07102	1057535844.3
7102021003	07102	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800	5586	0.1212555	07102	1568196237.5
7102021004	07102	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800	5245	0.1138534	07102	1472464959.8
7102021005	07102	760	2017	Constitución	280736.9	2017	7102	46068	12932986800	2544	0.0552227	07102	714194634.5
7102031001	07102	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800	4021	0.0872840	07102	1128843013.0
7102031002	07102	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800	4542	0.0985934	07102	1275106929.9
7102041001	07102	615	2017	Constitución	280736.9	2017	7102	46068	12932986800	2083	0.0452158	07102	584774930.7
7102071001	07102	367	2017	Constitución	280736.9	2017	7102	46068	12932986800	1093	0.0237258	07102	306845414.9
7102991999	07102	23	2017	Constitución	280736.9	2017	7102	46068	12932986800	75	0.0016280	07102	21055266.3
7103011001	07103	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120	3368	0.3564776	07103	949289450.8
7103991999	07103	10	2017	Curepto	281855.5	2017	7103	9448	2662971120	21	0.0022227	07103	5918966.3
7104011001	07104	895	2017	Empedrado	209235.2	2017	7104	4142	866652110	2983	0.7201835	07104	624148537.7
7104021001	07104	9	2017	Empedrado	209235.2	2017	7104	4142	866652110	22	0.0053114	07104	4603173.9
7104991999	07104	2	2017	Empedrado	209235.2	2017	7104	4142	866652110	3	0.0007243	07104	627705.5
7105011001	07105	1273	2017	Maule	245019.7	2017	7105	49721	12182624190	4782	0.0961767	07105	1171684175.2
7105051001	07105	221	2017	Maule	245019.7	2017	7105	49721	12182624190	620	0.0124696	07105	151912210.1
7105051002	07105	987	2017	Maule	245019.7	2017	7105	49721	12182624190	3942	0.0792824	07105	965867632.5

7 Análisis de regresión

Aplicaremos un análisis de regresión donde:

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

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

7.1 Diagrama de dispersión loess

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

7.2 Outliers

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

7.3 Modelo lineal

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

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

## 
## Call:
## lm(formula = multi_pob ~ (Freq.x), data = h_y_m_comuna_corr_01)
## 
## Residuals:
##        Min         1Q     Median         3Q        Max 
## -259374630  -49482525    3181871   39577197  531969693 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -9231951   10156753  -0.909    0.364    
## Freq.x        955494      11908  80.238   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 99800000 on 313 degrees of freedom
## Multiple R-squared:  0.9536, Adjusted R-squared:  0.9535 
## F-statistic:  6438 on 1 and 313 DF,  p-value: < 2.2e-16

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

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

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

8 Modelos alternativos

### 8.1 Modelo cuadrático

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

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

kbl(modelos_bind) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

	modelo	dato	sintaxis	latex
3	logarítmico	0.607594150779086	linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = ''beta_0 + ''beta_1 ln X \]
6	log-raíz	0.832319620889849	linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = e^{''beta_0 + ''beta_1 ''sqrt{X}} \]
7	raíz-log	0.837319916586905	linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = {''beta_0}^2 + 2 ''beta_0 ''beta_1 ''ln{X}+ ''beta_1^2 ln^2X \]
4	raíz cuadrada	0.87488729407859	linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = ''beta_0 + ''beta_1 ''sqrt {X} \]
1	cuadrático	0.953489035956054	linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = ''beta_0 + ''beta_1 X^2 \]
2	cúbico	0.953489035956054	linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = ''beta_0 + ''beta_1 X^3 \]
5	raíz-raíz	0.973569200632626	linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = {''beta_0}^2 + 2 ''beta_0 ''beta_1 ''sqrt{X}+ ''beta_1^2 X \]
8	log-log	0.983903732898933	linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr_01)	\[ ''hat Y = e^{''beta_0+''beta_1 ln{X}} \]

h_y_m_comuna_corr <- h_y_m_comuna_corr_01
metodo <- 8


switch (metodo,
        case = linearMod <- lm( multi_pob~(Freq.x^2) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~(Freq.x^3) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~log(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( multi_pob~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( sqrt(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( log(multi_pob)~sqrt(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( sqrt(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr),
        case = linearMod <- lm( log(multi_pob)~log(Freq.x) , data=h_y_m_comuna_corr)
)
summary(linearMod)

## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.97827 -0.08984 -0.00657  0.10866  0.73246 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 13.452623   0.046887   286.9   <2e-16 ***
## log(Freq.x)  1.044784   0.007541   138.5   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1933 on 313 degrees of freedom
## Multiple R-squared:  0.984,  Adjusted R-squared:  0.9839 
## F-statistic: 1.919e+04 on 1 and 313 DF,  p-value: < 2.2e-16

aa <- linearMod$coefficients[1]
aa

## (Intercept) 
##    13.45262

bb <- linearMod$coefficients[2]
bb

## log(Freq.x) 
##    1.044784

9 Modelo log-log (log-log)

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

9.1 Diagrama de dispersión sobre log-log

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

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

9.2 Modelo log-log

Observemos nuevamente el resultado sobre log-log.

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

## 
## Call:
## lm(formula = log(multi_pob) ~ log(Freq.x), data = h_y_m_comuna_corr)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.97827 -0.08984 -0.00657  0.10866  0.73246 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 13.452623   0.046887   286.9   <2e-16 ***
## log(Freq.x)  1.044784   0.007541   138.5   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.1933 on 313 degrees of freedom
## Multiple R-squared:  0.984,  Adjusted R-squared:  0.9839 
## F-statistic: 1.919e+04 on 1 and 313 DF,  p-value: < 2.2e-16

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

9.3 Análisis de residuos

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

9.4 Ecuación del modelo

\[ \hat Y = e^{13.45262 +1.044784 \cdot ln{X}} \]

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

Esta nueva variable se llamará: est_ing

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

r3_100 <- h_y_m_comuna_corr[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

zona	código.x	Freq.x	anio	comuna.x	promedio_i	año	comuna.y	personas	Ingresos_expandidos	Freq.y	p_poblacional	código.y	multi_pob	est_ing
7101011001	07101	417	2017	Talca	307377.4	2017	7101	220357	67732753814	1693	0.0076830	07101	520389877.4	380081029
7101011002	07101	687	2017	Talca	307377.4	2017	7101	220357	67732753814	1889	0.0085725	07101	580635840.7	640334628
7101021001	07101	408	2017	Talca	307377.4	2017	7101	220357	67732753814	951	0.0043157	07101	292315873.2	371514640
7101021002	07101	335	2017	Talca	307377.4	2017	7101	220357	67732753814	912	0.0041387	07101	280328156.0	302361417
7101021003	07101	262	2017	Talca	307377.4	2017	7101	220357	67732753814	857	0.0038891	07101	263422401.0	233885041
7101031001	07101	610	2017	Talca	307377.4	2017	7101	220357	67732753814	1465	0.0066483	07101	450307838.4	565546115
7101031002	07101	415	2017	Talca	307377.4	2017	7101	220357	67732753814	1428	0.0064804	07101	438934875.9	378176665
7101031003	07101	318	2017	Talca	307377.4	2017	7101	220357	67732753814	1037	0.0047060	07101	318750326.5	286349067
7101041001	07101	261	2017	Talca	307377.4	2017	7101	220357	67732753814	1306	0.0059267	07101	401434837.5	232952452
7101041002	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1682	0.0076331	07101	517008726.4	465214416
7101041003	07101	557	2017	Talca	307377.4	2017	7101	220357	67732753814	1627	0.0073835	07101	500102971.3	514310686
7101041004	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1684	0.0076421	07101	517623481.1	465214416
7101041005	07101	397	2017	Talca	307377.4	2017	7101	220357	67732753814	1196	0.0054276	07101	367623327.4	361056114
7101051001	07101	969	2017	Talca	307377.4	2017	7101	220357	67732753814	2870	0.0130243	07101	882173034.9	917198409
7101051002	07101	336	2017	Talca	307377.4	2017	7101	220357	67732753814	1049	0.0047605	07101	322438854.9	303304472
7101061001	07101	644	2017	Talca	307377.4	2017	7101	220357	67732753814	1756	0.0079689	07101	539754651.3	598520446
7101061002	07101	383	2017	Talca	307377.4	2017	7101	220357	67732753814	1371	0.0062217	07101	421414366.1	347764070
7101061003	07101	826	2017	Talca	307377.4	2017	7101	220357	67732753814	3032	0.0137595	07101	931968167.9	776272262
7101061004	07101	834	2017	Talca	307377.4	2017	7101	220357	67732753814	3215	0.0145900	07101	988218225.5	784129039
7101061005	07101	842	2017	Talca	307377.4	2017	7101	220357	67732753814	2884	0.0130879	07101	886476318.0	791989193
7101061006	07101	994	2017	Talca	307377.4	2017	7101	220357	67732753814	3367	0.0152798	07101	1034939584.8	941935855
7101071001	07101	676	2017	Talca	307377.4	2017	7101	220357	67732753814	2136	0.0096934	07101	656558049.7	629626515
7101071002	07101	1117	2017	Talca	307377.4	2017	7101	220357	67732753814	3210	0.0145673	07101	986681338.7	1064038108
7101071003	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1773	0.0080460	07101	544980066.5	496958599
7101071004	07101	559	2017	Talca	307377.4	2017	7101	220357	67732753814	1769	0.0080279	07101	543750557.0	516240262
7101071005	07101	609	2017	Talca	307377.4	2017	7101	220357	67732753814	1956	0.0088765	07101	601230124.1	564577505
7101071006	07101	854	2017	Talca	307377.4	2017	7101	220357	67732753814	3097	0.0140545	07101	951947696.5	803785685
7101071007	07101	1028	2017	Talca	307377.4	2017	7101	220357	67732753814	3759	0.0170587	07101	1155431511.5	975623399
7101071008	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2586	0.0117355	07101	794877863.5	843160412
7101071009	07101	1185	2017	Talca	307377.4	2017	7101	220357	67732753814	3775	0.0171313	07101	1160349549.4	1131805372
7101071010	07101	684	2017	Talca	307377.4	2017	7101	220357	67732753814	2262	0.0102652	07101	695287597.5	637413467
7101081001	07101	331	2017	Talca	307377.4	2017	7101	220357	67732753814	1170	0.0053096	07101	359631516.0	298590460
7101081002	07101	587	2017	Talca	307377.4	2017	7101	220357	67732753814	1841	0.0083546	07101	565881727.3	543286315
7101081003	07101	468	2017	Talca	307377.4	2017	7101	220357	67732753814	1281	0.0058133	07101	393750403.4	428775648
7101081004	07101	664	2017	Talca	307377.4	2017	7101	220357	67732753814	1562	0.0070885	07101	480123442.7	617953838
7101081005	07101	893	2017	Talca	307377.4	2017	7101	220357	67732753814	3349	0.0151981	07101	1029406792.3	842175067
7101081006	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1555	0.0070567	07101	477971801.1	496958599
7101091001	07101	40	2017	Talca	307377.4	2017	7101	220357	67732753814	179	0.0008123	07101	55020548.2	32825136
7101101001	07101	378	2017	Talca	307377.4	2017	7101	220357	67732753814	1347	0.0061128	07101	414037309.4	343022142
7101111001	07101	1458	2017	Talca	307377.4	2017	7101	220357	67732753814	5385	0.0244376	07101	1655227105.5	1405540116
7101111002	07101	1044	2017	Talca	307377.4	2017	7101	220357	67732753814	3446	0.0156383	07101	1059222396.6	991493738
7101111003	07101	927	2017	Talca	307377.4	2017	7101	220357	67732753814	3435	0.0155883	07101	1055841245.6	875704181
7101111004	07101	102	2017	Talca	307377.4	2017	7101	220357	67732753814	347	0.0015747	07101	106659945.3	87287739
7101111005	07101	1674	2017	Talca	307377.4	2017	7101	220357	67732753814	5190	0.0235527	07101	1595288519.5	1623783510
7101111006	07101	2449	2017	Talca	307377.4	2017	7101	220357	67732753814	7789	0.0353472	07101	2394162288.7	2416358008
7101111007	07101	509	2017	Talca	307377.4	2017	7101	220357	67732753814	1912	0.0086768	07101	587705520.1	468096509
7101111008	07101	1501	2017	Talca	307377.4	2017	7101	220357	67732753814	5035	0.0228493	07101	1547645028.1	1448877708
7101111009	07101	1179	2017	Talca	307377.4	2017	7101	220357	67732753814	3737	0.0169588	07101	1148669209.5	1125818749
7101111010	07101	566	2017	Talca	307377.4	2017	7101	220357	67732753814	1695	0.0076921	07101	521004632.1	522996204
7101121001	07101	1182	2017	Talca	307377.4	2017	7101	220357	67732753814	4067	0.0184564	07101	1250103739.7	1128811890
7101121002	07101	540	2017	Talca	307377.4	2017	7101	220357	67732753814	1855	0.0084182	07101	570185010.3	497921931
7101121003	07101	914	2017	Talca	307377.4	2017	7101	220357	67732753814	3773	0.0171222	07101	1159734794.6	862877610
7101121004	07101	1224	2017	Talca	307377.4	2017	7101	220357	67732753814	4426	0.0200856	07101	1360452213.4	1170751222
7101121005	07101	97	2017	Talca	307377.4	2017	7101	220357	67732753814	344	0.0015611	07101	105737813.2	82822291
7101121006	07101	1899	2017	Talca	307377.4	2017	7101	220357	67732753814	6821	0.0309543	07101	2096621000.3	1852466827
7101121007	07101	715	2017	Talca	307377.4	2017	7101	220357	67732753814	2524	0.0114541	07101	775820466.9	667626037
7101121008	07101	634	2017	Talca	307377.4	2017	7101	220357	67732753814	2216	0.0100564	07101	681148238.8	588813829
7101121009	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2511	0.0113951	07101	771824561.2	843160412
7101121010	07101	643	2017	Talca	307377.4	2017	7101	220357	67732753814	2406	0.0109186	07101	739549938.0	597549478
7101131001	07101	514	2017	Talca	307377.4	2017	7101	220357	67732753814	1494	0.0067799	07101	459221781.9	472901686
7101131002	07101	1408	2017	Talca	307377.4	2017	7101	220357	67732753814	4280	0.0194230	07101	1315575118.2	1355219617
7101131003	07101	1030	2017	Talca	307377.4	2017	7101	220357	67732753814	3547	0.0160966	07101	1090267510.4	977606590
7101131004	07101	1087	2017	Talca	307377.4	2017	7101	220357	67732753814	3431	0.0155702	07101	1054611736.1	1034198832
7101131005	07101	923	2017	Talca	307377.4	2017	7101	220357	67732753814	2875	0.0130470	07101	883709921.7	871756680
7101131006	07101	607	2017	Talca	307377.4	2017	7101	220357	67732753814	2089	0.0094801	07101	642111313.5	562640500
7101131007	07101	859	2017	Talca	307377.4	2017	7101	220357	67732753814	2855	0.0129562	07101	877562374.4	808703088
7101131008	07101	1075	2017	Talca	307377.4	2017	7101	220357	67732753814	3485	0.0158152	07101	1071210113.8	1022273389
7101141001	07101	811	2017	Talca	307377.4	2017	7101	220357	67732753814	3226	0.0146399	07101	991599376.5	761550013
7101141002	07101	835	2017	Talca	307377.4	2017	7101	220357	67732753814	2620	0.0118898	07101	805328693.9	785111375
7101141003	07101	753	2017	Talca	307377.4	2017	7101	220357	67732753814	2174	0.0098658	07101	668238389.5	704740688
7101141004	07101	1284	2017	Talca	307377.4	2017	7101	220357	67732753814	5692	0.0258308	07101	1749591956.3	1230775947
7101141005	07101	1119	2017	Talca	307377.4	2017	7101	220357	67732753814	4770	0.0216467	07101	1466190026.6	1066028681
7101141006	07101	1514	2017	Talca	307377.4	2017	7101	220357	67732753814	4860	0.0220551	07101	1493853989.4	1461990794
7101141007	07101	1413	2017	Talca	307377.4	2017	7101	220357	67732753814	4908	0.0222729	07101	1508608102.9	1360248112
7101141008	07101	1505	2017	Talca	307377.4	2017	7101	220357	67732753814	4971	0.0225588	07101	1527972876.8	1452911964
7101141009	07101	1110	2017	Talca	307377.4	2017	7101	220357	67732753814	3856	0.0174989	07101	1185247115.9	1057072364
7101141010	07101	658	2017	Talca	307377.4	2017	7101	220357	67732753814	2182	0.0099021	07101	670697408.4	612121031
7101151001	07101	79	2017	Talca	307377.4	2017	7101	220357	67732753814	282	0.0012797	07101	86680416.7	66835982
7101161001	07101	341	2017	Talca	307377.4	2017	7101	220357	67732753814	1644	0.0074606	07101	505328386.5	308021626
7101991999	07101	306	2017	Talca	307377.4	2017	7101	220357	67732753814	883	0.0040071	07101	271414212.5	275069178
7102011001	07102	535	2017	Constitución	280736.9	2017	7102	46068	12932986800	1744	0.0378571	07102	489605126.8	493106072
7102011002	07102	658	2017	Constitución	280736.9	2017	7102	46068	12932986800	2069	0.0449119	07102	580844614.3	612121031
7102021001	07102	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800	4504	0.0977685	07102	1264438928.3	1443835565
7102021002	07102	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800	3767	0.0817704	07102	1057535844.3	1070010303
7102021003	07102	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800	5586	0.1212555	07102	1568196237.5	1600481530
7102021004	07102	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800	5245	0.1138534	07102	1472464959.8	1536729378
7102021005	07102	760	2017	Constitución	280736.9	2017	7102	46068	12932986800	2544	0.0552227	07102	714194634.5	711586880
7102031001	07102	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800	4021	0.0872840	07102	1128843013.0	1303975122
7102031002	07102	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800	4542	0.0985934	07102	1275106929.9	1368295362
7102041001	07102	615	2017	Constitución	280736.9	2017	7102	46068	12932986800	2083	0.0452158	07102	584774930.7	570390227
7102071001	07102	367	2017	Constitución	280736.9	2017	7102	46068	12932986800	1093	0.0237258	07102	306845414.9	332599835
7102991999	07102	23	2017	Constitución	280736.9	2017	7102	46068	12932986800	75	0.0016280	07102	21055266.3	18412439
7103011001	07103	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120	3368	0.3564776	07103	949289450.8	1060057443
7103991999	07103	10	2017	Curepto	281855.5	2017	7103	9448	2662971120	21	0.0022227	07103	5918966.3	7712299
7104011001	07104	895	2017	Empedrado	209235.2	2017	7104	4142	866652110	2983	0.7201835	07104	624148537.7	844145806
7104021001	07104	9	2017	Empedrado	209235.2	2017	7104	4142	866652110	22	0.0053114	07104	4603173.9	6908395
7104991999	07104	2	2017	Empedrado	209235.2	2017	7104	4142	866652110	3	0.0007243	07104	627705.5	1435196
7105011001	07105	1273	2017	Maule	245019.7	2017	7105	49721	12182624190	4782	0.0961767	07105	1171684175.2	1219761831
7105051001	07105	221	2017	Maule	245019.7	2017	7105	49721	12182624190	620	0.0124696	07105	151912210.1	195786833
7105051002	07105	987	2017	Maule	245019.7	2017	7105	49721	12182624190	3942	0.0792824	07105	965867632.5	935006530

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

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

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

r3_100 <- h_y_m_comuna_corr[c(1:100),]
kbl(r3_100) %>%
  kable_styling(bootstrap_options = c("striped", "hover")) %>%
  kable_paper() %>%
  scroll_box(width = "100%", height = "300px")

zona	código.x	Freq.x	anio	comuna.x	promedio_i	año	comuna.y	personas	Ingresos_expandidos	Freq.y	p_poblacional	código.y	multi_pob	est_ing	ing_medio_zona
7101011001	07101	417	2017	Talca	307377.4	2017	7101	220357	67732753814	1693	0.0076830	07101	520389877.4	380081029	224501.5
7101011002	07101	687	2017	Talca	307377.4	2017	7101	220357	67732753814	1889	0.0085725	07101	580635840.7	640334628	338980.7
7101021001	07101	408	2017	Talca	307377.4	2017	7101	220357	67732753814	951	0.0043157	07101	292315873.2	371514640	390656.8
7101021002	07101	335	2017	Talca	307377.4	2017	7101	220357	67732753814	912	0.0041387	07101	280328156.0	302361417	331536.6
7101021003	07101	262	2017	Talca	307377.4	2017	7101	220357	67732753814	857	0.0038891	07101	263422401.0	233885041	272911.4
7101031001	07101	610	2017	Talca	307377.4	2017	7101	220357	67732753814	1465	0.0066483	07101	450307838.4	565546115	386038.3
7101031002	07101	415	2017	Talca	307377.4	2017	7101	220357	67732753814	1428	0.0064804	07101	438934875.9	378176665	264829.6
7101031003	07101	318	2017	Talca	307377.4	2017	7101	220357	67732753814	1037	0.0047060	07101	318750326.5	286349067	276132.2
7101041001	07101	261	2017	Talca	307377.4	2017	7101	220357	67732753814	1306	0.0059267	07101	401434837.5	232952452	178370.9
7101041002	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1682	0.0076331	07101	517008726.4	465214416	276584.1
7101041003	07101	557	2017	Talca	307377.4	2017	7101	220357	67732753814	1627	0.0073835	07101	500102971.3	514310686	316109.8
7101041004	07101	506	2017	Talca	307377.4	2017	7101	220357	67732753814	1684	0.0076421	07101	517623481.1	465214416	276255.6
7101041005	07101	397	2017	Talca	307377.4	2017	7101	220357	67732753814	1196	0.0054276	07101	367623327.4	361056114	301886.4
7101051001	07101	969	2017	Talca	307377.4	2017	7101	220357	67732753814	2870	0.0130243	07101	882173034.9	917198409	319581.3
7101051002	07101	336	2017	Talca	307377.4	2017	7101	220357	67732753814	1049	0.0047605	07101	322438854.9	303304472	289136.8
7101061001	07101	644	2017	Talca	307377.4	2017	7101	220357	67732753814	1756	0.0079689	07101	539754651.3	598520446	340843.1
7101061002	07101	383	2017	Talca	307377.4	2017	7101	220357	67732753814	1371	0.0062217	07101	421414366.1	347764070	253657.2
7101061003	07101	826	2017	Talca	307377.4	2017	7101	220357	67732753814	3032	0.0137595	07101	931968167.9	776272262	256026.5
7101061004	07101	834	2017	Talca	307377.4	2017	7101	220357	67732753814	3215	0.0145900	07101	988218225.5	784129039	243897.1
7101061005	07101	842	2017	Talca	307377.4	2017	7101	220357	67732753814	2884	0.0130879	07101	886476318.0	791989193	274614.8
7101061006	07101	994	2017	Talca	307377.4	2017	7101	220357	67732753814	3367	0.0152798	07101	1034939584.8	941935855	279755.2
7101071001	07101	676	2017	Talca	307377.4	2017	7101	220357	67732753814	2136	0.0096934	07101	656558049.7	629626515	294769.0
7101071002	07101	1117	2017	Talca	307377.4	2017	7101	220357	67732753814	3210	0.0145673	07101	986681338.7	1064038108	331476.0
7101071003	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1773	0.0080460	07101	544980066.5	496958599	280292.5
7101071004	07101	559	2017	Talca	307377.4	2017	7101	220357	67732753814	1769	0.0080279	07101	543750557.0	516240262	291826.0
7101071005	07101	609	2017	Talca	307377.4	2017	7101	220357	67732753814	1956	0.0088765	07101	601230124.1	564577505	288638.8
7101071006	07101	854	2017	Talca	307377.4	2017	7101	220357	67732753814	3097	0.0140545	07101	951947696.5	803785685	259536.9
7101071007	07101	1028	2017	Talca	307377.4	2017	7101	220357	67732753814	3759	0.0170587	07101	1155431511.5	975623399	259543.3
7101071008	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2586	0.0117355	07101	794877863.5	843160412	326048.1
7101071009	07101	1185	2017	Talca	307377.4	2017	7101	220357	67732753814	3775	0.0171313	07101	1160349549.4	1131805372	299816.0
7101071010	07101	684	2017	Talca	307377.4	2017	7101	220357	67732753814	2262	0.0102652	07101	695287597.5	637413467	281792.0
7101081001	07101	331	2017	Talca	307377.4	2017	7101	220357	67732753814	1170	0.0053096	07101	359631516.0	298590460	255205.5
7101081002	07101	587	2017	Talca	307377.4	2017	7101	220357	67732753814	1841	0.0083546	07101	565881727.3	543286315	295103.9
7101081003	07101	468	2017	Talca	307377.4	2017	7101	220357	67732753814	1281	0.0058133	07101	393750403.4	428775648	334719.5
7101081004	07101	664	2017	Talca	307377.4	2017	7101	220357	67732753814	1562	0.0070885	07101	480123442.7	617953838	395617.1
7101081005	07101	893	2017	Talca	307377.4	2017	7101	220357	67732753814	3349	0.0151981	07101	1029406792.3	842175067	251470.6
7101081006	07101	539	2017	Talca	307377.4	2017	7101	220357	67732753814	1555	0.0070567	07101	477971801.1	496958599	319587.5
7101091001	07101	40	2017	Talca	307377.4	2017	7101	220357	67732753814	179	0.0008123	07101	55020548.2	32825136	183380.6
7101101001	07101	378	2017	Talca	307377.4	2017	7101	220357	67732753814	1347	0.0061128	07101	414037309.4	343022142	254656.4
7101111001	07101	1458	2017	Talca	307377.4	2017	7101	220357	67732753814	5385	0.0244376	07101	1655227105.5	1405540116	261010.2
7101111002	07101	1044	2017	Talca	307377.4	2017	7101	220357	67732753814	3446	0.0156383	07101	1059222396.6	991493738	287723.1
7101111003	07101	927	2017	Talca	307377.4	2017	7101	220357	67732753814	3435	0.0155883	07101	1055841245.6	875704181	254935.7
7101111004	07101	102	2017	Talca	307377.4	2017	7101	220357	67732753814	347	0.0015747	07101	106659945.3	87287739	251549.7
7101111005	07101	1674	2017	Talca	307377.4	2017	7101	220357	67732753814	5190	0.0235527	07101	1595288519.5	1623783510	312867.7
7101111006	07101	2449	2017	Talca	307377.4	2017	7101	220357	67732753814	7789	0.0353472	07101	2394162288.7	2416358008	310227.0
7101111007	07101	509	2017	Talca	307377.4	2017	7101	220357	67732753814	1912	0.0086768	07101	587705520.1	468096509	244820.4
7101111008	07101	1501	2017	Talca	307377.4	2017	7101	220357	67732753814	5035	0.0228493	07101	1547645028.1	1448877708	287761.2
7101111009	07101	1179	2017	Talca	307377.4	2017	7101	220357	67732753814	3737	0.0169588	07101	1148669209.5	1125818749	301262.7
7101111010	07101	566	2017	Talca	307377.4	2017	7101	220357	67732753814	1695	0.0076921	07101	521004632.1	522996204	308552.3
7101121001	07101	1182	2017	Talca	307377.4	2017	7101	220357	67732753814	4067	0.0184564	07101	1250103739.7	1128811890	277553.9
7101121002	07101	540	2017	Talca	307377.4	2017	7101	220357	67732753814	1855	0.0084182	07101	570185010.3	497921931	268421.5
7101121003	07101	914	2017	Talca	307377.4	2017	7101	220357	67732753814	3773	0.0171222	07101	1159734794.6	862877610	228698.0
7101121004	07101	1224	2017	Talca	307377.4	2017	7101	220357	67732753814	4426	0.0200856	07101	1360452213.4	1170751222	264516.8
7101121005	07101	97	2017	Talca	307377.4	2017	7101	220357	67732753814	344	0.0015611	07101	105737813.2	82822291	240762.5
7101121006	07101	1899	2017	Talca	307377.4	2017	7101	220357	67732753814	6821	0.0309543	07101	2096621000.3	1852466827	271582.9
7101121007	07101	715	2017	Talca	307377.4	2017	7101	220357	67732753814	2524	0.0114541	07101	775820466.9	667626037	264511.1
7101121008	07101	634	2017	Talca	307377.4	2017	7101	220357	67732753814	2216	0.0100564	07101	681148238.8	588813829	265710.2
7101121009	07101	894	2017	Talca	307377.4	2017	7101	220357	67732753814	2511	0.0113951	07101	771824561.2	843160412	335786.7
7101121010	07101	643	2017	Talca	307377.4	2017	7101	220357	67732753814	2406	0.0109186	07101	739549938.0	597549478	248358.1
7101131001	07101	514	2017	Talca	307377.4	2017	7101	220357	67732753814	1494	0.0067799	07101	459221781.9	472901686	316533.9
7101131002	07101	1408	2017	Talca	307377.4	2017	7101	220357	67732753814	4280	0.0194230	07101	1315575118.2	1355219617	316640.1
7101131003	07101	1030	2017	Talca	307377.4	2017	7101	220357	67732753814	3547	0.0160966	07101	1090267510.4	977606590	275615.1
7101131004	07101	1087	2017	Talca	307377.4	2017	7101	220357	67732753814	3431	0.0155702	07101	1054611736.1	1034198832	301427.8
7101131005	07101	923	2017	Talca	307377.4	2017	7101	220357	67732753814	2875	0.0130470	07101	883709921.7	871756680	303219.7
7101131006	07101	607	2017	Talca	307377.4	2017	7101	220357	67732753814	2089	0.0094801	07101	642111313.5	562640500	269334.8
7101131007	07101	859	2017	Talca	307377.4	2017	7101	220357	67732753814	2855	0.0129562	07101	877562374.4	808703088	283258.5
7101131008	07101	1075	2017	Talca	307377.4	2017	7101	220357	67732753814	3485	0.0158152	07101	1071210113.8	1022273389	293335.3
7101141001	07101	811	2017	Talca	307377.4	2017	7101	220357	67732753814	3226	0.0146399	07101	991599376.5	761550013	236066.3
7101141002	07101	835	2017	Talca	307377.4	2017	7101	220357	67732753814	2620	0.0118898	07101	805328693.9	785111375	299660.8
7101141003	07101	753	2017	Talca	307377.4	2017	7101	220357	67732753814	2174	0.0098658	07101	668238389.5	704740688	324167.7
7101141004	07101	1284	2017	Talca	307377.4	2017	7101	220357	67732753814	5692	0.0258308	07101	1749591956.3	1230775947	216229.1
7101141005	07101	1119	2017	Talca	307377.4	2017	7101	220357	67732753814	4770	0.0216467	07101	1466190026.6	1066028681	223486.1
7101141006	07101	1514	2017	Talca	307377.4	2017	7101	220357	67732753814	4860	0.0220551	07101	1493853989.4	1461990794	300821.2
7101141007	07101	1413	2017	Talca	307377.4	2017	7101	220357	67732753814	4908	0.0222729	07101	1508608102.9	1360248112	277149.2
7101141008	07101	1505	2017	Talca	307377.4	2017	7101	220357	67732753814	4971	0.0225588	07101	1527972876.8	1452911964	292277.6
7101141009	07101	1110	2017	Talca	307377.4	2017	7101	220357	67732753814	3856	0.0174989	07101	1185247115.9	1057072364	274137.0
7101141010	07101	658	2017	Talca	307377.4	2017	7101	220357	67732753814	2182	0.0099021	07101	670697408.4	612121031	280532.1
7101151001	07101	79	2017	Talca	307377.4	2017	7101	220357	67732753814	282	0.0012797	07101	86680416.7	66835982	237007.0
7101161001	07101	341	2017	Talca	307377.4	2017	7101	220357	67732753814	1644	0.0074606	07101	505328386.5	308021626	187361.1
7101991999	07101	306	2017	Talca	307377.4	2017	7101	220357	67732753814	883	0.0040071	07101	271414212.5	275069178	311516.6
7102011001	07102	535	2017	Constitución	280736.9	2017	7102	46068	12932986800	1744	0.0378571	07102	489605126.8	493106072	282744.3
7102011002	07102	658	2017	Constitución	280736.9	2017	7102	46068	12932986800	2069	0.0449119	07102	580844614.3	612121031	295853.6
7102021001	07102	1496	2017	Constitución	280736.9	2017	7102	46068	12932986800	4504	0.0977685	07102	1264438928.3	1443835565	320567.4
7102021002	07102	1123	2017	Constitución	280736.9	2017	7102	46068	12932986800	3767	0.0817704	07102	1057535844.3	1070010303	284048.4
7102021003	07102	1651	2017	Constitución	280736.9	2017	7102	46068	12932986800	5586	0.1212555	07102	1568196237.5	1600481530	286516.6
7102021004	07102	1588	2017	Constitución	280736.9	2017	7102	46068	12932986800	5245	0.1138534	07102	1472464959.8	1536729378	292989.4
7102021005	07102	760	2017	Constitución	280736.9	2017	7102	46068	12932986800	2544	0.0552227	07102	714194634.5	711586880	279711.8
7102031001	07102	1357	2017	Constitución	280736.9	2017	7102	46068	12932986800	4021	0.0872840	07102	1128843013.0	1303975122	324291.3
7102031002	07102	1421	2017	Constitución	280736.9	2017	7102	46068	12932986800	4542	0.0985934	07102	1275106929.9	1368295362	301253.9
7102041001	07102	615	2017	Constitución	280736.9	2017	7102	46068	12932986800	2083	0.0452158	07102	584774930.7	570390227	273831.1
7102071001	07102	367	2017	Constitución	280736.9	2017	7102	46068	12932986800	1093	0.0237258	07102	306845414.9	332599835	304299.9
7102991999	07102	23	2017	Constitución	280736.9	2017	7102	46068	12932986800	75	0.0016280	07102	21055266.3	18412439	245499.2
7103011001	07103	1113	2017	Curepto	281855.5	2017	7103	9448	2662971120	3368	0.3564776	07103	949289450.8	1060057443	314743.9
7103991999	07103	10	2017	Curepto	281855.5	2017	7103	9448	2662971120	21	0.0022227	07103	5918966.3	7712299	367252.3
7104011001	07104	895	2017	Empedrado	209235.2	2017	7104	4142	866652110	2983	0.7201835	07104	624148537.7	844145806	282985.5
7104021001	07104	9	2017	Empedrado	209235.2	2017	7104	4142	866652110	22	0.0053114	07104	4603173.9	6908395	314017.9
7104991999	07104	2	2017	Empedrado	209235.2	2017	7104	4142	866652110	3	0.0007243	07104	627705.5	1435196	478398.6
7105011001	07105	1273	2017	Maule	245019.7	2017	7105	49721	12182624190	4782	0.0961767	07105	1171684175.2	1219761831	255073.6
7105051001	07105	221	2017	Maule	245019.7	2017	7105	49721	12182624190	620	0.0124696	07105	151912210.1	195786833	315785.2
7105051002	07105	987	2017	Maule	245019.7	2017	7105	49721	12182624190	3942	0.0792824	07105	965867632.5	935006530	237190.9

Guardamos:

saveRDS(h_y_m_comuna_corr, "P03C/region_07_P03C_u.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