library(haven)
EHPM_2019 <- read_sav("C:/Users/johan/OneDrive/Escritorio/prueba/EHPM 2019.sav")

# Cargar librerias 
library(openxlsx)
library(dplyr)
library(tidyverse)
library(knitr)
library(kableExtra)

# Selección de variables de interés 
educacion2019 <- EHPM_2019 %>% select(r204, r206a, r208, r219)

# Exportar a excel
write.xlsx(educacion2019, file = "C:/Users/johan/OneDrive/Escritorio/prueba/educacion2019.xlsx")

nivel_edu2019 <- educacion2019 %>%
  filter(!is.na(r204)) %>% # Excluir filas donde r204 es NA
  group_by(r204) %>%
  summarise(total = n()) %>%
  mutate(nivel = case_when(
    r204 == 0 ~ "Educación inicial",
    r204 == 1 ~ "Parvularia (1° a 3°)",
    r204 == 2 ~ "Básica (1° a 9°)",
    r204 == 3 ~ "Media (10° a 13°)",
    r204 == 4 ~ "Superior universitario (1° a 15°)",
    r204 == 4.1 ~ "Curso de nivelación (1°)",
    r204 == 5 ~ "Superior no universitario (1° a 3°)",
    r204 == 6 ~ "Educación especial (ciclos I, II, III, IV)",
    r204 == 7 ~ "Otros",
    TRUE ~ "No clasificado"
  ))

g2019 <- nivel_edu2019 %>%  select(nivel, total)

if (!"total" %in% g2019$nivel) {
  g2019 <- g2019 %>%
    bind_rows(tibble(nivel = "total", total = sum(g2019$total)))
}

g2019%>%
  kable(caption = "Tabla I: Frecuencia de cada nivel educativo durante año 2019") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2019")

repe_edu2019 <- educacion2019 %>%
  filter(!is.na(r206a)) %>% # Excluir filas donde r206a es NA
  group_by(r206a) %>%
  summarise(total = n()) %>%
  mutate(Repite = case_when(
    r206a == 1 ~ "Sí",
    r206a == 2 ~ "No",
    TRUE ~ "No clasificado"
  ))

gg2019 <- repe_edu2019 %>% select(Repite, total)

if (!"total" %in% gg2019$Repite) {
  gg2019 <- gg2019 %>%
    bind_rows(tibble(Repite = "total", total = sum(gg2019$total)))
}

gg2019%>%
  kable(caption = "Tabla II: Frecuencia de repetir curso durante año 2019") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2019")

razon_edu2019 <- educacion2019 %>%
  filter(!is.na(r208)) %>% # Excluir filas donde r204 es NA
  group_by(r208) %>%
  summarise(total = n()) %>%
  mutate("Razón de abandono" = case_when(
    r208 == 1 ~ "Necesita trabajar",
    r208 == 2 ~ "Problemas en el hogar",
    r208 == 3 ~ "Muy caro",
    r208 == 4 ~ "Por enfermedad",
    r208 == 5 ~ "Discapacidad",
    r208 == 6 ~ "Padre y/o madre no quieren que estudie",
    r208 == 7 ~ "Por la edad",
    r208 == 8 ~ "La escuela está muy lejos",
    r208 == 9 ~ "No le interesa estudiar",
    r208 == 10 ~ "Quehaceres domésticos",
    r208 == 11 ~ "Repite mucho",
    r208 == 12 ~ "Docentes no enseñan",
    r208 == 13 ~ "La escuela está inhabilitada",
    r208 == 14 ~ "Por las maras",
    r208 == 15 ~ "Otros___________________",
    TRUE ~ "No clasificado"
  ))

ggg2019 <- razon_edu2019 %>% select("Razón de abandono", total)

if (!"total" %in% ggg2019$"Razón de abandono") {
  ggg2019 <- ggg2019 %>%
    bind_rows(tibble("Razón de abandono" = "total", total = sum(ggg2019$total)))
}

ggg2019%>%
  kable(caption = "Tabla III: Frecuencia de cada razón de abandono durante el año 2019") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2019")

