library(openxlsx)
library(rmarkdown)
library(tidyverse)
library(haven)
library(foreign)
library(survey) Una ruta hacia una carpeta donde almacenaremos los excel que se elaborarán posteriormente, que contarán con dataset de información construida.
ruta <- "C:/Users/Trabajo/Desktop/RDATA" choose.files()## character(0)#participacion2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\PARTICIPACION - 759-Modulo84\\enaho01-2022-800b.dta")
#participacion2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\PARTICIPACION - 759-Modulo84\\enaho01-2022-800b.dta")
#participacion2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\PARTICIPACION - 759-Modulo84\\enaho01-2022-800b.dta")
#participacion2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\PARTICIPACION - 759-Modulo84\\enaho01-2022-800b.dta")
participacion2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\PARTICIPACION - 759-Modulo84\\enaho01-2022-800b.dta")
#caracter2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\CMH - 759-Modulo02\\enaho01-2022-200.dta")
#caracter2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\CMH - 759-Modulo02\\enaho01-2022-200.dta")
#caracter2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\CMH - 759-Modulo02\\enaho01-2022-200.dta")
#caracter2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\CMH - 759-Modulo02\\enaho01-2022-200.dta")
caracter2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\CMH - 759-Modulo02\\enaho01-2022-200.dta")
#sumaria2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SUMARIA - 759-Modulo34\\sumaria-2022.dta")
#sumaria2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SUMARIA - 759-Modulo34\\sumaria-2022.dta")
#sumaria2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SUMARIA - 759-Modulo34\\sumaria-2022.dta")
#sumaria2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SUMARIA - 759-Modulo34\\sumaria-2022.dta")
sumaria2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SUMARIA - 759-Modulo34\\sumaria-2022.dta")
#educacion2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EDUACIÓN - 759-Modulo03\\enaho01a-2022-300.dta")
#educacion2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EDUACIÓN - 759-Modulo03\\enaho01a-2022-300.dta")
#educacion2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EDUACIÓN - 759-Modulo03\\enaho01a-2022-300.dta")
#educacion2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EDUACIÓN - 759-Modulo03\\enaho01a-2022-300.dta")
educacion2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EDUACIÓN - 759-Modulo03\\enaho01a-2022-300.dta")
#salud2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SALUD - 759-Modulo04\\enaho01a-2022-400.dta")
#salud2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SALUD - 759-Modulo04\\enaho01a-2022-400.dta")
#salud2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SALUD - 759-Modulo04\\enaho01a-2022-400.dta")
#salud2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SALUD - 759-Modulo04\\enaho01a-2022-400.dta")
salud2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\SALUD - 759-Modulo04\\enaho01a-2022-400.dta")
#empleo2018 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EMPLEO - 759-Modulo05\\enaho01a-2022-500.dta")
#empleo2019 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EMPLEO - 759-Modulo05\\enaho01a-2022-500.dta")
#empleo2020 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EMPLEO - 759-Modulo05\\enaho01a-2022-500.dta")
#empleo2021 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EMPLEO - 759-Modulo05\\enaho01a-2022-500.dta")
empleo2022 <- read_dta("C:\\Users\\Trabajo\\Desktop\\ENAHO\\EMPLEO - 759-Modulo05\\enaho01a-2022-500.dta")
participacion2022 <- subset(participacion2022, select=c("conglome","hogar","vivienda","codperso","p803"))
sumaria2022 <- subset(sumaria2022, select= c("conglome","hogar","vivienda","pobreza"))
empleo2022 <- subset(empleo2022, select = c("codperso","conglome","vivienda","hogar","p558c"))
educacion2022 <- subset(educacion2022, select = c("codperso","conglome","vivienda","hogar","p300a"))
salud2022 <- subset(salud2022, select = c("codperso","conglome","vivienda","hogar","p401h1","p401h2","p401h3","p401h4","p401h5","p401h6"))Crear variable ‘dup’ Lo que nos permitirá rastrear las filas repetidas ya que pueden haber personas que participen más de una forma
participacion2022 <- participacion2022 %>%
arrange(conglome, vivienda, hogar, codperso) %>%
group_by(conglome, vivienda, hogar, codperso) %>%
mutate(dup = row_number())
table(participacion2022$dup)##
## 1 2 3 4 5 6 7
## 18037 2995 457 69 13 2 1participacion2022 <- participacion2022 %>% filter(dup == 1)enaho_inicial <- left_join(caracter2022, participacion2022, by=c("conglome","vivienda","hogar","codperso"))
enaho_inicial2 <- left_join(enaho_inicial, salud2022, by=c("codperso","conglome","vivienda","hogar"))
enaho_inicial3 <- left_join(enaho_inicial2, educacion2022, by=c("codperso","conglome","vivienda","hogar"))
enaho_inicial4 <- left_join(enaho_inicial3, empleo2022, by=c("codperso","conglome","vivienda","hogar"))
enaho_inicial5 <- left_join(enaho_inicial4, sumaria2022, by=c("conglome","vivienda","hogar"))
enaho <- enaho_inicial5
#ÚLTIMO FILTRO: En años anteriores la variable es "tipodecuestionario".
enaho <- filter(enaho, ticuest01 !=1)
enaho <- filter(enaho, p204==1 & p203 != 8 & p203 != 9)Aquellas variables que nos sirven para realizar las desagregaciones posteriores.
enaho <- enaho %>% mutate(sexo = ifelse(p207==1,"hombre",
ifelse(p207==2,"Mujer",NA)))
table(enaho$sexo, useNA = "alw")##
## hombre Mujer <NA>
## 55419 58559 0Hay que convertir la variable ubigeo.x (el “.x” es producto de la unión de bases de datos) a numérico para no encontrar problemas al momento de recodificar
enaho$ubigeonum <- as.numeric(enaho$ubigeo)
enaho <- enaho %>%
mutate(regiones2 =
ifelse(ubigeonum >= 010101 & ubigeonum <= 010707, "Amazonas",
ifelse(ubigeonum >= 020101 & ubigeonum <= 022008, "Ancash",
ifelse(ubigeonum >= 030101 & ubigeonum <= 030714, "Apurimac",
ifelse(ubigeonum >= 040101 & ubigeonum <= 040811, "Arequipa",
ifelse(ubigeonum >= 050101 & ubigeonum <= 051108, "Ayacucho",
ifelse(ubigeonum >= 060101 & ubigeonum <= 061311, "Cajamarca",
ifelse(ubigeonum >= 070101 & ubigeonum <= 070107, "Callao",
ifelse(ubigeonum >= 080101 & ubigeonum <= 081307, "Cusco",
ifelse(ubigeonum >= 090101 & ubigeonum <= 090723, "Huancavelica",
ifelse(ubigeonum >= 100101 & ubigeonum <= 101108, "Huanuco",
ifelse(ubigeonum >= 110101 & ubigeonum <= 110508, "Ica",
ifelse(ubigeonum >= 120101 & ubigeonum <= 120909, "Junin",
ifelse(ubigeonum >= 130101 & ubigeonum <= 131203, "La Libertad",
ifelse(ubigeonum >= 140101 & ubigeonum <= 140312, "Lambayeque",
ifelse(ubigeonum >= 150101 & ubigeonum <= 150143, "Lima Metropolitana",
ifelse(ubigeonum >= 150201 & ubigeonum <= 151033, "Lima Region",
ifelse(ubigeonum >= 160101 & ubigeonum <= 160804, "Loreto",
ifelse(ubigeonum >= 170101 & ubigeonum <= 170303, "Madre de Dios",
ifelse(ubigeonum >= 180101 & ubigeonum <= 180303, "Moquegua",
ifelse(ubigeonum >= 190101 & ubigeonum <= 190308, "Pasco",
ifelse(ubigeonum >= 200101 & ubigeonum <= 200806, "Piura",
ifelse(ubigeonum >= 210101 & ubigeonum <= 211307, "Puno",
ifelse(ubigeonum >= 220101 & ubigeonum <= 221005, "San Martín",
ifelse(ubigeonum >= 230101 & ubigeonum <= 230408, "Tacna",
ifelse(ubigeonum >= 240101 & ubigeonum <= 240304, "Tumbes",
ifelse(ubigeonum >= 250101 & ubigeonum <= 250401,"Ucayali",NA)))))))))))))))))))))))))))
table(enaho$regiones2, useNA = "alw")##
## Amazonas Ancash Apurimac Arequipa
## 4303 4615 2901 4973
## Ayacucho Cajamarca Callao Cusco
## 3420 4584 3554 3838
## Huancavelica Huanuco Ica Junin
## 3350 4075 4932 4683
## La Libertad Lambayeque Lima Metropolitana Lima Region
## 5334 5112 10416 4710
## Loreto Madre de Dios Moquegua Pasco
## 5986 1980 2754 2887
## Piura Puno San Martín Tacna
## 6002 3265 4800 3972
## Tumbes Ucayali <NA>
## 2996 4536 0enaho <- enaho %>%
mutate(regnat = ifelse(dominio>=1 & dominio<=3 | dominio==8,"Costa",
ifelse(dominio>=4 & dominio<=6,"Sierra",
ifelse(dominio==7,"Selva",NA))))
table(enaho$regnat, useNA = "alw")##
## Costa Selva Sierra <NA>
## 49300 26356 38322 0enaho <- enaho %>%
mutate(area = ifelse((dominio==8 |
(dominio>=1 & dominio<=7) &
(estrato>=1 & estrato<=5)), "Urbano",
ifelse(((dominio>=1 & dominio<=7) &
(estrato>=6 & estrato<=8)), "Rural", NA)))
table(enaho$area, useNA = "alw")##
## Rural Urbano <NA>
## 38771 75207 0enaho <- enaho %>%
mutate(pobreza3 = ifelse(pobreza==1, "Pobre extremo",
ifelse(pobreza==2, "Pobre no extremo",
ifelse(pobreza==3, "No pobre", NA
))))
table(enaho$pobreza3, useNA = "alw")##
## No pobre Pobre extremo Pobre no extremo <NA>
## 84562 6653 22763 0enaho <- enaho %>%
mutate(lengua = ifelse(p300a==4, "Castellano",
ifelse(p300a==1 | p300a==2 | p300a==3, "Originaria", NA)))
enaho$lengua <- as.factor(enaho$lengua)
table(enaho$lengua, useNA = "alw")##
## Castellano Originaria <NA>
## 88412 19065 6501enaho <- enaho %>%
mutate(discapacidad =ifelse(p401h1==1|p401h2==1|p401h3==1|
p401h4==1|p401h5==1|p401h6==1,1,0))
table(enaho$discapacidad, useNA = "alw")##
## 0 1 <NA>
## 108237 5678 63enaho <- enaho %>%
mutate(defiet2 = case_when(
p558c == 1 ~ "Quechua",
p558c == 2 ~ "Aimara",
p558c == 3 ~ "Nativo o indigena de la Amazonia",
p558c == 4 ~ "Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente",
p558c == 5 ~ "Blanco",
p558c == 6 ~ "Mestizo",
p558c == 7 ~ "otro",
p558c == 8 ~ "No sabe/No responde",
p558c == 9 ~ "Nativo o indigena de la Amazonia",
TRUE ~ NA_character_
))
enaho$defiet2 <- as.factor(enaho$defiet2)
table(enaho$defiet2, useNA = "alw")##
## Aimara
## 2875
## Blanco
## 3683
## Mestizo
## 44363
## Nativo o indigena de la Amazonia
## 1985
## Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente
## 6007
## No sabe/No responde
## 3477
## otro
## 3610
## Quechua
## 21051
## <NA>
## 26927enaho <- enaho %>%
mutate(participacion = ifelse(p803>=1 & p803<=20,"Participa", ifelse(is.na(p803),"no participa","X")))
enaho$participacion <- replace(is.na(enaho$participacion),0,NA)
enaho$participacion <- factor(enaho$participacion, levels = c(FALSE, TRUE), labels = c("Participa", "No participa"))
table(enaho$participacion, useNA = "alw")##
## Participa No participa <NA>
## 18037 95941 0enaho <- enaho %>% mutate(Juventud = ifelse(p208a>=15 & p208a<=29,"jovenes","no jovenes"))
table(enaho$Juventud)##
## jovenes no jovenes
## 25136 88842table(enaho$participacion, enaho$Juventud)##
## jovenes no jovenes
## Participa 1929 16108
## No participa 23207 72734encuesta = svydesign(data=enaho, id=~conglome, strata=~estrato,
weights=~facpob07)tabla0 <- svyby(~participacion, ~Juventud, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla0ic0 <-as.table(confint(tabla0)) #INTERVALOS DE CONFIANZA
ic0## 2.5 % 97.5 %
## jovenes:participacionParticipa 0.05383798 0.06197437
## no jovenes:participacionParticipa 0.13897198 0.14812880
## jovenes:participacionNo participa 0.93802563 0.94616202
## no jovenes:participacionNo participa 0.85187120 0.86102802cv0 <-cv(tabla0) #COEFICIENTE DE VARIACIÓN
cv0workbook0 <- createWorkbook()
addWorksheet(workbook0, sheetName = "Tabla 0")
addWorksheet(workbook0, sheetName = "IC 0")
addWorksheet(workbook0, sheetName = "CV 0")
writeData(workbook0, sheet = "Tabla 0", x = tabla0, colNames = TRUE)
writeData(workbook0, sheet = "IC 0", x = ic0, colNames = TRUE)
writeData(workbook0, sheet = "CV 0", x = cv0, colNames = TRUE)
saveWorkbook(workbook0, "datos0.xlsx")tabla1 <- svyby(~participacion, ~Juventud+sexo, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla1ic1 <-as.table(confint(tabla1)) #INTERVALOS DE CONFIANZA
ic1## 2.5 % 97.5 %
## jovenes.hombre:participacionParticipa 0.03404464 0.04263607
## no jovenes.hombre:participacionParticipa 0.14790743 0.15946187
## jovenes.Mujer:participacionParticipa 0.07183936 0.08456453
## no jovenes.Mujer:participacionParticipa 0.12906000 0.13946167
## jovenes.hombre:participacionNo participa 0.95736393 0.96595536
## no jovenes.hombre:participacionNo participa 0.84053813 0.85209257
## jovenes.Mujer:participacionNo participa 0.91543547 0.92816064
## no jovenes.Mujer:participacionNo participa 0.86053833 0.87094000cv1 <-cv(tabla1) #COEFICIENTE DE VARIACIÓN
cv1workbook1 <- createWorkbook()
addWorksheet(workbook1, sheetName = "Tabla 1")
addWorksheet(workbook1, sheetName = "IC 1")
addWorksheet(workbook1, sheetName = "CV 1")
writeData(workbook1, sheet = "Tabla 1", x = tabla1, colNames = TRUE)
writeData(workbook1, sheet = "IC 1", x = ic1, colNames = TRUE)
writeData(workbook1, sheet = "CV 1", x = cv1, colNames = TRUE)
saveWorkbook(workbook1, "datos1.xlsx")tabla2 <- svyby(~participacion, ~Juventud+area, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla2ic2 <-as.table(confint(tabla2)) #INTERVALOS DE CONFIANZA
ic2## 2.5 % 97.5 %
## jovenes.Rural:participacionParticipa 0.13104103 0.15430804
## no jovenes.Rural:participacionParticipa 0.30520864 0.32498169
## jovenes.Urbano:participacionParticipa 0.03616037 0.04437153
## no jovenes.Urbano:participacionParticipa 0.09381284 0.10339408
## jovenes.Rural:participacionNo participa 0.84569196 0.86895897
## no jovenes.Rural:participacionNo participa 0.67501831 0.69479136
## jovenes.Urbano:participacionNo participa 0.95562847 0.96383963
## no jovenes.Urbano:participacionNo participa 0.89660592 0.90618716cv2 <-cv(tabla2) #COEFICIENTE DE VARIACIÓN
cv2workbook2 <- createWorkbook()
addWorksheet(workbook2, sheetName = "Tabla 2")
addWorksheet(workbook2, sheetName = "IC 2")
addWorksheet(workbook2, sheetName = "CV 2")
writeData(workbook2, sheet = "Tabla 2", x = tabla2, colNames = TRUE)
writeData(workbook2, sheet = "IC 2", x = ic2, colNames = TRUE)
writeData(workbook2, sheet = "CV 2", x = cv2, colNames = TRUE)
saveWorkbook(workbook2, "datos2.xlsx")tabla3 <- svyby(~participacion, ~Juventud+regnat, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla3ic3 <-as.table(confint(tabla3)) #INTERVALOS DE CONFIANZA
ic3## 2.5 % 97.5 %
## jovenes.Costa:participacionParticipa 0.02623198 0.03505729
## no jovenes.Costa:participacionParticipa 0.05925477 0.06797604
## jovenes.Selva:participacionParticipa 0.09050756 0.11287565
## no jovenes.Selva:participacionParticipa 0.16244850 0.18117348
## jovenes.Sierra:participacionParticipa 0.08396200 0.10256236
## no jovenes.Sierra:participacionParticipa 0.26596454 0.28502316
## jovenes.Costa:participacionNo participa 0.96494271 0.97376802
## no jovenes.Costa:participacionNo participa 0.93202396 0.94074523
## jovenes.Selva:participacionNo participa 0.88712435 0.90949244
## no jovenes.Selva:participacionNo participa 0.81882652 0.83755150
## jovenes.Sierra:participacionNo participa 0.89743764 0.91603800
## no jovenes.Sierra:participacionNo participa 0.71497684 0.73403546cv3 <-cv(tabla3) #COEFICIENTE DE VARIACIÓN
cv3workbook3 <- createWorkbook()
addWorksheet(workbook3, sheetName = "Tabla 3")
addWorksheet(workbook3, sheetName = "IC 3")
addWorksheet(workbook3, sheetName = "CV 3")
writeData(workbook3, sheet = "Tabla 3", x = tabla3, colNames = TRUE)
writeData(workbook3, sheet = "IC 3", x = ic3, colNames = TRUE)
writeData(workbook3, sheet = "CV 3", x = cv3, colNames = TRUE)
saveWorkbook(workbook3, "datos3.xlsx")tabla4 <- svyby(~participacion, ~Juventud+regiones2, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla4ic4 <-as.table(confint(tabla4)) #INTERVALOS DE CONFIANZA
ic4## 2.5 % 97.5 %
## jovenes.Amazonas:participacionParticipa 0.115806529 0.185795099
## no jovenes.Amazonas:participacionParticipa 0.258766557 0.314760060
## jovenes.Ancash:participacionParticipa 0.064418089 0.118417565
## no jovenes.Ancash:participacionParticipa 0.157035388 0.207374487
## jovenes.Apurimac:participacionParticipa 0.026590017 0.072416138
## no jovenes.Apurimac:participacionParticipa 0.285823961 0.342416956
## jovenes.Arequipa:participacionParticipa 0.043782157 0.086156468
## no jovenes.Arequipa:participacionParticipa 0.133053916 0.165329412
## jovenes.Ayacucho:participacionParticipa 0.100487721 0.160361685
## no jovenes.Ayacucho:participacionParticipa 0.305027390 0.353096605
## jovenes.Cajamarca:participacionParticipa 0.055581333 0.100441876
## no jovenes.Cajamarca:participacionParticipa 0.109404687 0.145840588
## jovenes.Callao:participacionParticipa 0.003175891 0.018664999
## no jovenes.Callao:participacionParticipa 0.020195656 0.045987645
## jovenes.Cusco:participacionParticipa 0.019069919 0.046526697
## no jovenes.Cusco:participacionParticipa 0.254283038 0.303355997
## jovenes.Huancavelica:participacionParticipa 0.102785682 0.162058597
## no jovenes.Huancavelica:participacionParticipa 0.392075482 0.465154073
## jovenes.Huanuco:participacionParticipa 0.157082038 0.245997322
## no jovenes.Huanuco:participacionParticipa 0.380760444 0.425681275
## jovenes.Ica:participacionParticipa 0.031358047 0.064977565
## no jovenes.Ica:participacionParticipa 0.062682015 0.091075332
## jovenes.Junin:participacionParticipa 0.053317339 0.095311921
## no jovenes.Junin:participacionParticipa 0.210771503 0.255325620
## jovenes.La Libertad:participacionParticipa 0.025516434 0.050858104
## no jovenes.La Libertad:participacionParticipa 0.078194276 0.104989861
## jovenes.Lambayeque:participacionParticipa 0.037712414 0.074893557
## no jovenes.Lambayeque:participacionParticipa 0.068867523 0.097208636
## jovenes.Lima Metropolitana:participacionParticipa 0.016387759 0.029950086
## no jovenes.Lima Metropolitana:participacionParticipa 0.044258775 0.057888046
## jovenes.Lima Region:participacionParticipa 0.010453189 0.039847635
## no jovenes.Lima Region:participacionParticipa 0.063543791 0.088974235
## jovenes.Loreto:participacionParticipa 0.087843419 0.134299819
## no jovenes.Loreto:participacionParticipa 0.140601165 0.170309635
## jovenes.Madre de Dios:participacionParticipa 0.036090335 0.093629517
## no jovenes.Madre de Dios:participacionParticipa 0.118219431 0.175611635
## jovenes.Moquegua:participacionParticipa 0.060515124 0.109223205
## no jovenes.Moquegua:participacionParticipa 0.248194087 0.310999410
## jovenes.Pasco:participacionParticipa 0.133616404 0.222550156
## no jovenes.Pasco:participacionParticipa 0.317174716 0.375065570
## jovenes.Piura:participacionParticipa 0.054366772 0.089293704
## no jovenes.Piura:participacionParticipa 0.112960789 0.147039505
## jovenes.Puno:participacionParticipa 0.059165693 0.114919893
## no jovenes.Puno:participacionParticipa 0.312232072 0.377370052
## jovenes.San Martín:participacionParticipa 0.045165265 0.081142651
## no jovenes.San Martín:participacionParticipa 0.079329675 0.106100168
## jovenes.Tacna:participacionParticipa 0.016182308 0.056958085
## no jovenes.Tacna:participacionParticipa 0.087846308 0.137116539
## jovenes.Tumbes:participacionParticipa 0.016538943 0.064837822
## no jovenes.Tumbes:participacionParticipa 0.065861358 0.099505666
## jovenes.Ucayali:participacionParticipa 0.068453161 0.117732995
## no jovenes.Ucayali:participacionParticipa 0.136032305 0.170647880
## jovenes.Amazonas:participacionNo participa 0.814204901 0.884193471
## no jovenes.Amazonas:participacionNo participa 0.685239940 0.741233443
## jovenes.Ancash:participacionNo participa 0.881582435 0.935581911
## no jovenes.Ancash:participacionNo participa 0.792625513 0.842964612
## jovenes.Apurimac:participacionNo participa 0.927583862 0.973409983
## no jovenes.Apurimac:participacionNo participa 0.657583044 0.714176039
## jovenes.Arequipa:participacionNo participa 0.913843532 0.956217843
## no jovenes.Arequipa:participacionNo participa 0.834670588 0.866946084
## jovenes.Ayacucho:participacionNo participa 0.839638315 0.899512279
## no jovenes.Ayacucho:participacionNo participa 0.646903395 0.694972610
## jovenes.Cajamarca:participacionNo participa 0.899558124 0.944418667
## no jovenes.Cajamarca:participacionNo participa 0.854159412 0.890595313
## jovenes.Callao:participacionNo participa 0.981335001 0.996824109
## no jovenes.Callao:participacionNo participa 0.954012355 0.979804344
## jovenes.Cusco:participacionNo participa 0.953473303 0.980930081
## no jovenes.Cusco:participacionNo participa 0.696644003 0.745716962
## jovenes.Huancavelica:participacionNo participa 0.837941403 0.897214318
## no jovenes.Huancavelica:participacionNo participa 0.534845927 0.607924518
## jovenes.Huanuco:participacionNo participa 0.754002678 0.842917962
## no jovenes.Huanuco:participacionNo participa 0.574318725 0.619239556
## jovenes.Ica:participacionNo participa 0.935022435 0.968641953
## no jovenes.Ica:participacionNo participa 0.908924668 0.937317985
## jovenes.Junin:participacionNo participa 0.904688079 0.946682661
## no jovenes.Junin:participacionNo participa 0.744674380 0.789228497
## jovenes.La Libertad:participacionNo participa 0.949141896 0.974483566
## no jovenes.La Libertad:participacionNo participa 0.895010139 0.921805724
## jovenes.Lambayeque:participacionNo participa 0.925106443 0.962287586
## no jovenes.Lambayeque:participacionNo participa 0.902791364 0.931132477
## jovenes.Lima Metropolitana:participacionNo participa 0.970049914 0.983612241
## no jovenes.Lima Metropolitana:participacionNo participa 0.942111954 0.955741225
## jovenes.Lima Region:participacionNo participa 0.960152365 0.989546811
## no jovenes.Lima Region:participacionNo participa 0.911025765 0.936456209
## jovenes.Loreto:participacionNo participa 0.865700181 0.912156581
## no jovenes.Loreto:participacionNo participa 0.829690365 0.859398835
## jovenes.Madre de Dios:participacionNo participa 0.906370483 0.963909665
## no jovenes.Madre de Dios:participacionNo participa 0.824388365 0.881780569
## jovenes.Moquegua:participacionNo participa 0.890776795 0.939484876
## no jovenes.Moquegua:participacionNo participa 0.689000590 0.751805913
## jovenes.Pasco:participacionNo participa 0.777449844 0.866383596
## no jovenes.Pasco:participacionNo participa 0.624934430 0.682825284
## jovenes.Piura:participacionNo participa 0.910706296 0.945633228
## no jovenes.Piura:participacionNo participa 0.852960495 0.887039211
## jovenes.Puno:participacionNo participa 0.885080107 0.940834307
## no jovenes.Puno:participacionNo participa 0.622629948 0.687767928
## jovenes.San Martín:participacionNo participa 0.918857349 0.954834735
## no jovenes.San Martín:participacionNo participa 0.893899832 0.920670325
## jovenes.Tacna:participacionNo participa 0.943041915 0.983817692
## no jovenes.Tacna:participacionNo participa 0.862883461 0.912153692
## jovenes.Tumbes:participacionNo participa 0.935162178 0.983461057
## no jovenes.Tumbes:participacionNo participa 0.900494334 0.934138642
## jovenes.Ucayali:participacionNo participa 0.882267005 0.931546839
## no jovenes.Ucayali:participacionNo participa 0.829352120 0.863967695cv4 <-cv(tabla4) #COEFICIENTE DE VARIACIÓN
cv4workbook4 <- createWorkbook()
addWorksheet(workbook4, sheetName = "Tabla 4")
addWorksheet(workbook4, sheetName = "IC 4")
addWorksheet(workbook4, sheetName = "CV 4")
writeData(workbook4, sheet = "Tabla 4", x = tabla4, colNames = TRUE)
writeData(workbook4, sheet = "IC 4", x = ic4, colNames = TRUE)
writeData(workbook4, sheet = "CV 4", x = cv4, colNames = TRUE)
saveWorkbook(workbook4, "datos4.xlsx")tabla5 <- svyby(~participacion, ~Juventud+pobreza3, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla5ic5 <-as.table(confint(tabla5)) #INTERVALOS DE CONFIANZA
ic5## 2.5 % 97.5 %
## jovenes.No pobre:participacionParticipa 0.04761141 0.05643316
## no jovenes.No pobre:participacionParticipa 0.13804698 0.14774270
## jovenes.Pobre extremo:participacionParticipa 0.11719984 0.16750239
## no jovenes.Pobre extremo:participacionParticipa 0.15985931 0.18943338
## jovenes.Pobre no extremo:participacionParticipa 0.05326823 0.07051882
## no jovenes.Pobre no extremo:participacionParticipa 0.13002624 0.14680546
## jovenes.No pobre:participacionNo participa 0.94356684 0.95238859
## no jovenes.No pobre:participacionNo participa 0.85225730 0.86195302
## jovenes.Pobre extremo:participacionNo participa 0.83249761 0.88280016
## no jovenes.Pobre extremo:participacionNo participa 0.81056662 0.84014069
## jovenes.Pobre no extremo:participacionNo participa 0.92948118 0.94673177
## no jovenes.Pobre no extremo:participacionNo participa 0.85319454 0.86997376cv5 <-cv(tabla5) #COEFICIENTE DE VARIACIÓN
cv5workbook5 <- createWorkbook()
addWorksheet(workbook5, sheetName = "Tabla 5")
addWorksheet(workbook5, sheetName = "IC 5")
addWorksheet(workbook5, sheetName = "CV 5")
writeData(workbook5, sheet = "Tabla 5", x = tabla5, colNames = TRUE)
writeData(workbook5, sheet = "IC 5", x = ic5, colNames = TRUE)
writeData(workbook5, sheet = "CV 5", x = cv5, colNames = TRUE)
saveWorkbook(workbook5, "datos5.xlsx")tabla6 <- svyby(~participacion, ~Juventud+discapacidad, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla6ic6 <-as.table(confint(tabla6)) #INTERVALOS DE CONFIANZA
ic6## 2.5 % 97.5 %
## jovenes.0:participacionParticipa 0.05357655 0.06182728
## no jovenes.0:participacionParticipa 0.13629386 0.14535305
## jovenes.1:participacionParticipa 0.03847253 0.09675796
## no jovenes.1:participacionParticipa 0.17545759 0.20882859
## jovenes.0:participacionNo participa 0.93817272 0.94642345
## no jovenes.0:participacionNo participa 0.85464695 0.86370614
## jovenes.1:participacionNo participa 0.90324204 0.96152747
## no jovenes.1:participacionNo participa 0.79117141 0.82454241cv6<-cv(tabla6) #COEFICIENTE DE VARIACIÓN
cv6workbook6 <- createWorkbook()
addWorksheet(workbook6, sheetName = "Tabla 6")
addWorksheet(workbook6, sheetName = "IC 6")
addWorksheet(workbook6, sheetName = "CV 6")
writeData(workbook6, sheet = "Tabla 6", x = tabla6, colNames = TRUE)
writeData(workbook6, sheet = "IC 6", x = ic6, colNames = TRUE)
writeData(workbook6, sheet = "CV 6", x = cv6, colNames = TRUE)
saveWorkbook(workbook6, "datos6.xlsx") tabla7 <- svyby(~participacion, ~Juventud+defiet2, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla7 ic7 <-as.table(confint(tabla7)) #INTERVALOS DE CONFIANZA
ic7 ## 2.5 %
## jovenes.Aimara:participacionParticipa 0.05973868
## no jovenes.Aimara:participacionParticipa 0.31061569
## jovenes.Blanco:participacionParticipa 0.04914724
## no jovenes.Blanco:participacionParticipa 0.11965139
## jovenes.Mestizo:participacionParticipa 0.03449326
## no jovenes.Mestizo:participacionParticipa 0.12263897
## jovenes.Nativo o indigena de la Amazonia:participacionParticipa 0.13914683
## no jovenes.Nativo o indigena de la Amazonia:participacionParticipa 0.35379767
## jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionParticipa 0.05228936
## no jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionParticipa 0.16955430
## jovenes.No sabe/No responde:participacionParticipa 0.02435586
## no jovenes.No sabe/No responde:participacionParticipa 0.10204877
## jovenes.otro:participacionParticipa 0.03926058
## no jovenes.otro:participacionParticipa 0.15688502
## jovenes.Quechua:participacionParticipa 0.08958439
## no jovenes.Quechua:participacionParticipa 0.34859801
## jovenes.Aimara:participacionNo participa 0.87124284
## no jovenes.Aimara:participacionNo participa 0.59319446
## jovenes.Blanco:participacionNo participa 0.91386123
## no jovenes.Blanco:participacionNo participa 0.84666012
## jovenes.Mestizo:participacionNo participa 0.95671602
## no jovenes.Mestizo:participacionNo participa 0.86449671
## jovenes.Nativo o indigena de la Amazonia:participacionNo participa 0.77729476
## no jovenes.Nativo o indigena de la Amazonia:participacionNo participa 0.54538318
## jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionNo participa 0.91864269
## no jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionNo participa 0.79689038
## jovenes.No sabe/No responde:participacionNo participa 0.94845120
## no jovenes.No sabe/No responde:participacionNo participa 0.85983847
## jovenes.otro:participacionNo participa 0.92495537
## no jovenes.otro:participacionNo participa 0.79870650
## jovenes.Quechua:participacionNo participa 0.88706529
## no jovenes.Quechua:participacionNo participa 0.62419012
## 97.5 %
## jovenes.Aimara:participacionParticipa 0.12875716
## no jovenes.Aimara:participacionParticipa 0.40680554
## jovenes.Blanco:participacionParticipa 0.08613877
## no jovenes.Blanco:participacionParticipa 0.15333988
## jovenes.Mestizo:participacionParticipa 0.04328398
## no jovenes.Mestizo:participacionParticipa 0.13550329
## jovenes.Nativo o indigena de la Amazonia:participacionParticipa 0.22270524
## no jovenes.Nativo o indigena de la Amazonia:participacionParticipa 0.45461682
## jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionParticipa 0.08135731
## no jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionParticipa 0.20310962
## jovenes.No sabe/No responde:participacionParticipa 0.05154880
## no jovenes.No sabe/No responde:participacionParticipa 0.14016153
## jovenes.otro:participacionParticipa 0.07504463
## no jovenes.otro:participacionParticipa 0.20129350
## jovenes.Quechua:participacionParticipa 0.11293471
## no jovenes.Quechua:participacionParticipa 0.37580988
## jovenes.Aimara:participacionNo participa 0.94026132
## no jovenes.Aimara:participacionNo participa 0.68938431
## jovenes.Blanco:participacionNo participa 0.95085276
## no jovenes.Blanco:participacionNo participa 0.88034861
## jovenes.Mestizo:participacionNo participa 0.96550674
## no jovenes.Mestizo:participacionNo participa 0.87736103
## jovenes.Nativo o indigena de la Amazonia:participacionNo participa 0.86085317
## no jovenes.Nativo o indigena de la Amazonia:participacionNo participa 0.64620233
## jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionNo participa 0.94771064
## no jovenes.Negro/Moreno/Zambo/Mulato/Pueblo Afro peruano o Afrodescendiente:participacionNo participa 0.83044570
## jovenes.No sabe/No responde:participacionNo participa 0.97564414
## no jovenes.No sabe/No responde:participacionNo participa 0.89795123
## jovenes.otro:participacionNo participa 0.96073942
## no jovenes.otro:participacionNo participa 0.84311498
## jovenes.Quechua:participacionNo participa 0.91041561
## no jovenes.Quechua:participacionNo participa 0.65140199cv7 <-cv(tabla7) #COEFICIENTE DE VARIACIÓN
cv7 workbook7 <- createWorkbook()
addWorksheet(workbook7, sheetName = "Tabla 7")
addWorksheet(workbook7, sheetName = "IC 7")
addWorksheet(workbook7, sheetName = "CV 7")
writeData(workbook7, sheet = "Tabla 7", x = tabla7, colNames = TRUE)
writeData(workbook7, sheet = "IC 7", x = ic7, colNames = TRUE)
writeData(workbook7, sheet = "CV 7", x = cv7, colNames = TRUE)
saveWorkbook(workbook7, "datos7.xlsx")tabla8 <- svyby(~participacion, ~Juventud+lengua, encuesta, svymean, deff=F,na.rm=T) #PROMEDIO
tabla8 ic8 <-as.table(confint(tabla8)) #INTERVALOS DE CONFIANZA
ic8 ## 2.5 % 97.5 %
## jovenes.Castellano:participacionParticipa 0.04356575 0.05138710
## no jovenes.Castellano:participacionParticipa 0.09745765 0.10566093
## jovenes.Originaria:participacionParticipa 0.13512360 0.17615460
## no jovenes.Originaria:participacionParticipa 0.35196310 0.37849637
## jovenes.Castellano:participacionNo participa 0.94861290 0.95643425
## no jovenes.Castellano:participacionNo participa 0.89433907 0.90254235
## jovenes.Originaria:participacionNo participa 0.82384540 0.86487640
## no jovenes.Originaria:participacionNo participa 0.62150363 0.64803690cv8 <-cv(tabla8) #COEFICIENTE DE VARIACIÓN
cv8 workbook8 <- createWorkbook()
addWorksheet(workbook8, sheetName = "Tabla 8")
addWorksheet(workbook8, sheetName = "IC 8")
addWorksheet(workbook8, sheetName = "CV 8")
writeData(workbook8, sheet = "Tabla 8", x = tabla8, colNames = TRUE)
writeData(workbook8, sheet = "IC 8", x = ic8, colNames = TRUE)
writeData(workbook8, sheet = "CV 8", x = cv8, colNames = TRUE)
saveWorkbook(workbook8, "datos8.xlsx")#save(enaho,file=paste(ruta,"BASEDEDATOSIndicador1ENDES.RData",sep = "/"))
#BORRAMOS TODO MENOS "RUTA"
#rm(list=setdiff(ls(), c("ruta")))
#VOLVEMOS A CARGAR NUESTRA BD LIMPIA
#load(paste(ruta,"BASEDEDATOSIndicador1ENDES.RData",sep="/"))