Estimación del ingreso por trabajo en los municipios y las delegaciones de México utilizando técnicas de estimación para áreas pequeñas

 

Vol.6 Núm.3                                                                                 Estimación del ingreso por trabajo…

En el muestreo de poblaciones finitas —y en particular en el sistema estadístico nacional— existe una demanda creciente de estimaciones precisas sobre promedios o indicadores de interés en áreas cada vez más pequeñas (entidades federativas, municipios o localidades, o bien, en subgrupos o pequeños dominios de la población total, como las subclases de alguna actividad económica). Estos cálculos se consideran un subproducto de los trabajos de muestreo en áreas grandes, donde se diseña una muestra para obtener valores del total o de la media de una característica de interés con una precisión prefijada. Esto implica que el número de observaciones muestrales en un área pequeña es reducido o, incluso, nulo, por lo que utilizar estimadores basados en un diseño muestral conduce a grandes errores en los resultados y, de hecho, es imposible obtenerlos en áreas no muestreadas. Para resolver esta problemática, es necesario, por un lado, aumentar el tamaño de muestra (lo cual implica elevar los costos) y, por el otro, aplicar la técnica estadística estimación para áreas pequeñas (EAP).

Por tradición, los procedimientos de inferencia sobre características de poblaciones finitas se han basado en el diseño de la muestra; sin embargo, en la actualidad son cada vez más numerosos los de aproximación basada en modelos (Cassel et al., 1977), en la cual se le da más peso al modelo de la distribución de la característica en estudio en la población que al proceso de diseño de la muestra y se supone la aceptación de algunos riesgos (Hansen y Madow, 1983), como el no estar seguros de que el modelo que se adopta es el correcto; aun así, deben ser aceptados, ya que sólo es posible considerar la inferencia basada en el diseño libre de supuestos cuando la muestra es grande, lo cual no ocurre con la estimación para áreas pequeñas. Al recurrir a este tipo de inferencia se pueden solventar problemas no abordables por la vía del diseño, como la no respuesta o la estimación en áreas pequeñas (Särndal, 1984).

Harter (1993) recopiló diferentes métodos para la EAP para obtener valoraciones de la población en ciudades y municipios en periodos intercensales; en fecha más reciente, Rao (2003) hizo una compilación de aquéllos y otros métodos en su libro Small Area Estimation.

La EAP relaciona mediante un modelo a la variable de interés obtenida de una encuesta con la información emanada tanto de eventos censales como de registros administrativos y geográficos, contenida en las denominadas variables auxiliares, de las cuales se conoce, en unos casos, el valor para cada elemento de la población y en otros, sólo la información agregada (promedios, totales o proporciones) de cada dominio o área pequeña, de donde se tienen dos grandes tipos de modelos de EAP: los de nivel de unidad y los de nivel de área. Para los primeros, es necesario asociar la información auxiliar contenida en los elementos de la población (censos o registros administrativos) y su unidad muestral correspondiente, lo cual no siempre es posible, mientras que los segundos nada más necesitan información agregada de cada dominio o área pequeña como regresores del modelo y los estimadores directos obtenidos del muestreo como variables de respuesta.

En este trabajo se optó por utilizar la EAP basada en modelos de área con el propósito de obtener una estimación del promedio por municipio¹ de la variable ingreso por trabajo en la vivienda, la cual se eligió por dos razones:

Para conseguir ese propósito, este documento se organiza de la siguiente manera: primero se muestra brevemente la ENIGH 2010, centrándose en su diseño de muestreo, los estimadores que utiliza y la variable ingreso por trabajo; se continúa con la presentación de la muestra del CPV 2010, haciendo hincapié en su diseño de muestreo y la misma variable; después, se comenta de forma concisa la estimación basada en el modelo y sus ventajas cuando hay escasez de muestra; se expone, además, un resumen de la base teórica del modelo a nivel de área, que es el que aquí se utiliza para realizar las estimaciones en las áreas pequeñas; se continúa con una explicación sobre la construcción del modelo, la selección de las variables que mejor se ajustan al modelo planteado (el ajuste se basa en el cumplimiento de pruebas estadísticas al modelo, mismas que se muestran en su formulación matemática); se prosigue con la obtención del promedio municipal del ingreso mensual por trabajo en la vivienda, obteniéndolo primero con la formulación de la estimación directa y, después, aplicando el modelo a nivel área, para hacer una comparación gráfica de los resultados; a manera de validación, se cotejan también gráficamente las EAP con las de la muestra del CPV 2010; por último, se exponen las conclusiones logradas con los resultados y el método de estimación aplicado.

ENIGH 2010

