library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(readxl)

datos <- read_excel("Prueba_Saber_11.xlsx")

datos <- datos %>%
  mutate(
    GENERO = as.factor(ESTU_GENERO),
    SI_TIENE_ETNIA = as.factor(ESTU_TIENEETNIA),
    DEPARTAMENTO_DE_RESIDENCIA = as.factor(ESTU_DEPTO_RESIDE),
    ESTRATO_DE_LA_VIVIENDA = as.factor(FAMI_ESTRATOVIVIENDA),
    TIPO_DE_COLEGIO = as.factor(COLE_NATURALEZA),
    CALENDARIO = as.factor(COLE_CALENDARIO)
  )
head(datos)
## # A tibble: 6 × 31
##   ESTU_GENERO PERIODO ESTU_TIENEETNIA ESTU_DEPTO_RESIDE ESTU_COD_RESIDE_DEPTO
##   <chr>         <dbl> <chr>           <chr>                             <dbl>
## 1 M             20194 No              "MAGDALENA"                          47
## 2 M             20194 No              "BOGOTÃ\u0081"                       11
## 3 M             20194 No              "BOLIVAR"                            13
## 4 M             20194 Si              "BOGOTÃ\u0081"                       11
## 5 M             20194 No              "BOGOTÃ\u0081"                       11
## 6 M             20194 No              "ATLANTICO"                           8
## # ℹ 26 more variables: FAMI_ESTRATOVIVIENDA <chr>, COLE_NATURALEZA <chr>,
## #   COLE_CALENDARIO <chr>, PUNT_LECTURA_CRITICA <dbl>,
## #   PERCENTIL_LECTURA_CRITICA <dbl>, DESEMP_LECTURA_CRITICA <dbl>,
## #   PUNT_MATEMATICAS <dbl>, PERCENTIL_MATEMATICAS <dbl>,
## #   DESEMP_MATEMATICAS <dbl>, PUNT_C_NATURALES <dbl>,
## #   PERCENTIL_C_NATURALES <dbl>, DESEMP_C_NATURALES <dbl>,
## #   PUNT_SOCIALES_CIUDADANAS <dbl>, PERCENTIL_SOCIALES_CIUDADANAS <dbl>, …
str(datos)
## tibble [429 × 31] (S3: tbl_df/tbl/data.frame)
##  $ ESTU_GENERO                  : chr [1:429] "M" "M" "M" "M" ...
##  $ PERIODO                      : num [1:429] 20194 20194 20194 20194 20194 ...
##  $ ESTU_TIENEETNIA              : chr [1:429] "No" "No" "No" "Si" ...
##  $ ESTU_DEPTO_RESIDE            : chr [1:429] "MAGDALENA" "BOGOTÃ\u0081" "BOLIVAR" "BOGOTÃ\u0081" ...
##  $ ESTU_COD_RESIDE_DEPTO        : num [1:429] 47 11 13 11 11 8 76 68 25 70 ...
##  $ FAMI_ESTRATOVIVIENDA         : chr [1:429] "Estrato 3" "Estrato 3" "Estrato 1" "Estrato 3" ...
##  $ COLE_NATURALEZA              : chr [1:429] "OFICIAL" "NO OFICIAL" "OFICIAL" "OFICIAL" ...
##  $ COLE_CALENDARIO              : chr [1:429] "A" "A" "A" "A" ...
##  $ PUNT_LECTURA_CRITICA         : num [1:429] 47 60 66 62 63 49 76 57 62 68 ...
##  $ PERCENTIL_LECTURA_CRITICA    : num [1:429] 33 76 91 81 85 37 100 67 83 94 ...
##  $ DESEMP_LECTURA_CRITICA       : num [1:429] 2 3 4 3 3 2 4 3 3 4 ...
##  $ PUNT_MATEMATICAS             : num [1:429] 48 65 57 54 57 29 70 65 62 66 ...
##  $ PERCENTIL_MATEMATICAS        : num [1:429] 42 88 70 60 70 4 96 88 82 91 ...
##  $ DESEMP_MATEMATICAS           : num [1:429] 2 3 3 3 3 1 3 3 3 3 ...
##  $ PUNT_C_NATURALES             : num [1:429] 37 54 41 61 55 41 70 63 66 63 ...
##  $ PERCENTIL_C_NATURALES        : num [1:429] 15 71 29 87 74 26 98 90 95 91 ...
##  $ DESEMP_C_NATURALES           : num [1:429] 1 2 2 3 2 2 3 3 3 3 ...
##  $ PUNT_SOCIALES_CIUDADANAS     : num [1:429] 30 59 74 73 57 41 68 66 39 77 ...
##  $ PERCENTIL_SOCIALES_CIUDADANAS: num [1:429] 8 83 99 99 79 39 96 95 33 100 ...
##  $ DESEMP_SOCIALES_CIUDADANAS   : num [1:429] 1 3 4 4 3 2 3 3 1 4 ...
##  $ PUNT_INGLES                  : num [1:429] 54 63 64 53 52 35 72 60 63 51 ...
##  $ PERCENTIL_INGLES             : num [1:429] 67 88 89 67 63 14 96 83 88 61 ...
##  $ DESEMP_INGLES                : chr [1:429] "A1" "A2" "A2" "A1" ...
##  $ PUNT_GLOBAL                  : num [1:429] 208 299 299 309 288 198 355 313 288 336 ...
##  $ PERCENTIL_GLOBAL             : num [1:429] 25 83 84 88 78 19 99 89 78 95 ...
##  $ GENERO                       : Factor w/ 2 levels "F","M": 2 2 2 2 2 2 2 2 2 2 ...
##  $ SI_TIENE_ETNIA               : Factor w/ 2 levels "No","Si": 1 1 1 2 1 1 2 1 1 1 ...
##  $ DEPARTAMENTO_DE_RESIDENCIA   : Factor w/ 30 levels "ANTIOQUIA","ARAUCA",..: 17 4 5 4 4 3 28 25 14 26 ...
##  $ ESTRATO_DE_LA_VIVIENDA       : Factor w/ 7 levels "-","Estrato 1",..: 4 4 2 4 4 1 5 4 4 4 ...
##  $ TIPO_DE_COLEGIO              : Factor w/ 2 levels "NO OFICIAL","OFICIAL": 2 1 2 2 2 1 1 2 2 1 ...
##  $ CALENDARIO                   : Factor w/ 2 levels "A","OTRO": 1 1 1 1 1 1 1 1 1 1 ...
summary(datos)
##  ESTU_GENERO           PERIODO      ESTU_TIENEETNIA    ESTU_DEPTO_RESIDE 
##  Length:429         Min.   :20194   Length:429         Length:429        
##  Class :character   1st Qu.:20194   Class :character   Class :character  
##  Mode  :character   Median :20194   Mode  :character   Mode  :character  
##                     Mean   :20194                                        
##                     3rd Qu.:20194                                        
##                     Max.   :20194                                        
##                                                                          
##  ESTU_COD_RESIDE_DEPTO FAMI_ESTRATOVIVIENDA COLE_NATURALEZA   
##  Min.   : 5.00         Length:429           Length:429        
##  1st Qu.:11.00         Class :character     Class :character  
##  Median :20.00         Mode  :character     Mode  :character  
##  Mean   :29.44                                                
##  3rd Qu.:47.00                                                
##  Max.   :99.00                                                
##                                                               
##  COLE_CALENDARIO    PUNT_LECTURA_CRITICA PERCENTIL_LECTURA_CRITICA
##  Length:429         Min.   : 0.00        Min.   :  1.0            
##  Class :character   1st Qu.:45.00        1st Qu.: 26.0            
##  Mode  :character   Median :54.00        Median : 54.0            
##                     Mean   :52.78        Mean   : 52.4            
##                     3rd Qu.:61.00        3rd Qu.: 80.0            
##                     Max.   :79.00        Max.   :100.0            
##                                                                   
##  DESEMP_LECTURA_CRITICA PUNT_MATEMATICAS PERCENTIL_MATEMATICAS
##  Min.   :1.000          Min.   : 18.00   Min.   :  1.00       
##  1st Qu.:2.000          1st Qu.: 44.00   1st Qu.: 29.00       
##  Median :3.000          Median : 52.00   Median : 55.00       
##  Mean   :2.648          Mean   : 51.77   Mean   : 53.33       
##  3rd Qu.:3.000          3rd Qu.: 61.00   3rd Qu.: 80.00       
##  Max.   :4.000          Max.   :100.00   Max.   :100.00       
##                                                               
##  DESEMP_MATEMATICAS PUNT_C_NATURALES PERCENTIL_C_NATURALES DESEMP_C_NATURALES
##  Min.   :1.000      Min.   :22.00    Min.   :  1.00        Min.   :1.000     
##  1st Qu.:2.000      1st Qu.:41.00    1st Qu.: 27.00        1st Qu.:2.000     
##  Median :3.000      Median :49.00    Median : 54.00        Median :2.000     
##  Mean   :2.517      Mean   :49.12    Mean   : 52.74        Mean   :2.096     
##  3rd Qu.:3.000      3rd Qu.:58.00    3rd Qu.: 80.00        3rd Qu.:3.000     
##  Max.   :4.000      Max.   :79.00    Max.   :100.00        Max.   :4.000     
##                                                                              
##  PUNT_SOCIALES_CIUDADANAS PERCENTIL_SOCIALES_CIUDADANAS
##  Min.   :20.00            Min.   :  1.00               
##  1st Qu.:36.00            1st Qu.: 23.00               
##  Median :46.00            Median : 53.00               
##  Mean   :46.76            Mean   : 51.23               
##  3rd Qu.:56.00            3rd Qu.: 78.00               
##  Max.   :81.00            Max.   :100.00               
##                                                        
##  DESEMP_SOCIALES_CIUDADANAS  PUNT_INGLES     PERCENTIL_INGLES
##  Min.   :1.000              Min.   : 23.00   Min.   :  1.00  
##  1st Qu.:1.000              1st Qu.: 37.00   1st Qu.: 21.00  
##  Median :2.000              Median : 48.00   Median : 51.00  
##  Mean   :1.946              Mean   : 48.46   Mean   : 49.92  
##  3rd Qu.:3.000              3rd Qu.: 58.00   3rd Qu.: 79.00  
##  Max.   :4.000              Max.   :100.00   Max.   :100.00  
##                                                              
##  DESEMP_INGLES       PUNT_GLOBAL    PERCENTIL_GLOBAL GENERO  SI_TIENE_ETNIA
##  Length:429         Min.   : 80.0   Min.   :  1.00   F: 79   No:376        
##  Class :character   1st Qu.:208.0   1st Qu.: 25.00   M:350   Si: 53        
##  Mode  :character   Median :248.0   Median : 53.00                         
##                     Mean   :249.9   Mean   : 52.41                         
##                     3rd Qu.:292.0   3rd Qu.: 80.00                         
##                     Max.   :399.0   Max.   :100.00                         
##                                                                            
##  DEPARTAMENTO_DE_RESIDENCIA ESTRATO_DE_LA_VIVIENDA   TIPO_DE_COLEGIO CALENDARIO
##  ATLANTICO   : 80           -          : 38        NO OFICIAL:108    A   :428  
##  BOLIVAR     : 53           Estrato 1  :139        OFICIAL   :321    OTRO:  1  
##  CORDOBA     : 39           Estrato 2  :109                                    
##  BOGOTÃ\u0081: 32           Estrato 3  : 80                                    
##  MAGDALENA   : 26           Estrato 4  : 31                                    
##  ANTIOQUIA   : 25           Estrato 5  : 12                                    
##  (Other)     :174           Sin Estrato: 20
ggplot(datos, aes(x = PUNT_GLOBAL)) +
  geom_histogram(binwidth = 10, fill = "blue", color = "black") +
  labs(title = "Distribución del Puntaje Global", x = "Puntaje Global", y = "Frecuencia")

