casen_2006 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2006_c.rds")
casen_2006 <- mutate_if(casen_2006, is.factor, as.character)
casen_2009 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2009_c.rds")
casen_2009 <- mutate_if(casen_2009, is.factor, as.character)
casen_2011 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2011_c.rds")
casen_2011 <- mutate_if(casen_2011, is.factor, as.character)
casen_2013 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2013_c.rds")
casen_2013 <- mutate_if(casen_2013, is.factor, as.character)
casen_2015 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2015_c.rds")
casen_2015 <- mutate_if(casen_2015, is.factor, as.character)
casen_2017 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2017_c.rds")
casen_2017 <- mutate_if(casen_2017, is.factor, as.character)
casen_2020 <- readRDS(file = "C:/Users/enamo/Desktop/Shiny-R/Casen_en_pandemia_2020/casen/casen_2020_c.rds")
casen_2020 <- mutate_if(casen_2020, is.factor, as.character)

casen_2011$z[casen_2011$z == "Urbana"] <- "Urbano"

receptaculo <- data.frame(
  Variable = c("Rural","Urbano")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"Z","ZONA","z","zona","zona","zona","zona")
################ si solo son 2 categorias no hay que modificar esta seccion
tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
anios <- seq(2006,2020,1)
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")# <---- aqui se establecen las cabeceras que puden variar dependiendo los años en que estan las variables
tablav <- mutate_if(tablav, is.character, as.numeric)

dt1 <- data.frame()
dt2 <- data.frame()

linearMod_r <- lm(Rural ~ año, data=tablav)
linearMod_u <- lm(Urbano ~ año, data=tablav)


for (n in anios) {
  r <- linearMod_r$coefficients[1]+linearMod_r$coefficients[2]*n
  dt1 <- rbind(dt1,r)
  
  u <- linearMod_u$coefficients[1]+linearMod_u$coefficients[2]*n
  dt2 <- rbind(dt2,u)

  dt3 <- cbind(dt1,dt2)
}

dt4 <- as.data.frame(t(dt3), row.names = F)
colnames(dt4) <- anios
################
row.names(dt4) <- c("Rural","Urbano")
datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'ruralidad'),
          list(extend='pdf',
            filename= 'ruralidad')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c(paste(anios)) ,mark = "", digits=0)

casen_2006$CORTE[casen_2006$CORTE == "Pobre no indigente"] <- "Pobre"
casen_2006$CORTE[casen_2006$CORTE == "Indigente"] <- "Pobre extremo"

casen_2009$CORTE[casen_2009$CORTE == "Pobre no Indigente"] <- "Pobre"
casen_2009$CORTE[casen_2009$CORTE == "Indigente"] <- "Pobre extremo"

casen_2011$corte[casen_2011$corte == "Pobreza extrema"] <- "Pobre extremo"

casen_2013$pobreza_MN[casen_2013$pobreza_MN == "Pobres no extremos"] <- "Pobre"
casen_2013$pobreza_MN[casen_2013$pobreza_MN == "No pobres"] <- "No pobre"
casen_2013$pobreza_MN[casen_2013$pobreza_MN == "Pobres extremos"] <- "Pobre extremo"

casen_2015$pobreza[casen_2015$pobreza == "Pobres no extremos"] <- "Pobre"
casen_2015$pobreza[casen_2015$pobreza == "No pobres"] <- "No pobre"
casen_2015$pobreza[casen_2015$pobreza == "Pobres extremos"] <- "Pobre extremo"

casen_2017$pobreza[casen_2017$pobreza == "Pobres no extremos"] <- "Pobre"
casen_2017$pobreza[casen_2017$pobreza == "No pobres"] <- "No pobre"
casen_2017$pobreza[casen_2017$pobreza == "Pobres extremos"] <- "Pobre extremo"

casen_2020$pobreza[casen_2020$pobreza == "Pobres no extremos"] <- "Pobre"
casen_2020$pobreza[casen_2020$pobreza == "No pobres"] <- "No pobre"
casen_2020$pobreza[casen_2020$pobreza == "Pobres extremos"] <- "Pobre extremo"

receptaculo <- data.frame(
  Variable = c("No pobre","Pobre","Pobre extremo")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"CORTE","CORTE","corte","pobreza_MN","pobreza","pobreza","pobreza")

tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")
tablav <- mutate_if(tablav, is.character, as.numeric)

dt1 <- data.frame()
dt2 <- data.frame()
dt3 <- data.frame()