Tiene sus antecedentes en varios sondeos realizados por diferentes dependencias públicas, como la Secretaría de Industria y Comercio (SIC), el Banco de México, la Secretaría del Trabajo y Previsión Social (STPS) o la Secretaría de Programación y Presupuesto (SPP), pero fue a partir de 1984 que se integró como tal y ha sido levantada de manera formal por el Instituto Nacional de Estadística y Geografía (INEGI). Desde 1992 se ha realizado con una periodicidad bienal, con excepción del 2005, ya que fue un levantamiento ex profeso.

La ENIGH 2010 se levantó del 21 de agosto al 28 de noviembre con el objetivo de proporcionar información sobre la distribución, monto y estructura del ingreso y gasto de los hogares, así como de ofrecer datos tanto de las características sociodemográficas y ocupacionales de los integrantes del hogar como de la infraestructura de la vivienda y el equipamiento del hogar.
Los resultados de este levantamiento se encuentran a nivel nacional y para los ámbitos rural y urbano. Además, se tiene información para las entidades que, en su momento, convinieron con el INEGI una ampliación de la muestra para este operativo (Chiapas, Guanajuato, Distrito Federal, estado de México y Yucatán).

El marco de muestreo utilizado fue el de propósitos múltiples del INEGI, constituido con la información demográfica y cartográfica obtenida a partir del levantamiento del XII Censo General de Población y Vivienda 2000. El marco es probabilístico, estratificado,² unietápico y por conglomerados; estos últimos se consideran unidades primarias de muestreo (UPM),³ pues es en ellos donde, en una segunda etapa, se seleccionan las viviendas que forman parte de la muestra (INEGI, c2011), por lo que la elección de la muestra para la ENIGH 2010 se realizó en dos etapas: en la primera se escogieron las UPM y en la segunda, las viviendas objeto de entrevista de la Encuesta.
Por lo tanto, el esquema de muestreo de la ENIGH 2010 se considera como complejo, siendo éste probabilístico, estratificado, bietápico y por conglomerados, donde la unidad es la vivienda.
El estimador directo que se utiliza para la media en un dominio S (entidad, ciudad, municipio, etc.) es un derivado del estimador de Horvitz- Thompson, que tiene la forma:

donde y_i es el valor i-ésimo de la variable de interés; n_s, el número de elementos en el dominio S (por simplicidad, las ecuaciones de esta sección se refieren al dominio S) y π_i;la probabilidad de selección de y_i, dado por:

donde:

Al inverso de la probabilidad de selección w_i=1/π_i se le conoce como factor de expansión; para la ENIGH 2010, este factor tiene ajustes por no respuesta y por encuadre con proyecciones de población, siendo ésta la razón por la que no es un estimador de Horvitz-Thompson puro.

La varianza estimada directa del estimador de la media se obtiene con la siguiente expresión:

donde:

La variable a estimar es la de ingreso por trabajo;⁴ en la ENIGH 2010 se consideró que un integrante del hogar percibe ingreso por trabajo sólo si tiene o ha tenido participación directa en actividades reconocidas como económicas. Por sus fuentes,⁵ los ingresos del trabajo pueden provenir de: las remuneraciones por el trabajo subordinado, los ingresos por el trabajo independiente y otros ingresos provenientes del trabajo. La referencia temporal de todos ellos es del mes inmediato anterior al levantamiento y de los cinco meses anteriores (INEGI, 2011c).

De los 2 456 municipios existentes al 2010, únicamente en 600 hubo muestra de la ENIGH 2010 con viviendas que reportaron ingresos por trabajo, con una o más UPM por municipio.

Muestra del CPV 2010

Se elaboró un cuestionario ampliado para el levantamiento de la muestra censal⁶ 2010, que fue aplicado del 31 de mayo al 25 de junio de ese año en alrededor de 2.9 millones de viviendas en el país. Dados los requerimientos de información, se decidió incluir con certeza en la muestra las viviendas habitadas de los municipios con menos de 1 100 de éstas, así como las de los 125 con menor índice de desarrollo humano (IDH), sin importar su tamaño. Bajo este criterio, en 766 municipios se censaron la totalidad de viviendas con el cuestionario ampliado; para los 794 que tenían entre 1 101 y 4 mil viviendas habitadas el tamaño de la muestra fue de 800; en los municipios sin localidades de 50 mil y más habitantes, la muestra fue de 1 100 viviendas y para los 193 restantes el tamaño fue variable, pero mayor a 2 mil. Los tamaños de muestra fijados garantizan estimaciones municipales aceptables para proporciones cercanas a 0.01 o mayores.