no_edu2019 <- educacion2019 %>%
  filter(!is.na(r219)) %>% # Excluir filas donde r204 es NA
  group_by(r219) %>%
  summarise(total = n()) %>%
  mutate("Razón de no estudio" = case_when(
    r219 == 1 ~ "Necesita trabajar",
    r219 == 2 ~ "Causas del hogar",
    r219 == 3 ~ "Muy caro",
    r219 == 4 ~ "Por enfermedad",
    r219 == 5 ~ "Discapacidad",
    r219 == 6 ~ "Padre y/o madre no quieren",
    r219 == 7 ~ "Por la edad",
    r219 == 8 ~ "Finalizó sus estudios",
    r219 == 9 ~ "No hay escuela cercana",
    r219 == 10 ~ "No le interesa",
    r219 == 11 ~ "Repite mucho",
    r219 == 12 ~ "Trabajo doméstico y cuidado",
    r219 == 13 ~ "Centro de enseñanza inhabilitado",
    r219 == 14 ~ "No hay escuela para adultos",
    r219 == 15 ~ "Por embarazo ",
    r219 == 16 ~ "Por maternidad",
    r219 == 17 ~ "otros",
    r219 == 18 ~ "Inseguridad",
    TRUE ~ "No clasificado"
  ))

gggg2019 <- no_edu2019 %>% select("Razón de no estudio", total)

if (!"total" %in% gggg2019$"Razón de no estudio") {
  gggg2019 <- gggg2019 %>%
    bind_rows(tibble("Razón de no estudio" = "total", total = sum(gggg2019$total)))
}

gggg2019%>%
  kable(caption = "Tabla IV: Frecuencia de cada razón de no estudio durante el año 2019") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2019")

library(haven)
EHPM_2020 <- read_sav("C:/Users/johan/OneDrive/Escritorio/prueba/EHPM 2020.sav")

# Selección de variables de interés 
educacion2020 <- EHPM_2020 %>% select(r204, r205, r208, r219)

# Exportar a excel
write.xlsx(educacion2020, file = "C:/Users/johan/OneDrive/Escritorio/prueba/educacion2020.xlsx")

nivel_edu2020 <- educacion2020 %>%
  filter(!is.na(r204)) %>% # Excluir filas donde r204 es NA
  group_by(r204) %>%
  summarise(total = n()) %>%
  mutate(nivel = case_when(
    r204 == 0 ~ "Educación inicial",
    r204 == 1 ~ "Parvularia (1° a 3°)",
    r204 == 2 ~ "Básica (1° a 9°)",
    r204 == 3 ~ "Media (10° a 13°)",
    r204 == 4 ~ "Superior universitario (1° a 15°)",
    r204 == 4.1 ~ "Curso de nivelación (1°)",
    r204 == 5 ~ "Superior no universitario (1° a 3°)",
    r204 == 6 ~ "Educación especial (ciclos I, II, III, IV)",
    r204 == 7 ~ "Otros",
    TRUE ~ "No clasificado"
  ))

g2020 <- nivel_edu2020 %>% select(nivel, total)

if (!"total" %in% g2020$nivel) {
  g2020 <- g2020 %>%
    bind_rows(tibble(nivel = "total", total = sum(g2020$total)))
}

g2020 %>%
  kable(caption = "Tabla V: Frecuencia de cada nivel educativo durante año 2020") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2020")

repe_edu2020 <- educacion2020 %>%
  filter(!is.na(r205)) %>% # Excluir filas donde r205 es NA
  group_by(r205) %>%
  summarise(total = n()) %>%
  mutate(Repite = case_when(
    r205 == 1 ~ "Sí",
    r205 == 2 ~ "No",
    TRUE ~ "No clasificado"
  ))

gg2020 <- repe_edu2020 %>% select(Repite, total)

if (!"total" %in% gg2020$Repite) {
  gg2020 <- gg2020 %>%
    bind_rows(tibble(Repite = "total", total = sum(gg2020$total)))
}

gg2020 %>%
  kable(caption = "Tabla VI: Frecuencia de repetir curso durante año 2020") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2020")