linearMod_1 <- lm(`No pobre` ~ año, data=tablav)
linearMod_2 <- lm(Pobre ~ año, data=tablav)
linearMod_3 <- lm(`Pobre extremo` ~ año, data=tablav)

for (n in seq(2006,2020,1)) {
  p1 <- linearMod_1$coefficients[1]+linearMod_1$coefficients[2]*n
  dt1 <- rbind(dt1,p1)
  
  p2 <- linearMod_2$coefficients[1]+linearMod_2$coefficients[2]*n
  dt2 <- rbind(dt2,p2)
  
  p3 <- linearMod_3$coefficients[1]+linearMod_3$coefficients[2]*n
  dt3 <- rbind(dt3,p3)

  dt4 <- cbind(dt1,dt2,dt3)
}

dt4 <- as.data.frame(t(dt4), row.names = F)
colnames(dt4) <- seq(2006,2020,1)
row.names(dt4) <- c("Pobre","No pobre","Pobre extremo")
datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'ruralidad'),
          list(extend='pdf',
            filename= 'ruralidad')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c(paste(seq(2006,2020,1))) ,mark = "", digits=0)

casen_2006$E1[casen_2006$E1 == "No sabe /Sin dato"] <- NA

casen_2011$e1[casen_2011$e1 == "Sí, lee y escribe"] <- "Sí"
casen_2011$e1[casen_2011$e1 == "No, sólo lee"] <- "No"
casen_2011$e1[casen_2011$e1 == "No, ninguno"] <- "No"
casen_2011$e1[casen_2011$e1 == "No, sólo escribe"] <- "No"

casen_2013$e1[casen_2013$e1 == "Sí, lee y escribe"] <- "Sí"
casen_2013$e1[casen_2013$e1 == "No, ninguno"] <- "No"
casen_2013$e1[casen_2013$e1 == "No, sólo lee"] <- "No"
casen_2013$e1[casen_2013$e1 == "No, sólo escribe"] <- "No"
casen_2013$e1[casen_2013$e1 == "NS/NR"] <- NA

casen_2015$e1[casen_2015$e1 == "Sí, lee y escribe"] <- "Sí"
casen_2015$e1[casen_2015$e1 == "No, ninguno"] <- "No"
casen_2015$e1[casen_2015$e1 == "No, sólo lee"] <- "No"
casen_2015$e1[casen_2015$e1 == "No, sólo escribe"] <- "No"

casen_2017$e1[casen_2017$e1 == "Sí, lee y escribe"] <- "Sí"
casen_2017$e1[casen_2017$e1 == "No, sólo lee"] <- "No"
casen_2017$e1[casen_2017$e1 == "No, ninguno"] <- "No"
casen_2017$e1[casen_2017$e1 == "No sabe/responde"] <- NA
casen_2017$e1[casen_2017$e1 == "No, sólo escribe"] <- "No"

casen_2020$e1[casen_2020$e1 == 1] <- "Sí"
casen_2020$e1[casen_2020$e1 == 0] <- "No"

receptaculo <- data.frame(
  Variable = c("Sí","No")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"E1","E1","e1","e1","e1","e1","e1")

tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")
tablav <- mutate_if(tablav, is.character, as.numeric)

dt1 <- data.frame()
dt2 <- data.frame() 

linearMod_1 <- lm(Sí ~ año, data=tablav)
linearMod_2 <- lm(No ~ año, data=tablav) 

for (n in seq(2006,2020,1)) {
  p1 <- linearMod_1$coefficients[1]+linearMod_1$coefficients[2]*n
  dt1 <- rbind(dt1,p1)
  
  p2 <- linearMod_2$coefficients[1]+linearMod_2$coefficients[2]*n
  dt2 <- rbind(dt2,p2)
 

  dt3 <- cbind(dt1,dt2)
}

dt4 <- as.data.frame(t(dt3), row.names = F)
colnames(dt4) <- seq(2006,2020,1)
row.names(dt4) <- c("Sí","No")
datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'ruralidad'),
          list(extend='pdf',
            filename= 'ruralidad')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c(paste(seq(2006,2020,1))) ,mark = "", digits=0)

receptaculo <- data.frame(
  Variable = c("Sí","No")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"O1","O1","o1","o1","o1","o1","o1")

tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")
tablav <- mutate_if(tablav, is.character, as.numeric)

dt1 <- data.frame()
dt2 <- data.frame() 

linearMod_1 <- lm(Sí ~ año, data=tablav)
linearMod_2 <- lm(No ~ año, data=tablav) 