La variable estimada a nivel municipal es la de ingresos por trabajo, definida para la muestra del CPV 2010 como la percepción monetaria que la población ocupada obtiene o recibe del(los) trabajo(s) que desempeñó en la semana de referencia. Se consideran los ingresos por concepto de ganancia, comisión, sueldo, salario, jornal, propina o cualquier otro devengado de su participación en alguna actividad económica. Los ingresos están calculados de forma mensual (INEGI, 2011c).

Estimación basada en el modelo

Permite relacionar, mediante el uso de información auxiliar, a las áreas con escasez o inexistencia de muestra con aquéllas próximas de mayor información muestral para, de esta forma, incrementar la precisión del estimador; además, el uso de la estimación basada en modelos tiene las siguientes ventajas:

Para la estimación basada en modelos en el caso de las áreas pequeñas se tienen, como ya se dijo antes, dos grandes tipos de modelos según la disponibilidad de la información auxiliar: a nivel de área y de unidad. En este trabajo se utilizó sólo el primero para una variable cuantitativa, razón por la que este modelo se presenta más adelante con cierto nivel de detalle.

Una premisa importante en estos modelos es que el comportamiento de la variable bajo estudio en las áreas pequeñas presenta ligeras diferencias entre sí y respecto al comportamiento de la misma en el área mayor que las contiene; esto se refleja considerando efectos aleatorios de área en el modelo lineal que se plantea, con lo que se introduce a los modelos lineales mixtos, llamados así porque abarcan en su expresión efectos fijos y aleatorios.

Resumen teórico del modelo a nivel de área

Se basa en el modelo de Fay-Herriot (1979), donde se tienen p variables como información auxiliar (promedios poblacionales por área) en el vector y se supone que están relacionadas con la media en la subpoblación de la variable de interés , o una cierta función de ésta, a través del modelo lineal con efectos aleatorios:

donde β es el vector de parámetros de regresión y ʋ_α , el efecto aleatorio que se supone independiente e idénticamente distribuido (iid) con media 0 y varianza σ²_ʋ por lo general, se supone la normalidad de ʋ_α . En caso de que no todas las áreas sean seleccionadas en la muestra, se continúa bajo el supuesto de que las muestreadas (m) obedecen al modelo de población.

Por otro lado, sea el estimador directo de la media de la variable Y en el área pequeña α-ésima, esto supone que el tamaño de la muestra en el área (n_α) es mayor o igual a 1.
Ahora, se tiene que:

donde e_α son los errores muestrales que son independientes con media 0 y la varianza Ψ_α se supone conocida; se puede relajar este supuesto reemplazando Ψ_α por estimadores suavizados basados en las varianzas calculadas de los datos a nivel unidad (Rao, 2003).
Sustituyendo (3) en (4) se obtiene el modelo lineal mixto:

Se observa que este modelo involucra tanto errores del diseño muestral e_α como los del modelo ʋ_α y se supone que e_α y ʋ_α son independientes.

Bajo el modelo (5), el mejor estimador lineal insesgado (BLUP, por sus siglas en inglés) de θ_α, el cual minimiza el error cuadrático medio (MSE, por sus siglas en inglés) es (Rao & Molina, 2012):

donde:

El estimador puede expresarse como:

Se observa que el estimador BLUP de θ_α se expresa como un promedio ponderado del estimador directo y el estimador de regresión sintética donde el ponderador mide la incertidumbre en el modelizado del predictor para cada área pequeña α (Azula et al., 2004).
Como depende de σ²_ʋ a través de y se puede utilizar el BLUP empírico EBLUP, reemplazando σ²_ʋ por su estimador :

En lo referente a la estimación de θ para las k áreas no muestreadas sólo se aplica el estimador de regresión sintética:

Una versión más general de (6) que no se limita a la linealidad de θ_α es el mejor estimador (B) o también llamado estimador de Bayes, el que bajo normalidad se expresa como:

es insesgado ya que:

También, el mejor estimador empírico bajo normalidad de θ_α coincide con el EBLUP, esto es:

Bajo el supuesto de normalidad de los efectos aleatorios se pueden estimar los componentes de la varianza por máxima verosimilitud —Maximum Likelihood (ML)— o máxima verosimilitud restringida —Restricted Maximum Likelihood (REML)—,⁷ esta última reduce el sesgo de la estimación ML, ya que no depende de β . La función log-verosímil del modelo mixto en cuestión es:

donde . Derivando parcialmente esta función con respecto a β y σ²_ʋ tomando a Ψ como constante conocida e igualando a cero se tienen las ecuaciones para la estimación de β y σ²_ʋ y (Rao, 2003):

Es claro que la estimación ML de β es justamente su estimación por mínimos cuadrados generalizados con el parámetro , la ecuación de la varianza σ²_ʋno tiene una solución analítica, por lo que se debe resolver de forma numérica.