razon_edu2020 <- educacion2020 %>%
  filter(!is.na(r208)) %>% # Excluir filas donde r208 es NA
  group_by(r208) %>%
  summarise(total = n()) %>%
  mutate("Razón de abandono" = case_when(
    r208 == 1 ~ "Necesita trabajar",
    r208 == 2 ~ "Problemas en el hogar",
    r208 == 3 ~ "Muy caro",
    r208 == 4 ~ "Por enfermedad",
    r208 == 5 ~ "Discapacidad",
    r208 == 6 ~ "Padre y/o madre no quieren que estudie",
    r208 == 7 ~ "Por la edad",
    r208 == 8 ~ "La escuela está muy lejos",
    r208 == 9 ~ "No le interesa estudiar",
    r208 == 10 ~ "Quehaceres domésticos",
    r208 == 11 ~ "Repite mucho",
    r208 == 12 ~ "Docentes no enseñan",
    r208 == 13 ~ "La escuela está inhabilitada",
    r208 == 14 ~ "Por las maras",
    r208 == 15 ~ "Otros___________________",
    TRUE ~ "No clasificado"
  ))

ggg2020 <- razon_edu2020 %>% select("Razón de abandono", total)

if (!"total" %in% ggg2020$"Razón de abandono") {
  ggg2020 <- ggg2020 %>%
    bind_rows(tibble("Razón de abandono" = "total", total = sum(ggg2020$total)))
}

ggg2020 %>%
  kable(caption = "Tabla VII: Frecuencia de cada razón de abandono durante el año 2020") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2020")

no_edu2020 <- educacion2020 %>%
  filter(!is.na(r219)) %>% # Excluir filas donde r219 es NA
  group_by(r219) %>%
  summarise(total = n()) %>%
  mutate("Razón de no estudio" = case_when(
    r219 == 1 ~ "Necesita trabajar",
    r219 == 2 ~ "Causas del hogar",
    r219 == 3 ~ "Muy caro",
    r219 == 4 ~ "Por enfermedad",
    r219 == 5 ~ "Discapacidad",
    r219 == 6 ~ "Padre y/o madre no quieren",
    r219 == 7 ~ "Por la edad",
    r219 == 8 ~ "Finalizó sus estudios",
    r219 == 9 ~ "No hay escuela cercana",
    r219 == 10 ~ "No le interesa",
    r219 == 11 ~ "Repite mucho",
    r219 == 12 ~ "Trabajo doméstico y cuidado",
    r219 == 13 ~ "Centro de enseñanza inhabilitado",
    r219 == 14 ~ "No hay escuela para adultos",
    r219 == 15 ~ "Por embarazo ",
    r219 == 16 ~ "Por maternidad",
    r219 == 17 ~ "otros",
    r219 == 18 ~ "Inseguridad", 
    TRUE ~ "No clasificado"
  ))

gggg2020 <- no_edu2020 %>% select("Razón de no estudio", total)

if (!"total" %in% gggg2020$"Razón de no estudio") {
  gggg2020 <- gggg2020 %>%
    bind_rows(tibble("Razón de no estudio" = "total", total = sum(gggg2020$total)))
}

gggg2020 %>%
  kable(caption = "Tabla VIII: Frecuencia de cada razón de no estudio durante el año 2020") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2020")

library(haven)
EHPM_2021 <- read_sav("C:/Users/johan/OneDrive/Escritorio/prueba/EHPM 2021.sav")

# Selección de variables de interés 
educacion2021 <- EHPM_2021 %>% select(r204, r205, r208, r219)

# Exportar a excel
write.xlsx(educacion2021, file = "C:/Users/johan/OneDrive/Escritorio/prueba/educacion2021.xlsx")

nivel_edu2021 <- educacion2021 %>%
  filter(!is.na(r204)) %>% # Excluir filas donde r204 es NA
  group_by(r204) %>%
  summarise(total = n()) %>%
  mutate(nivel = case_when(
    r204 == 0 ~ "Educación inicial",
    r204 == 1 ~ "Parvularia (1° a 3°)",
    r204 == 2 ~ "Básica (1° a 9°)",
    r204 == 3 ~ "Media (10° a 13°)",
    r204 == 4 ~ "Superior universitario (1° a 15°)",
    r204 == 4.1 ~ "Curso de nivelación (1°)",
    r204 == 5 ~ "Superior no universitario (1° a 3°)",
    r204 == 6 ~ "Educación especial (ciclos I, II, III, IV)",
    r204 == 7 ~ "Otros",
    TRUE ~ "No clasificado"
  ))