ggplot(datos, aes(x = GENERO, y = PUNT_GLOBAL, fill = GENERO)) +
  geom_boxplot() +
  labs(title = "Comparación del Puntaje Global por Género", x = "Género", y = "Puntaje Global")

ggplot(datos, aes(x = TIPO_DE_COLEGIO, y = PUNT_GLOBAL, fill = TIPO_DE_COLEGIO)) +
  geom_boxplot() +
  labs(title = "Comparación del Puntaje Global por Tipo de Colegio", x = "Tipo de Colegio", y = "Puntaje Global")

modelo <- lm(
  PUNT_GLOBAL ~ GENERO + SI_TIENE_ETNIA + DEPARTAMENTO_DE_RESIDENCIA +
  ESTRATO_DE_LA_VIVIENDA + TIPO_DE_COLEGIO + CALENDARIO +
  PUNT_LECTURA_CRITICA + PUNT_MATEMATICAS +
  PUNT_C_NATURALES + PUNT_SOCIALES_CIUDADANAS +
  PUNT_INGLES, data = datos
)

summary(modelo)
## 
## Call:
## lm(formula = PUNT_GLOBAL ~ GENERO + SI_TIENE_ETNIA + DEPARTAMENTO_DE_RESIDENCIA + 
##     ESTRATO_DE_LA_VIVIENDA + TIPO_DE_COLEGIO + CALENDARIO + PUNT_LECTURA_CRITICA + 
##     PUNT_MATEMATICAS + PUNT_C_NATURALES + PUNT_SOCIALES_CIUDADANAS + 
##     PUNT_INGLES, data = datos)
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -0.62162 -0.22231  0.00225  0.21677  0.59010 
## 
## Coefficients:
##                                            Estimate Std. Error t value Pr(>|t|)
## (Intercept)                                0.021885   0.110296   0.198   0.8428
## GENEROM                                   -0.044944   0.038398  -1.170   0.2425
## SI_TIENE_ETNIASi                           0.045415   0.055878   0.813   0.4169
## DEPARTAMENTO_DE_RESIDENCIAARAUCA           0.227955   0.178729   1.275   0.2029
## DEPARTAMENTO_DE_RESIDENCIAATLANTICO        0.008486   0.068051   0.125   0.9008
## DEPARTAMENTO_DE_RESIDENCIABOGOTÃ\u0081     0.015534   0.078867   0.197   0.8440
## DEPARTAMENTO_DE_RESIDENCIABOLIVAR          0.042515   0.071219   0.597   0.5509
## DEPARTAMENTO_DE_RESIDENCIABOYACA           0.272890   0.213531   1.278   0.2020
## DEPARTAMENTO_DE_RESIDENCIACALDAS          -0.441704   0.296387  -1.490   0.1370
## DEPARTAMENTO_DE_RESIDENCIACAQUETA         -0.027181   0.295852  -0.092   0.9268
## DEPARTAMENTO_DE_RESIDENCIACASANARE         0.010481   0.177152   0.059   0.9529
## DEPARTAMENTO_DE_RESIDENCIACAUCA           -0.111432   0.111833  -0.996   0.3197
## DEPARTAMENTO_DE_RESIDENCIACESAR            0.012019   0.086932   0.138   0.8901
## DEPARTAMENTO_DE_RESIDENCIACHOCO           -0.094032   0.147990  -0.635   0.5255
## DEPARTAMENTO_DE_RESIDENCIACORDOBA          0.014389   0.076243   0.189   0.8504
## DEPARTAMENTO_DE_RESIDENCIACUNDINAMARCA    -0.094221   0.096944  -0.972   0.3317
## DEPARTAMENTO_DE_RESIDENCIAHUILA           -0.081137   0.110345  -0.735   0.4626
## DEPARTAMENTO_DE_RESIDENCIALA GUAJIRA      -0.084938   0.096352  -0.882   0.3786
## DEPARTAMENTO_DE_RESIDENCIAMAGDALENA       -0.057816   0.081336  -0.711   0.4776
## DEPARTAMENTO_DE_RESIDENCIAMETA            -0.041184   0.141525  -0.291   0.7712
## DEPARTAMENTO_DE_RESIDENCIANARIÑO          0.167258   0.176995   0.945   0.3453
## DEPARTAMENTO_DE_RESIDENCIANORTE SANTANDER  0.007917   0.106012   0.075   0.9405
## DEPARTAMENTO_DE_RESIDENCIAPUTUMAYO         0.392058   0.297024   1.320   0.1876
## DEPARTAMENTO_DE_RESIDENCIAQUINDIO          0.313253   0.212674   1.473   0.1416
## DEPARTAMENTO_DE_RESIDENCIARISARALDA        0.091268   0.155878   0.586   0.5585
## DEPARTAMENTO_DE_RESIDENCIASAN ANDRES      -0.380079   0.178441  -2.130   0.0338
## DEPARTAMENTO_DE_RESIDENCIASANTANDER       -0.075150   0.088387  -0.850   0.3957
## DEPARTAMENTO_DE_RESIDENCIASUCRE           -0.038064   0.091860  -0.414   0.6788
## DEPARTAMENTO_DE_RESIDENCIATOLIMA           0.136548   0.213117   0.641   0.5221
## DEPARTAMENTO_DE_RESIDENCIAVALLE            0.084582   0.095282   0.888   0.3753
## DEPARTAMENTO_DE_RESIDENCIAVAUPES          -0.373155   0.305234  -1.223   0.2223
## DEPARTAMENTO_DE_RESIDENCIAVICHADA          0.007455   0.217837   0.034   0.9727
## ESTRATO_DE_LA_VIVIENDAEstrato 1           -0.026459   0.054981  -0.481   0.6306
## ESTRATO_DE_LA_VIVIENDAEstrato 2           -0.010904   0.057835  -0.189   0.8506
## ESTRATO_DE_LA_VIVIENDAEstrato 3            0.037722   0.060762   0.621   0.5351
## ESTRATO_DE_LA_VIVIENDAEstrato 4           -0.010206   0.074021  -0.138   0.8904
## ESTRATO_DE_LA_VIVIENDAEstrato 5            0.114532   0.100717   1.137   0.2562
## ESTRATO_DE_LA_VIVIENDASin Estrato         -0.005961   0.082893  -0.072   0.9427
## TIPO_DE_COLEGIOOFICIAL                     0.021368   0.038106   0.561   0.5753
## CALENDARIOOTRO                             0.202782   0.296910   0.683   0.4950
## PUNT_LECTURA_CRITICA                       1.151797   0.002275 506.332   <2e-16
## PUNT_MATEMATICAS                           1.154963   0.002142 539.100   <2e-16
## PUNT_C_NATURALES                           1.152857   0.002512 458.927   <2e-16
## PUNT_SOCIALES_CIUDADANAS                   1.156430   0.001949 593.283   <2e-16
## PUNT_INGLES                                0.384005   0.001778 216.002   <2e-16
##                                              
## (Intercept)                                  
## GENEROM                                      
## SI_TIENE_ETNIASi                             
## DEPARTAMENTO_DE_RESIDENCIAARAUCA             
## DEPARTAMENTO_DE_RESIDENCIAATLANTICO          
## DEPARTAMENTO_DE_RESIDENCIABOGOTÃ\u0081       
## DEPARTAMENTO_DE_RESIDENCIABOLIVAR            
## DEPARTAMENTO_DE_RESIDENCIABOYACA             
## DEPARTAMENTO_DE_RESIDENCIACALDAS             
## DEPARTAMENTO_DE_RESIDENCIACAQUETA            
## DEPARTAMENTO_DE_RESIDENCIACASANARE           
## DEPARTAMENTO_DE_RESIDENCIACAUCA              
## DEPARTAMENTO_DE_RESIDENCIACESAR              
## DEPARTAMENTO_DE_RESIDENCIACHOCO              
## DEPARTAMENTO_DE_RESIDENCIACORDOBA            
## DEPARTAMENTO_DE_RESIDENCIACUNDINAMARCA       
## DEPARTAMENTO_DE_RESIDENCIAHUILA              
## DEPARTAMENTO_DE_RESIDENCIALA GUAJIRA         
## DEPARTAMENTO_DE_RESIDENCIAMAGDALENA          
## DEPARTAMENTO_DE_RESIDENCIAMETA               
## DEPARTAMENTO_DE_RESIDENCIANARIÑO            
## DEPARTAMENTO_DE_RESIDENCIANORTE SANTANDER    
## DEPARTAMENTO_DE_RESIDENCIAPUTUMAYO           
## DEPARTAMENTO_DE_RESIDENCIAQUINDIO            
## DEPARTAMENTO_DE_RESIDENCIARISARALDA          
## DEPARTAMENTO_DE_RESIDENCIASAN ANDRES      *  
## DEPARTAMENTO_DE_RESIDENCIASANTANDER          
## DEPARTAMENTO_DE_RESIDENCIASUCRE              
## DEPARTAMENTO_DE_RESIDENCIATOLIMA             
## DEPARTAMENTO_DE_RESIDENCIAVALLE              
## DEPARTAMENTO_DE_RESIDENCIAVAUPES             
## DEPARTAMENTO_DE_RESIDENCIAVICHADA            
## ESTRATO_DE_LA_VIVIENDAEstrato 1              
## ESTRATO_DE_LA_VIVIENDAEstrato 2              
## ESTRATO_DE_LA_VIVIENDAEstrato 3              
## ESTRATO_DE_LA_VIVIENDAEstrato 4              
## ESTRATO_DE_LA_VIVIENDAEstrato 5              
## ESTRATO_DE_LA_VIVIENDASin Estrato            
## TIPO_DE_COLEGIOOFICIAL                       
## CALENDARIOOTRO                               
## PUNT_LECTURA_CRITICA                      ***
## PUNT_MATEMATICAS                          ***
## PUNT_C_NATURALES                          ***
## PUNT_SOCIALES_CIUDADANAS                  ***
## PUNT_INGLES                               ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.2853 on 384 degrees of freedom
## Multiple R-squared:      1,  Adjusted R-squared:      1 
## F-statistic: 3.566e+05 on 44 and 384 DF,  p-value: < 2.2e-16
coef(modelo)
##                               (Intercept) 
##                               0.021885079 
##                                   GENEROM 
##                              -0.044944158 
##                          SI_TIENE_ETNIASi 
##                               0.045414948 
##          DEPARTAMENTO_DE_RESIDENCIAARAUCA 
##                               0.227954523 
##       DEPARTAMENTO_DE_RESIDENCIAATLANTICO 
##                               0.008486486 
##    DEPARTAMENTO_DE_RESIDENCIABOGOTÃ\u0081 
##                               0.015534404 
##         DEPARTAMENTO_DE_RESIDENCIABOLIVAR 
##                               0.042515321 
##          DEPARTAMENTO_DE_RESIDENCIABOYACA 
##                               0.272889973 
##          DEPARTAMENTO_DE_RESIDENCIACALDAS 
##                              -0.441704115 
##         DEPARTAMENTO_DE_RESIDENCIACAQUETA 
##                              -0.027181377 
##        DEPARTAMENTO_DE_RESIDENCIACASANARE 
##                               0.010481037 
##           DEPARTAMENTO_DE_RESIDENCIACAUCA 
##                              -0.111431968 
##           DEPARTAMENTO_DE_RESIDENCIACESAR 
##                               0.012018748 
##           DEPARTAMENTO_DE_RESIDENCIACHOCO 
##                              -0.094032077 
##         DEPARTAMENTO_DE_RESIDENCIACORDOBA 
##                               0.014388922 
##    DEPARTAMENTO_DE_RESIDENCIACUNDINAMARCA 
##                              -0.094221026 
##           DEPARTAMENTO_DE_RESIDENCIAHUILA 
##                              -0.081137465 
##      DEPARTAMENTO_DE_RESIDENCIALA GUAJIRA 
##                              -0.084937979 
##       DEPARTAMENTO_DE_RESIDENCIAMAGDALENA 
##                              -0.057815807 
##            DEPARTAMENTO_DE_RESIDENCIAMETA 
##                              -0.041183588 
##         DEPARTAMENTO_DE_RESIDENCIANARIÑO 
##                               0.167258202 
## DEPARTAMENTO_DE_RESIDENCIANORTE SANTANDER 
##                               0.007916685 
##        DEPARTAMENTO_DE_RESIDENCIAPUTUMAYO 
##                               0.392058406 
##         DEPARTAMENTO_DE_RESIDENCIAQUINDIO 
##                               0.313252736 
##       DEPARTAMENTO_DE_RESIDENCIARISARALDA 
##                               0.091268359 
##      DEPARTAMENTO_DE_RESIDENCIASAN ANDRES 
##                              -0.380078763 
##       DEPARTAMENTO_DE_RESIDENCIASANTANDER 
##                              -0.075150005 
##           DEPARTAMENTO_DE_RESIDENCIASUCRE 
##                              -0.038064020 
##          DEPARTAMENTO_DE_RESIDENCIATOLIMA 
##                               0.136548365 
##           DEPARTAMENTO_DE_RESIDENCIAVALLE 
##                               0.084582108 
##          DEPARTAMENTO_DE_RESIDENCIAVAUPES 
##                              -0.373154898 
##         DEPARTAMENTO_DE_RESIDENCIAVICHADA 
##                               0.007454583 
##           ESTRATO_DE_LA_VIVIENDAEstrato 1 
##                              -0.026458930 
##           ESTRATO_DE_LA_VIVIENDAEstrato 2 
##                              -0.010904246 
##           ESTRATO_DE_LA_VIVIENDAEstrato 3 
##                               0.037722346 
##           ESTRATO_DE_LA_VIVIENDAEstrato 4 
##                              -0.010205917 
##           ESTRATO_DE_LA_VIVIENDAEstrato 5 
##                               0.114531983 
##         ESTRATO_DE_LA_VIVIENDASin Estrato 
##                              -0.005961418 
##                    TIPO_DE_COLEGIOOFICIAL 
##                               0.021367816 
##                            CALENDARIOOTRO 
##                               0.202782224 
##                      PUNT_LECTURA_CRITICA 
##                               1.151796808 
##                          PUNT_MATEMATICAS 
##                               1.154963096 
##                          PUNT_C_NATURALES 
##                               1.152856876 
##                  PUNT_SOCIALES_CIUDADANAS 
##                               1.156430208 
##                               PUNT_INGLES 
##                               0.384005169
par(mfrow = c(2, 2))
plot(modelo)
## Warning: not plotting observations with leverage one:
##   149, 219, 241, 279, 297