La ecuación de estimación de σ²_ʋ por REML es más simple, ya que no depende de β (Rao, 2003):

con P = V^-1 – V^-1X(X^TV^-1X)^-1X^TV^-1, sin embargo, tampoco tiene una solución analítica, sólo por métodos numéricos. Para esto, se utiliza el algoritmo de Fisher scoring (Rao, 2003), mismo que se describe con los siguientes pasos:

Las funciones del paso 5 son la primera y segunda derivadas de la función de ML o de REML, según sea el caso. Para ML se tienen las siguientes expresiones:

para REML se tiene:

donde P = V^-1 – V^-1x̄(x̄^TV^-1x̄)^-1 x̄^TV^-1 . Ahora, bajo supuestos de normalidad de e_α y, ʋ_α, el error cuadrático medio de con respecto al modelo (Rao, 2003) es:

donde:

es la varianza asintótica de σ²_ʋ, su forma depende del método de estimación usado para σ²_ʋ .
En los sumandos de (8), g1α representa el error debido al efecto aleatorio, g2α representa el error por la estimación de y β y g3α, el que ocurre debido a la estimación de la varianza σ²_ʋ .

Cuando se obtiene por REML, el estimador insesgado del MSE es:

Si se obtiene por ML se debe añadir un término extra de sesgo (Rao, 2003).

El error cuadrático medio para las áreas no muestreadas es:

donde es la matriz de varianzas y covarianzas de los coeficientes β estimados.

Construcción del modelo

Si se tienen p variables explicativas, se puede conformar un total de 2^p-1 modelos a elegir. A medida que el número de variables aumenta, la cantidad de cálculos necesarios se incrementa muy rápido. Cuando el número de variables es grande (p>40), el trabajo de cómputo es enorme; una forma de aligerar la carga se obtiene al aplicar el método de selección hacia adelante (forward). El procedimiento inicia eligiendo de las p variables aquella que tiene la mayor correlación con la que se va a explicar para después ir añadiendo una de las restantes —la que mejor contribuya al modelo— y, así, hasta llegar a un máximo determinado o una condición a cumplir o bien, hasta que nuevas variables no aporten estadísticamente más al modelo. Es pertinente repetirlo iniciado con otra variable también representativa y comparar los modelos resultantes.

Las condiciones a cumplir en la elección del modelo se dividen en tres apartados: los criterios estadísticos utilizados primordialmente para comparar dos modelos y seleccionar el mejor, la significancia estadística de los efectos fijos y el cumplimiento de los supuestos estadísticos del modelo.

Primer apartado

Segundo apartado

Tercer apartado

Aplicación de la EAP con modelo a nivel de área para estimar el promedio municipal del ingreso por trabajo

El objetivo de la aplicación de la EAP es obtener las estimaciones del promedio del ingreso por trabajo de todos los municipios del país utilizando como insumos:

Una vez obtenidas las estimaciones del modelo a nivel de área y las directas, se realizará un comparativo de ambas con las cifras obtenidas de la muestra del CPV 2010.

Previo al cálculo de las estimaciones se realizaron dos importantes ajustes: el primero es para hacer comparable la información de los ingresos respecto a la diferencia temporal entre la ENIGH 2010 y la muestra del CPV 2010, el método utilizado es el que sugieren Cortés y Rubalcava (1994), donde se toma cada rubro en la ENIGH 2010 que forma parte del ingreso por trabajo para cada uno de los seis meses captados y se deflacta a la mitad del periodo de referencia del CPV 2010, es decir, a junio del 2010; el segundo es respecto a la muestra a nivel municipal para asegurar que la expansión de la muestra del número de viviendas sea el contabilizado por el CPV 2010 para cada municipio; se hace utilizando la siguiente expresión:

donde:

La estimación directa del promedio municipal del ingreso por trabajo por vivienda de los municipios en los que hay muestra en la ENIGH 2010 se calcula con la ecuación (1) tomando a w’_αi=1/π_i y al dominio S como el municipio (área pequeña) α, la varianza respectiva se calcula con la ecuación (2). Estos cálculos se realizan utilizando el paquete estadístico R con la librería survey; en los resultados se aprecia que los 600 municipios con muestra están ordenados de menor a mayor según su número de viviendas censales (ver gráfica 1); su correspondiente intervalo de confianza se obtiene con la siguiente expresión:

Se debe aclarar que para varios municipios no fue posible calcular la varianza debido a que sólo una UPM fue muestreada en ellos; por esta razón, en la gráfica 1 varios promedios estimados no incluyen su intervalo de confianza.