g2021 <- nivel_edu2021 %>% select(nivel, total)

if (!"total" %in% g2021$nivel) {
  g2021 <- g2021 %>%
    bind_rows(tibble(nivel = "total", total = sum(g2021$total)))
}

g2021 %>%
  kable(caption = "Tabla IX: Frecuencia de cada nivel educativo durante año 2021") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2021")

repe_edu2021 <- educacion2021 %>%
  filter(!is.na(r205)) %>% # Excluir filas donde r205 es NA
  group_by(r205) %>%
  summarise(total = n()) %>%
  mutate(Repite = case_when(
    r205 == 1 ~ "Sí",
    r205 == 2 ~ "No",
    TRUE ~ "No clasificado"
  ))

gg2021 <- repe_edu2021 %>% select(Repite, total)

if (!"total" %in% gg2021$Repite) {
  gg2021 <- gg2021 %>%
    bind_rows(tibble(Repite = "total", total = sum(gg2021$total)))
}

gg2021 %>%
  kable(caption = "Tabla X: Frecuencia de repetir curso durante año 2021") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2021")

razon_edu2021 <- educacion2021 %>%
  filter(!is.na(r208)) %>% # Excluir filas donde r208 es NA
  group_by(r208) %>%
  summarise(total = n()) %>%
  mutate("Razón de abandono" = case_when(
    r208 == 1 ~ "Necesita trabajar",
    r208 == 2 ~ "Problemas en el hogar",
    r208 == 3 ~ "Muy caro",
    r208 == 4 ~ "Por enfermedad",
    r208 == 5 ~ "Discapacidad",
    r208 == 6 ~ "Padre y/o madre no quieren que estudie",
    r208 == 7 ~ "Por la edad",
    r208 == 8 ~ "La escuela está muy lejos",
    r208 == 9 ~ "No le interesa estudiar",
    r208 == 10 ~ "Quehaceres domésticos",
    r208 == 11 ~ "Repite mucho",
    r208 == 12 ~ "Docentes no enseñan",
    r208 == 13 ~ "La escuela está inhabilitada",
    r208 == 14 ~ "Por las maras",
    r208 == 15 ~ "Otros___________________",
    TRUE ~ "No clasificado"
  ))

ggg2021 <- razon_edu2021 %>% select("Razón de abandono", total)

if (!"total" %in% ggg2021$"Razón de abandono") {
  ggg2021 <- ggg2021 %>%
    bind_rows(tibble("Razón de abandono" = "total", total = sum(ggg2021$total)))
}

ggg2021 %>%
  kable(caption = "Tabla XI: Frecuencia de cada razón de abandono durante el año 2021") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2021")

no_edu2021 <- educacion2021 %>%
  filter(!is.na(r219)) %>% # Excluir filas donde r219 es NA
  group_by(r219) %>%
  summarise(total = n()) %>%
  mutate("Razón de no estudio" = case_when(
    r219 == 1 ~ "Necesita trabajar",
    r219 == 2 ~ "Causas del hogar",
    r219 == 3 ~ "Muy caro",
    r219 == 4 ~ "Por enfermedad",
    r219 == 5 ~ "Discapacidad",
    r219 == 6 ~ "Padre y/o madre no quieren",
    r219 == 7 ~ "Por la edad",
    r219 == 8 ~ "Finalizó sus estudios",
    r219 == 9 ~ "No hay escuela cercana",
    r219 == 10 ~ "No le interesa",
    r219 == 11 ~ "Repite mucho",
    r219 == 12 ~ "Trabajo doméstico y cuidado",
    r219 == 13 ~ "Centro de enseñanza inhabilitado",
    r219 == 14 ~ "No hay escuela para adultos",
    r219 == 15 ~ "Por embarazo",
    r219 == 16 ~ "Por maternidad",
    r219 == 17 ~ "otros",
    r219 == 18 ~ "Inseguridad",  # Cambié r219 == 13 a r219 == 18 porque parece ser un error.
    TRUE ~ "No clasificado"
  ))

