class: center, middle, inverse, title-slide .title[ # Análisis de la anemia en niños de 6 a 35 meses de vida ] .author[ ### Rodrigo Carrasco Pastor ] .institute[ ### Universidad Nacional Mayor de San Marcos ] --- class: middle,center ## La anemia ### Definición La anemia es una enfermedad en la cual los glóbulos rojos del cuerpo disminuyen por debajo del nivel normal para la edad del niño. Puede causar que el niño luzca pálido y se sienta irritable, cansado o débil. --- class: middle,center ## Evolución de la anemia en niños  .footnote[Fuente: ENDES-INEI] --- class: middle,center ### Tasa de la anemia por regiones para el año 2021 <!-- --> .footnote[ La region con mayor incidencia de anemia es Puno (70%) y la region con menos incidencia es Tacna (25,3%)] --- class: middle,center ### Mapa de la anemia por regiones para el año 2021 Los shapefile se obtuvieron de la pagina [**GEO GPS PERU**](https://www.geogpsperu.com/2014/03/base-de-datos-peru-shapefile-shp-minam.html) <!-- --> --- class: middle,center ### Anemia por quintiles de ingreso para el año 2021 <!-- --> .footnote[ La mayor incidencia de anemia se encuentra en el quintil inferior, es decir, en el quintil más bajo aproximadamente el 50% de los niños de 6 a 35 meses tienen anemia] --- class: middle,center ### Análisis econométrico Para analizar la incidencia de la anemia se va a realizar un modelo de respuesta discreta **logit**. #### Modelo logit Es un modelo de respuesta discreta que se usa para modelizar el comportamiento de una variable dependiente (respuesta) de naturaleza discreta (categórica) a partir de un conjunto de variables independientes (predictoras) que pueden ser continuas o discretas --- class: middle ### Formulación del modelo Un modelo logit viene expresado a través de: `$$\ln \left[\frac{p\left(x_{i}\right)}{1-p\left(x_{i}\right)}\right]=a+\sum_{j} b_{j} x_{i j}, \quad(i=1, \ldots, N),(j=1, \ldots, R)$$` Donde J es el número de variables independientes, n el número de observaciones en cada una, `\(x_{i}\)` = `\(x_{i 1},\ldots,x_{i R}\)` el vector que contiene las observaciones de cada variable para el individuo i-ésimo y el cociente `\(\left[\frac{p\left(x_{i}\right)}{1-p\left(x_{i}\right)}\right]\)` representa la ventaja de respuesta Y=1 para los valores observados de las variables independientes .footnote[ Para la formulación del modelo se tomó como referencia el libro [**Modelo de respuesta discreta en R y aplicación con datos reales**](https://docplayer.es/17444675-Alejandro-salazar-couso-master-oficial-en-estadistica-aplicada-universidad-de-granada.html)] --- class: middle ### Aplicación Para la aplicacion del modelo se utilizó los módulos 1638 (peso y talla-anemia) y 1630 (caracteristicas de la vivienda) del [**INEI-ENDES**](http://iinei.inei.gob.pe/microdatos/) ```r colnames(anemia2) ``` ``` ## [1] "ID1" "HHID" "HC0" "HV101" ## [5] "años_estudio_JH" "edad_JH" "indice_riqueza" "dpto" ## [9] "longitudx" "latitudy" "peso" "aguaVivienda" ## [13] "generoJH" "servHig" "anemia" "_merge" ``` --- ### Resumen de las variables ```r summary(anemia2) ``` ``` ## ID1 HHID HC0 HV101 años_estudio_JH ## Min. :2021 Length:11112 Min. :1 Min. :1 Min. : 0.00 ## 1st Qu.:2021 Class :character 1st Qu.:1 1st Qu.:1 1st Qu.: 8.00 ## Median :2021 Mode :character Median :1 Median :1 Median :11.00 ## Mean :2021 Mean :1 Mean :1 Mean :10.71 ## 3rd Qu.:2021 3rd Qu.:1 3rd Qu.:1 3rd Qu.:13.00 ## Max. :2021 Max. :1 Max. :1 Max. :98.00 ## edad_JH indice_riqueza dpto longitudx ## Min. :15.00 Min. :1.000 Min. : 1.0 Min. :-81.31 ## 1st Qu.:30.00 1st Qu.:1.000 1st Qu.: 7.0 1st Qu.:-77.68 ## Median :36.00 Median :2.000 Median :14.0 Median :-76.19 ## Mean :38.78 Mean :2.444 Mean :13.3 Mean :-75.58 ## 3rd Qu.:45.00 3rd Qu.:3.000 3rd Qu.:19.0 3rd Qu.:-73.43 ## Max. :97.00 Max. :5.000 Max. :25.0 Max. :-43.09 ## latitudy peso aguaVivienda generoJH ## Min. :-18.328 Min. :0.01525 Min. :0.0000 Min. :0.0000 ## 1st Qu.:-13.293 1st Qu.:0.13874 1st Qu.:0.0000 1st Qu.:0.0000 ## Median :-11.798 Median :0.31090 Median :1.0000 Median :0.0000 ## Mean :-10.614 Mean :0.35978 Mean :0.6614 Mean :0.2694 ## 3rd Qu.: -7.068 3rd Qu.:0.47680 3rd Qu.:1.0000 3rd Qu.:1.0000 ## Max. : -2.134 Max. :1.45293 Max. :1.0000 Max. :1.0000 ## servHig anemia _merge ## Min. :0.0000 Min. :0.0000 Min. :3 ## 1st Qu.:0.0000 1st Qu.:0.0000 1st Qu.:3 ## Median :1.0000 Median :0.0000 Median :3 ## Mean :0.5427 Mean :0.4198 Mean :3 ## 3rd Qu.:1.0000 3rd Qu.:1.0000 3rd Qu.:3 ## Max. :1.0000 Max. :1.0000 Max. :3 ``` --- ### Modelo logit de anemia #### Parámatros del modelo ```r mod_log <- glm(anemia~años_estudio_JH+edad_JH+generoJH+indice_riqueza+servHig+aguaVivienda, family = binomial(link = "logit"), data = anemia2) ```  Vemos que las variables género del jefe del hogar (generoJH1) y si la vivienda tiene agua (aguaVivienda1) son no significativas --- ### Modelo logit de anemia #### Modelo final ```r mod_log2 <- glm(anemia~años_estudio_JH+edad_JH+indice_riqueza+servHig, family = binomial(link = "logit"), data = anemia2) ```  En el modelo final todas las variables son estadisticamente significativas --- ### Interpretación del modelo ```r exp(coef(summary(mod_log2))) ``` ``` ## Estimate Std. Error z value Pr(>|z|) ## (Intercept) 1.7143026 1.080199 1.082322e+03 1.000000 ## años_estudio_JH 0.9929431 1.002827 8.138104e-02 1.012194 ## edad_JH 0.9902492 1.001660 2.722294e-03 1.000000 ## indice_riqueza2 0.8088646 1.054993 1.901797e-02 1.000074 ## indice_riqueza3 0.6287989 1.069288 9.826559e-04 1.000000 ## indice_riqueza4 0.4711116 1.082112 7.208868e-05 1.000000 ## indice_riqueza5 0.3811097 1.100055 4.042577e-05 1.000000 ## servHig1 0.8580459 1.053659 5.344938e-02 1.003406 ``` -- a) Para el caso de la variable servHig1 se considera la inversa porque el odd ratio es menor a 1 (1/0.858). Entonces indica que la incidencia de la anemia aumenta 1.16 veces mas para los hogares que no tienen servicios higiénicos con respecto a los hogares que si tienen servicios higienicos. -- b)Para la variable indice_riqueza5 que indica el quintil más alto de la distribucion del ingreso se interpreta que la incidencia de la anemia es 2.65 veces superior en los hogares del primer quintil(quintil base) en comparación con el quintil más alto. --- class: middle,center ### GRACIAS !!!