Base de datos y libreria
setwd("C:/Users/MariaLourdes/Desktop/QUINTO SEMESTRE/ECONOMETRIA I/PROYECTO FINAL")
library("foreign", lib.loc="C:/Program Files/R/R-3.2.2/library")
graduandos<-read.spss("GRADUANDOS2015.sav", to.data.frame = TRUE)
Variables categoricas volverlas Factor()
- Arreglo de todas las variables
graduandos$GENERO<-as.factor(graduandos$GENERO)
graduandos$AREA<-as.factor(graduandos$AREA)
graduandos$TEC_RECIBIDO_CURSO_COMPUTACION<-as.factor(graduandos$TEC_RECIBIDO_CURSO_COMPUTACION)
graduandos$FM_ASISTIO_ESCUELA_PAPA<-as.factor(graduandos$FM_ASISTIO_ESCUELA_PAPA)
graduandos$FM_ASISTIO_ESCUELA_MAMA<-as.factor(graduandos$FM_ASISTIO_ESCUELA_MAMA)
graduandos$ED_TRABAJA_ACTUALMENTE<-as.factor(graduandos$ED_TRABAJA_ACTUALMENTE)
graduandos$CC_ELECTRICIDAD<-as.factor(graduandos$CC_ELECTRICIDAD)
graduandos$CC_SERVICIO_INTERNET<-as.factor(graduandos$CC_SERVICIO_INTERNET)
graduandos$CC_FAMILIA_AUTOMOVIL<-as.factor(graduandos$CC_FAMILIA_AUTOMOVIL)
graduandos$TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOS<-as.factor(graduandos$TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOS)
graduandos$MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada<-as.factor(graduandos$MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada)
graduandos$LEC_PERIODOS_LECTURA_SEMANA_Recodificada<-as.factor(graduandos$LEC_PERIODOS_LECTURA_SEMANA_Recodificada)
graduandos$ED_REPITIO_ALGUN_GRADO_PRIMARIA<-as.factor(graduandos$ED_REPITIO_ALGUN_GRADO_PRIMARIA)
graduandos$IE_IDENTIFICACION_ETNICA_Recodificada<-as.factor(graduandos$IE_IDENTIFICACION_ETNICA_Recodificada)
graduandos$IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA<-as.factor(graduandos$IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA)
graduandos$EDAD_Recodificada<-as.factor(graduandos$EDAD_Recodificada)
regresion logistica primera
mylogit1<-glm(as.factor(LOGRO_MATEMÁTICAS)~GENERO+AREA+TEC_RECIBIDO_CURSO_COMPUTACION+CC_ELECTRICIDAD+FM_ASISTIO_ESCUELA_PAPA+FM_ASISTIO_ESCUELA_MAMA+ED_TRABAJA_ACTUALMENTE+IE_IDENTIFICACION_ETNICA_Recodificada+CC_SERVICIO_INTERNET+CC_FAMILIA_AUTOMOVIL+ED_REPITIO_ALGUN_GRADO_PRIMARIA+IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA+MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada+EDAD_Recodificada, data=graduandos,family=binomial(link = "probit"))
summary(mylogit1)
##
## Call:
## glm(formula = as.factor(LOGRO_MATEMÁTICAS) ~ GENERO + AREA +
## TEC_RECIBIDO_CURSO_COMPUTACION + CC_ELECTRICIDAD + FM_ASISTIO_ESCUELA_PAPA +
## FM_ASISTIO_ESCUELA_MAMA + ED_TRABAJA_ACTUALMENTE + IE_IDENTIFICACION_ETNICA_Recodificada +
## CC_SERVICIO_INTERNET + CC_FAMILIA_AUTOMOVIL + ED_REPITIO_ALGUN_GRADO_PRIMARIA +
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA + MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada +
## EDAD_Recodificada, family = binomial(link = "probit"), data = graduandos)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.3754 -0.4411 -0.2510 -0.1225 4.2372
##
## Coefficients:
## Estimate
## (Intercept) -3.358588
## GENEROMASCULINO 0.463674
## AREAURBANA 0.173504
## TEC_RECIBIDO_CURSO_COMPUTACIONSI 0.026492
## CC_ELECTRICIDADSI 0.407582
## FM_ASISTIO_ESCUELA_PAPASI 0.231932
## FM_ASISTIO_ESCUELA_PAPANO SABE 0.192867
## FM_ASISTIO_ESCUELA_MAMASI 0.187810
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.066583
## ED_TRABAJA_ACTUALMENTESI -0.197680
## IE_IDENTIFICACION_ETNICA_RecodificadaMAYA 0.253269
## IE_IDENTIFICACION_ETNICA_RecodificadaLADINO 0.505517
## IE_IDENTIFICACION_ETNICA_RecodificadaGARÍFUNA 0.075602
## IE_IDENTIFICACION_ETNICA_RecodificadaXINKA -0.100274
## IE_IDENTIFICACION_ETNICA_RecodificadaEXTRANJERO 0.482264
## IE_IDENTIFICACION_ETNICA_RecodificadaDOBLE O MÁS RESPUESTAS -0.007129
## CC_SERVICIO_INTERNETSI 0.414153
## CC_FAMILIA_AUTOMOVILSI 0.085663
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI -0.540373
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 0.431400
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO -0.040376
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.100768
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 0.168209
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 0.218119
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 0.412244
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 0.545917
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 0.695474
## EDAD_Recodificada17 AÑOS 0.019037
## EDAD_Recodificada18 AÑOS -0.130046
## EDAD_Recodificada19 AÑOS -0.291279
## EDAD_Recodificada20 AÑOS -0.414918
## EDAD_RecodificadaMAYOR DE 20 AÑOS -0.356941
## Std. Error
## (Intercept) 0.113147
## GENEROMASCULINO 0.011487
## AREAURBANA 0.020188
## TEC_RECIBIDO_CURSO_COMPUTACIONSI 0.012394
## CC_ELECTRICIDADSI 0.065074
## FM_ASISTIO_ESCUELA_PAPASI 0.025344
## FM_ASISTIO_ESCUELA_PAPANO SABE 0.033882
## FM_ASISTIO_ESCUELA_MAMASI 0.020082
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.039981
## ED_TRABAJA_ACTUALMENTESI 0.014376
## IE_IDENTIFICACION_ETNICA_RecodificadaMAYA 0.085960
## IE_IDENTIFICACION_ETNICA_RecodificadaLADINO 0.085180
## IE_IDENTIFICACION_ETNICA_RecodificadaGARÍFUNA 0.127355
## IE_IDENTIFICACION_ETNICA_RecodificadaXINKA 0.169688
## IE_IDENTIFICACION_ETNICA_RecodificadaEXTRANJERO 0.107527
## IE_IDENTIFICACION_ETNICA_RecodificadaDOBLE O MÁS RESPUESTAS 0.280399
## CC_SERVICIO_INTERNETSI 0.012309
## CC_FAMILIA_AUTOMOVILSI 0.012280
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI 0.018176
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 0.012198
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO 0.024548
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.026147
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 0.025413
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 0.021890
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 0.019547
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 0.025628
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 0.026164
## EDAD_Recodificada17 AÑOS 0.018360
## EDAD_Recodificada18 AÑOS 0.019267
## EDAD_Recodificada19 AÑOS 0.024549
## EDAD_Recodificada20 AÑOS 0.033714
## EDAD_RecodificadaMAYOR DE 20 AÑOS 0.028727
## z value
## (Intercept) -29.683
## GENEROMASCULINO 40.367
## AREAURBANA 8.594
## TEC_RECIBIDO_CURSO_COMPUTACIONSI 2.137
## CC_ELECTRICIDADSI 6.263
## FM_ASISTIO_ESCUELA_PAPASI 9.151
## FM_ASISTIO_ESCUELA_PAPANO SABE 5.692
## FM_ASISTIO_ESCUELA_MAMASI 9.352
## FM_ASISTIO_ESCUELA_MAMANO SABE 1.665
## ED_TRABAJA_ACTUALMENTESI -13.751
## IE_IDENTIFICACION_ETNICA_RecodificadaMAYA 2.946
## IE_IDENTIFICACION_ETNICA_RecodificadaLADINO 5.935
## IE_IDENTIFICACION_ETNICA_RecodificadaGARÍFUNA 0.594
## IE_IDENTIFICACION_ETNICA_RecodificadaXINKA -0.591
## IE_IDENTIFICACION_ETNICA_RecodificadaEXTRANJERO 4.485
## IE_IDENTIFICACION_ETNICA_RecodificadaDOBLE O MÁS RESPUESTAS -0.025
## CC_SERVICIO_INTERNETSI 33.646
## CC_FAMILIA_AUTOMOVILSI 6.976
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI -29.730
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 35.368
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO -1.645
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 3.854
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 6.619
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 9.964
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 21.090
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 21.301
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 26.581
## EDAD_Recodificada17 AÑOS 1.037
## EDAD_Recodificada18 AÑOS -6.750
## EDAD_Recodificada19 AÑOS -11.865
## EDAD_Recodificada20 AÑOS -12.307
## EDAD_RecodificadaMAYOR DE 20 AÑOS -12.425
## Pr(>|z|)
## (Intercept) < 2e-16 ***
## GENEROMASCULINO < 2e-16 ***
## AREAURBANA < 2e-16 ***
## TEC_RECIBIDO_CURSO_COMPUTACIONSI 0.032561 *
## CC_ELECTRICIDADSI 3.77e-10 ***
## FM_ASISTIO_ESCUELA_PAPASI < 2e-16 ***
## FM_ASISTIO_ESCUELA_PAPANO SABE 1.25e-08 ***
## FM_ASISTIO_ESCUELA_MAMASI < 2e-16 ***
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.095838 .
## ED_TRABAJA_ACTUALMENTESI < 2e-16 ***
## IE_IDENTIFICACION_ETNICA_RecodificadaMAYA 0.003215 **
## IE_IDENTIFICACION_ETNICA_RecodificadaLADINO 2.94e-09 ***
## IE_IDENTIFICACION_ETNICA_RecodificadaGARÍFUNA 0.552762
## IE_IDENTIFICACION_ETNICA_RecodificadaXINKA 0.554565
## IE_IDENTIFICACION_ETNICA_RecodificadaEXTRANJERO 7.29e-06 ***
## IE_IDENTIFICACION_ETNICA_RecodificadaDOBLE O MÁS RESPUESTAS 0.979715
## CC_SERVICIO_INTERNETSI < 2e-16 ***
## CC_FAMILIA_AUTOMOVILSI 3.04e-12 ***
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI < 2e-16 ***
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO 0.100009
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.000116 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 3.61e-11 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS < 2e-16 ***
## EDAD_Recodificada17 AÑOS 0.299806
## EDAD_Recodificada18 AÑOS 1.48e-11 ***
## EDAD_Recodificada19 AÑOS < 2e-16 ***
## EDAD_Recodificada20 AÑOS < 2e-16 ***
## EDAD_RecodificadaMAYOR DE 20 AÑOS < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 82896 on 137497 degrees of freedom
## Residual deviance: 65229 on 137466 degrees of freedom
## (12154 observations deleted due to missingness)
## AIC: 65293
##
## Number of Fisher Scoring iterations: 7
Regresion logistica final
mylogit<-glm(as.factor(LOGRO_MATEMÁTICAS)~GENERO+AREA+FM_ASISTIO_ESCUELA_PAPA+FM_ASISTIO_ESCUELA_MAMA+ED_TRABAJA_ACTUALMENTE+CC_SERVICIO_INTERNET+TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOS+ED_REPITIO_ALGUN_GRADO_PRIMARIA+IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA+MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada+EDAD_Recodificada, data=graduandos,family=binomial(link = "probit"))
summary(mylogit)
##
## Call:
## glm(formula = as.factor(LOGRO_MATEMÁTICAS) ~ GENERO + AREA +
## FM_ASISTIO_ESCUELA_PAPA + FM_ASISTIO_ESCUELA_MAMA + ED_TRABAJA_ACTUALMENTE +
## CC_SERVICIO_INTERNET + TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOS +
## ED_REPITIO_ALGUN_GRADO_PRIMARIA + IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA +
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada + EDAD_Recodificada,
## family = binomial(link = "probit"), data = graduandos)
##
## Deviance Residuals:
## Min 1Q Median 3Q Max
## -1.3271 -0.4426 -0.2502 -0.1184 4.0436
##
## Coefficients:
## Estimate
## (Intercept) -2.78385
## GENEROMASCULINO 0.45955
## AREAURBANA 0.16849
## FM_ASISTIO_ESCUELA_PAPASI 0.25185
## FM_ASISTIO_ESCUELA_PAPANO SABE 0.21142
## FM_ASISTIO_ESCUELA_MAMASI 0.26190
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.11020
## ED_TRABAJA_ACTUALMENTESI -0.20976
## CC_SERVICIO_INTERNETSI 0.47597
## TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOSSI 0.29289
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI -0.54905
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 0.35185
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO -0.02992
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.09315
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 0.15615
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 0.20393
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 0.40778
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 0.54973
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 0.70889
## EDAD_Recodificada17 AÑOS 0.02840
## EDAD_Recodificada18 AÑOS -0.12853
## EDAD_Recodificada19 AÑOS -0.29795
## EDAD_Recodificada20 AÑOS -0.41396
## EDAD_RecodificadaMAYOR DE 20 AÑOS -0.33767
## Std. Error
## (Intercept) 0.04314
## GENEROMASCULINO 0.01137
## AREAURBANA 0.02005
## FM_ASISTIO_ESCUELA_PAPASI 0.02519
## FM_ASISTIO_ESCUELA_PAPANO SABE 0.03372
## FM_ASISTIO_ESCUELA_MAMASI 0.01975
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.03976
## ED_TRABAJA_ACTUALMENTESI 0.01424
## CC_SERVICIO_INTERNETSI 0.01146
## TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOSSI 0.02036
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI 0.01803
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 0.01114
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO 0.02439
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.02599
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 0.02520
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 0.02175
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 0.01943
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 0.02541
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 0.02595
## EDAD_Recodificada17 AÑOS 0.01818
## EDAD_Recodificada18 AÑOS 0.01907
## EDAD_Recodificada19 AÑOS 0.02436
## EDAD_Recodificada20 AÑOS 0.03334
## EDAD_RecodificadaMAYOR DE 20 AÑOS 0.02849
## z value
## (Intercept) -64.528
## GENEROMASCULINO 40.422
## AREAURBANA 8.403
## FM_ASISTIO_ESCUELA_PAPASI 9.999
## FM_ASISTIO_ESCUELA_PAPANO SABE 6.270
## FM_ASISTIO_ESCUELA_MAMASI 13.262
## FM_ASISTIO_ESCUELA_MAMANO SABE 2.772
## ED_TRABAJA_ACTUALMENTESI -14.729
## CC_SERVICIO_INTERNETSI 41.534
## TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOSSI 14.387
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI -30.454
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI 31.590
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO -1.226
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 3.584
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 6.195
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS 9.374
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS 20.984
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS 21.635
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS 27.313
## EDAD_Recodificada17 AÑOS 1.562
## EDAD_Recodificada18 AÑOS -6.742
## EDAD_Recodificada19 AÑOS -12.234
## EDAD_Recodificada20 AÑOS -12.416
## EDAD_RecodificadaMAYOR DE 20 AÑOS -11.850
## Pr(>|z|)
## (Intercept) < 2e-16 ***
## GENEROMASCULINO < 2e-16 ***
## AREAURBANA < 2e-16 ***
## FM_ASISTIO_ESCUELA_PAPASI < 2e-16 ***
## FM_ASISTIO_ESCUELA_PAPANO SABE 3.61e-10 ***
## FM_ASISTIO_ESCUELA_MAMASI < 2e-16 ***
## FM_ASISTIO_ESCUELA_MAMANO SABE 0.005572 **
## ED_TRABAJA_ACTUALMENTESI < 2e-16 ***
## CC_SERVICIO_INTERNETSI < 2e-16 ***
## TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOSSI < 2e-16 ***
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI < 2e-16 ***
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO 0.220030
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS 0.000338 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS 5.82e-10 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS < 2e-16 ***
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS < 2e-16 ***
## EDAD_Recodificada17 AÑOS 0.118307
## EDAD_Recodificada18 AÑOS 1.56e-11 ***
## EDAD_Recodificada19 AÑOS < 2e-16 ***
## EDAD_Recodificada20 AÑOS < 2e-16 ***
## EDAD_RecodificadaMAYOR DE 20 AÑOS < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## (Dispersion parameter for binomial family taken to be 1)
##
## Null deviance: 84013 on 140014 degrees of freedom
## Residual deviance: 66316 on 139991 degrees of freedom
## (9637 observations deleted due to missingness)
## AIC: 66364
##
## Number of Fisher Scoring iterations: 7
Tablas
#genero
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$GENERO),1)
##
## FEMENINO MASCULINO
## NO LOGRO 0.5062375 0.4937625
## LOGRO 0.3540865 0.6459135
#area
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$AREA),1)
##
## RURAL URBANA
## NO LOGRO 0.10430635 0.89569365
## LOGRO 0.06885452 0.93114548
#Asistio PAPA
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$FM_ASISTIO_ESCUELA_PAPA),1)
##
## NO SI NO SABE
## NO LOGRO 0.14337765 0.78452215 0.07210020
## LOGRO 0.03230636 0.92041604 0.04727760
#Asistio MAMA
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$FM_ASISTIO_ESCUELA_MAMA),1)
##
## NO SI NO SABE
## NO LOGRO 0.22365500 0.74059801 0.03574699
## LOGRO 0.06176332 0.91977055 0.01846613
#Trabaja actualmente
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$ED_TRABAJA_ACTUALMENTE),1)
##
## NO SI
## NO LOGRO 0.5961772 0.4038228
## LOGRO 0.8317500 0.1682500
#Si tiene servicio de internet
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$CC_SERVICIO_INTERNET),1)
##
## NO SI
## NO LOGRO 0.7085671 0.2914329
## LOGRO 0.3270984 0.6729016
#Repitio algun grado
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$ED_REPITIO_ALGUN_GRADO_PRIMARIA),1)
##
## NO SI
## NO LOGRO 0.67949224 0.32050776
## LOGRO 0.93507105 0.06492895
#Habla otro idioma
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA),1)
##
## NO SI
## NO LOGRO 0.6918513 0.3081487
## LOGRO 0.5202109 0.4797891
#Periodos de mate
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada),1)
##
## NINGUNO 1 PERÍODO 2 PERÍODOS 3 PERÍODOS 4 PERÍODOS
## NO LOGRO 0.17077143 0.23616285 0.10343719 0.08270885 0.13184919
## LOGRO 0.07207349 0.06006124 0.05370181 0.07144540 0.14485358
##
## 5 PERÍODOS 6 PERÍODOS MÁS DE 6 PERÍODOS
## NO LOGRO 0.20607096 0.03890033 0.03009919
## LOGRO 0.36342938 0.11219282 0.12224229
#Edad del estudiante
prop.table(table(graduandos$LOGRO_MATEMÁTICAS,graduandos$EDAD_Recodificada),1)
##
## 16 AÑOS Ó MENOS 17 AÑOS 18 AÑOS 19 AÑOS 20 AÑOS
## NO LOGRO 0.05866454 0.20955171 0.27639028 0.16119418 0.09144832
## LOGRO 0.14077098 0.44830023 0.28059983 0.07058177 0.02190469
##
## MAYOR DE 20 AÑOS
## NO LOGRO 0.20275097
## LOGRO 0.03784251
Marginal Effects
#marginal effects
coef(mylogit)
## (Intercept)
## -2.78385228
## GENEROMASCULINO
## 0.45955186
## AREAURBANA
## 0.16848726
## FM_ASISTIO_ESCUELA_PAPASI
## 0.25184731
## FM_ASISTIO_ESCUELA_PAPANO SABE
## 0.21141741
## FM_ASISTIO_ESCUELA_MAMASI
## 0.26189958
## FM_ASISTIO_ESCUELA_MAMANO SABE
## 0.11020239
## ED_TRABAJA_ACTUALMENTESI
## -0.20976400
## CC_SERVICIO_INTERNETSI
## 0.47597175
## TEC_POSEEN_COMPUTADORAS_ESTABLECIMIENTO_USO_ALUMNOSSI
## 0.29288691
## ED_REPITIO_ALGUN_GRADO_PRIMARIASI
## -0.54905131
## IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMASI
## 0.35185414
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada1 PERÍODO
## -0.02991512
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada2 PERÍODOS
## 0.09315167
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada3 PERÍODOS
## 0.15614998
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada4 PERÍODOS
## 0.20392617
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada5 PERÍODOS
## 0.40777950
## MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada6 PERÍODOS
## 0.54973225
## MAT_PERIODOS_MATEMATICAS_SEMANA_RecodificadaMÁS DE 6 PERÍODOS
## 0.70888648
## EDAD_Recodificada17 AÑOS
## 0.02839801
## EDAD_Recodificada18 AÑOS
## -0.12853283
## EDAD_Recodificada19 AÑOS
## -0.29795261
## EDAD_Recodificada20 AÑOS
## -0.41395605
## EDAD_RecodificadaMAYOR DE 20 AÑOS
## -0.33766542
require(effects)
## Loading required package: effects
## Warning: package 'effects' was built under R version 3.2.5
#Plot marginal effects
mod.cowles <- mylogit
eff.cowles <- allEffects(mod.cowles)
# the following are equivalent:
par(mfrow=c(2,2))
plot(eff.cowles, 'GENERO', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'AREA', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'FM_ASISTIO_ESCUELA_PAPA', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'FM_ASISTIO_ESCUELA_MAMA', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'ED_TRABAJA_ACTUALMENTE', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'CC_SERVICIO_INTERNET', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'ED_REPITIO_ALGUN_GRADO_PRIMARIA', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada', ylab="Prob(Logro matematicas)")
plot(eff.cowles, 'EDAD_Recodificada', ylab="Prob(Logro matematicas)")
significancia de las variables
Pseudo R-squared
#calculo R2
probit0<-update(mylogit, formula=as.factor(LOGRO_MATEMÁTICAS)~1)
McFadden<-1- as.vector(logLik(mylogit)/logLik(probit0))
McFadden
## [1] 0.2388059
#Valores entre 0.2 a 0.4 estan considerados altamente satisfactorios.
Multicolinealidad
#Explicar por que no se aplica
Heteroscedasticidad
require(glmx)
ah1<-hetglm(as.factor(LOGRO_MATEMÁTICAS)~GENERO+AREA+FM_ASISTIO_ESCUELA_PAPA+FM_ASISTIO_ESCUELA_MAMA+ED_TRABAJA_ACTUALMENTE+CC_SERVICIO_INTERNET+ED_REPITIO_ALGUN_GRADO_PRIMARIA+IE_DOMINA_OTRO_COMO_SEGUNDO_IDIOMA+MAT_PERIODOS_MATEMATICAS_SEMANA_Recodificada+ EDAD_Recodificada, data=graduandos)
#Agregar edad
summary(ah1)
#Encontrar como se interpreta
- Interpretar bien
- Log-likelihood: -3.264e+04 on 45 Df
- LR test for homoskedasticity: 1341 on 22 Df, p-value: < 2.2e-16
- Dispersion: 1
- Number of iterations in nlminb optimization: 9