gggg2021 <- no_edu2021 %>% select("Razón de no estudio", total)

if (!"total" %in% gggg2021$"Razón de no estudio") {
  gggg2021 <- gggg2021 %>%
    bind_rows(tibble("Razón de no estudio" = "total", total = sum(gggg2021$total)))
}

gggg2021 %>%
  kable(caption = "Tabla XII: Frecuencia de cada razón de no estudio durante el año 2021") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2021")

library(haven)
EHPM_2022 <- read_sav("C:/Users/johan/OneDrive/Escritorio/prueba/EHPM 2022.sav")

# Selección de variables de interés 
educacion2022 <- EHPM_2022 %>% select(r204, r205, r208, r219)

# Exportar a excel
write.xlsx(educacion2022, file = "C:/Users/johan/OneDrive/Escritorio/prueba/educacion2022.xlsx")

nivel_edu2022 <- educacion2022 %>%
  filter(!is.na(r204)) %>% # Excluir filas donde r204 es NA
  group_by(r204) %>%
  summarise(total = n()) %>%
  mutate(nivel = case_when(
    r204 == 0 ~ "Educación inicial",
    r204 == 1 ~ "Parvularia (1° a 3°)",
    r204 == 2 ~ "Básica (1° a 9°)",
    r204 == 3 ~ "Media (10° a 13°)",
    r204 == 4 ~ "Superior universitario (1° a 15°)",
    r204 == 4.1 ~ "Curso de nivelación (1°)",
    r204 == 5 ~ "Superior no universitario (1° a 3°)",
    r204 == 6 ~ "Educación especial (ciclos I, II, III, IV)",
    r204 == 7 ~ "Otros",
    TRUE ~ "No clasificado"
  ))

g2022 <- nivel_edu2022 %>% select(nivel, total)

if (!"total" %in% g2022$nivel) {
  g2022 <- g2022 %>%
    bind_rows(tibble(nivel = "total", total = sum(g2022$total)))
}

g2022 %>%
  kable(caption = "Tabla XIII: Frecuencia de cada nivel educativo durante año 2022") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2022")

repe_edu2022 <- educacion2022 %>%
  filter(!is.na(r205)) %>% # Excluir filas donde r205 es NA
  group_by(r205) %>%
  summarise(total = n()) %>%
  mutate(Repite = case_when(
    r205 == 1 ~ "Sí",
    r205 == 2 ~ "No",
    TRUE ~ "No clasificado"
  ))

gg2022 <- repe_edu2022 %>% select(Repite, total)

if (!"total" %in% gg2022$Repite) {
  gg2022 <- gg2022 %>%
    bind_rows(tibble(Repite = "total", total = sum(gg2022$total)))
}

gg2022 %>%
  kable(caption = "Tabla XIV: Frecuencia de repetir curso durante año 2022") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2022")

razon_edu2022 <- educacion2022 %>%
  filter(!is.na(r208)) %>% # Excluir filas donde r208 es NA
  group_by(r208) %>%
  summarise(total = n()) %>%
  mutate("Razón de abandono" = case_when(
    r208 == 1 ~ "Necesita trabajar",
    r208 == 2 ~ "Problemas en el hogar",
    r208 == 3 ~ "Muy caro",
    r208 == 4 ~ "Por enfermedad",
    r208 == 5 ~ "Discapacidad",
    r208 == 6 ~ "Padre y/o madre no quieren que estudie",
    r208 == 7 ~ "Por la edad",
    r208 == 8 ~ "La escuela está muy lejos",
    r208 == 9 ~ "No le interesa estudiar",
    r208 == 10 ~ "Quehaceres domésticos",
    r208 == 11 ~ "Repite mucho",
    r208 == 12 ~ "Docentes no enseñan",
    r208 == 13 ~ "La escuela está inhabilitada",
    r208 == 14 ~ "Por las maras",
    r208 == 15 ~ "Otros___________________",
    TRUE ~ "No clasificado"
  ))