En el análisis exploratorio de la variable a estimar y de su varianza se determinó que existe una relación entre éstas debido a que su correlación es cercana a 0.5 y a que, para cumplir con el supuesto de homogeneidad de la varianza y la normalidad de la variable dependiente, era necesario realizar una transformación por lo que, partiendo de la expresión (3), el modelo propuesto ahora es:

Uno de los supuestos fundamentales del modelo a nivel de área es que se conoce la varianza Ψ_α para cada área pequeña, pero por lo general no es así, entonces se sustituye por su estimador suavizado ; éste se obtiene partiendo de la estimación directa de la varianza de los datos a nivel unidad , pero como premisa es sabido que los tamaños de la muestra en los municipios son pequeños, por lo que no es un buen estimador, ya que puede tomar valores en extremo grandes o increíblemente pequeños. Para solucionar esta situación, las estimaciones de se suavizan sustituyendo los valores extremos por los obtenidos del siguiente modelo lineal para la varianza:

Para el ajuste de los parámetros del modelo (12) se toma la varianza de los municipios en muestra, excluyendo los 333 que sólo cuentan con una UPM y cuya varianza no es posible calcular, además de excluir también 13 municipios cuyo coeficiente de variación es mayor a 0.35.⁸
Así, utilizando las predicciones de este modelo lineal se pretende obtener valores más realistas de para estos 13 municipios y para los 333 con una sola UPM.
Es importante comentar que los municipios fueron excluidos sólo para modelar la varianza; no lo son para las demás actividades del ajuste del modelo a nivel área.
Ahora, para la estimación en áreas pequeñas, es necesario establecer el modelo a nivel de área, que obedece a la ecuación (6) con la transformación del modelo (11) y que cumpla las condiciones señaladas en los tres apartados descritos en la sección Construcción del modelo; estas condiciones se verifican utilizando funciones de las librerías stats, nnet, MASS, survival, car y nlme del paquete estadístico R, así como a través de la elaboración de programas de cómputo propios.
Las variables regresoras contenidas en X_α se seleccionan bajo los criterios del primer apartado, resultando que las que mejor explican el promedio municipal del ingreso por trabajo en la vivienda para la ENIGH 2010 son cinco; sus mnemónicos y descripciones son los siguientes:

Los resultados de la validación de las condiciones del segundo apartado se presentan en el cuadro 1.

El cumplimiento de las condiciones del tercer apartado se aplica para verificar la normalidad de los residuales del modelo, la ausencia de multicolinealidad en los valores de X, la homocedasticidad de los residuales y la normalidad de los efectos aleatorios, siendo este último el supuesto subyacente para la estimación de los componentes de la varianza.

Para verificar la normalidad de los residuales del modelo, se utiliza la prueba de normalidad de Shapiro-Wilks, aplicando la función resid (modelo, type = normalized) de R, la cual usa la descomposición de Cholesky, que elimina la dependencia entre los residuales dada la estructura de correlación (Valencia, M., 2010); el resultado es:

Shapiro-Wilk normality test
data: resnorm
W = 0.9948, p-value = 0.04058.

Por ello, no se rechaza la hipótesis de normalidad a un nivel de significancia de 4 por ciento.
Los resultados de la aplicación de las pruebas de multicolinealidad VIF y número de condición se observan en el cuadro 2, ambas indican la ausencia de problemas severos de multicolinealidad, ya que los VIF obtenidos son siempre menores a 10 y el número de condición es menor que 30.

Para verificar la homocedasticidad en el modelo se aplica la prueba de Breush-Pagan, utilizando la función ncvtest de R, con el resultado siguiente:

Non-constant Variance Score Test
Variance formula: ~ fitted.values
Chisquare = 0.7538, Df = 1, p = 0.3852.

Esto indica que la varianza del error puede considerarse como constante.
La prueba de Shapiro-Wilks es la empleada para comprobar la normalidad de los efectos aleatorios, dando el resultado siguiente:

Shapiro-Wilk normality test
data: eblup1a$randeff
W = 0.998, p-value = 0.7216.

El valor p observado indica que hay evidencia de normalidad en los efectos aleatorios ʋ.
Una vez que el modelo (6) con las variables seleccionadas cumple con las condiciones establecidas, se obtienen las estimaciones del promedio del ingreso por trabajo en las viviendas para los municipios con muestra en la ENIGH 2010, mediante la aplicación de la librería sae del paquete estadístico R. Como resultado se obtienen las estimaciones de los promedios que se presentan en la gráfica 2, así como la

varianza de los efectos aleatorios, con un valor estimado de = 112.467.

