Abrir las bases de datos.

load("datos.rdata")
load("segregacion2017.rdata")

Unir bases de datos.

EC2016 <- select(EC2016, COMUNA_15R, ev.policial)
BBDD <- merge(EC2016, indice ,by="COMUNA_15R")

Colapsar a nivel comunal.

BBDD <- BBDD %>% 
  group_by(COMUNA_15R) %>% 
  summarise(segregacion = mean(IDC),
            acumulacion = mean(IDS),
            media.ev.policial = mean(ev.policial))

Correlación y diagrama de dispersión.

Segregación

cor(BBDD$segregacion, BBDD$media.ev.policial, method = "spearman") # 0.261
## [1] 0.26064
grafico <- ggplot(BBDD, aes(segregacion, media.ev.policial)) + 
  geom_point() +
  geom_point(aes(colour = cut_number(media.ev.policial,4))) + 
  theme_classic() +
  labs(title = "Segregación urbana y evaluación de las policías - nivel comunal",
       caption = "Fuente: Elaboración propia en base a Censo 2017 y Encuestas comunales de seguridad ciudadana 2016",
       x = "Índice Duncan de disimilitud",
       y = "Media de evaluación de las policías",
       colour = "Cuartiles - ev. policías")
grafico

Concentración desventajas sociales

cor(BBDD$acumulacion, BBDD$media.ev.policial, method = "spearman") # -0.241
## [1] -0.2406385
grafico2 <- ggplot(BBDD, aes(acumulacion, media.ev.policial)) + 
  geom_point() +
  geom_point(aes(colour = cut_number(media.ev.policial,4))) + 
  theme_classic() +
  labs(title = "Acumulación de desventajas sociales y evaluación de las policías - nivel comunal",
       caption = "Fuente: Elaboración propia en base a Censo 2017 y Encuestas comunales de seguridad ciudadana 2016",
       x = "Proporción de jefes de hogar sin escolaridad obligatoria",
       y = "Media de evaluación de las policías",
       colour = "Cuartiles - ev. policías")
grafico2

cor(BBDD$acumulacion, BBDD$segregacion, method = "spearman") #-0.411
## [1] -0.410719
grafico3 <- ggplot(BBDD, aes(acumulacion, segregacion)) + 
  geom_point() +
  geom_point(aes(colour = cut_number(media.ev.policial,4))) + 
  theme_classic() +
  labs(title = "Acumulación de desventajas sociales y segregación urbana - nivel comunal",
       caption = "Fuente: Elaboración propia en base a Censo 2017 y Encuestas comunales de seguridad ciudadana 2016",
       x = "Proporción de jefes de hogar sin escolaridad obligatoria",
       y = "Índice Duncan de disimilitud",
       colour = "Cuartiles - ev. policías")
grafico3