for (n in seq(2006,2020,1)) {
  p1 <- linearMod_1$coefficients[1]+linearMod_1$coefficients[2]*n
  dt1 <- rbind(dt1,p1)
  
  p2 <- linearMod_2$coefficients[1]+linearMod_2$coefficients[2]*n
  dt2 <- rbind(dt2,p2)
 

  dt3 <- cbind(dt1,dt2)
}

dt4 <- as.data.frame(t(dt3), row.names = F)
colnames(dt4) <- seq(2006,2020,1)
row.names(dt4) <- c("Sí","No")
datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'ruralidad'),
          list(extend='pdf',
            filename= 'ruralidad')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c(paste(seq(2006,2020,1))) ,mark = "", digits=0)

receptaculo <- data.frame(
  Variable = c("Sí","No")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"O3","O3","o3","o3","o3","o3","o3")

tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")
tablav <- mutate_if(tablav, is.character, as.numeric)

dt1 <- data.frame()
dt2 <- data.frame() 

linearMod_1 <- lm(Sí ~ año, data=tablav)
linearMod_2 <- lm(No ~ año, data=tablav) 

for (n in seq(2006,2020,1)) {
  p1 <- linearMod_1$coefficients[1]+linearMod_1$coefficients[2]*n
  dt1 <- rbind(dt1,p1)
  
  p2 <- linearMod_2$coefficients[1]+linearMod_2$coefficients[2]*n
  dt2 <- rbind(dt2,p2)
 

  dt3 <- cbind(dt1,dt2)
}

dt4 <- as.data.frame(t(dt3), row.names = F)
colnames(dt4) <- seq(2006,2020,1)
row.names(dt4) <- c("Sí","No")
datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'ruralidad'),
          list(extend='pdf',
            filename= 'ruralidad')),
          text = 'Download')), scrollX = TRUE))%>%
    formatRound(columns=c(paste(seq(2006,2020,1))) ,mark = "", digits=0)

 receptaculo <- data.frame(
  Variable = c("No pobre","Pobre","Pobre extremo")
)
for (i in 1:7) {
  
  switch(i,
        case = casen <- casen_2006,
        case = casen <- casen_2009,
        case = casen <- casen_2011,
        case = casen <- casen_2013,
        case = casen <- casen_2015,
        case = casen <- casen_2017,
        case = casen <- casen_2020
)

expan <-switch(i,"EXPC","EXPC","expc_full","expc","expc_todas","expc","expc")
var1  <-switch(i,"CORTE","CORTE","corte","pobreza_MN","pobreza","pobreza","pobreza")

tabla_matp <-xtabs(casen[,(expan)]~casen[,(var1)], data = casen)
tabla_matp <- as.data.frame(tabla_matp)
receptaculo <- cbind(receptaculo,tabla_matp[,2])
}

tablav <- as.data.frame(t(receptaculo), row.names = F)
colnames(tablav) <- tablav[1,]
tablav <- tablav[-1,]
tablav$año <- c("2006","2009","2011","2013","2015","2017","2020")
tablav <- mutate_if(tablav, is.character, as.numeric)
 
dt1 <- data.frame()
 dt2 <- data.frame(v="V")
 dt3 <- data.frame()
 
for (n in seq(2006,2020,1)) {

 funcion_lm <- function(m){
   linearMod <- lm(tablav[,m] ~ tablav[,ncol(tablav)], data=tablav)
   valores <- linearMod$coefficients[1]+linearMod$coefficients[2]*n
   dt1 <- rbind(dt1,valores)
   dt2 <<- cbind(dt2,dt1)
 
 }
 for (m in 1:3) {
   funcion_lm(m)
 
 }}
 dt2<-dt2[,-1]
 dt3 <- as.data.frame(matrix(unlist(dt2, use.names=FALSE),ncol=3, byrow=TRUE))

 dt4 <- as.data.frame(t(dt3), row.names = F)
 colnames(dt4) <- seq(2006,2020,1)
 row.names(dt4) <- c("No pobre","Pobre","Pobre extremo")
 datatable(dt4, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
           options = list(dom = 'Bfrtip',
           buttons = list('colvis', list(extend = 'collection',
           buttons = list(
           list(extend='copy'),
           list(extend='excel',
             filename = 'ruralidad'),
           list(extend='pdf',
             filename= 'ruralidad')),
           text = 'Download')), scrollX = TRUE))%>%
     formatRound(columns=c(paste(seq(2006,2020,1))) ,mark = "", digits=0)