ENCOVI - Principales índices e indicadores de desigualdad

Data Sources: ENCOVI

Attaching package: <U+393C><U+3E31>dplyr<U+393C><U+3E32>

The following objects are masked from <U+393C><U+3E31>package:stats<U+393C><U+3E32>:

    filter, lag

The following objects are masked from <U+393C><U+3E31>package:base<U+393C><U+3E32>:

    intersect, setdiff, setequal, union

indicadores_desigualdad<- read_excel("indicadores_desigualdad.xlsx") 

names(indicadores_desigualdad) <- make.names(names(indicadores_desigualdad))

indicadores_desigualdad<- indicadores_desigualdad[complete.cases(indicadores_desigualdad), ]

str(indicadores_desigualdad)

Classes ‘tbl_df’, ‘tbl’ and 'data.frame':   22 obs. of  4 variables:
 $ Departamento    : chr  "Guatemala" "El Progreso" "Sacatepéquez" "Chimaltenango" ...
 $ Gini            : num  0.469 0.421 0.479 0.505 0.424 ...
 $ Atkinson..e...1.: num  0.325 0.267 0.337 0.364 0.266 ...
 $ Theil           : num  0.474 0.299 0.711 0.465 0.344 ...

head(indicadores_desigualdad)

summary(indicadores_desigualdad)

 Departamento            Gini        Atkinson..e...1.     Theil       
 Length:22          Min.   :0.3951   Min.   :0.2384   Min.   :0.2777  
 Class :character   1st Qu.:0.4538   1st Qu.:0.3153   1st Qu.:0.3760  
 Mode  :character   Median :0.4805   Median :0.3419   Median :0.4674  
                    Mean   :0.4862   Mean   :0.3535   Mean   :0.4924  
                    3rd Qu.:0.5069   3rd Qu.:0.3812   3rd Qu.:0.5318  
                    Max.   :0.6124   Max.   :0.5225   Max.   :0.9684  

indicadores_desigualdad <- indicadores_desigualdad %>% gather(indicadores.desigualdad, Penetracion, -Departamento)

indicadores_desigualdad$indicadores.desigualdad <- factor(indicadores_desigualdad$indicadores.desigualdad)
indicadores_desigualdad$Penetracion <- as.numeric(indicadores_desigualdad$Penetracion) / 100

str(indicadores_desigualdad)

Classes ‘tbl_df’, ‘tbl’ and 'data.frame':   66 obs. of  3 variables:
 $ Departamento           : chr  "Guatemala" "El Progreso" "Sacatepéquez" "Chimaltenango" ...
 $ indicadores.desigualdad: Factor w/ 3 levels "Atkinson..e...1.",..: 2 2 2 2 2 2 2 2 2 2 ...
 $ Penetracion            : num  0.00469 0.00421 0.00479 0.00505 0.00424 ...

head(indicadores_desigualdad)

summary(indicadores_desigualdad)

 Departamento           indicadores.desigualdad  Penetracion      
 Length:66          Atkinson..e...1.:22         Min.   :0.002384  
 Class :character   Gini            :22         1st Qu.:0.003575  
 Mode  :character   Theil           :22         Median :0.004406  
                                                Mean   :0.004440  
                                                3rd Qu.:0.004940  
                                                Max.   :0.009684  

ggplot(indicadores_desigualdad, aes(indicadores.desigualdad, Penetracion)) +
    geom_point() + theme(axis.text.x = element_text(angle = 90))

ggplot(indicadores_desigualdad, aes(indicadores.desigualdad, Penetracion)) +
    geom_boxplot() + theme(axis.text.x = element_text(angle = 90))

indicadores_desigualdad %>% filter(indicadores.desigualdad== "Gini") %>% arrange(desc(Penetracion)) %>% ggplot(aes(Departamento, Penetracion)) +
    geom_col(fill = "orange", color = "black") +
    labs(title = "Penetracion Gini") +
    theme(axis.text.x = element_text(angle = 90))

indicadores_desigualdad %>% ggplot(aes(Departamento, Penetracion)) +
    geom_col(aes(fill = indicadores.desigualdad, group = indicadores.desigualdad), position = "stack") +
    theme(axis.text.x = element_text(angle = 90))

indicadores_desigualdad %>% ggplot(aes(Departamento, Penetracion)) +
    geom_line(aes(color = indicadores.desigualdad, group = indicadores.desigualdad)) +
    theme(axis.text.x = element_text(angle = 90))

sc_pread <- spread(indicadores_desigualdad, indicadores.desigualdad, Penetracion)

sc_cor <- cor(select(sc_pread, -Departamento))