1. Lectura del marco.
library(survey)
## Warning: package 'survey' was built under R version 4.5.2
## Loading required package: grid
## Loading required package: Matrix
## Loading required package: survival
##
## Attaching package: 'survey'
## The following object is masked from 'package:graphics':
##
## dotchart
library(sampling)
##
## Attaching package: 'sampling'
## The following objects are masked from 'package:survival':
##
## cluster, strata
base <- readRDS("/Users/usermac/Documents/especializacion/muestreo estadistico/practicas_en_r/Copia de Marco.rds")
#base<-readRDS(file = "Marco.rds")
str(base)
## 'data.frame': 546212 obs. of 82 variables:
## $ ESTU_TIPODOCUMENTO : chr "TI" "TI" "TI" "TI" ...
## $ ESTU_NACIONALIDAD : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_GENERO : chr "M" "M" "M" "M" ...
## $ ESTU_FECHANACIMIENTO : chr "02/18/2002 12:00:00 AM" "06/28/2003 12:00:00 AM" "06/02/2001 12:00:00 AM" "04/26/2002 12:00:00 AM" ...
## $ PERIODO : int 20194 20194 20194 20194 20194 20194 20194 20194 20194 20194 ...
## $ ESTU_CONSECUTIVO : chr "SB11201940031558" "SB11201940303339" "SB11201940255017" "SB11201940154799" ...
## $ ESTU_ESTUDIANTE : chr "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" ...
## $ ESTU_TIENEETNIA : chr "No" "No" "No" "Si" ...
## $ ESTU_PAIS_RESIDE : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_ETNIA : chr "-" "-" "-" "Ninguno" ...
## $ ESTU_DEPTO_RESIDE : chr "MAGDALENA" "BOGOTÁ" "BOLIVAR" "BOGOTÁ" ...
## $ ESTU_COD_RESIDE_DEPTO : chr "47" "11" "13" "11" ...
## $ ESTU_MCPIO_RESIDE : chr "SANTA ANA" "BOGOTÁ D.C." "CARTAGENA DE INDIAS" "BOGOTÁ D.C." ...
## $ ESTU_COD_RESIDE_MCPIO : chr "47707" "11001" "13001" "11001" ...
## $ FAMI_ESTRATOVIVIENDA : chr "Estrato 3" "Estrato 3" "Estrato 1" "Estrato 3" ...
## $ FAMI_PERSONASHOGAR : chr "5 a 6" "9 o más" "5 a 6" "3 a 4" ...
## $ FAMI_CUARTOSHOGAR : chr "Cinco" "Tres" "Dos" "Dos" ...
## $ FAMI_EDUCACIONPADRE : chr "Educación profesional completa" "Técnica o tecnológica completa" "Secundaria (Bachillerato) completa" "Primaria incompleta" ...
## $ FAMI_EDUCACIONMADRE : chr "Secundaria (Bachillerato) completa" "Técnica o tecnológica completa" "Secundaria (Bachillerato) completa" "Secundaria (Bachillerato) completa" ...
## $ FAMI_TRABAJOLABORPADRE : chr "No aplica" "Trabaja por cuenta propia (por ejemplo plomero, electricista)" "Es vendedor o trabaja en atención al público" "No aplica" ...
## $ FAMI_TRABAJOLABORMADRE : chr "Trabaja en el hogar, no trabaja o estudia" "Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc" "Es vendedor o trabaja en atención al público" "Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc" ...
## $ FAMI_TIENEINTERNET : chr "Si" "Si" "No" "No" ...
## $ FAMI_TIENESERVICIOTV : chr "Si" "Si" "Si" "No" ...
## $ FAMI_TIENECOMPUTADOR : chr "No" "Si" "No" "Si" ...
## $ FAMI_TIENELAVADORA : chr "Si" "Si" "Si" "No" ...
## $ FAMI_TIENEHORNOMICROOGAS : chr "Si" "Si" "No" "No" ...
## $ FAMI_TIENEAUTOMOVIL : chr "No" "Si" "No" "No" ...
## $ FAMI_TIENEMOTOCICLETA : chr "No" "No" "No" "No" ...
## $ FAMI_TIENECONSOLAVIDEOJUEGOS : chr "No" "No" "No" "No" ...
## $ FAMI_NUMLIBROS : chr "MÁS DE 100 LIBROS" "MÁS DE 100 LIBROS" "26 A 100 LIBROS" "11 A 25 LIBROS" ...
## $ FAMI_COMELECHEDERIVADOS : chr "1 o 2 veces por semana" "Todos o casi todos los días" "Todos o casi todos los días" "Todos o casi todos los días" ...
## $ FAMI_COMECARNEPESCADOHUEVO : chr "3 a 5 veces por semana" "3 a 5 veces por semana" "Todos o casi todos los días" "Todos o casi todos los días" ...
## $ FAMI_COMECEREALFRUTOSLEGUMBRE: chr "Nunca o rara vez comemos eso" "1 o 2 veces por semana" "3 a 5 veces por semana" "3 a 5 veces por semana" ...
## $ FAMI_SITUACIONECONOMICA : chr "Peor" "Mejor" "Igual" "Igual" ...
## $ ESTU_DEDICACIONLECTURADIARIA : chr "Entre 30 y 60 minutos" "Entre 30 y 60 minutos" "Entre 30 y 60 minutos" "30 minutos o menos" ...
## $ ESTU_DEDICACIONINTERNET : chr "Entre 30 y 60 minutos" "Entre 30 y 60 minutos" "Más de 3 horas" "Entre 30 y 60 minutos" ...
## $ ESTU_HORASSEMANATRABAJA : chr "Menos de 10 horas" "Menos de 10 horas" "0" "Más de 30 horas" ...
## $ ESTU_TIPOREMUNERACION : chr "Si, en efectivo" "Si, en efectivo" "Si, en efectivo" "Si, en efectivo" ...
## $ COLE_CODIGO_ICFES : int 39008 128611 35493 138032 217406 706101 16659 709261 9233 179630 ...
## $ COLE_COD_DANE_ESTABLECIMIENTO: num 1.48e+11 3.11e+11 1.13e+11 1.11e+11 1.11e+11 ...
## $ COLE_NOMBRE_ESTABLECIMIENTO : chr "INSTITUCION EDUCATIVA DEPARTAMENTAL ANTONIO BRUJES CARMONA" "GIMN SAN ANGELO" "INSTITUCION EDUCATIVA PROMOCION SOCIAL DE C/GENA." "COLEGIO VIRGINIA GUTIERREZ DE PINEDA (IED)" ...
## $ COLE_GENERO : chr "MIXTO" "MIXTO" "MIXTO" "MIXTO" ...
## $ COLE_NATURALEZA : chr "OFICIAL" "NO OFICIAL" "OFICIAL" "OFICIAL" ...
## $ COLE_CALENDARIO : chr "A" "A" "A" "A" ...
## $ COLE_BILINGUE : chr "N" "-" "N" "-" ...
## $ COLE_CARACTER : chr "TÉCNICO/ACADÉMICO" "ACADÉMICO" "TÉCNICO/ACADÉMICO" "ACADÉMICO" ...
## $ COLE_COD_DANE_SEDE : num 1.48e+11 3.11e+11 1.13e+11 1.11e+11 1.11e+11 ...
## $ COLE_NOMBRE_SEDE : chr "COL DPTAL ANTONIO BRUGES CARMONA" "GIMN SAN ANGELO" "INSTITUCION EDUCATIVA PROMOCION SOCIAL DE C/GENA." "COL VIRGINIA GUTIERREZ DE PINEDA (INS EDUC DIST)" ...
## $ COLE_SEDE_PRINCIPAL : chr "S" "S" "S" "S" ...
## $ COLE_AREA_UBICACION : chr "URBANO" "URBANO" "URBANO" "URBANO" ...
## $ COLE_JORNADA : chr "MAÑANA" "COMPLETA" "MAÑANA" "MAÑANA" ...
## $ COLE_COD_MCPIO_UBICACION : int 47707 11001 13001 11001 11001 8758 76001 68001 25754 70001 ...
## $ COLE_MCPIO_UBICACION : chr "SANTA ANA" "BOGOTÁ D.C." "CARTAGENA DE INDIAS" "BOGOTÁ D.C." ...
## $ COLE_COD_DEPTO_UBICACION : int 47 11 13 11 11 8 76 68 25 70 ...
## $ COLE_DEPTO_UBICACION : chr "MAGDALENA" "BOGOTÁ" "BOLIVAR" "BOGOTÁ" ...
## $ ESTU_PRIVADO_LIBERTAD : chr "N" "N" "N" "N" ...
## $ ESTU_COD_MCPIO_PRESENTACION : int 47707 11001 13001 11001 11001 8758 76001 68001 25754 70001 ...
## $ ESTU_MCPIO_PRESENTACION : chr "SANTA ANA" "BOGOTÁ D.C." "CARTAGENA DE INDIAS" "BOGOTÁ D.C." ...
## $ ESTU_DEPTO_PRESENTACION : chr "MAGDALENA" "BOGOTÁ" "BOLIVAR" "BOGOTÁ" ...
## $ ESTU_COD_DEPTO_PRESENTACION : int 47 11 13 11 11 8 76 68 25 70 ...
## $ PUNT_LECTURA_CRITICA : int 47 60 66 62 63 49 76 57 62 68 ...
## $ PERCENTIL_LECTURA_CRITICA : int 33 76 91 81 85 37 100 67 83 94 ...
## $ DESEMP_LECTURA_CRITICA : int 2 3 4 3 3 2 4 3 3 4 ...
## $ PUNT_MATEMATICAS : int 48 65 57 54 57 29 70 65 62 66 ...
## $ PERCENTIL_MATEMATICAS : int 42 88 70 60 70 4 96 88 82 91 ...
## $ DESEMP_MATEMATICAS : int 2 3 3 3 3 1 3 3 3 3 ...
## $ PUNT_C_NATURALES : int 37 54 41 61 55 41 70 63 66 63 ...
## $ PERCENTIL_C_NATURALES : int 15 71 29 87 74 26 98 90 95 91 ...
## $ DESEMP_C_NATURALES : int 1 2 2 3 2 2 3 3 3 3 ...
## $ PUNT_SOCIALES_CIUDADANAS : int 30 59 74 73 57 41 68 66 39 77 ...
## $ PERCENTIL_SOCIALES_CIUDADANAS: int 8 83 99 99 79 39 96 95 33 100 ...
## $ DESEMP_SOCIALES_CIUDADANAS : int 1 3 4 4 3 2 3 3 1 4 ...
## $ PUNT_INGLES : int 54 63 64 53 52 35 72 60 63 51 ...
## $ PERCENTIL_INGLES : int 67 88 89 67 63 14 96 83 88 61 ...
## $ DESEMP_INGLES : chr "A1" "A2" "A2" "A1" ...
## $ PUNT_GLOBAL : int 208 299 299 309 288 198 355 313 288 336 ...
## $ PERCENTIL_GLOBAL : chr "25" "83" "84" "88" ...
## $ ESTU_INSE_INDIVIDUAL : chr "50.9421557404314" "59.767864431576" "50.7642383644749" "45.7791353123516" ...
## $ ESTU_NSE_INDIVIDUAL : chr "2" "3" "2" "2" ...
## $ ESTU_NSE_ESTABLECIMIENTO : int 2 4 3 3 3 2 3 3 2 3 ...
## $ ESTU_ESTADOINVESTIGACION : chr "PUBLICAR" "PUBLICAR" "PUBLICAR" "PUBLICAR" ...
## $ ESTU_GENERACION.E : chr "NO" "NO" "GENERACION E - GRATUIDAD" "GENERACION E - GRATUIDAD" ...
2. Base de Funza.
funza <- base[base$COLE_COD_MCPIO_UBICACION == 25286,]
head(funza)
## ESTU_TIPODOCUMENTO ESTU_NACIONALIDAD ESTU_GENERO ESTU_FECHANACIMIENTO
## 41 PE VENEZUELA M 03/03/2002 12:00:00 AM
## 1937 TI COLOMBIA F 04/28/2001 12:00:00 AM
## 2943 TI COLOMBIA F 01/18/2002 12:00:00 AM
## 3223 CC COLOMBIA F 09/25/2000 12:00:00 AM
## 4267 TI COLOMBIA M 03/27/2003 12:00:00 AM
## 4582 TI COLOMBIA F 05/02/2002 12:00:00 AM
## PERIODO ESTU_CONSECUTIVO ESTU_ESTUDIANTE ESTU_TIENEETNIA ESTU_PAIS_RESIDE
## 41 20194 SB11201940498839 ESTUDIANTE No VENEZUELA
## 1937 20194 SB11201940504543 ESTUDIANTE No COLOMBIA
## 2943 20194 SB11201940357444 ESTUDIANTE No COLOMBIA
## 3223 20194 SB11201940494987 ESTUDIANTE No COLOMBIA
## 4267 20194 SB11201940039213 ESTUDIANTE No COLOMBIA
## 4582 20194 SB11201940040952 ESTUDIANTE No COLOMBIA
## ESTU_ETNIA ESTU_DEPTO_RESIDE ESTU_COD_RESIDE_DEPTO ESTU_MCPIO_RESIDE
## 41 - CUNDINAMARCA 25 FUNZA
## 1937 - CUNDINAMARCA 25 FUNZA
## 2943 - CUNDINAMARCA 25 FUNZA
## 3223 - CUNDINAMARCA 25 FUNZA
## 4267 - CUNDINAMARCA 25 FUNZA
## 4582 - CUNDINAMARCA 25 FUNZA
## ESTU_COD_RESIDE_MCPIO FAMI_ESTRATOVIVIENDA FAMI_PERSONASHOGAR
## 41 25286 Estrato 3 3 a 4
## 1937 25286 Estrato 2 3 a 4
## 2943 25286 Estrato 2 3 a 4
## 3223 25286 Estrato 3 3 a 4
## 4267 25286 Estrato 2 7 a 8
## 4582 25286 Estrato 2 3 a 4
## FAMI_CUARTOSHOGAR FAMI_EDUCACIONPADRE
## 41 Cuatro Secundaria (Bachillerato) incompleta
## 1937 Dos Secundaria (Bachillerato) incompleta
## 2943 Cuatro Primaria completa
## 3223 Dos Técnica o tecnológica completa
## 4267 Cinco Educación profesional completa
## 4582 Tres Primaria incompleta
## FAMI_EDUCACIONMADRE
## 41 Técnica o tecnológica completa
## 1937 Secundaria (Bachillerato) incompleta
## 2943 Primaria completa
## 3223 Secundaria (Bachillerato) completa
## 4267 Postgrado
## 4582 Primaria incompleta
## FAMI_TRABAJOLABORPADRE
## 41 No sabe
## 1937 Es operario de máquinas o conduce vehículos (taxita, chofer)
## 2943 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## 3223 No sabe
## 4267 Es dueño de un negocio grande, tiene un cargo de nivel directivo o gerencial
## 4582 Es agricultor, pesquero o jornalero
## FAMI_TRABAJOLABORMADRE
## 41 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## 1937 Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc
## 2943 Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc
## 3223 Trabaja en el hogar, no trabaja o estudia
## 4267 Trabaja como profesional (por ejemplo médico, abogado, ingeniero)
## 4582 Trabaja en el hogar, no trabaja o estudia
## FAMI_TIENEINTERNET FAMI_TIENESERVICIOTV FAMI_TIENECOMPUTADOR
## 41 No No No
## 1937 Si Si Si
## 2943 Si No Si
## 3223 Si Si Si
## 4267 Si Si Si
## 4582 No Si No
## FAMI_TIENELAVADORA FAMI_TIENEHORNOMICROOGAS FAMI_TIENEAUTOMOVIL
## 41 No No No
## 1937 Si No Si
## 2943 Si Si No
## 3223 Si Si No
## 4267 Si Si Si
## 4582 Si No No
## FAMI_TIENEMOTOCICLETA FAMI_TIENECONSOLAVIDEOJUEGOS FAMI_NUMLIBROS
## 41 No No -
## 1937 No No 0 A 10 LIBROS
## 2943 Si No MÁS DE 100 LIBROS
## 3223 No Si 26 A 100 LIBROS
## 4267 Si Si 26 A 100 LIBROS
## 4582 Si No 11 A 25 LIBROS
## FAMI_COMELECHEDERIVADOS FAMI_COMECARNEPESCADOHUEVO
## 41 3 a 5 veces por semana 3 a 5 veces por semana
## 1937 Todos o casi todos los días Todos o casi todos los días
## 2943 Todos o casi todos los días 1 o 2 veces por semana
## 3223 Nunca o rara vez comemos eso Todos o casi todos los días
## 4267 Todos o casi todos los días Todos o casi todos los días
## 4582 1 o 2 veces por semana 1 o 2 veces por semana
## FAMI_COMECEREALFRUTOSLEGUMBRE FAMI_SITUACIONECONOMICA
## 41 1 o 2 veces por semana Igual
## 1937 1 o 2 veces por semana Mejor
## 2943 Todos o casi todos los días Mejor
## 3223 1 o 2 veces por semana Mejor
## 4267 3 a 5 veces por semana Mejor
## 4582 1 o 2 veces por semana Igual
## ESTU_DEDICACIONLECTURADIARIA ESTU_DEDICACIONINTERNET
## 41 30 minutos o menos Entre 30 y 60 minutos
## 1937 30 minutos o menos Entre 30 y 60 minutos
## 2943 Entre 1 y 2 horas 30 minutos o menos
## 3223 30 minutos o menos Entre 30 y 60 minutos
## 4267 Entre 30 y 60 minutos Entre 1 y 3 horas
## 4582 Entre 30 y 60 minutos Entre 30 y 60 minutos
## ESTU_HORASSEMANATRABAJA ESTU_TIPOREMUNERACION COLE_CODIGO_ICFES
## 41 Entre 21 y 30 horas No 234781
## 1937 Menos de 10 horas Si, en efectivo 141887
## 2943 Menos de 10 horas No 173872
## 3223 Menos de 10 horas Si, en efectivo 234781
## 4267 0 Si, en efectivo 150599
## 4582 0 No 717744
## COLE_COD_DANE_ESTABLECIMIENTO COLE_NOMBRE_ESTABLECIMIENTO COLE_GENERO
## 41 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## 1937 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## 2943 3.25286e+11 FUNDACIÓN COLOMBIA APRENDE MIXTO
## 3223 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## 4267 3.25286e+11 COLEGIO INGLES NUEVA INGLATERRA MIXTO
## 4582 1.25286e+11 I.E.D. FUNZA MIXTO
## COLE_NATURALEZA COLE_CALENDARIO COLE_BILINGUE COLE_CARACTER
## 41 OFICIAL A N ACADÉMICO
## 1937 OFICIAL A N ACADÉMICO
## 2943 NO OFICIAL A N ACADÉMICO
## 3223 OFICIAL A N ACADÉMICO
## 4267 NO OFICIAL A - TÉCNICO/ACADÉMICO
## 4582 OFICIAL A N ACADÉMICO
## COLE_COD_DANE_SEDE COLE_NOMBRE_SEDE
## 41 125286000123 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## 1937 125286000123 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## 2943 325286001161 FUNDACIÓN COLOMBIA APRENDE - SEDE PRINCIPAL
## 3223 125286000123 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## 4267 325286000840 COL INGLES NUEVA INGLATERRA
## 4582 125286800006 FURATENA
## COLE_SEDE_PRINCIPAL COLE_AREA_UBICACION COLE_JORNADA
## 41 N URBANO MAÑANA
## 1937 N URBANO NOCHE
## 2943 S URBANO NOCHE
## 3223 N URBANO MAÑANA
## 4267 S URBANO COMPLETA
## 4582 N URBANO MAÑANA
## COLE_COD_MCPIO_UBICACION COLE_MCPIO_UBICACION COLE_COD_DEPTO_UBICACION
## 41 25286 FUNZA 25
## 1937 25286 FUNZA 25
## 2943 25286 FUNZA 25
## 3223 25286 FUNZA 25
## 4267 25286 FUNZA 25
## 4582 25286 FUNZA 25
## COLE_DEPTO_UBICACION ESTU_PRIVADO_LIBERTAD ESTU_COD_MCPIO_PRESENTACION
## 41 CUNDINAMARCA N 25286
## 1937 CUNDINAMARCA N 25286
## 2943 CUNDINAMARCA N 25286
## 3223 CUNDINAMARCA N 25286
## 4267 CUNDINAMARCA N 25286
## 4582 CUNDINAMARCA N 25286
## ESTU_MCPIO_PRESENTACION ESTU_DEPTO_PRESENTACION
## 41 FUNZA CUNDINAMARCA
## 1937 FUNZA CUNDINAMARCA
## 2943 FUNZA CUNDINAMARCA
## 3223 FUNZA CUNDINAMARCA
## 4267 FUNZA CUNDINAMARCA
## 4582 FUNZA CUNDINAMARCA
## ESTU_COD_DEPTO_PRESENTACION PUNT_LECTURA_CRITICA PERCENTIL_LECTURA_CRITICA
## 41 25 57 67
## 1937 25 37 9
## 2943 25 35 6
## 3223 25 58 69
## 4267 25 53 52
## 4582 25 72 98
## DESEMP_LECTURA_CRITICA PUNT_MATEMATICAS PERCENTIL_MATEMATICAS
## 41 3 67 92
## 1937 2 37 14
## 2943 1 32 6
## 3223 3 41 22
## 4267 3 68 92
## 4582 4 60 77
## DESEMP_MATEMATICAS PUNT_C_NATURALES PERCENTIL_C_NATURALES
## 41 3 66 94
## 1937 2 40 24
## 2943 1 35 12
## 3223 2 52 65
## 4267 3 63 91
## 4582 3 62 89
## DESEMP_C_NATURALES PUNT_SOCIALES_CIUDADANAS PERCENTIL_SOCIALES_CIUDADANAS
## 41 3 52 67
## 1937 1 28 4
## 2943 1 27 4
## 3223 2 49 60
## 4267 3 60 86
## 4582 3 69 97
## DESEMP_SOCIALES_CIUDADANAS PUNT_INGLES PERCENTIL_INGLES DESEMP_INGLES
## 41 2 67 93 A2
## 1937 1 44 39 A-
## 2943 1 33 11 A-
## 3223 2 54 68 A1
## 4267 3 52 61 A1
## 4582 3 40 28 A-
## PUNT_GLOBAL PERCENTIL_GLOBAL ESTU_INSE_INDIVIDUAL ESTU_NSE_INDIVIDUAL
## 41 305 86 42.9347982318362 2
## 1937 181 10 52.9804696713036 3
## 2943 162 4 53.3997281953332 3
## 3223 252 56 56.3996702122112 3
## 4267 302 85 74.4007858105476 4
## 4582 319 91 38.1705705683999 1
## ESTU_NSE_ESTABLECIMIENTO ESTU_ESTADOINVESTIGACION ESTU_GENERACION.E
## 41 2 PUBLICAR GENERACION E - GRATUIDAD
## 1937 2 PUBLICAR GENERACION E - GRATUIDAD
## 2943 3 PUBLICAR GENERACION E - GRATUIDAD
## 3223 2 PUBLICAR GENERACION E - GRATUIDAD
## 4267 3 PUBLICAR NO
## 4582 3 PUBLICAR NO
nrow(funza)
## [1] 1489
funza$academico=ifelse(funza$COLE_CARACTER=="ACADÉMICO",1,0)
attach(funza)
4. Selección de la muestra.
Coordinado negativo.
s.mas=function(base, n, seed){
N=nrow(base)
set.seed(seed)
base$u = runif(nrow(base))
base=base[with(base, order(base$u)),]
base=base[1:n, ]
base$pik = n/N
return(base)
}
n=428
seed=123
muestra = s.mas(funza, n, seed)
str(muestra)
## 'data.frame': 428 obs. of 85 variables:
## $ ESTU_TIPODOCUMENTO : chr "CC" "TI" "TI" "TI" ...
## $ ESTU_NACIONALIDAD : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_GENERO : chr "M" "M" "F" "M" ...
## $ ESTU_FECHANACIMIENTO : chr "10/19/2000 12:00:00 AM" "08/01/2002 12:00:00 AM" "07/24/2002 12:00:00 AM" "05/06/2002 12:00:00 AM" ...
## $ PERIODO : int 20194 20194 20194 20194 20194 20194 20194 20194 20194 20194 ...
## $ ESTU_CONSECUTIVO : chr "SB11201940505953" "SB11201940387065" "SB11201940156899" "SB11201940497663" ...
## $ ESTU_ESTUDIANTE : chr "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" ...
## $ ESTU_TIENEETNIA : chr "No" "No" "No" "No" ...
## $ ESTU_PAIS_RESIDE : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_ETNIA : chr "-" "-" "-" "-" ...
## $ ESTU_DEPTO_RESIDE : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_COD_RESIDE_DEPTO : chr "25" "25" "25" "25" ...
## $ ESTU_MCPIO_RESIDE : chr "FUNZA" "MOSQUERA" "FUNZA" "FUNZA" ...
## $ ESTU_COD_RESIDE_MCPIO : chr "25286" "25473" "25286" "25286" ...
## $ FAMI_ESTRATOVIVIENDA : chr "-" "Estrato 3" "Estrato 2" "Estrato 2" ...
## $ FAMI_PERSONASHOGAR : chr "" "3 a 4" "3 a 4" "9 o más" ...
## $ FAMI_CUARTOSHOGAR : chr "-" "Tres" "Dos" "Seis o mas" ...
## $ FAMI_EDUCACIONPADRE : chr "" "Educación profesional incompleta" "Secundaria (Bachillerato) completa" "Técnica o tecnológica completa" ...
## $ FAMI_EDUCACIONMADRE : chr "-" "Educación profesional completa" "Secundaria (Bachillerato) completa" "Técnica o tecnológica completa" ...
## $ FAMI_TRABAJOLABORPADRE : chr "" "Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)" "Es operario de máquinas o conduce vehículos (taxita, chofer)" "Trabaja como profesional (por ejemplo médico, abogado, ingeniero)" ...
## $ FAMI_TRABAJOLABORMADRE : chr "-" "Es dueño de un negocio grande, tiene un cargo de nivel directivo o gerencial" "Trabaja en el hogar, no trabaja o estudia" "Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)" ...
## $ FAMI_TIENEINTERNET : chr "" "Si" "Si" "Si" ...
## $ FAMI_TIENESERVICIOTV : chr "-" "Si" "Si" "Si" ...
## $ FAMI_TIENECOMPUTADOR : chr "" "Si" "Si" "Si" ...
## $ FAMI_TIENELAVADORA : chr "-" "Si" "Si" "Si" ...
## $ FAMI_TIENEHORNOMICROOGAS : chr "" "Si" "Si" "Si" ...
## $ FAMI_TIENEAUTOMOVIL : chr "-" "Si" "No" "Si" ...
## $ FAMI_TIENEMOTOCICLETA : chr "" "No" "No" "Si" ...
## $ FAMI_TIENECONSOLAVIDEOJUEGOS : chr "-" "Si" "No" "No" ...
## $ FAMI_NUMLIBROS : chr "" "11 A 25 LIBROS" "0 A 10 LIBROS" "26 A 100 LIBROS" ...
## $ FAMI_COMELECHEDERIVADOS : chr "-" "Todos o casi todos los días" "Todos o casi todos los días" "Todos o casi todos los días" ...
## $ FAMI_COMECARNEPESCADOHUEVO : chr "" "Todos o casi todos los días" "Todos o casi todos los días" "Todos o casi todos los días" ...
## $ FAMI_COMECEREALFRUTOSLEGUMBRE: chr "-" "3 a 5 veces por semana" "1 o 2 veces por semana" "Todos o casi todos los días" ...
## $ FAMI_SITUACIONECONOMICA : chr "" "Igual" "Igual" "Igual" ...
## $ ESTU_DEDICACIONLECTURADIARIA : chr "-" "No leo por entretenimiento" "30 minutos o menos" "30 minutos o menos" ...
## $ ESTU_DEDICACIONINTERNET : chr "" "Entre 30 y 60 minutos" "Entre 30 y 60 minutos" "Entre 1 y 3 horas" ...
## $ ESTU_HORASSEMANATRABAJA : chr "-" "0" "0" "0" ...
## $ ESTU_TIPOREMUNERACION : chr "" "No" "No" "No" ...
## $ COLE_CODIGO_ICFES : int 141887 122200 667253 234781 44370 41327 150599 125682 44370 44370 ...
## $ COLE_COD_DANE_ESTABLECIMIENTO: num 1.25e+11 3.25e+11 1.25e+11 1.25e+11 1.25e+11 ...
## $ COLE_NOMBRE_ESTABLECIMIENTO : chr "I.E.D. MIGUEL ANTONIO CARO" "COL CRISTIANO INTEGRAL" "I.E.D. BICENTENARIO" "I.E.D. MIGUEL ANTONIO CARO" ...
## $ COLE_GENERO : chr "MIXTO" "MIXTO" "MIXTO" "MIXTO" ...
## $ COLE_NATURALEZA : chr "OFICIAL" "NO OFICIAL" "OFICIAL" "OFICIAL" ...
## $ COLE_CALENDARIO : chr "A" "A" "A" "A" ...
## $ COLE_BILINGUE : chr "N" "N" "N" "N" ...
## $ COLE_CARACTER : chr "ACADÉMICO" "ACADÉMICO" "TÉCNICO/ACADÉMICO" "ACADÉMICO" ...
## $ COLE_COD_DANE_SEDE : num 1.25e+11 3.25e+11 1.25e+11 1.25e+11 1.25e+11 ...
## $ COLE_NOMBRE_SEDE : chr "CONCENTRACION URBANA MIGUEL ANTONIO CARO" "COL CRISTIANO INTEGRAL" "I.E.D. BICENTENARIO - SEDE PRINCIPAL" "CONCENTRACION URBANA MIGUEL ANTONIO CARO" ...
## $ COLE_SEDE_PRINCIPAL : chr "N" "S" "S" "N" ...
## $ COLE_AREA_UBICACION : chr "URBANO" "URBANO" "URBANO" "URBANO" ...
## $ COLE_JORNADA : chr "NOCHE" "COMPLETA" "UNICA" "MAÑANA" ...
## $ COLE_COD_MCPIO_UBICACION : int 25286 25286 25286 25286 25286 25286 25286 25286 25286 25286 ...
## $ COLE_MCPIO_UBICACION : chr "FUNZA" "FUNZA" "FUNZA" "FUNZA" ...
## $ COLE_COD_DEPTO_UBICACION : int 25 25 25 25 25 25 25 25 25 25 ...
## $ COLE_DEPTO_UBICACION : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_PRIVADO_LIBERTAD : chr "N" "N" "N" "N" ...
## $ ESTU_COD_MCPIO_PRESENTACION : int 25286 25473 25286 25286 25286 25286 25286 25286 25286 25286 ...
## $ ESTU_MCPIO_PRESENTACION : chr "FUNZA" "MOSQUERA" "FUNZA" "FUNZA" ...
## $ ESTU_DEPTO_PRESENTACION : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_COD_DEPTO_PRESENTACION : int 25 25 25 25 25 25 25 25 25 25 ...
## $ PUNT_LECTURA_CRITICA : int 59 70 39 65 43 52 68 47 55 53 ...
## $ PERCENTIL_LECTURA_CRITICA : int 73 96 12 88 20 50 93 33 60 52 ...
## $ DESEMP_LECTURA_CRITICA : int 3 4 2 3 2 3 4 2 3 3 ...
## $ PUNT_MATEMATICAS : int 60 69 52 71 49 37 56 36 59 42 ...
## $ PERCENTIL_MATEMATICAS : int 78 94 56 96 46 14 66 11 76 24 ...
## $ DESEMP_MATEMATICAS : int 3 3 3 4 2 2 3 2 3 2 ...
## $ PUNT_C_NATURALES : int 42 64 43 66 41 48 51 37 53 42 ...
## $ PERCENTIL_C_NATURALES : int 31 92 35 95 28 51 62 17 69 31 ...
## $ DESEMP_C_NATURALES : int 2 3 2 3 2 2 2 1 2 2 ...
## $ PUNT_SOCIALES_CIUDADANAS : int 48 69 42 54 39 41 64 39 39 48 ...
## $ PERCENTIL_SOCIALES_CIUDADANAS: int 58 97 42 73 33 38 91 33 31 58 ...
## $ DESEMP_SOCIALES_CIUDADANAS : int 2 3 2 2 1 2 3 1 1 2 ...
## $ PUNT_INGLES : int 51 60 49 56 42 55 67 36 40 45 ...
## $ PERCENTIL_INGLES : int 58 83 53 75 35 70 93 17 27 41 ...
## $ DESEMP_INGLES : chr "A1" "A2" "A1" "A1" ...
## $ PUNT_GLOBAL : int 261 337 222 317 215 227 302 197 253 231 ...
## $ PERCENTIL_GLOBAL : chr "62" "96" "35" "90" ...
## $ ESTU_INSE_INDIVIDUAL : chr "-" "68.1938073220759" "54.5456046491048" "65.4885945948435" ...
## $ ESTU_NSE_INDIVIDUAL : chr "" "4" "3" "4" ...
## $ ESTU_NSE_ESTABLECIMIENTO : int 2 3 3 2 3 3 3 3 3 3 ...
## $ ESTU_ESTADOINVESTIGACION : chr "PUBLICAR" "PUBLICAR" "PUBLICAR" "PUBLICAR" ...
## $ ESTU_GENERACION.E : chr "NO" "NO" "NO" "NO" ...
## $ academico : num 1 1 0 1 1 1 0 0 1 1 ...
## $ u : num 0.000465 0.000625 0.001156 0.001192 0.002479 ...
## $ pik : num 0.287 0.287 0.287 0.287 0.287 ...
head(muestra, n=4)
## ESTU_TIPODOCUMENTO ESTU_NACIONALIDAD ESTU_GENERO ESTU_FECHANACIMIENTO
## 177595 CC COLOMBIA M 10/19/2000 12:00:00 AM
## 31496 TI COLOMBIA M 08/01/2002 12:00:00 AM
## 443512 TI COLOMBIA F 07/24/2002 12:00:00 AM
## 267426 TI COLOMBIA M 05/06/2002 12:00:00 AM
## PERIODO ESTU_CONSECUTIVO ESTU_ESTUDIANTE ESTU_TIENEETNIA
## 177595 20194 SB11201940505953 ESTUDIANTE No
## 31496 20194 SB11201940387065 ESTUDIANTE No
## 443512 20194 SB11201940156899 ESTUDIANTE No
## 267426 20194 SB11201940497663 ESTUDIANTE No
## ESTU_PAIS_RESIDE ESTU_ETNIA ESTU_DEPTO_RESIDE ESTU_COD_RESIDE_DEPTO
## 177595 COLOMBIA - CUNDINAMARCA 25
## 31496 COLOMBIA - CUNDINAMARCA 25
## 443512 COLOMBIA - CUNDINAMARCA 25
## 267426 COLOMBIA - CUNDINAMARCA 25
## ESTU_MCPIO_RESIDE ESTU_COD_RESIDE_MCPIO FAMI_ESTRATOVIVIENDA
## 177595 FUNZA 25286 -
## 31496 MOSQUERA 25473 Estrato 3
## 443512 FUNZA 25286 Estrato 2
## 267426 FUNZA 25286 Estrato 2
## FAMI_PERSONASHOGAR FAMI_CUARTOSHOGAR FAMI_EDUCACIONPADRE
## 177595 -
## 31496 3 a 4 Tres Educación profesional incompleta
## 443512 3 a 4 Dos Secundaria (Bachillerato) completa
## 267426 9 o más Seis o mas Técnica o tecnológica completa
## FAMI_EDUCACIONMADRE
## 177595 -
## 31496 Educación profesional completa
## 443512 Secundaria (Bachillerato) completa
## 267426 Técnica o tecnológica completa
## FAMI_TRABAJOLABORPADRE
## 177595
## 31496 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## 443512 Es operario de máquinas o conduce vehículos (taxita, chofer)
## 267426 Trabaja como profesional (por ejemplo médico, abogado, ingeniero)
## FAMI_TRABAJOLABORMADRE
## 177595 -
## 31496 Es dueño de un negocio grande, tiene un cargo de nivel directivo o gerencial
## 443512 Trabaja en el hogar, no trabaja o estudia
## 267426 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## FAMI_TIENEINTERNET FAMI_TIENESERVICIOTV FAMI_TIENECOMPUTADOR
## 177595 -
## 31496 Si Si Si
## 443512 Si Si Si
## 267426 Si Si Si
## FAMI_TIENELAVADORA FAMI_TIENEHORNOMICROOGAS FAMI_TIENEAUTOMOVIL
## 177595 - -
## 31496 Si Si Si
## 443512 Si Si No
## 267426 Si Si Si
## FAMI_TIENEMOTOCICLETA FAMI_TIENECONSOLAVIDEOJUEGOS FAMI_NUMLIBROS
## 177595 -
## 31496 No Si 11 A 25 LIBROS
## 443512 No No 0 A 10 LIBROS
## 267426 Si No 26 A 100 LIBROS
## FAMI_COMELECHEDERIVADOS FAMI_COMECARNEPESCADOHUEVO
## 177595 -
## 31496 Todos o casi todos los días Todos o casi todos los días
## 443512 Todos o casi todos los días Todos o casi todos los días
## 267426 Todos o casi todos los días Todos o casi todos los días
## FAMI_COMECEREALFRUTOSLEGUMBRE FAMI_SITUACIONECONOMICA
## 177595 -
## 31496 3 a 5 veces por semana Igual
## 443512 1 o 2 veces por semana Igual
## 267426 Todos o casi todos los días Igual
## ESTU_DEDICACIONLECTURADIARIA ESTU_DEDICACIONINTERNET
## 177595 -
## 31496 No leo por entretenimiento Entre 30 y 60 minutos
## 443512 30 minutos o menos Entre 30 y 60 minutos
## 267426 30 minutos o menos Entre 1 y 3 horas
## ESTU_HORASSEMANATRABAJA ESTU_TIPOREMUNERACION COLE_CODIGO_ICFES
## 177595 - 141887
## 31496 0 No 122200
## 443512 0 No 667253
## 267426 0 No 234781
## COLE_COD_DANE_ESTABLECIMIENTO COLE_NOMBRE_ESTABLECIMIENTO COLE_GENERO
## 177595 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## 31496 3.25286e+11 COL CRISTIANO INTEGRAL MIXTO
## 443512 1.25286e+11 I.E.D. BICENTENARIO MIXTO
## 267426 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## COLE_NATURALEZA COLE_CALENDARIO COLE_BILINGUE COLE_CARACTER
## 177595 OFICIAL A N ACADÉMICO
## 31496 NO OFICIAL A N ACADÉMICO
## 443512 OFICIAL A N TÉCNICO/ACADÉMICO
## 267426 OFICIAL A N ACADÉMICO
## COLE_COD_DANE_SEDE COLE_NOMBRE_SEDE
## 177595 1.25286e+11 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## 31496 3.25286e+11 COL CRISTIANO INTEGRAL
## 443512 1.25286e+11 I.E.D. BICENTENARIO - SEDE PRINCIPAL
## 267426 1.25286e+11 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## COLE_SEDE_PRINCIPAL COLE_AREA_UBICACION COLE_JORNADA
## 177595 N URBANO NOCHE
## 31496 S URBANO COMPLETA
## 443512 S URBANO UNICA
## 267426 N URBANO MAÑANA
## COLE_COD_MCPIO_UBICACION COLE_MCPIO_UBICACION COLE_COD_DEPTO_UBICACION
## 177595 25286 FUNZA 25
## 31496 25286 FUNZA 25
## 443512 25286 FUNZA 25
## 267426 25286 FUNZA 25
## COLE_DEPTO_UBICACION ESTU_PRIVADO_LIBERTAD ESTU_COD_MCPIO_PRESENTACION
## 177595 CUNDINAMARCA N 25286
## 31496 CUNDINAMARCA N 25473
## 443512 CUNDINAMARCA N 25286
## 267426 CUNDINAMARCA N 25286
## ESTU_MCPIO_PRESENTACION ESTU_DEPTO_PRESENTACION
## 177595 FUNZA CUNDINAMARCA
## 31496 MOSQUERA CUNDINAMARCA
## 443512 FUNZA CUNDINAMARCA
## 267426 FUNZA CUNDINAMARCA
## ESTU_COD_DEPTO_PRESENTACION PUNT_LECTURA_CRITICA
## 177595 25 59
## 31496 25 70
## 443512 25 39
## 267426 25 65
## PERCENTIL_LECTURA_CRITICA DESEMP_LECTURA_CRITICA PUNT_MATEMATICAS
## 177595 73 3 60
## 31496 96 4 69
## 443512 12 2 52
## 267426 88 3 71
## PERCENTIL_MATEMATICAS DESEMP_MATEMATICAS PUNT_C_NATURALES
## 177595 78 3 42
## 31496 94 3 64
## 443512 56 3 43
## 267426 96 4 66
## PERCENTIL_C_NATURALES DESEMP_C_NATURALES PUNT_SOCIALES_CIUDADANAS
## 177595 31 2 48
## 31496 92 3 69
## 443512 35 2 42
## 267426 95 3 54
## PERCENTIL_SOCIALES_CIUDADANAS DESEMP_SOCIALES_CIUDADANAS PUNT_INGLES
## 177595 58 2 51
## 31496 97 3 60
## 443512 42 2 49
## 267426 73 2 56
## PERCENTIL_INGLES DESEMP_INGLES PUNT_GLOBAL PERCENTIL_GLOBAL
## 177595 58 A1 261 62
## 31496 83 A2 337 96
## 443512 53 A1 222 35
## 267426 75 A1 317 90
## ESTU_INSE_INDIVIDUAL ESTU_NSE_INDIVIDUAL ESTU_NSE_ESTABLECIMIENTO
## 177595 - 2
## 31496 68.1938073220759 4 3
## 443512 54.5456046491048 3 3
## 267426 65.4885945948435 4 2
## ESTU_ESTADOINVESTIGACION ESTU_GENERACION.E academico u
## 177595 PUBLICAR NO 1 0.0004653491
## 31496 PUBLICAR NO 1 0.0006247733
## 443512 PUBLICAR NO 0 0.0011558197
## 267426 PUBLICAR NO 1 0.0011916284
## pik
## 177595 0.2874412
## 31496 0.2874412
## 443512 0.2874412
## 267426 0.2874412
names(muestra)
## [1] "ESTU_TIPODOCUMENTO" "ESTU_NACIONALIDAD"
## [3] "ESTU_GENERO" "ESTU_FECHANACIMIENTO"
## [5] "PERIODO" "ESTU_CONSECUTIVO"
## [7] "ESTU_ESTUDIANTE" "ESTU_TIENEETNIA"
## [9] "ESTU_PAIS_RESIDE" "ESTU_ETNIA"
## [11] "ESTU_DEPTO_RESIDE" "ESTU_COD_RESIDE_DEPTO"
## [13] "ESTU_MCPIO_RESIDE" "ESTU_COD_RESIDE_MCPIO"
## [15] "FAMI_ESTRATOVIVIENDA" "FAMI_PERSONASHOGAR"
## [17] "FAMI_CUARTOSHOGAR" "FAMI_EDUCACIONPADRE"
## [19] "FAMI_EDUCACIONMADRE" "FAMI_TRABAJOLABORPADRE"
## [21] "FAMI_TRABAJOLABORMADRE" "FAMI_TIENEINTERNET"
## [23] "FAMI_TIENESERVICIOTV" "FAMI_TIENECOMPUTADOR"
## [25] "FAMI_TIENELAVADORA" "FAMI_TIENEHORNOMICROOGAS"
## [27] "FAMI_TIENEAUTOMOVIL" "FAMI_TIENEMOTOCICLETA"
## [29] "FAMI_TIENECONSOLAVIDEOJUEGOS" "FAMI_NUMLIBROS"
## [31] "FAMI_COMELECHEDERIVADOS" "FAMI_COMECARNEPESCADOHUEVO"
## [33] "FAMI_COMECEREALFRUTOSLEGUMBRE" "FAMI_SITUACIONECONOMICA"
## [35] "ESTU_DEDICACIONLECTURADIARIA" "ESTU_DEDICACIONINTERNET"
## [37] "ESTU_HORASSEMANATRABAJA" "ESTU_TIPOREMUNERACION"
## [39] "COLE_CODIGO_ICFES" "COLE_COD_DANE_ESTABLECIMIENTO"
## [41] "COLE_NOMBRE_ESTABLECIMIENTO" "COLE_GENERO"
## [43] "COLE_NATURALEZA" "COLE_CALENDARIO"
## [45] "COLE_BILINGUE" "COLE_CARACTER"
## [47] "COLE_COD_DANE_SEDE" "COLE_NOMBRE_SEDE"
## [49] "COLE_SEDE_PRINCIPAL" "COLE_AREA_UBICACION"
## [51] "COLE_JORNADA" "COLE_COD_MCPIO_UBICACION"
## [53] "COLE_MCPIO_UBICACION" "COLE_COD_DEPTO_UBICACION"
## [55] "COLE_DEPTO_UBICACION" "ESTU_PRIVADO_LIBERTAD"
## [57] "ESTU_COD_MCPIO_PRESENTACION" "ESTU_MCPIO_PRESENTACION"
## [59] "ESTU_DEPTO_PRESENTACION" "ESTU_COD_DEPTO_PRESENTACION"
## [61] "PUNT_LECTURA_CRITICA" "PERCENTIL_LECTURA_CRITICA"
## [63] "DESEMP_LECTURA_CRITICA" "PUNT_MATEMATICAS"
## [65] "PERCENTIL_MATEMATICAS" "DESEMP_MATEMATICAS"
## [67] "PUNT_C_NATURALES" "PERCENTIL_C_NATURALES"
## [69] "DESEMP_C_NATURALES" "PUNT_SOCIALES_CIUDADANAS"
## [71] "PERCENTIL_SOCIALES_CIUDADANAS" "DESEMP_SOCIALES_CIUDADANAS"
## [73] "PUNT_INGLES" "PERCENTIL_INGLES"
## [75] "DESEMP_INGLES" "PUNT_GLOBAL"
## [77] "PERCENTIL_GLOBAL" "ESTU_INSE_INDIVIDUAL"
## [79] "ESTU_NSE_INDIVIDUAL" "ESTU_NSE_ESTABLECIMIENTO"
## [81] "ESTU_ESTADOINVESTIGACION" "ESTU_GENERACION.E"
## [83] "academico" "u"
## [85] "pik"
Selección y rechazo
Fan_Muller=function(base, n, seed){
N=nrow(base)
j = 0
m = numeric(N)
set.seed(seed)
for (k in 1:N) if (runif(1)< (n-j)/(N-k+1)) {
j = j + 1
m[k] = 1
}
return(m)
}
n = 428
seed = 123
m = Fan_Muller(funza, n, seed)
muestra2=funza[which(m==1),]
str(muestra2)
## 'data.frame': 428 obs. of 83 variables:
## $ ESTU_TIPODOCUMENTO : chr "TI" "TI" "CC" "TI" ...
## $ ESTU_NACIONALIDAD : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_GENERO : chr "F" "M" "M" "F" ...
## $ ESTU_FECHANACIMIENTO : chr "05/02/2002 12:00:00 AM" "01/23/2003 12:00:00 AM" "03/08/2001 12:00:00 AM" "05/25/2003 12:00:00 AM" ...
## $ PERIODO : int 20194 20194 20194 20194 20194 20194 20194 20194 20194 20194 ...
## $ ESTU_CONSECUTIVO : chr "SB11201940040952" "SB11201940111850" "SB11201940028648" "SB11201940066729" ...
## $ ESTU_ESTUDIANTE : chr "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" "ESTUDIANTE" ...
## $ ESTU_TIENEETNIA : chr "No" "No" "No" "No" ...
## $ ESTU_PAIS_RESIDE : chr "COLOMBIA" "COLOMBIA" "COLOMBIA" "COLOMBIA" ...
## $ ESTU_ETNIA : chr "-" "-" "-" "-" ...
## $ ESTU_DEPTO_RESIDE : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_COD_RESIDE_DEPTO : chr "25" "25" "25" "25" ...
## $ ESTU_MCPIO_RESIDE : chr "FUNZA" "FUNZA" "FUNZA" "MOSQUERA" ...
## $ ESTU_COD_RESIDE_MCPIO : chr "25286" "25286" "25286" "25473" ...
## $ FAMI_ESTRATOVIVIENDA : chr "Estrato 2" "Estrato 4" "Estrato 3" "Estrato 2" ...
## $ FAMI_PERSONASHOGAR : chr "3 a 4" "7 a 8" "3 a 4" "3 a 4" ...
## $ FAMI_CUARTOSHOGAR : chr "Tres" "Cuatro" "Dos" "Tres" ...
## $ FAMI_EDUCACIONPADRE : chr "Primaria incompleta" "Secundaria (Bachillerato) completa" "Secundaria (Bachillerato) completa" "Técnica o tecnológica completa" ...
## $ FAMI_EDUCACIONMADRE : chr "Primaria incompleta" "Secundaria (Bachillerato) completa" "Técnica o tecnológica completa" "Secundaria (Bachillerato) completa" ...
## $ FAMI_TRABAJOLABORPADRE : chr "Es agricultor, pesquero o jornalero" "Trabaja como profesional (por ejemplo médico, abogado, ingeniero)" "Es operario de máquinas o conduce vehículos (taxita, chofer)" "Trabaja como profesional (por ejemplo médico, abogado, ingeniero)" ...
## $ FAMI_TRABAJOLABORMADRE : chr "Trabaja en el hogar, no trabaja o estudia" "Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc" "Trabaja en el hogar, no trabaja o estudia" "Trabaja en el hogar, no trabaja o estudia" ...
## $ FAMI_TIENEINTERNET : chr "No" "Si" "Si" "Si" ...
## $ FAMI_TIENESERVICIOTV : chr "Si" "Si" "Si" "Si" ...
## $ FAMI_TIENECOMPUTADOR : chr "No" "Si" "Si" "Si" ...
## $ FAMI_TIENELAVADORA : chr "Si" "Si" "Si" "Si" ...
## $ FAMI_TIENEHORNOMICROOGAS : chr "No" "Si" "No" "No" ...
## $ FAMI_TIENEAUTOMOVIL : chr "No" "Si" "No" "No" ...
## $ FAMI_TIENEMOTOCICLETA : chr "Si" "Si" "No" "Si" ...
## $ FAMI_TIENECONSOLAVIDEOJUEGOS : chr "No" "No" "No" "No" ...
## $ FAMI_NUMLIBROS : chr "11 A 25 LIBROS" "26 A 100 LIBROS" "-" "MÁS DE 100 LIBROS" ...
## $ FAMI_COMELECHEDERIVADOS : chr "1 o 2 veces por semana" "Todos o casi todos los días" "1 o 2 veces por semana" "Todos o casi todos los días" ...
## $ FAMI_COMECARNEPESCADOHUEVO : chr "1 o 2 veces por semana" "3 a 5 veces por semana" "3 a 5 veces por semana" "3 a 5 veces por semana" ...
## $ FAMI_COMECEREALFRUTOSLEGUMBRE: chr "1 o 2 veces por semana" "3 a 5 veces por semana" "1 o 2 veces por semana" "3 a 5 veces por semana" ...
## $ FAMI_SITUACIONECONOMICA : chr "Igual" "Igual" "Mejor" "Igual" ...
## $ ESTU_DEDICACIONLECTURADIARIA : chr "Entre 30 y 60 minutos" "30 minutos o menos" "Entre 30 y 60 minutos" "30 minutos o menos" ...
## $ ESTU_DEDICACIONINTERNET : chr "Entre 30 y 60 minutos" "Entre 1 y 3 horas" "Entre 1 y 3 horas" "Entre 1 y 3 horas" ...
## $ ESTU_HORASSEMANATRABAJA : chr "0" "0" "0" "0" ...
## $ ESTU_TIPOREMUNERACION : chr "No" "No" "No" "No" ...
## $ COLE_CODIGO_ICFES : int 717744 150599 613026 175422 24935 234781 24927 41327 159525 122200 ...
## $ COLE_COD_DANE_ESTABLECIMIENTO: num 1.25e+11 3.25e+11 3.25e+11 3.25e+11 3.25e+11 ...
## $ COLE_NOMBRE_ESTABLECIMIENTO : chr "I.E.D. FUNZA" "COLEGIO INGLES NUEVA INGLATERRA" "I E SANTA ANA" "COL DE FORMACION INTEGRAL SEMBRADORES DE PAZ" ...
## $ COLE_GENERO : chr "MIXTO" "MIXTO" "MIXTO" "MIXTO" ...
## $ COLE_NATURALEZA : chr "OFICIAL" "NO OFICIAL" "NO OFICIAL" "NO OFICIAL" ...
## $ COLE_CALENDARIO : chr "A" "A" "A" "A" ...
## $ COLE_BILINGUE : chr "N" "-" "N" "N" ...
## $ COLE_CARACTER : chr "ACADÉMICO" "TÉCNICO/ACADÉMICO" "ACADÉMICO" "ACADÉMICO" ...
## $ COLE_COD_DANE_SEDE : num 1.25e+11 3.25e+11 3.25e+11 3.25e+11 3.25e+11 ...
## $ COLE_NOMBRE_SEDE : chr "FURATENA" "COL INGLES NUEVA INGLATERRA" "COLEGIO SANTA ANA" "COL DE FORMACION INTEGRAL SEMBRADORES DE PAZ" ...
## $ COLE_SEDE_PRINCIPAL : chr "N" "S" "S" "S" ...
## $ COLE_AREA_UBICACION : chr "URBANO" "URBANO" "URBANO" "URBANO" ...
## $ COLE_JORNADA : chr "MAÑANA" "COMPLETA" "MAÑANA" "MAÑANA" ...
## $ COLE_COD_MCPIO_UBICACION : int 25286 25286 25286 25286 25286 25286 25286 25286 25286 25286 ...
## $ COLE_MCPIO_UBICACION : chr "FUNZA" "FUNZA" "FUNZA" "FUNZA" ...
## $ COLE_COD_DEPTO_UBICACION : int 25 25 25 25 25 25 25 25 25 25 ...
## $ COLE_DEPTO_UBICACION : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_PRIVADO_LIBERTAD : chr "N" "N" "N" "N" ...
## $ ESTU_COD_MCPIO_PRESENTACION : int 25286 25286 25286 25473 25473 25286 25286 25473 25286 25286 ...
## $ ESTU_MCPIO_PRESENTACION : chr "FUNZA" "FUNZA" "FUNZA" "MOSQUERA" ...
## $ ESTU_DEPTO_PRESENTACION : chr "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" "CUNDINAMARCA" ...
## $ ESTU_COD_DEPTO_PRESENTACION : int 25 25 25 25 25 25 25 25 25 25 ...
## $ PUNT_LECTURA_CRITICA : int 72 56 47 63 71 48 74 57 66 62 ...
## $ PERCENTIL_LECTURA_CRITICA : int 98 62 33 85 97 36 99 67 91 82 ...
## $ DESEMP_LECTURA_CRITICA : int 4 3 2 3 4 2 4 3 4 3 ...
## $ PUNT_MATEMATICAS : int 60 51 49 61 69 43 63 48 55 70 ...
## $ PERCENTIL_MATEMATICAS : int 77 52 43 80 95 28 84 42 65 95 ...
## $ DESEMP_MATEMATICAS : int 3 3 2 3 3 2 3 2 3 3 ...
## $ PUNT_C_NATURALES : int 62 52 42 60 68 43 66 47 46 52 ...
## $ PERCENTIL_C_NATURALES : int 89 64 30 85 97 34 95 47 45 63 ...
## $ DESEMP_C_NATURALES : int 3 2 2 3 3 2 3 2 2 2 ...
## $ PUNT_SOCIALES_CIUDADANAS : int 69 46 32 63 80 37 61 44 38 39 ...
## $ PERCENTIL_SOCIALES_CIUDADANAS: int 97 53 13 91 100 25 86 48 31 32 ...
## $ DESEMP_SOCIALES_CIUDADANAS : int 3 2 1 3 4 1 3 2 1 1 ...
## $ PUNT_INGLES : int 40 58 49 63 73 43 71 39 69 44 ...
## $ PERCENTIL_INGLES : int 28 78 55 88 97 36 95 26 94 39 ...
## $ DESEMP_INGLES : chr "A-" "A2" "A1" "A2" ...
## $ PUNT_GLOBAL : int 319 259 215 309 360 214 332 241 263 274 ...
## $ PERCENTIL_GLOBAL : chr "91" "61" "30" "88" ...
## $ ESTU_INSE_INDIVIDUAL : chr "38.1705705683999" "62.3941798982287" "52.0019623816" "58.2733846535184" ...
## $ ESTU_NSE_INDIVIDUAL : chr "1" "3" "3" "3" ...
## $ ESTU_NSE_ESTABLECIMIENTO : int 3 3 3 3 3 2 4 3 3 3 ...
## $ ESTU_ESTADOINVESTIGACION : chr "PUBLICAR" "PUBLICAR" "PUBLICAR" "PUBLICAR" ...
## $ ESTU_GENERACION.E : chr "NO" "NO" "NO" "GENERACION E - GRATUIDAD" ...
## $ academico : num 1 0 1 1 1 1 0 1 1 1 ...
head(muestra2, n=4)
## ESTU_TIPODOCUMENTO ESTU_NACIONALIDAD ESTU_GENERO ESTU_FECHANACIMIENTO
## 4582 TI COLOMBIA F 05/02/2002 12:00:00 AM
## 9063 TI COLOMBIA M 01/23/2003 12:00:00 AM
## 9598 CC COLOMBIA M 03/08/2001 12:00:00 AM
## 10246 TI COLOMBIA F 05/25/2003 12:00:00 AM
## PERIODO ESTU_CONSECUTIVO ESTU_ESTUDIANTE ESTU_TIENEETNIA ESTU_PAIS_RESIDE
## 4582 20194 SB11201940040952 ESTUDIANTE No COLOMBIA
## 9063 20194 SB11201940111850 ESTUDIANTE No COLOMBIA
## 9598 20194 SB11201940028648 ESTUDIANTE No COLOMBIA
## 10246 20194 SB11201940066729 ESTUDIANTE No COLOMBIA
## ESTU_ETNIA ESTU_DEPTO_RESIDE ESTU_COD_RESIDE_DEPTO ESTU_MCPIO_RESIDE
## 4582 - CUNDINAMARCA 25 FUNZA
## 9063 - CUNDINAMARCA 25 FUNZA
## 9598 - CUNDINAMARCA 25 FUNZA
## 10246 - CUNDINAMARCA 25 MOSQUERA
## ESTU_COD_RESIDE_MCPIO FAMI_ESTRATOVIVIENDA FAMI_PERSONASHOGAR
## 4582 25286 Estrato 2 3 a 4
## 9063 25286 Estrato 4 7 a 8
## 9598 25286 Estrato 3 3 a 4
## 10246 25473 Estrato 2 3 a 4
## FAMI_CUARTOSHOGAR FAMI_EDUCACIONPADRE
## 4582 Tres Primaria incompleta
## 9063 Cuatro Secundaria (Bachillerato) completa
## 9598 Dos Secundaria (Bachillerato) completa
## 10246 Tres Técnica o tecnológica completa
## FAMI_EDUCACIONMADRE
## 4582 Primaria incompleta
## 9063 Secundaria (Bachillerato) completa
## 9598 Técnica o tecnológica completa
## 10246 Secundaria (Bachillerato) completa
## FAMI_TRABAJOLABORPADRE
## 4582 Es agricultor, pesquero o jornalero
## 9063 Trabaja como profesional (por ejemplo médico, abogado, ingeniero)
## 9598 Es operario de máquinas o conduce vehículos (taxita, chofer)
## 10246 Trabaja como profesional (por ejemplo médico, abogado, ingeniero)
## FAMI_TRABAJOLABORMADRE
## 4582 Trabaja en el hogar, no trabaja o estudia
## 9063 Es dueño de un negocio pequeño (tiene pocos empleados o no tiene, por ejemplo tienda, papelería, etc
## 9598 Trabaja en el hogar, no trabaja o estudia
## 10246 Trabaja en el hogar, no trabaja o estudia
## FAMI_TIENEINTERNET FAMI_TIENESERVICIOTV FAMI_TIENECOMPUTADOR
## 4582 No Si No
## 9063 Si Si Si
## 9598 Si Si Si
## 10246 Si Si Si
## FAMI_TIENELAVADORA FAMI_TIENEHORNOMICROOGAS FAMI_TIENEAUTOMOVIL
## 4582 Si No No
## 9063 Si Si Si
## 9598 Si No No
## 10246 Si No No
## FAMI_TIENEMOTOCICLETA FAMI_TIENECONSOLAVIDEOJUEGOS FAMI_NUMLIBROS
## 4582 Si No 11 A 25 LIBROS
## 9063 Si No 26 A 100 LIBROS
## 9598 No No -
## 10246 Si No MÁS DE 100 LIBROS
## FAMI_COMELECHEDERIVADOS FAMI_COMECARNEPESCADOHUEVO
## 4582 1 o 2 veces por semana 1 o 2 veces por semana
## 9063 Todos o casi todos los días 3 a 5 veces por semana
## 9598 1 o 2 veces por semana 3 a 5 veces por semana
## 10246 Todos o casi todos los días 3 a 5 veces por semana
## FAMI_COMECEREALFRUTOSLEGUMBRE FAMI_SITUACIONECONOMICA
## 4582 1 o 2 veces por semana Igual
## 9063 3 a 5 veces por semana Igual
## 9598 1 o 2 veces por semana Mejor
## 10246 3 a 5 veces por semana Igual
## ESTU_DEDICACIONLECTURADIARIA ESTU_DEDICACIONINTERNET
## 4582 Entre 30 y 60 minutos Entre 30 y 60 minutos
## 9063 30 minutos o menos Entre 1 y 3 horas
## 9598 Entre 30 y 60 minutos Entre 1 y 3 horas
## 10246 30 minutos o menos Entre 1 y 3 horas
## ESTU_HORASSEMANATRABAJA ESTU_TIPOREMUNERACION COLE_CODIGO_ICFES
## 4582 0 No 717744
## 9063 0 No 150599
## 9598 0 No 613026
## 10246 0 No 175422
## COLE_COD_DANE_ESTABLECIMIENTO
## 4582 1.25286e+11
## 9063 3.25286e+11
## 9598 3.25286e+11
## 10246 3.25286e+11
## COLE_NOMBRE_ESTABLECIMIENTO COLE_GENERO COLE_NATURALEZA
## 4582 I.E.D. FUNZA MIXTO OFICIAL
## 9063 COLEGIO INGLES NUEVA INGLATERRA MIXTO NO OFICIAL
## 9598 I E SANTA ANA MIXTO NO OFICIAL
## 10246 COL DE FORMACION INTEGRAL SEMBRADORES DE PAZ MIXTO NO OFICIAL
## COLE_CALENDARIO COLE_BILINGUE COLE_CARACTER COLE_COD_DANE_SEDE
## 4582 A N ACADÉMICO 125286800006
## 9063 A - TÉCNICO/ACADÉMICO 325286000840
## 9598 A N ACADÉMICO 325286000149
## 10246 A N ACADÉMICO 325286000858
## COLE_NOMBRE_SEDE COLE_SEDE_PRINCIPAL
## 4582 FURATENA N
## 9063 COL INGLES NUEVA INGLATERRA S
## 9598 COLEGIO SANTA ANA S
## 10246 COL DE FORMACION INTEGRAL SEMBRADORES DE PAZ S
## COLE_AREA_UBICACION COLE_JORNADA COLE_COD_MCPIO_UBICACION
## 4582 URBANO MAÑANA 25286
## 9063 URBANO COMPLETA 25286
## 9598 URBANO MAÑANA 25286
## 10246 URBANO MAÑANA 25286
## COLE_MCPIO_UBICACION COLE_COD_DEPTO_UBICACION COLE_DEPTO_UBICACION
## 4582 FUNZA 25 CUNDINAMARCA
## 9063 FUNZA 25 CUNDINAMARCA
## 9598 FUNZA 25 CUNDINAMARCA
## 10246 FUNZA 25 CUNDINAMARCA
## ESTU_PRIVADO_LIBERTAD ESTU_COD_MCPIO_PRESENTACION ESTU_MCPIO_PRESENTACION
## 4582 N 25286 FUNZA
## 9063 N 25286 FUNZA
## 9598 N 25286 FUNZA
## 10246 N 25473 MOSQUERA
## ESTU_DEPTO_PRESENTACION ESTU_COD_DEPTO_PRESENTACION PUNT_LECTURA_CRITICA
## 4582 CUNDINAMARCA 25 72
## 9063 CUNDINAMARCA 25 56
## 9598 CUNDINAMARCA 25 47
## 10246 CUNDINAMARCA 25 63
## PERCENTIL_LECTURA_CRITICA DESEMP_LECTURA_CRITICA PUNT_MATEMATICAS
## 4582 98 4 60
## 9063 62 3 51
## 9598 33 2 49
## 10246 85 3 61
## PERCENTIL_MATEMATICAS DESEMP_MATEMATICAS PUNT_C_NATURALES
## 4582 77 3 62
## 9063 52 3 52
## 9598 43 2 42
## 10246 80 3 60
## PERCENTIL_C_NATURALES DESEMP_C_NATURALES PUNT_SOCIALES_CIUDADANAS
## 4582 89 3 69
## 9063 64 2 46
## 9598 30 2 32
## 10246 85 3 63
## PERCENTIL_SOCIALES_CIUDADANAS DESEMP_SOCIALES_CIUDADANAS PUNT_INGLES
## 4582 97 3 40
## 9063 53 2 58
## 9598 13 1 49
## 10246 91 3 63
## PERCENTIL_INGLES DESEMP_INGLES PUNT_GLOBAL PERCENTIL_GLOBAL
## 4582 28 A- 319 91
## 9063 78 A2 259 61
## 9598 55 A1 215 30
## 10246 88 A2 309 88
## ESTU_INSE_INDIVIDUAL ESTU_NSE_INDIVIDUAL ESTU_NSE_ESTABLECIMIENTO
## 4582 38.1705705683999 1 3
## 9063 62.3941798982287 3 3
## 9598 52.0019623816 3 3
## 10246 58.2733846535184 3 3
## ESTU_ESTADOINVESTIGACION ESTU_GENERACION.E academico
## 4582 PUBLICAR NO 1
## 9063 PUBLICAR NO 0
## 9598 PUBLICAR NO 1
## 10246 PUBLICAR GENERACION E - GRATUIDAD 1
5. El pi-estimador bajo el MAS en una etapa.
salida = function(est, alpha){
est = as.data.frame(est)
names(est)[2]="se"
est$cv = 100*(est$se/est[,1])
est$ic_low = est[,1]-qnorm(1-alpha/2)*est$se
est$ic_upp = est[,1]+qnorm(1-alpha/2)*est$se
return(round(est, 2))
}
N=nrow(funza)
muestra$ind = rep(1, nrow(muestra))
#factor de expansión.
muestra$Fexp=1/muestra$pik
head(muestra, n=4)
## ESTU_TIPODOCUMENTO ESTU_NACIONALIDAD ESTU_GENERO ESTU_FECHANACIMIENTO
## 177595 CC COLOMBIA M 10/19/2000 12:00:00 AM
## 31496 TI COLOMBIA M 08/01/2002 12:00:00 AM
## 443512 TI COLOMBIA F 07/24/2002 12:00:00 AM
## 267426 TI COLOMBIA M 05/06/2002 12:00:00 AM
## PERIODO ESTU_CONSECUTIVO ESTU_ESTUDIANTE ESTU_TIENEETNIA
## 177595 20194 SB11201940505953 ESTUDIANTE No
## 31496 20194 SB11201940387065 ESTUDIANTE No
## 443512 20194 SB11201940156899 ESTUDIANTE No
## 267426 20194 SB11201940497663 ESTUDIANTE No
## ESTU_PAIS_RESIDE ESTU_ETNIA ESTU_DEPTO_RESIDE ESTU_COD_RESIDE_DEPTO
## 177595 COLOMBIA - CUNDINAMARCA 25
## 31496 COLOMBIA - CUNDINAMARCA 25
## 443512 COLOMBIA - CUNDINAMARCA 25
## 267426 COLOMBIA - CUNDINAMARCA 25
## ESTU_MCPIO_RESIDE ESTU_COD_RESIDE_MCPIO FAMI_ESTRATOVIVIENDA
## 177595 FUNZA 25286 -
## 31496 MOSQUERA 25473 Estrato 3
## 443512 FUNZA 25286 Estrato 2
## 267426 FUNZA 25286 Estrato 2
## FAMI_PERSONASHOGAR FAMI_CUARTOSHOGAR FAMI_EDUCACIONPADRE
## 177595 -
## 31496 3 a 4 Tres Educación profesional incompleta
## 443512 3 a 4 Dos Secundaria (Bachillerato) completa
## 267426 9 o más Seis o mas Técnica o tecnológica completa
## FAMI_EDUCACIONMADRE
## 177595 -
## 31496 Educación profesional completa
## 443512 Secundaria (Bachillerato) completa
## 267426 Técnica o tecnológica completa
## FAMI_TRABAJOLABORPADRE
## 177595
## 31496 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## 443512 Es operario de máquinas o conduce vehículos (taxita, chofer)
## 267426 Trabaja como profesional (por ejemplo médico, abogado, ingeniero)
## FAMI_TRABAJOLABORMADRE
## 177595 -
## 31496 Es dueño de un negocio grande, tiene un cargo de nivel directivo o gerencial
## 443512 Trabaja en el hogar, no trabaja o estudia
## 267426 Tiene un trabajo de tipo auxiliar administrativo (por ejemplo, secretario o asistente)
## FAMI_TIENEINTERNET FAMI_TIENESERVICIOTV FAMI_TIENECOMPUTADOR
## 177595 -
## 31496 Si Si Si
## 443512 Si Si Si
## 267426 Si Si Si
## FAMI_TIENELAVADORA FAMI_TIENEHORNOMICROOGAS FAMI_TIENEAUTOMOVIL
## 177595 - -
## 31496 Si Si Si
## 443512 Si Si No
## 267426 Si Si Si
## FAMI_TIENEMOTOCICLETA FAMI_TIENECONSOLAVIDEOJUEGOS FAMI_NUMLIBROS
## 177595 -
## 31496 No Si 11 A 25 LIBROS
## 443512 No No 0 A 10 LIBROS
## 267426 Si No 26 A 100 LIBROS
## FAMI_COMELECHEDERIVADOS FAMI_COMECARNEPESCADOHUEVO
## 177595 -
## 31496 Todos o casi todos los días Todos o casi todos los días
## 443512 Todos o casi todos los días Todos o casi todos los días
## 267426 Todos o casi todos los días Todos o casi todos los días
## FAMI_COMECEREALFRUTOSLEGUMBRE FAMI_SITUACIONECONOMICA
## 177595 -
## 31496 3 a 5 veces por semana Igual
## 443512 1 o 2 veces por semana Igual
## 267426 Todos o casi todos los días Igual
## ESTU_DEDICACIONLECTURADIARIA ESTU_DEDICACIONINTERNET
## 177595 -
## 31496 No leo por entretenimiento Entre 30 y 60 minutos
## 443512 30 minutos o menos Entre 30 y 60 minutos
## 267426 30 minutos o menos Entre 1 y 3 horas
## ESTU_HORASSEMANATRABAJA ESTU_TIPOREMUNERACION COLE_CODIGO_ICFES
## 177595 - 141887
## 31496 0 No 122200
## 443512 0 No 667253
## 267426 0 No 234781
## COLE_COD_DANE_ESTABLECIMIENTO COLE_NOMBRE_ESTABLECIMIENTO COLE_GENERO
## 177595 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## 31496 3.25286e+11 COL CRISTIANO INTEGRAL MIXTO
## 443512 1.25286e+11 I.E.D. BICENTENARIO MIXTO
## 267426 1.25286e+11 I.E.D. MIGUEL ANTONIO CARO MIXTO
## COLE_NATURALEZA COLE_CALENDARIO COLE_BILINGUE COLE_CARACTER
## 177595 OFICIAL A N ACADÉMICO
## 31496 NO OFICIAL A N ACADÉMICO
## 443512 OFICIAL A N TÉCNICO/ACADÉMICO
## 267426 OFICIAL A N ACADÉMICO
## COLE_COD_DANE_SEDE COLE_NOMBRE_SEDE
## 177595 1.25286e+11 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## 31496 3.25286e+11 COL CRISTIANO INTEGRAL
## 443512 1.25286e+11 I.E.D. BICENTENARIO - SEDE PRINCIPAL
## 267426 1.25286e+11 CONCENTRACION URBANA MIGUEL ANTONIO CARO
## COLE_SEDE_PRINCIPAL COLE_AREA_UBICACION COLE_JORNADA
## 177595 N URBANO NOCHE
## 31496 S URBANO COMPLETA
## 443512 S URBANO UNICA
## 267426 N URBANO MAÑANA
## COLE_COD_MCPIO_UBICACION COLE_MCPIO_UBICACION COLE_COD_DEPTO_UBICACION
## 177595 25286 FUNZA 25
## 31496 25286 FUNZA 25
## 443512 25286 FUNZA 25
## 267426 25286 FUNZA 25
## COLE_DEPTO_UBICACION ESTU_PRIVADO_LIBERTAD ESTU_COD_MCPIO_PRESENTACION
## 177595 CUNDINAMARCA N 25286
## 31496 CUNDINAMARCA N 25473
## 443512 CUNDINAMARCA N 25286
## 267426 CUNDINAMARCA N 25286
## ESTU_MCPIO_PRESENTACION ESTU_DEPTO_PRESENTACION
## 177595 FUNZA CUNDINAMARCA
## 31496 MOSQUERA CUNDINAMARCA
## 443512 FUNZA CUNDINAMARCA
## 267426 FUNZA CUNDINAMARCA
## ESTU_COD_DEPTO_PRESENTACION PUNT_LECTURA_CRITICA
## 177595 25 59
## 31496 25 70
## 443512 25 39
## 267426 25 65
## PERCENTIL_LECTURA_CRITICA DESEMP_LECTURA_CRITICA PUNT_MATEMATICAS
## 177595 73 3 60
## 31496 96 4 69
## 443512 12 2 52
## 267426 88 3 71
## PERCENTIL_MATEMATICAS DESEMP_MATEMATICAS PUNT_C_NATURALES
## 177595 78 3 42
## 31496 94 3 64
## 443512 56 3 43
## 267426 96 4 66
## PERCENTIL_C_NATURALES DESEMP_C_NATURALES PUNT_SOCIALES_CIUDADANAS
## 177595 31 2 48
## 31496 92 3 69
## 443512 35 2 42
## 267426 95 3 54
## PERCENTIL_SOCIALES_CIUDADANAS DESEMP_SOCIALES_CIUDADANAS PUNT_INGLES
## 177595 58 2 51
## 31496 97 3 60
## 443512 42 2 49
## 267426 73 2 56
## PERCENTIL_INGLES DESEMP_INGLES PUNT_GLOBAL PERCENTIL_GLOBAL
## 177595 58 A1 261 62
## 31496 83 A2 337 96
## 443512 53 A1 222 35
## 267426 75 A1 317 90
## ESTU_INSE_INDIVIDUAL ESTU_NSE_INDIVIDUAL ESTU_NSE_ESTABLECIMIENTO
## 177595 - 2
## 31496 68.1938073220759 4 3
## 443512 54.5456046491048 3 3
## 267426 65.4885945948435 4 2
## ESTU_ESTADOINVESTIGACION ESTU_GENERACION.E academico u
## 177595 PUBLICAR NO 1 0.0004653491
## 31496 PUBLICAR NO 1 0.0006247733
## 443512 PUBLICAR NO 0 0.0011558197
## 267426 PUBLICAR NO 1 0.0011916284
## pik ind Fexp
## 177595 0.2874412 1 3.478972
## 31496 0.2874412 1 3.478972
## 443512 0.2874412 1 3.478972
## 267426 0.2874412 1 3.478972
attach(muestra)
## The following objects are masked from funza:
##
## academico, COLE_AREA_UBICACION, COLE_BILINGUE, COLE_CALENDARIO,
## COLE_CARACTER, COLE_COD_DANE_ESTABLECIMIENTO, COLE_COD_DANE_SEDE,
## COLE_COD_DEPTO_UBICACION, COLE_COD_MCPIO_UBICACION,
## COLE_CODIGO_ICFES, COLE_DEPTO_UBICACION, COLE_GENERO, COLE_JORNADA,
## COLE_MCPIO_UBICACION, COLE_NATURALEZA, COLE_NOMBRE_ESTABLECIMIENTO,
## COLE_NOMBRE_SEDE, COLE_SEDE_PRINCIPAL, DESEMP_C_NATURALES,
## DESEMP_INGLES, DESEMP_LECTURA_CRITICA, DESEMP_MATEMATICAS,
## DESEMP_SOCIALES_CIUDADANAS, ESTU_COD_DEPTO_PRESENTACION,
## ESTU_COD_MCPIO_PRESENTACION, ESTU_COD_RESIDE_DEPTO,
## ESTU_COD_RESIDE_MCPIO, ESTU_CONSECUTIVO, ESTU_DEDICACIONINTERNET,
## ESTU_DEDICACIONLECTURADIARIA, ESTU_DEPTO_PRESENTACION,
## ESTU_DEPTO_RESIDE, ESTU_ESTADOINVESTIGACION, ESTU_ESTUDIANTE,
## ESTU_ETNIA, ESTU_FECHANACIMIENTO, ESTU_GENERACION.E, ESTU_GENERO,
## ESTU_HORASSEMANATRABAJA, ESTU_INSE_INDIVIDUAL,
## ESTU_MCPIO_PRESENTACION, ESTU_MCPIO_RESIDE, ESTU_NACIONALIDAD,
## ESTU_NSE_ESTABLECIMIENTO, ESTU_NSE_INDIVIDUAL, ESTU_PAIS_RESIDE,
## ESTU_PRIVADO_LIBERTAD, ESTU_TIENEETNIA, ESTU_TIPODOCUMENTO,
## ESTU_TIPOREMUNERACION, FAMI_COMECARNEPESCADOHUEVO,
## FAMI_COMECEREALFRUTOSLEGUMBRE, FAMI_COMELECHEDERIVADOS,
## FAMI_CUARTOSHOGAR, FAMI_EDUCACIONMADRE, FAMI_EDUCACIONPADRE,
## FAMI_ESTRATOVIVIENDA, FAMI_NUMLIBROS, FAMI_PERSONASHOGAR,
## FAMI_SITUACIONECONOMICA, FAMI_TIENEAUTOMOVIL, FAMI_TIENECOMPUTADOR,
## FAMI_TIENECONSOLAVIDEOJUEGOS, FAMI_TIENEHORNOMICROOGAS,
## FAMI_TIENEINTERNET, FAMI_TIENELAVADORA, FAMI_TIENEMOTOCICLETA,
## FAMI_TIENESERVICIOTV, FAMI_TRABAJOLABORMADRE,
## FAMI_TRABAJOLABORPADRE, PERCENTIL_C_NATURALES, PERCENTIL_GLOBAL,
## PERCENTIL_INGLES, PERCENTIL_LECTURA_CRITICA, PERCENTIL_MATEMATICAS,
## PERCENTIL_SOCIALES_CIUDADANAS, PERIODO, PUNT_C_NATURALES,
## PUNT_GLOBAL, PUNT_INGLES, PUNT_LECTURA_CRITICA, PUNT_MATEMATICAS,
## PUNT_SOCIALES_CIUDADANAS
5.1. Especificación del diseño.
- id=~1: Indica que no hay conglomerados, cada
individuo se seleccionó directamente.
- fpc = ~rep(N,n): Es la corrección por población
finita. Le dice a survey que la población tiene tamaño N, para que
ajuste la varianza (porque al muestrear una proporción grande de la
población, la varianza se reduce).
- data = muestra: Son los datos.
- weights = ~Fexp: Los factores de expansión.
dsgn=svydesign(id=~1,fpc=~rep(N,n),data=muestra,weights=~Fexp)
5.2. Estimación del total.
(est =svytotal(~academico, dsgn, deff=TRUE))
## total SE DEff
## academico 1182.850 24.583 1
alpha = 0.05
(tabla=salida(est, alpha))
## total se deff cv ic_low ic_upp
## academico 1182.85 24.58 1 2.08 1134.67 1231.03
5.4. Estimación de proporción.
unique(muestra$FAMI_ESTRATOVIVIENDA)
## [1] "-" "Estrato 3" "Estrato 2" "Sin Estrato" "Estrato 1"
## [6] "Estrato 5" "Estrato 4" "Estrato 6"
muestra$FAMI_ESTRATOVIVIENDA <- as.factor(muestra$FAMI_ESTRATOVIVIENDA)
(est2=svymean(~FAMI_ESTRATOVIVIENDA,dsgn,deff=TRUE))
## mean SE DEff
## FAMI_ESTRATOVIVIENDA- 0.0420561 0.0081994 1
## FAMI_ESTRATOVIVIENDAEstrato 1 0.0700935 0.0104293 1
## FAMI_ESTRATOVIVIENDAEstrato 2 0.4112150 0.0201006 1
## FAMI_ESTRATOVIVIENDAEstrato 3 0.4158879 0.0201341 1
## FAMI_ESTRATOVIVIENDAEstrato 4 0.0467290 0.0086218 1
## FAMI_ESTRATOVIVIENDAEstrato 5 0.0046729 0.0027859 1
## FAMI_ESTRATOVIVIENDAEstrato 6 0.0023364 0.0019723 1
## FAMI_ESTRATOVIVIENDASin Estrato 0.0070093 0.0034081 1
alpha=0.05
(tabla2=salida(est2,alpha))
## mean se deff cv ic_low ic_upp
## FAMI_ESTRATOVIVIENDA- 0.04 0.01 1 19.50 0.03 0.06
## FAMI_ESTRATOVIVIENDAEstrato 1 0.07 0.01 1 14.88 0.05 0.09
## FAMI_ESTRATOVIVIENDAEstrato 2 0.41 0.02 1 4.89 0.37 0.45
## FAMI_ESTRATOVIVIENDAEstrato 3 0.42 0.02 1 4.84 0.38 0.46
## FAMI_ESTRATOVIVIENDAEstrato 4 0.05 0.01 1 18.45 0.03 0.06
## FAMI_ESTRATOVIVIENDAEstrato 5 0.00 0.00 1 59.62 0.00 0.01
## FAMI_ESTRATOVIVIENDAEstrato 6 0.00 0.00 1 84.41 0.00 0.01
## FAMI_ESTRATOVIVIENDASin Estrato 0.01 0.00 1 48.62 0.00 0.01
5.5. Estimación de más de una variable.
(est3=svymean(~PUNT_GLOBAL+PUNT_LECTURA_CRITICA+PUNT_MATEMATICAS+PUNT_C_NATURALES+
PUNT_SOCIALES_CIUDADANAS+PUNT_INGLES
,dsgn))
## mean SE
## PUNT_GLOBAL 260.103 1.9513
## PUNT_LECTURA_CRITICA 54.526 0.4086
## PUNT_MATEMATICAS 53.395 0.4715
## PUNT_C_NATURALES 50.836 0.3939
## PUNT_SOCIALES_CIUDADANAS 49.100 0.4682
## PUNT_INGLES 52.708 0.4696
alpha=0.05
(tabla3=salida(est3,alpha))
## mean se cv ic_low ic_upp
## PUNT_GLOBAL 260.10 1.95 0.75 256.28 263.93
## PUNT_LECTURA_CRITICA 54.53 0.41 0.75 53.72 55.33
## PUNT_MATEMATICAS 53.39 0.47 0.88 52.47 54.32
## PUNT_C_NATURALES 50.84 0.39 0.77 50.06 51.61
## PUNT_SOCIALES_CIUDADANAS 49.10 0.47 0.95 48.18 50.02
## PUNT_INGLES 52.71 0.47 0.89 51.79 53.63
6. Estimación del tamaño de un dominio usando MAS.
(est4=svyby(~ind,~FAMI_NUMLIBROS,dsgn,svytotal))
## FAMI_NUMLIBROS ind se
## 17.39486 6.535857
## - - 62.62150 12.208871
## 0 A 10 LIBROS 0 A 10 LIBROS 288.75467 24.048939
## 11 A 25 LIBROS 11 A 25 LIBROS 542.71963 29.274810
## 26 A 100 LIBROS 26 A 100 LIBROS 420.95561 27.391083
## MÁS DE 100 LIBROS MÁS DE 100 LIBROS 156.55374 18.657480
est4=as.data.frame(est4)
table(FAMI_NUMLIBROS)
## FAMI_NUMLIBROS
## - 0 A 10 LIBROS 11 A 25 LIBROS
## 5 18 83 156
## 26 A 100 LIBROS MÁS DE 100 LIBROS
## 121 45
(est4=est4[,-1])
## ind se
## 17.39486 6.535857
## - 62.62150 12.208871
## 0 A 10 LIBROS 288.75467 24.048939
## 11 A 25 LIBROS 542.71963 29.274810
## 26 A 100 LIBROS 420.95561 27.391083
## MÁS DE 100 LIBROS 156.55374 18.657480
alpha=0.05
(tabla4=salida(est4,alpha))
## ind se cv ic_low ic_upp
## 17.39 6.54 37.57 4.58 30.20
## - 62.62 12.21 19.50 38.69 86.55
## 0 A 10 LIBROS 288.75 24.05 8.33 241.62 335.89
## 11 A 25 LIBROS 542.72 29.27 5.39 485.34 600.10
## 26 A 100 LIBROS 420.96 27.39 6.51 367.27 474.64
## MÁS DE 100 LIBROS 156.55 18.66 11.92 119.99 193.12
6.1. Estimación del total de una variable asociada a un dominio bajo
el MAS
(est5=svyby(~FAMI_ESTRATOVIVIENDA,~FAMI_NUMLIBROS,dsgn,deff=TRUE,svytotal))
## FAMI_NUMLIBROS FAMI_ESTRATOVIVIENDA-
## 17.394860
## - - 41.747664
## 0 A 10 LIBROS 0 A 10 LIBROS 0.000000
## 11 A 25 LIBROS 11 A 25 LIBROS 3.478972
## 26 A 100 LIBROS 26 A 100 LIBROS 0.000000
## MÁS DE 100 LIBROS MÁS DE 100 LIBROS 0.000000
## FAMI_ESTRATOVIVIENDAEstrato 1 FAMI_ESTRATOVIVIENDAEstrato 2
## 0.000000 0.000000
## - 0.000000 3.478972
## 0 A 10 LIBROS 38.268692 149.595794
## 11 A 25 LIBROS 52.184579 243.528037
## 26 A 100 LIBROS 10.436916 170.469626
## MÁS DE 100 LIBROS 3.478972 45.226636
## FAMI_ESTRATOVIVIENDAEstrato 3 FAMI_ESTRATOVIVIENDAEstrato 4
## 0.00000 0.000000
## - 13.91589 3.478972
## 0 A 10 LIBROS 80.01636 13.915888
## 11 A 25 LIBROS 212.21729 27.831776
## 26 A 100 LIBROS 219.17523 17.394860
## MÁS DE 100 LIBROS 93.93224 6.957944
## FAMI_ESTRATOVIVIENDAEstrato 5 FAMI_ESTRATOVIVIENDAEstrato 6
## 0.000000 0.000000
## - 0.000000 0.000000
## 0 A 10 LIBROS 0.000000 0.000000
## 11 A 25 LIBROS 3.478972 0.000000
## 26 A 100 LIBROS 3.478972 0.000000
## MÁS DE 100 LIBROS 0.000000 3.478972
## FAMI_ESTRATOVIVIENDASin Estrato se.FAMI_ESTRATOVIVIENDA-
## 0.000000 6.535857
## - 0.000000 10.041177
## 0 A 10 LIBROS 6.957944 0.000000
## 11 A 25 LIBROS 0.000000 2.936711
## 26 A 100 LIBROS 0.000000 0.000000
## MÁS DE 100 LIBROS 3.478972 0.000000
## se.FAMI_ESTRATOVIVIENDAEstrato 1
## 0.000000
## - 0.000000
## 0 A 10 LIBROS 9.625243
## 11 A 25 LIBROS 11.185824
## 26 A 100 LIBROS 5.074607
## MÁS DE 100 LIBROS 2.936711
## se.FAMI_ESTRATOVIVIENDAEstrato 2
## 0.000000
## - 2.936711
## 0 A 10 LIBROS 18.285714
## 11 A 25 LIBROS 22.497687
## 26 A 100 LIBROS 19.367117
## MÁS DE 100 LIBROS 10.438619
## se.FAMI_ESTRATOVIVIENDAEstrato 3
## 0.000000
## - 5.852754
## 0 A 10 LIBROS 13.716356
## 11 A 25 LIBROS 21.264016
## 26 A 100 LIBROS 21.550832
## MÁS DE 100 LIBROS 14.787726
## se.FAMI_ESTRATOVIVIENDAEstrato 4
## 0.000000
## - 2.936711
## 0 A 10 LIBROS 5.852754
## 11 A 25 LIBROS 8.237909
## 26 A 100 LIBROS 6.535857
## MÁS DE 100 LIBROS 4.148271
## se.FAMI_ESTRATOVIVIENDAEstrato 5
## 0.000000
## - 0.000000
## 0 A 10 LIBROS 0.000000
## 11 A 25 LIBROS 2.936711
## 26 A 100 LIBROS 2.936711
## MÁS DE 100 LIBROS 0.000000
## se.FAMI_ESTRATOVIVIENDAEstrato 6
## 0.000000
## - 0.000000
## 0 A 10 LIBROS 0.000000
## 11 A 25 LIBROS 0.000000
## 26 A 100 LIBROS 0.000000
## MÁS DE 100 LIBROS 2.936711
## se.FAMI_ESTRATOVIVIENDASin Estrato DEff.FAMI_ESTRATOVIVIENDA-
## 0.000000 Inf
## - 0.000000 2.760343
## 0 A 10 LIBROS 4.148271 NaN
## 11 A 25 LIBROS 0.000000 1.000000
## 26 A 100 LIBROS 0.000000 NaN
## MÁS DE 100 LIBROS 2.936711 NaN
## DEff.FAMI_ESTRATOVIVIENDAEstrato 1
## NaN
## - NaN
## 0 A 10 LIBROS 1.112217
## 11 A 25 LIBROS 1.063248
## 26 A 100 LIBROS 1.012186
## MÁS DE 100 LIBROS 1.000000
## DEff.FAMI_ESTRATOVIVIENDAEstrato 2
## NaN
## - 1.000000
## 0 A 10 LIBROS 1.848361
## 11 A 25 LIBROS 1.511083
## 26 A 100 LIBROS 1.479313
## MÁS DE 100 LIBROS 1.336358
## DEff.FAMI_ESTRATOVIVIENDAEstrato 3
## NaN
## - 1.205754
## 0 A 10 LIBROS 1.296253
## 11 A 25 LIBROS 1.402317
## 26 A 100 LIBROS 1.768554
## MÁS DE 100 LIBROS 2.295602
## DEff.FAMI_ESTRATOVIVIENDAEstrato 4
## NaN
## - 1.000000
## 0 A 10 LIBROS 1.030682
## 11 A 25 LIBROS 1.030128
## 26 A 100 LIBROS 1.024792
## MÁS DE 100 LIBROS 1.020859
## DEff.FAMI_ESTRATOVIVIENDAEstrato 5
## NaN
## - NaN
## 0 A 10 LIBROS NaN
## 11 A 25 LIBROS 1
## 26 A 100 LIBROS 1
## MÁS DE 100 LIBROS NaN
## DEff.FAMI_ESTRATOVIVIENDAEstrato 6
## NaN
## - NaN
## 0 A 10 LIBROS NaN
## 11 A 25 LIBROS NaN
## 26 A 100 LIBROS NaN
## MÁS DE 100 LIBROS 1
## DEff.FAMI_ESTRATOVIVIENDASin Estrato
## NaN
## - NaN
## 0 A 10 LIBROS 1.009975
## 11 A 25 LIBROS NaN
## 26 A 100 LIBROS NaN
## MÁS DE 100 LIBROS 1.000000
est5=as.data.frame(est5)
est5=est5[,-1]
alpha=0.05
(tabla5=salida(est5,alpha))
## FAMI_ESTRATOVIVIENDA- se FAMI_ESTRATOVIVIENDAEstrato 2
## 17.39 0.00 0.00
## - 41.75 0.00 3.48
## 0 A 10 LIBROS 0.00 38.27 149.60
## 11 A 25 LIBROS 3.48 52.18 243.53
## 26 A 100 LIBROS 0.00 10.44 170.47
## MÁS DE 100 LIBROS 0.00 3.48 45.23
## FAMI_ESTRATOVIVIENDAEstrato 3 FAMI_ESTRATOVIVIENDAEstrato 4
## 0.00 0.00
## - 13.92 3.48
## 0 A 10 LIBROS 80.02 13.92
## 11 A 25 LIBROS 212.22 27.83
## 26 A 100 LIBROS 219.18 17.39
## MÁS DE 100 LIBROS 93.93 6.96
## FAMI_ESTRATOVIVIENDAEstrato 5 FAMI_ESTRATOVIVIENDAEstrato 6
## 0.00 0.00
## - 0.00 0.00
## 0 A 10 LIBROS 0.00 0.00
## 11 A 25 LIBROS 3.48 0.00
## 26 A 100 LIBROS 3.48 0.00
## MÁS DE 100 LIBROS 0.00 3.48
## FAMI_ESTRATOVIVIENDASin Estrato se.FAMI_ESTRATOVIVIENDA-
## 0.00 6.54
## - 0.00 10.04
## 0 A 10 LIBROS 6.96 0.00
## 11 A 25 LIBROS 0.00 2.94
## 26 A 100 LIBROS 0.00 0.00
## MÁS DE 100 LIBROS 3.48 0.00
## se.FAMI_ESTRATOVIVIENDAEstrato 1
## 0.00
## - 0.00
## 0 A 10 LIBROS 9.63
## 11 A 25 LIBROS 11.19
## 26 A 100 LIBROS 5.07
## MÁS DE 100 LIBROS 2.94
## se.FAMI_ESTRATOVIVIENDAEstrato 2
## 0.00
## - 2.94
## 0 A 10 LIBROS 18.29
## 11 A 25 LIBROS 22.50
## 26 A 100 LIBROS 19.37
## MÁS DE 100 LIBROS 10.44
## se.FAMI_ESTRATOVIVIENDAEstrato 3
## 0.00
## - 5.85
## 0 A 10 LIBROS 13.72
## 11 A 25 LIBROS 21.26
## 26 A 100 LIBROS 21.55
## MÁS DE 100 LIBROS 14.79
## se.FAMI_ESTRATOVIVIENDAEstrato 4
## 0.00
## - 2.94
## 0 A 10 LIBROS 5.85
## 11 A 25 LIBROS 8.24
## 26 A 100 LIBROS 6.54
## MÁS DE 100 LIBROS 4.15
## se.FAMI_ESTRATOVIVIENDAEstrato 5
## 0.00
## - 0.00
## 0 A 10 LIBROS 0.00
## 11 A 25 LIBROS 2.94
## 26 A 100 LIBROS 2.94
## MÁS DE 100 LIBROS 0.00
## se.FAMI_ESTRATOVIVIENDAEstrato 6
## 0.00
## - 0.00
## 0 A 10 LIBROS 0.00
## 11 A 25 LIBROS 0.00
## 26 A 100 LIBROS 0.00
## MÁS DE 100 LIBROS 2.94
## se.FAMI_ESTRATOVIVIENDASin Estrato DEff.FAMI_ESTRATOVIVIENDA-
## 0.00 Inf
## - 0.00 2.76
## 0 A 10 LIBROS 4.15 NaN
## 11 A 25 LIBROS 0.00 1.00
## 26 A 100 LIBROS 0.00 NaN
## MÁS DE 100 LIBROS 2.94 NaN
## DEff.FAMI_ESTRATOVIVIENDAEstrato 1
## NaN
## - NaN
## 0 A 10 LIBROS 1.11
## 11 A 25 LIBROS 1.06
## 26 A 100 LIBROS 1.01
## MÁS DE 100 LIBROS 1.00
## DEff.FAMI_ESTRATOVIVIENDAEstrato 2
## NaN
## - 1.00
## 0 A 10 LIBROS 1.85
## 11 A 25 LIBROS 1.51
## 26 A 100 LIBROS 1.48
## MÁS DE 100 LIBROS 1.34
## DEff.FAMI_ESTRATOVIVIENDAEstrato 3
## NaN
## - 1.21
## 0 A 10 LIBROS 1.30
## 11 A 25 LIBROS 1.40
## 26 A 100 LIBROS 1.77
## MÁS DE 100 LIBROS 2.30
## DEff.FAMI_ESTRATOVIVIENDAEstrato 4
## NaN
## - 1.00
## 0 A 10 LIBROS 1.03
## 11 A 25 LIBROS 1.03
## 26 A 100 LIBROS 1.02
## MÁS DE 100 LIBROS 1.02
## DEff.FAMI_ESTRATOVIVIENDAEstrato 5
## NaN
## - NaN
## 0 A 10 LIBROS NaN
## 11 A 25 LIBROS 1
## 26 A 100 LIBROS 1
## MÁS DE 100 LIBROS NaN
## DEff.FAMI_ESTRATOVIVIENDAEstrato 6
## NaN
## - NaN
## 0 A 10 LIBROS NaN
## 11 A 25 LIBROS NaN
## 26 A 100 LIBROS NaN
## MÁS DE 100 LIBROS 1
## DEff.FAMI_ESTRATOVIVIENDASin Estrato cv ic_low ic_upp
## NaN 0 17.39 17.39
## - NaN 0 41.75 41.75
## 0 A 10 LIBROS 1.01 Inf -75.01 75.01
## 11 A 25 LIBROS NaN 1500 -98.80 105.76
## 26 A 100 LIBROS NaN Inf -20.46 20.46
## MÁS DE 100 LIBROS 1.00 Inf -6.82 6.82
6.2. Estimación de la media de una variable asociada a un dominio
bajo el MAS.
(est6=svyby(~PUNT_MATEMATICAS,~FAMI_NUMLIBROS,dsgn,deff=T,svymean))
## FAMI_NUMLIBROS PUNT_MATEMATICAS se
## 53.80000 4.4166709
## - - 52.05556 2.1866033
## 0 A 10 LIBROS 0 A 10 LIBROS 49.95181 1.0283387
## 11 A 25 LIBROS 11 A 25 LIBROS 52.07692 0.7477841
## 26 A 100 LIBROS 26 A 100 LIBROS 56.68595 0.7968431
## MÁS DE 100 LIBROS MÁS DE 100 LIBROS 55.95556 1.7808682
## DEff.PUNT_MATEMATICAS
## 0.8018735
## - 0.9466563
## 0 A 10 LIBROS 0.9902655
## 11 A 25 LIBROS 0.9959167
## 26 A 100 LIBROS 0.9940581
## MÁS DE 100 LIBROS 0.9800677
est6=as.data.frame(est6)
est6=est6[,-1]
alpha=0.05
(tabla6=salida(est6,alpha))
## PUNT_MATEMATICAS se DEff.PUNT_MATEMATICAS cv ic_low
## 53.80 4.42 0.80 8.21 45.14
## - 52.06 2.19 0.95 4.20 47.77
## 0 A 10 LIBROS 49.95 1.03 0.99 2.06 47.94
## 11 A 25 LIBROS 52.08 0.75 1.00 1.44 50.61
## 26 A 100 LIBROS 56.69 0.80 0.99 1.41 55.12
## MÁS DE 100 LIBROS 55.96 1.78 0.98 3.18 52.47
## ic_upp
## 62.46
## - 56.34
## 0 A 10 LIBROS 51.97
## 11 A 25 LIBROS 53.54
## 26 A 100 LIBROS 58.25
## MÁS DE 100 LIBROS 59.45
(est6=svyby(~PUNT_MATEMATICAS,~FAMI_ESTRATOVIVIENDA,dsgn,deff=T,svymean))
## FAMI_ESTRATOVIVIENDA PUNT_MATEMATICAS se
## - - 52.05556 2.3601719
## Estrato 1 Estrato 1 51.13333 1.6597915
## Estrato 2 Estrato 2 51.80114 0.7055901
## Estrato 3 Estrato 3 55.98315 0.7111177
## Estrato 4 Estrato 4 50.90000 2.1545584
## Estrato 5 Estrato 5 55.00000 9.5614394
## Estrato 6 Estrato 6 76.00000 0.0000000
## Sin Estrato Sin Estrato 32.00000 1.5935732
## DEff.PUNT_MATEMATICAS
## - 0.9466563
## Estrato 1 0.9689305
## Estrato 2 0.9966468
## Estrato 3 0.9967108
## Estrato 4 0.9522248
## Estrato 5 0.5011710
## Estrato 6 NA
## Sin Estrato 0.6682279
est6=as.data.frame(est6)
est6=est6[,-1]
alpha=0.05
(tabla6=salida(est6,alpha))
## PUNT_MATEMATICAS se DEff.PUNT_MATEMATICAS cv ic_low ic_upp
## - 52.06 2.36 0.95 4.53 47.43 56.68
## Estrato 1 51.13 1.66 0.97 3.25 47.88 54.39
## Estrato 2 51.80 0.71 1.00 1.36 50.42 53.18
## Estrato 3 55.98 0.71 1.00 1.27 54.59 57.38
## Estrato 4 50.90 2.15 0.95 4.23 46.68 55.12
## Estrato 5 55.00 9.56 0.50 17.38 36.26 73.74
## Estrato 6 76.00 0.00 NA 0.00 76.00 76.00
## Sin Estrato 32.00 1.59 0.67 4.98 28.88 35.12
(est6=svyby(~academico,~FAMI_NUMLIBROS,dsgn,deff=T,svymean))
## FAMI_NUMLIBROS academico se DEff.academico
## 0.8000000 0.15117963 0.8018735
## - - 0.8333333 0.07423621 0.9466563
## 0 A 10 LIBROS 0 A 10 LIBROS 0.8554217 0.03262278 0.9902655
## 11 A 25 LIBROS 11 A 25 LIBROS 0.8141026 0.02632281 0.9959167
## 26 A 100 LIBROS 26 A 100 LIBROS 0.7438017 0.03353840 0.9940581
## MÁS DE 100 LIBROS MÁS DE 100 LIBROS 0.7333333 0.05571179 0.9800677
est6=as.data.frame(est6)
est6=est6[,-1]
alpha=0.05
(tabla6=salida(est6,alpha))
## academico se DEff.academico cv ic_low ic_upp
## 0.80 0.15 0.80 18.90 0.50 1.10
## - 0.83 0.07 0.95 8.91 0.69 0.98
## 0 A 10 LIBROS 0.86 0.03 0.99 3.81 0.79 0.92
## 11 A 25 LIBROS 0.81 0.03 1.00 3.23 0.76 0.87
## 26 A 100 LIBROS 0.74 0.03 0.99 4.51 0.68 0.81
## MÁS DE 100 LIBROS 0.73 0.06 0.98 7.60 0.62 0.84
6.3. Estimación de proporción de una variable asociada a un dominio
bajo el MAS.
salida.p=function(est){
for(i in 1:nrow(est)){
a=as.data.frame(t(est[i,]))
categoria=rownames(a)[1:(nrow(a)/2)]
dominio=rep(rownames(est)[i],(nrow(a)/2))
p=a[1:(nrow(a)/2),]
se=a[(nrow(a)/2+1):nrow(a),]
tab=data.frame(dominio,categoria,p,se)
ifelse(i==1,(tab2=tab),(tab2=rbind(tab2,tab)))
}
return(tab2)
}
(est7=svyby(~academico,~FAMI_NUMLIBROS,dsgn,deff=T,svymean))
## FAMI_NUMLIBROS academico se DEff.academico
## 0.8000000 0.15117963 0.8018735
## - - 0.8333333 0.07423621 0.9466563
## 0 A 10 LIBROS 0 A 10 LIBROS 0.8554217 0.03262278 0.9902655
## 11 A 25 LIBROS 11 A 25 LIBROS 0.8141026 0.02632281 0.9959167
## 26 A 100 LIBROS 26 A 100 LIBROS 0.7438017 0.03353840 0.9940581
## MÁS DE 100 LIBROS MÁS DE 100 LIBROS 0.7333333 0.05571179 0.9800677
est7=as.data.frame(est7)
est7=est7[,-1]
(tabla7=salida.p(est7))
## dominio categoria p se
## 1 academico 0.8000000 0.15117963
## 2 - academico 0.8333333 0.07423621
## 3 0 A 10 LIBROS academico 0.8554217 0.03262278
## 4 11 A 25 LIBROS academico 0.8141026 0.02632281
## 5 26 A 100 LIBROS academico 0.7438017 0.03353840
## 6 MÁS DE 100 LIBROS academico 0.7333333 0.05571179