La escala de la gráfica 2 es igual a la de la 1 con el propósito de facilitar la comparación visual entre ellas; al hacerlo se aprecia, a simple vista, que las estimaciones se condensan siguiendo el modelo planteado, además de que se homogeneizan los intervalos de confianza a una longitud razonable, mismos que se obtienen aplicando la ecuación:

donde ya está considerada la transformación a las unidades originales (pesos mexicanos). Para transformar el MSE se aplicó el método Delta que relaciona la varianza de Y con la de θ(Y) = √Y (para más detalles, ver Velasco, C., 2007).

La gráfica 3 se construye tomando como base la 2, a la cual se añaden los puntos en azul que corresponden a la estimación del promedio del ingreso por trabajo en la vivienda de la muestra del CPV 2010 en los municipios con muestra de la ENIGH 2010; se observa cómo la condensación citada en la gráfica 2 corresponde con las estimaciones de la muestra censal, aun cuando la ENIGH 2010 y la muestra del CPV 2010 fueron levantamientos y diseños muestrales por completo independientes.

El coeficiente de variación (CV) es una medida que puede proporcionar, en casos específicos,⁹ una buena idea de la precisión de las estimaciones. En la gráfica 4 se observan los CV de las estimaciones directas y del modelo de área de los municipios con muestra en la ENIGH 2010, mismos que están ordenados de menor a mayor tamaño de muestra. En la línea punteada se establece el valor del CV = 0.25 como umbral para indicar que los CV menores representan estimaciones razonablemente confiables; la gran mayoría de las estimaciones con el modelo de área (puntos rojos) están por debajo de este umbral. Los puntos azules con CV igual a cero representan a los municipios con sólo una UPM en muestra, para los que en realidad su CV está indeterminado.

Ahora, para la estimación del promedio municipal del ingreso por trabajo en la vivienda para los municipios que no tienen muestra en la ENIGH 2010, se utilizó la ecuación (7); la varianza respectiva se calcula mediante la ecuación (10) y los intervalos de confianza con la siguiente expresión:

donde está considerada la transformación a las unidades originales.
Los resultados se muestran en la gráfica 5, la cual presenta, también, la estimación de la misma variable, pero de la muestra del CPV 2010. Se aprecia que, en general, las estimaciones con el modelo son menores a las del Censo. Se conserva la misma escala de la gráfica 1 para facilitar su comparación visual.

Enseguida, se realiza el cálculo del MSE empírico de las estimaciones con el método directo, y de las del modelo a nivel de área, ambas respecto a la estimación de la muestra del CPV 2010; para el segundo caso se diferencia a los municipios con muestra de los que no la tuvieron. Dichos cálculos se confrontan en la gráfica 6, en la cual se puede apreciar que el MSE empírico entre el método directo y el del modelo a nivel de área se reduce aproximadamente 3.9 veces, lo cual es razonable en virtud de que el dominio de la muestra de la ENIGH 2010 no abarca un nivel municipal. El MSE empírico en las estimaciones de los municipios no muestreados es muy favorable, debido a que la mayoría de ellos son de ingresos bajos, por lo tanto, sus diferencias cuadráticas con los ingresos de la muestra censal son menores que las de los municipios con muestra; sin embargo, el modelo para estos municipios es conservador, ya que proporciona intervalos de confianza amplios (ver gráfica 5).

Otra forma de comparar las estimaciones se puede ver en el mapa con división municipal, donde se colorean las diferencias por rangos al restar de la estimación censal 2010 la del modelo EAP. En color negro se resaltan los municipios en los que esta diferencia absoluta es menor a mil pesos; con azul obscuro, las que están por abajo de la muestra censal entre mil y 3 mil pesos y con verde obscuro, entre mil y 3 mil pesos por arriba de la misma; en los demás colores (que son más claros) se aprecian los municipios con mayor diferencia. Es importante anotar que estos colores aparecen dispersos en todo el mapa, es decir, no hay una concentración que pueda indicar que el modelo sólo explica el comportamiento de alguna región en particular. Sin tomar en cuenta el negro, el predominio del verde obscuro indica que la estimación del modelo EAP es ligeramente menor a la de la muestra del CPV 2010.

Conclusiones

