tests/quantiles-chile.R

In bschneidr/fastsurvey: Analysis of Complex Survey Samples

library(survey)
datos <- readRDS("datos_ejemplo.rds")

design <- svydesign(id = ~id_directorio, strata = ~estrato, weights = ~f_pers, check.strata = TRUE, data = datos)
set.seed(234262762)
repdesign <- as.svrepdesign(design, type = "subbootstrap", replicates=20)
options(survey.lonely.psu="remove")

values<-datos$ing_t_p[datos$CL_GRUPO_OCU_08=="ISCO08_6"]

suppressWarnings({
f0<-coef(svyquantile(~ing_t_p, subset(design,CL_GRUPO_OCU_08=="ISCO08_6"),quantiles=c(0.5), qrule="math"))
f0.5<-coef(svyquantile(~ing_t_p, subset(design,CL_GRUPO_OCU_08=="ISCO08_6"),quantiles=c(0.5), qrule="school"))
})
all.equal(c(values[1],mean(values)), as.vector(c(f0,f0.5)))

suppressWarnings({
f0<-coef(svyquantile(~ing_t_p, subset(repdesign,CL_GRUPO_OCU_08=="ISCO08_6"),quantiles=c(0.5), qrule="math"))
f0.5<-coef(svyquantile(~ing_t_p, subset(repdesign,CL_GRUPO_OCU_08=="ISCO08_6"),quantiles=c(0.5), qrule="school"))
})
all.equal(c(values[1],mean(values)), as.vector(c(f0,f0.5)))