ggg2022 <- razon_edu2022 %>% select("Razón de abandono", total)

if (!"total" %in% ggg2022$"Razón de abandono") {
  ggg2022 <- ggg2022 %>%
    bind_rows(tibble("Razón de abandono" = "total", total = sum(ggg2022$total)))
}

ggg2022 %>%
  kable(caption = "Tabla XV: Frecuencia de cada razón de abandono durante el año 2022") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2022")

no_edu2022 <- educacion2022 %>%
  filter(!is.na(r219)) %>% # Excluir filas donde r219 es NA
  group_by(r219) %>%
  summarise(total = n()) %>%
  mutate("Razón de no estudio" = case_when(
    r219 == 1 ~ "Necesita trabajar",
    r219 == 2 ~ "Causas del hogar",
    r219 == 3 ~ "Muy caro",
    r219 == 4 ~ "Por enfermedad",
    r219 == 5 ~ "Discapacidad",
    r219 == 6 ~ "Padre y/o madre no quieren",
    r219 == 7 ~ "Por la edad",
    r219 == 8 ~ "Finalizó sus estudios",
    r219 == 9 ~ "No hay escuela cercana",
    r219 == 10 ~ "No le interesa",
    r219 == 11 ~ "Repite mucho",
    r219 == 12 ~ "Trabajo doméstico y cuidado",
    r219 == 13 ~ "Centro de enseñanza inhabilitado",
    r219 == 14 ~ "No hay escuela para adultos",
    r219 == 15 ~ "Por embarazo",
    r219 == 16 ~ "Por maternidad",
    r219 == 17 ~ "otros",
    r219 == 18 ~ "Inseguridad",
    TRUE ~ "No clasificado"
  ))

gggg2022 <- no_edu2022 %>% select("Razón de no estudio", total)

if (!"total" %in% gggg2022$"Razón de no estudio") {
  gggg2022 <- gggg2022 %>%
    bind_rows(tibble("Razón de no estudio" = "total", total = sum(gggg2022$total)))
}

gggg2022 %>%
  kable(caption = "Tabla XVI: Frecuencia de cada razón de no estudio durante el año 2022") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2022")

library(haven)
EHPM_2023 <- read_sav("C:/Users/johan/OneDrive/Escritorio/prueba/EHPM 2023.sav")

# Selección de variables de interés 
educacion2023 <- EHPM_2023 %>% select(r204, r205, r208, r219)

# Exportar a excel
write.xlsx(educacion2023, file = "C:/Users/johan/OneDrive/Escritorio/prueba/educacion2023.xlsx")

nivel_edu2023 <- educacion2023 %>%
  filter(!is.na(r204)) %>% # Excluir filas donde r204 es NA
  group_by(r204) %>%
  summarise(total = n()) %>%
  mutate(nivel = case_when(
    r204 == 0 ~ "Educación inicial",
    r204 == 1 ~ "Parvularia (1° a 3°)",
    r204 == 2 ~ "Básica (1° a 9°)",
    r204 == 3 ~ "Media (10° a 13°)",
    r204 == 4 ~ "Superior universitario (1° a 15°)",
    r204 == 4.1 ~ "Curso de nivelación (1°)",
    r204 == 5 ~ "Superior no universitario (1° a 3°)",
    r204 == 6 ~ "Educación especial (ciclos I, II, III, IV)",
    r204 == 7 ~ "Otros",
    TRUE ~ "No clasificado"
  ))

g2023 <- nivel_edu2023 %>% select(nivel, total)

if (!"total" %in% g2023$nivel) {
  g2023 <- g2023 %>%
    bind_rows(tibble(nivel = "total", total = sum(g2023$total)))
}

g2023 %>%
  kable(caption = "Tabla XVII: Frecuencia de cada nivel educativo durante año 2023") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2023")

repe_edu2023 <- educacion2023 %>%
  filter(!is.na(r205)) %>% # Excluir filas donde r205 es NA
  group_by(r205) %>%
  summarise(total = n()) %>%
  mutate(Repite = case_when(
    r205 == 1 ~ "Sí",
    r205 == 2 ~ "No",
    TRUE ~ "No clasificado"
  ))