En el ejercicio realizado se ha mostrado que la técnica de EAP con modelo a nivel de área mejora en todos los casos la eficiencia de la estimación directa y, en muchos, ésta es radical. Es claro que la disponibilidad de datos auxiliares y la selección de un buen modelo es fundamental para lograr el éxito deseado. La evidente tendencia de las estimaciones con el modelo a nivel de área hacia las estimaciones de la muestra del CPV 2010 indica una buena congruencia de los datos con la ENIGH 2010 del mismo año, sabiendo que ambos levantamientos fueron totalmente independientes.
Aún hay muchas situaciones teóricas por resolver en la gran gama de aplicaciones para la EAP, sin embargo, existen técnicas probadas para aplicar con modelos espaciales, temporales, espacio-temporales tanto para modelos a nivel de área como de unidad.
La aplicación de modelos para la estimación en áreas pequeñas en el ámbito de las estadísticas oficiales es cada vez más utilizada en los institutos nacionales de estadística debido al gran ahorro de recursos que pueden representar. En México estas técnicas aún son poco conocidas y mucho menos aplicadas, sin embargo, al emplearlas es importante asegurarse de que los métodos utilizados, las hipótesis que subyacen a los modelos, así como la calidad de los resultados sean descritos de forma clara a los usuarios.

 

Fuentes

Belsley, D. A.Conditioning Diagnostics, Collinearity and Weak Data in Regression. 1st ed. New York, NY; John Wiley & Sons, Inc.; 1991.

Breusch, T. S. y A. R. Pagan. “A Simple Test for Heteroscedasticity and Random Coefficient Variation”, en: Econometrica. 47. 1979, pp. 1287-1294.

Cassel, C. M., C. E. Särndal y J. H. Wretman. Foundations of Inference in Survey Sampling. John Wiley & Sons, Inc., 1977.

Chatterjee, S., A. S. Hadi y B. Price. Regression Analysis by Example. 3rd ed. New York, NY; John Wiley & Sons, Inc.; 2000.

Cortés, F. y R. Rubalcava. El ingreso de los hogares. Serie Monografías Censales. Vol. VII. Aguascalientes, México; INEGI-El Colegio de México-Instituto de Investigaciones Sociales, UNAM; 1995.

Eurarea Consortium. Enhancing Small Area Estimation Techniques to meet European Needs. Project Reference. Vol. 2. UK, Explanatory Appendices, 2004.

Fay, R. E. y R. A. Herriot. “Estimates of Income for Small Places: an Application of James-Stein procedures to census data”, en: Journal of the American Statistical Association. Vol. 74. 1979, pp. 269-277.

Fox J. y S. Weisberg. An{R} Companion to Applied Regression. Second edition. Thousand Oaks, CA, Sage, 2011.
Consultado en bit.ly/4bTv10G

Harter, R. M. Small Area Estimation using nested-error and for sugar models and auxiliary data. Ph. D. Theses. Iowa State University, Estados Unidos de América (EE.UU.), 1983.

Hasen, M. H., W. G. Madow y B. J. Tepping. “An Evaluation of Model-Dependent and Probability-Sampling Inferences in Sample Surveys”, en: Journal of the American Statistical Association. Vol. 78. 1983, pp. 776-793.

Instituto Nacional de Estadística y Geografía (INEGI). Diseño de la muestra censal 2010. México, INEGI, 2011c.

_______ Encuesta Nacional de Ingresos y Gastos de los Hogares (ENIGH) 2010. Diseño muestral. México, INEGI, 2011c.

_______ Marco conceptual del Censo de Población y Vivienda 2010. México, INEGI, 2011c.

_______ Nueva construcción de la Encuesta Nacional de Ingresos y Gastos de los Hogares (ENIGH) 2010. Nueva construcción de ingresos y gastos. México, INEGI, 2011c.

Kackar, R. N. y D. A. Harville. “Approximations for Standard Errors of Estimators of Fixed and Random Effects in Mixed Linear Models”, en: Journal of the American Statistical Association. 79. 1984, pp. 853-862.

León, E. “Métodos de estimación de componentes de varianza en poblaciones. Una reseña histórica”, en: Revista Computarizada de Producción Porcina. 11, 2004, pp. 28-29.

Lohr, S. L. Muestreo: diseño y análisis. Thomson International, 2000.

Lumley, T. Survey: analysis of complex survey samples. R package version 3.24. 2011.

Molina, I. y Y. Marhuenda. Sae: Small Area Estimation. R package version 1.0. 2012.

Molina I. y J. N. K. Rao. Taller de Aplicación de Técnicas de Estimación para Áreas Pequeñas a las Ciencias Sociales. México, DF; Universidad Iberoamericana; 2012.

Nel Pacheco, P. Verificación de supuestos. Consultado en bit.ly/4dg0ACZ el 25 de noviembre de 2013.

Neter, J., W. Wasserman y M. H. Kutner. Applied Linear Statistical Models. 3.a ed. Irwin, MA, 1990.

Pinheiro J. y D. Bates. Mixed Effects Models in S and S-PLUS. Corrected third printing. New York, Springer-Verlag, 2002.

Pinheiro, J.; D. Bates; S. DebRoy; D. Sarkar y the R Development Core Team. Nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.pp. 1-100. 2011.

