Tablas de contingencia sobre ingresos de Casen 2006-2020

por comunas y regiones

VE-CC-AJ

DataIntelligence
date: 31-08-2021

1 Las variables de ingreso a utilizar serán:

Tipo
Ingreso total Ingreso autónomo Ingreso del trabajo Ingreso de la ocupación principal
año
2017 ytotcor ytrabajocor ytrabajocor yoprcor
2015 ytotcor ytrabajocor ytrabajocor yoprcor
2013 ytotcor ytrabajocor ytrabajocor yoprcor
2011 ytrabaj yautaj ytrabaj yopraj
2009 ytrabaj yautaj ytrabaj yopraj
2006 ytrabaj yautaj ytrabaj yopraj


1.1 Lectura de bases de datos Casen

dataset_06 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2006_c.rds")
dataset_09 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2009_c.rds")
dataset_11 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2011_c.rds")
dataset_13 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2013_c.rds")
dataset_15 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2015_c.rds")
dataset_17 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2017_c.rds")
dataset_20 <<- readRDS("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2020.rds")

2 Generación de tablas de contingencia

for (var1 in 1:5) {
funcion1 <- function(n){

  comunales<-switch(n,"codigos_comunales_2006.rds","codigos_comunales_2009.rds","codigos_comunales_2011-2017.rds","codigos_comunales_2011-2017.rds","codigos_comunales_2011-2017.rds","codigos_comunales_2011-2017.rds","codigos_comunales_2011-2017.rds")
 xx<-switch(n,"2006","2009","2011","2013","2015","2017","2020")
 tanio <<- xx

 switch (var1,
        case = v1 <- switch(n,"S1","S1","s17","s14","s12","s12","s13"),
        case = v1 <- switch(n,"O1","O1","o1","o1","o1","o1","o1"),
        case = v1 <- switch(n,"O2","O2","o2","o2","o2","o2","o2"),
        case = v1 <- switch(n,"O3","O3","o3","o3","o3","o3","o3"),
        case = v1 <- switch(n,"O4","O4","o6","o6","o6","o6","o6")
)


if(xx==2006) {
eliminated <- dataset_06[!is.na(dataset_06$YTOTAJ),]
a <- eliminated$YTOTAJ
b <- eliminated$COMUNA
c <- eliminated[,c(v1)]
d <- eliminated$E1 #alfabetismo
e <- eliminated$T4 #etnia
f <- eliminated$SEXO
anio <- 2006
}

if(xx==2009) {
eliminated <- dataset_09[!is.na(dataset_09$YTOTAJ),]
a <- eliminated$YTOTAJ
b <- eliminated$COMUNA
c <- eliminated[,c(v1)]
d <- eliminated$E1 #alfabetismo
e <- eliminated$T5 #etnia
f <- eliminated$SEXO
anio <- 2009

}

if(xx==2011) {
eliminated <- dataset_11[!is.na(dataset_11$ytotaj),]
a <- eliminated$ytotaj
b <- eliminated$comuna
c <- eliminated[,c(v1)]
d <- eliminated$e1 #alfabetismo
e <- eliminated$r6 #etnia
f <- eliminated$sexo
anio <- 2011

}

if(xx==2013) {
eliminated <- dataset_13[!is.na(dataset_13$ytotcor),]
a <- eliminated$ytotcor
b <- eliminated$comuna
c <- eliminated[,c(v1)]
d <- eliminated$e1 #alfabetismo
e <- eliminated$r6 #etnia
f <- eliminated$sexo
anio <- 2013

}

if(xx==2015) {
eliminated <- dataset_15[!is.na(dataset_15$ytotcor),]
a <- eliminated$ytotcor
b <- eliminated$comuna
c <- eliminated[,c(v1)]
d <- eliminated$e1 #alfabetismo
e <- eliminated$r3 #etnia
f <- eliminated$sexo
anio <- 2015
}

if(xx==2017) {
eliminated <- dataset_17[!is.na(dataset_17$ytotcor),]
a <- eliminated$ytotcor
b <- eliminated$comuna
c <- eliminated[,c(v1)]
d <- eliminated$e1 #alfabetismo
e <- eliminated$r3 #etnia
f <- eliminated$sexo
anio <- 2017
}

if(xx==2020) {
eliminated <- dataset_20[!is.na(dataset_20$ytotcor),]
a <- eliminated$ytotcor
b <- eliminated$comuna
c <- eliminated[,c(v1)]
d <- eliminated$sexo #alfabetismo
e <- eliminated$r3 #etnia
f <- eliminated$sexo
anio <- 2020
}

promedios_grupales <-aggregate(a, by=list(b, c, d, e,f), FUN = mean , na.rm = TRUE)
promedios_grupales_sd <-aggregate(a, by=list(b, c, d, e,f), FUN = sd , na.rm = TRUE)
promedios_grupales_gini <-aggregate(a, by=list(b, c, d, e,f), FUN = gini , na.rm = TRUE)

promedios_grupales$sd = promedios_grupales_sd$sd.a
promedios_grupales$gini = promedios_grupales_gini$gini.a
d<-promedios_grupales
d$Año <- anio

names(d)[1] <- "Comuna"
names(d)[2] <- "Variable de estudio"
names(d)[3] <- "Alfabetismo"
names(d)[4] <- "Etnia"
names(d)[5] <- "Sexo"
names(d)[6] <- "Promedio del Ingreso"
names(d)[7] <- "Desviación standard"
names(d)[8] <- "Gini"
names(d)[9] <- "Año"

direc_cod_com <- paste0("C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/",comunales)
codigos_comunales <- readRDS(file = direc_cod_com)
names(codigos_comunales)[1] <- "Código"
names(codigos_comunales)[2] <- "Comuna"

tabla_df = merge( x = d, y = codigos_comunales, by = "Comuna", all.x = TRUE)
tabla_df <-tabla_df[!(tabla_df$Promedio == 0),]

tabla_df2 <<- tabla_df
# variable_in<-"ytotcor"
# ingreso_rds <- paste("e6a_",variable_in,"_",anio,"pruebas.rds", sep="")

}

data_df3 <- data.frame()
for (n in 1:7){
  funcion1(n)

  data_df3 <- rbind(data_df3,tabla_df2)

}
variable_e <- switch(var1,"s12","o1","o2","o3","o6")
assign(paste0("tabla_",variable_e),data_df3)
# print(data_df)
}