gg2023 <- repe_edu2023 %>% select(Repite, total)

if (!"total" %in% gg2023$Repite) {
  gg2023 <- gg2023 %>%
    bind_rows(tibble(Repite = "total", total = sum(gg2023$total)))
}

gg2023 %>%
  kable(caption = "Tabla XVIII: Frecuencia de repetir curso durante año 2023") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2023")

razon_edu2023 <- educacion2023 %>%
  filter(!is.na(r208)) %>% # Excluir filas donde r208 es NA
  group_by(r208) %>%
  summarise(total = n()) %>%
  mutate("Razón de abandono" = case_when(
    r208 == 1 ~ "Necesita trabajar",
    r208 == 2 ~ "Problemas en el hogar",
    r208 == 3 ~ "Muy caro",
    r208 == 4 ~ "Por enfermedad",
    r208 == 5 ~ "Discapacidad",
    r208 == 6 ~ "Padre y/o madre no quieren que estudie",
    r208 == 7 ~ "Por la edad",
    r208 == 8 ~ "La escuela está muy lejos",
    r208 == 9 ~ "No le interesa estudiar",
    r208 == 10 ~ "Quehaceres domésticos",
    r208 == 11 ~ "Repite mucho",
    r208 == 12 ~ "Docentes no enseñan",
    r208 == 13 ~ "La escuela está inhabilitada",
    r208 == 14 ~ "Por las maras",
    r208 == 15 ~ "Otros___________________",
    TRUE ~ "No clasificado"
  ))

ggg2023 <- razon_edu2023 %>% select("Razón de abandono", total)

if (!"total" %in% ggg2023$"Razón de abandono") {
  ggg2023 <- ggg2023 %>%
    bind_rows(tibble("Razón de abandono" = "total", total = sum(ggg2023$total)))
}

ggg2023 %>%
  kable(caption = "Tabla XIX: Frecuencia de cada razón de abandono durante el año 2023") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2023")

no_edu2023 <- educacion2023 %>%
  filter(!is.na(r219)) %>% # Excluir filas donde r219 es NA
  group_by(r219) %>%
  summarise(total = n()) %>%
  mutate("Razón de no estudio" = case_when(
    r219 == 1 ~ "Necesita trabajar",
    r219 == 2 ~ "Causas del hogar",
    r219 == 3 ~ "Muy caro",
    r219 == 4 ~ "Por enfermedad",
    r219 == 5 ~ "Discapacidad",
    r219 == 6 ~ "Padre y/o madre no quieren",
    r219 == 7 ~ "Por la edad",
    r219 == 8 ~ "Finalizó sus estudios",
    r219 == 9 ~ "No hay escuela cercana",
    r219 == 10 ~ "No le interesa",
    r219 == 11 ~ "Repite mucho",
    r219 == 12 ~ "Trabajo doméstico y cuidado",
    r219 == 13 ~ "Centro de enseñanza inhabilitado",
    r219 == 14 ~ "No hay escuela para adultos",
    r219 == 15 ~ "Por embarazo",
    r219 == 16 ~ "Por maternidad",
    r219 == 17 ~ "otros",
    r219 == 18 ~ "Inseguridad",
    TRUE ~ "No clasificado"
  ))

gggg2023 <- no_edu2023 %>% select("Razón de no estudio", total)

if (!"total" %in% gggg2023$"Razón de no estudio") {
  gggg2023 <- gggg2023 %>%
    bind_rows(tibble("Razón de no estudio" = "total", total = sum(gggg2023$total)))
}

gggg2023 %>%
  kable(caption = "Tabla XX: Frecuencia de cada razón de no estudio durante el año 2023") %>%
  kable_styling("hover", full_width = FALSE) %>%
  row_spec(0, bold = TRUE) %>%
  kable_classic(html_font = "Arial", font_size = 12) %>%
  footnote(general = "Elaboración propia en base a la Encuesta de Hogares de Propósitos Múltiples de 2023")