Prasad, N. G. N. y J. N. K. Rao. “The Estimation of Mean Squared Errors of Small Area Estimators”, en: Journal of American Statistical Association. 85. 1990, pp. 163-171.

R Development Core Team. R: A language and environment for statistical computing. Viena, Austria, R Foundation for Statistical Computing, 2011. Consultado en bit.ly/3Si0q5Q

Rachudel, W. J. Multicollinearity once again. Cambridge, Harvard Institute of Economic Research, 1971.

Rao, J. N. K. Small Area Estimation. New Jersey, EE.UU.; Wiley Interscience; 2003.

Roso, V. M.; F. S. Schenkel; S. P. Miller y L. R. Schaeffer. Estimation of genetic effects in the presence of multicollinearity in multibreed beef cattle evaluation. Consultado en bit.ly/4bSE3eb el 25 de abril de 2014.

Särndal, C. E. “Design-Consistent versus Model-Dependent Estimators for Small Domains”, en: Journal of the American Statistical Association. Vol. 79. 1984, pp. 624-631.

Särndal, C. E. y M. A. Hidiroglou. “Small Domain Estimation: A Conditional Analysis”, en: Journal of the American Statistical Association. Vol. 84. 1989, pp. 266-275.

Suárez M. Estimación del ingreso promedio por vivienda en los municipios del estado de Sonora. Tesis. CIMAT, 2010.

Székely M.; L. López-Calva; A. Meléndez; E. Rascón y L. Rodríguez. Poniendo a la pobreza de ingresos y a la desigualdad en el mapa de México. 2006. Consultado en bit.ly/3WwKH5j

Tarozzi A. y A. Deaton. “Using Census and Survey Data to Estimate Poverty and Inequality for Small Areas”, en: The Review of Economics and Statistics. Vol. 91, núm. 4. Noviembre del 2009, pp 773-792.

Terry Therneau y original Splus->R port by Thomas Lumley. Survival: Survival analysis, including penalised likelihood. R package version 2.36-5. Consultado en bit.ly/46cmW61. 2011.

Uriel E., Multicolinealidad. Universidad de Valencia. Consultado en bit.ly/3LCIC1D el 26 de abril de 2014.

Velasco, C. Curso de Estadística. Capítulo 5. Madrid, Universidad Carlos III. Consultado en bit.ly/3Wyh1ow

Valencia, M. Estimación en modelos lineales mixtos con datos continuos usando transformaciones y distribuciones no normales. Tesis. Universidad Nacional de Colombia, 2010.

Venables, W. N. y B. D. Ripley. Modern Applied Statistics with S. Fourth edition. New York, Springer, 2002.

¹ Para efectos de exposición de este documento, al utilizar las palabras municipio o municipios se referirá también a la(s) delegación(es) del Distrito Federal.

² La estratificación se realiza tomando en cuenta la división política del país, el ámbito urbano y rural, el tamaño de las localidades y las características sociodemográficas de los habitantes y de equipamiento de las viviendas (ENIGH 2010. Diseño muestral).

³ Son agrupaciones de 80 a 160 viviendas con características diferenciadas de acuerdo con el ámbito al que pertenecen (ENIGH 2010. Diseño muestral).

⁴ En el diseño conceptual de la ENIGH 2010 se denomina a la variable como ingreso del trabajo y el marco conceptual del CPV 2010 la nombra ingreso por trabajo; en este artículo se homologa con este último.

⁵ El detalle de estas fuentes de ingreso se puede ver en la nueva construcción de la ENIGH 2010.

⁶ Su diseño es estratificado por conglomerados.

⁷ Surge al plantear que el estimador ML suele producir estimaciones sesgadas de la varianza porque no tiene en cuenta los grados de libertad que se pierden al estimar la media. Para evitar este problema, se factoriza la verosimilitud completa en dos partes independientes, una de las cuales no contiene la media, asumiendo que por usar esta parte de la verosimilitud no se pierde información con respecto a la verosimilitud completa, así, la restringida corresponde en realidad con la verosimilitud asociada a una combinación lineal de las observaciones, cuya media es nula (León E., 2004).

⁸ Escárcega, Campeche; Namiquipa, Chihuahua; Nuevo Casas Grandes, Chihuahua; Cuajimalpa de Morelos, Distrito Federal; Lerdo, Durango; Chilapa de Álvarez, Guerrero; Huixquilucan, estado de México; Santa Lucía del Camino, Oaxaca; Matlapa, San Luís Potosí; Angostura, Sinaloa; Nuevo Laredo, Tamaulipas; Nativitas, Tlaxcala y Fresnillo, Zacatecas.

⁹ Medias positivas derivadas de valores positivos con distribución normal o aproximadamente normal.