3 s12 A que sistema previsional de salud pertenece

datatable(tabla_s12, extensions = 'Buttons', escape = FALSE, rownames = FALSE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'hitStats'),
          list(extend='pdf',
            filename= 'hitStats')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c("Promedio del Ingreso","Desviación standard","Gini") ,mark = "", digits=3)

4 o1 La semana pasada, ¿trabajó al menos una hora?

datatable(tabla_o1, extensions = 'Buttons', escape = FALSE, rownames = FALSE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'hitStats'),
          list(extend='pdf',
            filename= 'hitStats')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c("Promedio del Ingreso","Desviación standard","Gini") ,mark = "", digits=3)

5 o2 ¿Aunque no trabajó la semana pasada, ¿realizó alguna actividad?

datatable(tabla_o2, extensions = 'Buttons', escape = FALSE, rownames = FALSE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'hitStats'),
          list(extend='pdf',
            filename= 'hitStats')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c("Promedio del Ingreso","Desviación standard","Gini") ,mark = "", digits=3)

6 o3 ¿Aunque no trabajó la semana pasada, ¿tenía algún empleo del cual estuvo ausente temporalmente por licencia, huelga, enfermedad, vacaciones u otra razón?

datatable(tabla_o3, extensions = 'Buttons', escape = FALSE, rownames = FALSE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'hitStats'),
          list(extend='pdf',
            filename= 'hitStats')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c("Promedio del Ingreso","Desviación standard","Gini") ,mark = "", digits=3)

7 O6 ¿Buscó trabajo remunerado en las últimas cuatro semanas?

datatable(tabla_o6, extensions = 'Buttons', escape = FALSE, rownames = FALSE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'hitStats'),
          list(extend='pdf',
            filename= 'hitStats')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c("Promedio del Ingreso","Desviación standard","Gini") ,mark = "", digits=3)