Métodos para la suavización de indicadores de mortalidad: aplicación al análisis de desigualdades en mortalidad en ciudades del Estado español (Proyecto MEDEA)

Methods to smooth mortality indicators: application to analysis of inequalities in mortality in Spanish cities (the MEDEA Project)

M. Antònia Barceló^a,b, Marc Saez^a,b, Gemma Cano-Serral^b,c, Miguel Ángel Martínez-Beneito^d, José Miguel Martínez^e, Carme Borrell^b,c,e, Ricardo Ocaña-Riola^f, Imanol Montoya^g, Montse Calvo^h, Gonzalo López-Abente^b,i, Maica Rodríguez-Sanz^b,c, Silvia Toro^f, José Tomás Alcalá^j, Carme Saurina^a,b, Pablo Sánchez-Villegas^f, Adolfo Figueiras^b,k

^aGrup de Recerca en Estadística, Economia Aplicada i Salut, GRECS, Universitat de Girona.
^bCIBER en Epidemiología y Salud Pública (CIBERESP).
^cAgència de Salut Pública de Barcelona, Barcelona, España ]]> ^dDirecció General de Salut Pública, Conselleria de Sanitat, València, España
^eUniversitat Pompeu Fabra, Barcelona, España
^fEscuela Andaluza de Salud Pública, Granada, España
^gUnidad de Investigación de Atención Primaria de Bizkaia, Osakidetza/Servicio Vasco de Salud, Vitoria, España
^hDepartamento de Sanidad, Gobierno Vasco, Vitoria, España
ⁱCentro Nacional de Epidemiología, Instituto de Salud Carlos III, Madrid, España
^jUniversidad de Zaragoza, Zaragoza, España
^kDepartamento de Medicina Preventiva, Universidad de Santiago de Compostela, Santiago de Compostela, A Coruña, España.

Dirección para correspondencia

RESUMEN

Aunque la experiencia en el estudio de las desigualdades en la mortalidad en las ciudades españolas es amplia, quedan grandes núcleos urbanos que no han sido investigados utilizando la sección censal como unidad de análisis territorial. En este contexto se sitúa el proyecto coordinado «Desigualdades socioeconómicas y medioambientales en la mortalidad en ciudades de España. Proyecto MEDEA», en el cual participan 10 grupos de investigadores de Andalucía, Aragón, Cataluña, Galicia, Madrid, Comunitat Valenciana y País Vasco. Cabe señalar cuatro particularidades: a) se utiliza como área geográfica básica la sección censal; b) se emplean métodos estadísticos que tienen en cuenta la estructura geográfica de la región de estudio para la estimación de riesgos; c) se aprovechan las oportunidades que ofrecen 3 fuentes de datos complementarias (información sobre contaminación atmosférica, información sobre contaminación industrial y registros de mortalidad), y d) se emprende un análisis coordinado de gran alcance, favorecido por la implantación de la redes temáticas de investigación. El objetivo de este trabajo es explicar los métodos para la suavización de indicadores de mortalidad en el proyecto MEDEA. El artículo se centra en la metodología y los resultados del modelo de mapa de enfermedades de Besag, York y Mollié (BYM). Aunque en el proyecto se han suavizado, mediante el modelo BYM, las razones de mortalidad estandarizadas (RME) correspondientes a 17 grandes grupos de causas de defunción y 28 causas específicas, aquí se aplica esta metodología a la mortalidad por cáncer de tráquea, de bronquios y de pulmón en ambos sexos en la ciudad de Barcelona durante el período 1996-2003. Como resultado se aprecia un diferente patrón geográfico en las RME suavizadas en ambos sexos. En los hombres se observan unas RME mayores que la unidad en los barrios con mayor privación socioeconómica. En las mujeres este patrón se observa en las zonas con un mayor nivel socioeconómico.

Palabras clave: Proyecto MEDEA. RME. Sección censal. Modelo BYM.

ABSTRACT

Although there is some experience in the study of mortality inequalities in Spanish cities, there are large urban centers that have not yet been investigated using the census tract as the unit of territorial analysis. The coordinated project «Socioeconomic and environmental inequalities in mortality in Spanish cities. The MEDEA project» was designed to fill this gap, with the participation of 10 groups of researchers in Andalusia, Aragon, Catalonia, Galicia, Madrid, Valencia, and the Basque Country. The MEDEA project has four distinguishing features: a) the census tract is used as the basic geographical area; b) statistical methods that include the geographical structure of the region under study are employed for risk estimation; c) data are drawn from three complementary data sources (information on air pollution, information on industrial pollution, and the records of mortality registrars), and d) a coordinated, large-scale analysis, favored by the implantation of coordinated research networks, is carried out. The main objective of the present study was to explain the methods for smoothing mortality indicators in the context of the MEDEA project. This study focusses on the methodology and the results of the Besag, York and Mollié model (BYM) in disease mapping. In the MEDEA project, standardized mortality ratios (SMR), corresponding to 17 large groups of causes of death and 28 specific causes, were smoothed by means of the BYM model; however, in the present study this methodology was applied to mortality due to cancer of the trachea, bronchi and lung in men and women in the city of Barcelona from 1996 to 2003. As a result of smoothing, a different geographical pattern for SMR in both genders was observed. In men, a SMR higher than unity was found in highly deprived areas. In contrast, in women, this pattern was observed in more affluent areas.

Key words: MEDEA project. SMR. Census tract. BYM model.

Introducción

En estos últimos años abundan los estudios que han considerado el área geográfica como un determinante de la salud. De hecho, la salud está afectada no sólo por determinantes individuales sino también por otros ecológicos o contextuales^1-3. Así, por ejemplo, es ampliamente conocido que la pobreza afecta a la salud de las personas de forma individual, pero el análisis del contexto permite averiguar si la pobreza del área de residencia también es un determinante de la salud. Además, la salud de un grupo es más que la suma de la salud de los miembros individuales, y la perspectiva brindada por estos estudios permite un abordaje más cercano a la salud pública que los estudios individuales⁴.

Durante la década de 1990 se llevaron a cabo algunos estudios ecológicos sobre desigualdades en mortalidad en ciudades españolas, como Barcelona^7-10, Granada¹¹, Málaga¹² y Valencia⁹. Posteriormente, se realizaron investigaciones en otras ciudades españolas que, al igual que las anteriores, fueron desarrolladas tomando como unidad territorial de análisis grandes divisiones administrativas o sanitarias, fundamentalmente distritos y zonas básicas de salud^13,14.

Actualmente, además de ser la unidad de desagregación más pequeña disponible con datos oficiales, la sección censal está reconocida como el área geográfica óptima para el estudio de la variabilidad espacial de resultados en salud dentro de las ciudades¹⁵. En ellas la población es más homogénea en cuanto a hábitos de vida y condiciones socioeconómicas, el entorno presenta características medioambientales similares, la información necesaria para llevar a cabo estudios ecológicos se encuentra habitualmente disponible en organismos oficiales y hay pocos problemas de confidencialidad asociados. Todas estas propiedades, entre otras, han fomentado la aceptación de la sección censal como la unidad geográfica más apropiada para este tipo de estudios^15,16. Este hecho ha propiciado la reciente aparición de estudios de desigualdades en salud en algunas ciudades de Cataluña y el País Vasco que, a diferencia de las investigaciones realizadas en años anteriores, utilizan la sección censal como unidad de análisis territorial^17,18. En la misma línea de trabajo se sitúa el proyecto coordinado «Atlas de la mortalidad en las capitales de provincia de Andalucía y Cataluña» (AMCAC), el cual ha permitido además estudiar de forma conjunta la distribución geográfica de la mortalidad y su relación con factores socioeconómicos en las 12 ciudades que constituyen las capitales de provincia de Andalucía y Cataluña^19-21. Como particularidad, además de utilizar la sección censal como unidad geográfica, el proyecto AMCAC usó una metodología común basada en modelos jerárquicos bayesianos que hasta el momento no se había aplicado en España para el estudio de las desigualdades en salud en el ámbito urbano.

A pesar de que la experiencia en el estudio de las desigualdades en la mortalidad en las ciudades españolas es amplia, aún quedan grandes núcleos urbanos que no han sido investigados utilizando la sección censal como unidad de análisis territorial. En este contexto se sitúa el proyecto coordinado «Desigualdades socioeconómicas y medioambientales en la mortalidad en ciudades de España. Proyecto MEDEA», financiado por el Fondo de Investigaciones Sanitarias del Instituto Carlos III (años 2005-2007). En dicho proyecto, 10 grupos de investigadores pertenecientes a 7 comunidades autónomas (Andalucía, Aragón, Cataluña, Galicia, Madrid, Comunitat Valenciana y País Vasco) se plantearon los siguientes objetivos:

- Describir las desigualdades en indicadores socioeconómicos y de mortalidad en las secciones censales de grandes ciudades del Estado español (Madrid, Barcelona, Valencia, Sevilla, Zaragoza, Málaga, Bilbao, Córdoba, Alicante y Vigo) durante los años 1996-2003.

- Analizar la relación entre las desigualdades en mortalidad y los indicadores socioeconómicos en las áreas estudiadas de cada una de las ciudades seleccionadas.

- Comparar las desigualdades en mortalidad y según el nivel socioeconómico entre las ciudades seleccionadas, teniendo en cuenta el sexo, los grupos de edad y las causas de defunción.

- Estudiar las desigualdades medioambientales en salud, así como sus relaciones con las desigualdades socioeconómicas. Tratar de conocer la influencia de factores medioambientales, como la contaminación ambiental e industrial, en los patrones de mortalidad y su posible asociación con las desigualdades socioeconómicas.

Cabe señalar 4 particularidades del Proyecto MEDEA. En primer lugar, se utiliza como área geográfica básica la sección censal. En segundo lugar, se emplean métodos estadísticos que tienen en cuenta la estructura geográfica de la región de estudio para la estimación de riesgos. Hasta la fecha, en el Estado español no se han utilizado unidades geográficas inframunicipales, excepto en el proyecto AMCAC y en el Atlas de mortalidad en áreas pequeñas de Cataluña²², donde se analizó Barcelona por áreas básicas de salud. Por otra parte, el uso de métodos que tienen en cuenta la estructura espacial sólo tiene como precedente el proyecto AMCAC. Aunque son relativamente recientes, estos métodos se han convertido en referente en la literatura especializada²³. De hecho, comienza a ser habitual su aplicación en estudios de la variabilidad geográfica^24-30. En tercer lugar, se aprovechan las oportunidades que ofrecen 3 fuentes de datos complementarias, a saber, la información sobre contaminación atmosférica, la información sobre contaminación industrial aportada por el Registro Europeo de Emisiones Contaminantes (EPER)³¹ y el desarrollo de los registros de mortalidad. Y por último se emprende, por primera vez en nuestro país, un análisis coordinado de gran alcance favorecido por la implantación de la redes temáticas de investigación.

Métodos

Epidemiología espacial

El estudio de las desigualdades en salud, teniendo en cuenta la perspectiva geográfica, se enmarca en el contexto de la epidemiología espacial. El término epidemiología espacial o geográfica se emplea para describir estudios sobre las causas y la prevención de las enfermedades utilizando diferentes perspectivas de análisis en las cuales la localización geográfica de los eventos es uno de los componentes fundamentales^32,33. Uno de sus principales objetivos es mostrar qué parte de la variación espacial de la distribución de la frecuencia de una enfermedad no está explicada por la distribución espacial de factores explicativos, ni se puede atribuir al azar. De hecho, a menudo, interesa encontrar pistas sobre algún factor de riesgo desconocido de esta enfermedad.

El enfoque de análisis adoptado en el proyecto MEDEA se enmarca en las dos primeras perspectivas en que, siguiendo a Elliott y Watenberg³⁴, podría dividirse la epidemiología espacial: mapas de enfermedades y estudios de asociación geográfica.

Análisis de la mortalidad en áreas pequeñas

No hay una definición exacta de lo que se considera un área geográfica pequeña, ya que esta definición dependería del contexto de estudio y del número de casos de la enfermedad analizada. Algunos autores han definido un área pequeña como cualquier región cuyo tamaño muestral en un estudio es demasiado pequeño como para producir estimaciones con un nivel de precisión aceptable; por este motivo, es necesario hacer uso de áreas relacionadas para realizar estimaciones más precisas³⁵. Otros autores consideran que las áreas pequeñas son las que presentan menos de 20 casos de enfermedad y con una población que debería ser al menos de 10.000 habitantes³⁶. La ventaja de las áreas pequeñas es que suelen ser más homogéneas en cuanto a la exposición, y una limitación es el reducido número de eventos.

Los mapas de enfermedades proporcionan un resumen visual de información geográfica compleja y permiten identificar patrones geográficos de la enfermedad que de otro modo podrían pasar inadvertidos en las presentaciones tabulares^21,37. De hecho, se utilizan con propósitos descriptivos; con el objetivo de generar hipótesis etiológicas; para la vigilancia de la salud pública, a fin de detectar áreas con un aparente mayor riesgo, y como ayuda en la definición de políticas de salud y de asignación de recursos. Son también útiles para localizar conglomerados, o agrupaciones específicas de enfermedades, así como para situar los resultados de otros estudios de focos contaminantes^21,33,37-39.

Los mapas de enfermedades más simples representan los casos o muertes observadas en cada área geográfica. Sin embargo, cualquier interpretación de la estructura geográfica de los casos está limitada por la falta de información sobre la distribución espacial de la población que podría estar «en riesgo» y, consecuentemente, ha generado esos casos observados. Por tanto, se prefiere la representación de tasas que permitan incorporar el efecto de la población en riesgo, en vez de representar los casos brutos. Sin embargo, el uso directo de las tasas brutas no permite la comparación entre distintas áreas, ya que las diferencias observadas entre ellas pueden deberse a factores de riesgo que no hayan sido tenidos en cuenta, como la edad. Una medida de la mortalidad que tiene en cuenta la estructura de edad es la tasa estandarizada por edad. Hay 2 métodos para la estandarización por edad, que se conocen como estandarización directa e indirecta. En la representación de los mapas de mortalidad se prefiere la utilización del método indirecto, que consiste en la comparación de las muertes observadas en un área con las esperadas si los riesgos para cada grupo de edad fueran los mismos que en cierta población de referencia. La razón observada/esperada se llama razón de mortalidad estandarizada (RME), que no es más que un estimador del riesgo relativo del área, es decir, del riesgo de muerte en relación con el grupo considerado de referencia.

Además, la variabilidad de los casos observados suele ser bastante mayor que la esperada en una distribución de Poisson, produciéndose lo que se denomina extravariabilidad. De hecho, cuando se dispone de datos espaciales es importante distinguir 2 fuentes de extravariación. En primer lugar, la fuente más importante suele ser la denominada «dependencia espacial» y es consecuencia de la correlación de la unidad espacial con unidades espaciales vecinas, generalmente las contiguas geográficamente. De este modo, las RME de áreas contiguas, o cercanas, son más similares que las RME de áreas distantes espacialmente. Parte de esta dependencia no es realmente una dependencia estructural, sino que se debe, principalmente, a la existencia de variables no controladas, es decir, no incluidas en el análisis. Respecto a la segunda fuente, debe asumirse la existencia de extravariación independiente e incorrelacionada espacialmente, denominada «heterogeneidad incorrelacionada» o «no espacial», debida a variables no observadas sin estructura espacial que podrían influir en el riesgo^15,39.

Métodos de suavización de las razones de mortalidad estandarizada

Con el fin de solucionar los problemas derivados de la utilización directa de las RME se han propuesto varias alternativas para «suavizarlas», es decir, reducir la extravariación. Se parte del supuesto de que los casos observados se distribuyen como una Poisson:

O_i ~ Poisson (μ_iE_i)

Entonces, el estimador máximo verosímil del riesgo relativo en cada área, û , no es más que la RME.

Los métodos de suavización, en resumen, implican la formulación de modelos complejos cuyos parámetros necesitan ser estimados. La primera aproximación (al menos desde el punto de vista cronológico) consiste en el método de estimación denominado bayesiano empírico^40-42. En segundo lugar, se han utilizado modelos jerárquicos bayesianos completos^23,43,44. En un contexto completamente bayesiano se asignan a todos los parámetros del modelo distribuciones a priori. Alternativamente, en el contexto bayesiano empírico, los parámetros del modelo pueden estimarse maximizando directamente su distribución a posteriori.

El método bayesiano empírico, además de suponer una distribución de probabilidad para los casos observados, habitualmente una Poisson, asume que el riesgo relativo también es una variable aleatoria. Se puede suponer que éste se distribuye según una determinada distribución de probabilidad (habitualmente una gamma⁴⁰, aunque también pueden ser otras como una log-normal), o bien se puede dejar libre, y no suponer ninguna distribución determinada. El modelo resultante, a pesar de su nombre, se estima utilizando un enfoque frecuentista, por máxima verosimilitud no completa, habitualmente mediante el algoritmo EM.

Hace ya 20 años, Clayton y Kaldor⁴⁰ mostraron que ambos métodos, el empírico y el bayesiano completo, aunque proporcionan diferentes valores para los estimadores del riesgo relativo, en términos absolutos, coinciden en el rango y en el orden del riesgo, en términos relativos. Además, utilizando una gran variedad de escenarios con datos simulados, así como aplicaciones reales, Lawson et al⁴² muestran que ambos métodos son los que proporcionan los resultados más consistentes respecto a otras alternativas²³. El problema más importante de los métodos bayesianos emp&i acute;ricos, en relación con los modelos completamente bayesianos, es que no son capaces de modelizar, al menos fácilmente, la dependencia (extravariabilidad) espacial, ya que no tienen en cuenta la incertidumbre inicial de los parámetros del modelo⁴³. Sin embargo, los estimadores bayesianos empíricos pueden ser útiles a la hora de realizar una valoración inicial de la distribución espacial de los datos. Esto se debe a que en comparación con los completamente bayesianos su estimación es computacionalmente menos exigente. De todos modos, en la actualidad, los modelos de efectos aleatorios, como el BYM, son los más utilizados²³.

Un aspecto que debe tenerse muy en cuenta, ya sea mediante la utilización de un método u otro, es el grado de suavización que se consiga³⁴. Este grado determina un compromiso entre una alta sensibilidad (áreas de exceso de riesgo correctamente identificadas) y una alta especificidad (áreas que no tienen exceso de riesgo correctamente identificadas)⁴⁵. Un resultado sensible pero no específico generará falsos positivos, mientras que uno específico pero no sensible no identificará las áreas con un riesgo alto. Tras un exhaustivo estudio de simulación, Richardson et al⁴⁵ mostraron que, con un riesgo relativo en torno a 2, se requieren como mínimo 20 casos esperados (cuando el riesgo relativo es de 2-3, «sólo» se requieren 5 casos esperados) para que el método de suavización proporcione resultados «suficientemente sensibles» (en el sentido de detectar correctamente áreas con un riesgo alto). Si en vez de utilizar las razones suavizadas se utiliza la probabilidad a posteriori, PrP (es decir, que las razones sean mayores que la unidad), alternativa preferida por estos autores, el resultado será «suficientemente específico» cuando la probabilidad sea mayor del 80%.

Estudios ecológicos de asociación geográfica

El objetivo en este tipo de estudios es examinar, desde un punto de vista geográfico (y ecológico) la relación entre las variables respuesta de salud y las exposiciones tanto a factores medioambientales de tipo ecológico (contaminación atmosférica, del agua y del suelo) como de tipo ocupacional, las condiciones socioeconómicas y demográficas, o las variables relacionadas con el estilo de vida (como el hábito tabáquico y la dieta).

Suele ser habitual, sobre todo con datos espaciales, que la extravariabilidad no desaparezca por completo. Puede ocurrir, por ejemplo, que se trate de un problema de escala. Así, los datos pueden ser «escasos» (sparseness), es decir, puede que no se observe ningún caso en (relativamente) bastantes áreas geográficas. Sin embargo, parte de la extravariabilidad puede atribuirse a la falta de control de la confusión. De hecho, si se han omitido variables confusoras del riesgo esperado (privación, otras variables socioeconómicas), éstas deberían incluirse en el análisis.

Precisamente, en el proyecto MEDEA se incluyen, como variables explicativas, varios indicadores socioeconómicos, así como variables medioambientales. En este sentido, Domínguez-Berjón et al⁴⁶ describen la metodología seguida en la construcción de un índice de privación por sección censal, que permite identificar las secciones con situaciones más desfavorables, y analizan la relación de este índice con la mortalidad general.

Métodos de análisis en el proyecto MEDEA

Con el fin de cumplir los objetivos del proyecto MEDEA, el análisis se ha realizado estratificando por sexo y por ciudad, según 17 grandes grupos de causas de defunción y 28 causas específicas (3 de ellas opcionales) consensuadas por los miembros del proyecto (tabla 1).

Las fuentes de información utilizadas fueron el INE y los Registros de Mortalidad de las ciudades correspondientes. Como población de estudio se consideraron todas las personas residentes en estas ciudades. El período de estudio considerado fue 1996-2003.

En el caso del proyecto MEDEA, el modelo BYM se especifica como un modelo lineal generalizado mixto (GLMM) con variable respuesta de Poisson y considerando como offset los casos esperados³³:

O_i ~ Poisson (μ_iE_i)

Ln(μ_i) = α + Ln (E_i) + v_i + S_i

donde O_i denota los casos observados en la sección censal i, E_i son los casos esperados en esa sección censal, μ_i es el riesgo relativo en la sección i, v_i es el efecto aleatorio no espacial y S_i el efecto aleatorio espacial.

El efecto aleatorio no espacial, también denominado heterogeneidad, se supone distribuido normalmente con media cero y varianza constante. Para el efecto aleatorio, que recoge la variabilidad espacial, se utiliza un modelo condicional autorregresivo CAR^47,48. Esta aproximación, la más habitual y la computacionalmente más sencilla, aproxima la dependencia espacial como un promedio del efecto espacial de sus áreas vecinas³³.

En el proyecto MEDEA, las áreas consideradas son las secciones censales de cada una de las ciudades, y por vecindad se entiende contigüidad.

La estimación bayesiana completa requiere que se proporcione la función de verosimilitud, es decir, que se especifique el modelo, BYM en el caso del proyecto MEDEA, y las distribuciones a priori de los parámetros de interés. En concreto, distribuciones a priori para los efectos aleatorios (normal intercambiable para la heterogeneidad y CAR para la dependencia espacial), para los parámetros del modelo, si los hubiera (p. ej.,el asociado a alguna variable explicativa), así como para la desviación estándar (DE) de los efectos aleatorios.

Una vez especificado el modelo (la verosimilitud) y las distribuciones a priori, se obtiene la distribución a posteriori. Cuando se deben estimar varios nodos estocásticos, como en el caso del proyecto MEDEA (RME, DE de los efectos aleatorios, etc.), deben utilizarse integrales múltiples. Al no poderse obtener los resultados de una manera analítica, debe recurrirse a métodos de simulación. En particular, se utiliza el muestreo dependiente a partir de cadenas de Markov, MCMC^50,51. Hay varios algoritmos para diseñar cadenas de Markov. El más conocido es el denominado Gibbs sampling ^52-54, que es uno de los que utiliza el programa WinBUGS⁵⁵.

Como consecuencia de su forma iterativa, la convergencia de las simulaciones resulta crucial a la hora de validar los resultados obtenidos mediante estos métodos bayesianos completos. Por este motivo, en el proyecto MEDEA, para cada nodo estimado, se utilizan 3 cadenas y se ejecuta un gran número de simulaciones (50.000 como mínimo), de las que se «desprecian» el 10% de las iniciales (procedimiento denominado burn-in), y cada cierto número de iteraciones se guardan los nodos estimados, descartándose el resto de iteraciones con el fin de evitar la autocorrelación que se da en las cadenas (procedimiento denominado thinning).

Una vez obtenidos los resultados del modelo, la convergencia se comprueba en detalle mediante el estadístico de Brooks-Gelman-Rubin⁵⁶ (denominado R-hat en la librería R2WinBUGS en el entorno de software libre R) y el tamaño muestral efectivo de las cadenas^57,58 (n.eff en la librería R2WinBUGS en el entorno de software libre R). Como criterios de convergencia se acordó que el R-hat fuese menor que 1,1 y el n.eff mayor que 100 en las estimaciones de las DE de los efectos aleatorios espacial y heterogéneo, ya que son los parámetros más problemáticos en cuanto a la convergencia. Si no se cumplían dichos criterios se aumentaba el número de simulaciones hasta que se produjese la convergencia.

Una vez han convergido las simulaciones de las desviaciones de los efectos heterogéneo y espacial, se obtiene, como resultado, una muestra de la distribución posterior para las RME suavizadas, la PrP (probabilidad de que la RME suavizada sea mayor que la unidad), así como diversos histogramas, con el fin de comprobar la suavización y la distribución de frecuencias de las RME suavizadas.

Resultados

Aunque en el proyecto MEDEA se ha aplicado la metodología descrita a la totalidad de las causas contenidas en la tabla 1 y para todas las ciudades que forman parte del proyecto, en este artículo la metodología anteriormente expuesta se aplica, a modo de ejemplo, a la mortalidad por cáncer de tráquea, de bronquios y de pulmón (códigos de la Clasificación Internacional de Enfermedades, novena revisión [CIE-9] 162, y décima [CIE-10] C33, C34) en ambos sexos en la ciudad de Barcelona durante el período 1996-2003.

Como se ha mencionado, la unidad de análisis utilizada es la sección censal. El proyecto MEDEA utiliza el seccionado censal correspondiente al año 2001. En ese año, Barcelona contaba con 1.491 secciones censales, con una población media de 1.009 habitantes por sección (mediana de 993 habitantes); la sección menos poblada contaba con 91 habitantes y la más poblada con 7.003.

El 3,76% de las secciones de Barcelona en el caso de los hombres y el 53,92% en el caso de las mujeres no presentaron ninguna defunción por cáncer de pulmón en el período considerado. El número medio de defunciones en las secciones que presentaron al menos un fallecido fue de 4,11 en el caso de los hombres y de 1,42 en las mujeres (mediana de 4 y 1, respectivamente), con una DE de 2,57 en el caso de los hombres y de 0,65 en las mujeres.

Estratificando por sexo, las RME fueron suavizadas mediante el modelo BYM (especificado en la fórmula anterior), utilizando los métodos de análisis descritos anteriormente. Se utilizaron, además de los entornos de software libre R (versión 2.4.1)⁵⁹ y Win-BUGS⁵⁵, las librerías R2WinBUGS⁶⁰ (que permite conectar ambos entornos), maptools⁶¹ (herramientas para lectura y manipulación de objetos espaciales) y RColorBrewer⁶² (paletas para facilitar la coloración de mapas). La sintaxis y los códigos utilizados pueden ser consultados en la página web del grupo GRECS⁶³.

En el ajuste se utilizaron 3 cadenas, 200.000 simulaciones, con un burn-in de 20.000 y un thinning igual a 30. Los modelos cumplieron los criterios de convergencia (R-hat de las DE del efecto aleatorio espacial igual a 1,002282 en hombres y 1,000615 en mujeres, R-hat de las DE de la heterogeneidad igual a 1,094476 en hombres y 1,048364 en mujeres, y n.eff mayor que 100 en todos los casos).

En las figuras 2 y 3 se muestran los resultados de la suavización obtenida aplicando los métodos expuestos anteriormente, una vez comprobada la convergencia de las simulaciones. En la parte superior izquierda se muestra el mapa de la razón de mortalidad estandarizada suavizada; en la parte superior derecha, el mapa de la probabilidad a posteriori (probabilidad de que el riesgo relativo fuese mayor que 1). La gama de colores rojo muestra las secciones con un exceso de riesgo, y la de verdes las que presentan un defecto. Se consideró exceso de riesgo las secciones cuya PrP fuese mayor que 0,8 y defecto de riesgo aquellas cuya PrP fuese menor que 0,2. En la parte inferior, a la izquierda, se encuentra la función de densidad de probabilidad de la RME suavizada (con el mismo rango de colores que el mapa), y a la derecha la función de densidad de probabilidad de las RME crudas y de la razón estandarizada suavizada.

El resultado de la suavización puede observarse en las funciones de densidad de probabilidad de las RME crudas y suavizadas. En este sentido, obsérvese el caso de las mujeres, en las cuales el modelo BYM suaviza unas RME excesivamente inestables (compárese la figura 1 con la figura 3). Nótese asimismo el diferente patrón geográfico presente en las RME suavizadas en ambos sexos, patrones que, como muestran las PrP, se corresponden a agrupaciones de áreas con exceso o defecto de riesgo. En el caso de los hombres se observan unas RME mayores que la unidad en los barrios con mayor privación socioeconómica. Contrariamente, en las mujeres este patrón se observa en las zonas con un mayor nivel socioeconómico.

De hecho, este fenómeno podría explicarse, en gran parte, por el consumo de tabaco. En las décadas de 1980 y 1990 empezaron a fumar más los hombres de clases sociales menos favorecidas y, posteriormente, fueron las mujeres de clases sociales más privilegiadas las que iniciaron el hábito con más frecuencia^64-67.

Discusión

Cabría aquí señalar algunas limitaciones. En primer lugar, independientemente de la metodología empleada en los estudios ecológicos de desigualdades en salud, es importante destacar que las técnicas aplicadas son capaces de describir la distribución geográfic a de la mortalidad, pero no permiten explicar por qué el riesgo es más elevado en algunas secciones censales. Para ello, será necesario realizar estudios con datos individuales en los que se analicen los factores de riesgo relacionados con causas específicas de mortalidad.

De hecho, en los estudios ecológicos no se puede obtener una relación entre la exposición y el factor de riesgo estudiado, puesto que no se conoce si los individuos más expuestos a un determinado factor en cada área geográfica son los que han muerto. Cuando se extrapolan al ámbito individual los resultados obtenidos en estudios poblacionales aparece la llamada «falacia ecológica»⁶⁸. En este caso, los mapas obtenidos acerca de los factores socioeconómicos o ambientales nos ofrecerán una descripción de las secciones censales en las distintas zonas geográficas; observando dichos mapas se puede determinar cómo son las condiciones ambientales y sociales, así como la mortalidad, pero no se puede concluir que las personas que residen en una sección con una tasa de mortalidad alta y malas condiciones ambientales están más expuestas que las que viven en una sección con otras características. Ambas secciones son comparables a la población de referencia utilizada (España, 2001), pero no entre sí. Por otro lado, al realizar análisis con datos agregados, los riesgos para la salud estimados proporcionan una información realmente válida^69,70.

Cabe señalar también, como limitación del estudio, el denominado «efecto frontera». Este efecto se produce al ignorar la distribución del riesgo fuera del área de estudio, pero que puede influir en la distribución del riesgo dentro del área estudiada. Este efecto puede ser particularmente acusado en grandes ciudades en las que el cinturón metropolitano se entremezcla con otras ciudades circundantes.

Desde un punto de vista más metodológico, la sección censal no es más que una división meramente administrativa. Como ventajas de su utilización cabe mencionar que se trata de áreas de, más o menos, la misma población (lo que facilita la estandarización de las incertidumbres en la estimación de tasas); sobre todo, son áreas para las que se dispone de información procedente de los Censos de Población y Viviendas. Sin embargo, su uso no está libre de inconvenientes. El área de éstas es muy variable, por propia definición, lo que puede causar problemas en el análisis espacial y en la representación de mapas. Además, los límites de las secciones no permanecen constantes en el tiempo. Por último, suelen contener pocos casos (a veces muy pocos) de la enfermedad de interés⁷¹.

De hecho, un «exceso de ceros» en la distribución de los casos entre secciones censales, es decir, si hubiese muchas secciones sin ningún caso, podría implicar que la distribución de probabilidad no fuese, en realidad, una Poisson, sino una mixtura de distribuciones. Este problema podría limitar la aplicación de modelos como el BYM. Sin embargo, este aspecto presenta cierta controversia. Recientemente, Best et al²³ comparaban la mayoría de los modelos que se han aportado para obtener los riesgos relativos en estudios de mapas de enfermedades, incluido el BYM. Los resultados sugieren que el modelo BYM proporciona los mejores resultados, al menos por lo que se refiere al análisis de una única enfermedad. Por tanto, estos autores aconsejan el uso de este modelo para el análisis espacial en áreas pequeñas.

A pesar de sus limitaciones, el estudio de las desigualdades territoriales en mortalidad constituye un elemento importante para ampliar los conocimientos sobre la salud de la población. Los métodos de suavización resultan particularmente indicados en entornos urbanos, especialmente para secciones censales, ya que éstas son, salvo excepciones, pequeñas en extensión además de serlo en población. Por ello, los posibles factores de riesgo excederán habitualmente el ámbito geográfico de las secciones, por lo que resulta particularmente aconsejable la utilización de términos de efectos aleatorios con estructura espacial que basa las estimaciones del riesgo en las de las secciones vecinas.

Las conclusiones que puedan extraerse de estos estudios supondrán una aportación relevante para la planificación de recursos sanitarios, la implantación de programas de prevención y la disminución de las desigualdades en salud.

Respecto a las líneas de futuro del proyecto MEDEA, además de profundizar en diferentes aspectos medioambientales (comparación de desigualdades medioambientales en salud y su relación con factores socioeconómicos), se pretende realizar un atlas presentando las desigualdades en mortalidad en áreas pequeñas.

Agradecimientos

Estudio parcialmente financiado por el proyecto «Mortalidad en áreas pequeñas Españolas y Desigualdades socio-Económicas y Ambientales (MEDEA)»: PI04/0399 (GRECS, Universitat de Girona), PI04/2013 (ASPB, Barcelona), PI04/0388 (Departamento de Sanidad y Servicio Vasco de Salud, Gobierno Vasco), PI04/0041 (CNE, ISCIII), PI04/0170 (Conselleria de Sanitat, València); por la Red de Centros de Epidemiología y Salud Pública (FISS C03/09); y por el CIBER de Epidemiología y Salud Pública (CIBERESP).

Bibliografía

1. Krieger N. Epidemiology and the web of causation: has anyone seen the spider? Social Sci Med. 1994;38:887-903.        [ Links ]

2. Susser M, Susser E. Choosing a future for epidemiology (II). From black box to Chinese boxes and eco-epidemiology. Am J Public Health. 1996;86:674-7.        [ Links ]

3. Diez-Roux AV. Investigating neighborhood and area effects on health. Am J Public Health. 2001;91:1783-9.        [ Links ]

4. Susser M. The logic in ecological (I). The logic of analysis. Am J Public Health. 1994;84:825-9.        [ Links ]

5. Grupo de trabajo de la Sociedad Española de Epidemiología. La medición de la clase social en ciencias de la salud. Barcelona: SG Editores; 1995.        [ Links ]

6. Borrell C, Pasarín MI. Desigualdad en salud y territorio urbano. Gac Sanit. 2004;18:1-4.        [ Links ]

7. Costa J. Desigualtats en la mortalitat als barris de Barcelona: la seva relació amb l'atur, l'analfabetisme i la categoria professional [tesis doctoral]. Barcelona: Universitat Autònoma de Barcelona; 1991.        [ Links ]

8. Borrell C, Arias A. Desigualdades en mortalidad en los barrios de Barcelona. Gac Sanit. 1993;7:205-220.        [ Links ]

9. Arias A, Ribagliato M, Palumbo MA, Bellver R, Ashton J, Colomer C, et al. Desigualdades en salud en Barcelona y Valencia. Med Clin (Barc). 1993;100:281-7.        [ Links ]

10. Pasarín MI, Borrell C, Plasencia A. ¿Dos patrones de desigualdades en mortalidad en Barcelona? Gac Sanit. 1999; 13:431-40.        [ Links ]

11. Ruiz-Mariscal I, March-Cerdá J. Desigualdades en salud en los municipios periféricos de la aglomeración urbana de Granada. Aten Primaria. 1998;22:269-78.        [ Links ]

12. Martín FJ, March JC. Desigualdades sociales en salud en la ciudad de Málaga. Gac Sanit. 1992;6:198-206.        [ Links ]

13. Ruiz-Ramos M, Sánchez J, Garrucho G, Viciana F. Desigualdades en mortalidad en la ciudad de Sevilla. Gac Sanit. 2004;18:16-23.        [ Links ]

14. García-Gil C, Cruz-Rojo C, Álvarez-Girón M, Solano-Pares A. Health inequalities in Seville, Spain: use of indicators of social deprivation and mortality in small areas. Public Health. 2004;118:11-20.        [ Links ]

15. Lawson AB. Statistical methods in spatial epidemiology. 2nd ed. New Jersey: John Wiley & Sons; 2006.        [ Links ]

16. The Public Health Disparities Geocoding Project Monograph. Boston: Harvard School of Public Health 2004 [citado 20 Jul 2007]. Disponible en: http://www.hsph.harvard.edu/thegeocodingproject        [ Links ]

17. Domínguez-Berjón MF, Borrell C. Mortalidad y privación socioeconómica en las secciones censales y los distritos de Barcelona. Gac Sanit. 2005;19:363-9.        [ Links ]

18. Esnaola S, Aldasoro E, Ruiz R, Audicana C, Pérez Y, Calvo M. Desigualdades socioeconómicas en la mortalidad en la Comunidad Autónoma del País Vasco. Gac Sanit. 2006;20: 16-24.        [ Links ]

19. Ocaña-Riola R, Saez M, Sánchez-Cantalejo C, Barceló MA, Fernández A, Saurina C. Protocolo de Investigación del Atlas de Mortalidad de las Capitales de provincia de Andalucía y Cataluña (Proyecto AMCAC). Rev Esp Salud Pub. 2005;79: 613-20.        [ Links ]

20. Ocaña-Riola R, Saurina C, Fernández-Ajuria A, Lertxundi A, Sánchez-Cantalejo C, Sáez M, et al. Area deprivation and mortality in the provincial capital cities in Andalusia and Catalonia (Spain). J Epidemiol Comm Health. 2008;62:147-52.        [ Links ]

21. Ocaña-Riola R, Sánchez-Cantalejo C, Fernández-Ajuria A, Ruiz-Ramos M, Mayoral JM, Méndez C, et al. Atlas de Mortalidad de las Capitales de Provincia de Andalucía (1992-2002). Granada: EASP; 2007.        [ Links ]

22. Benach J, Martínez JM, Yasui Y, Borrell C, Pasarín MI, Español E, et al. Atlas de Mortalidad en Áreas Pequeñas en Cataluña (1984-1998). Barcelona: UPF, Fundació Jaume Bofill, Editorial Mediterrània; 2004.        [ Links ]

23. Best N, Richardson S, Thomson A. A comparison of Bayesian spatial models for disease mapping. Stat Met Med Res. 2005;14:35-9.        [ Links ]

24. Torrella A, Abellán JJ, Martínez-Beneito MA, Melchor I, Guallar A, Vanaclocha H, et al. Atlas de la mortalidad por cáncer en la Comunidad Valenciana 1987-1996 [citado 20 Jul 2007]. Disponible en: http://www.sp.san.gva.es/isum/Direct.jsp?ISUM_Shortcut=atlas_cancer        [ Links ]

25. Martínez-Beneito MA, López-Quílez A, Amador A, Melchor I, Botella P, Abellán C, et al. Atlas de la mortalidad de la Comunidad Valenciana 1991-2000. València: Generalitat Valenciana, Conselleria de Salut; 2005.        [ Links ]

26. Lope V, Pollán M, Pérez-Gómez B, Aragonés N, Ramis R, Gómez-Barroso D, et al. Municipal mortality due to thyroid cancer in Spain. BMC Public Health. 2006;6:302.        [ Links ]

27. López-Abente G, Ramis R, Pollán M, Aragonés N, Pérez-Gómez B, Gómez-Barroso D, et al. Atlas municipal de mortalidad por cáncer en España, 1989-1998. Madrid: Instituto de Salud Carlos III; 2006.        [ Links ]

28. Abellana R, Ascaso C. Métodos estadísticos para analizar el riesgo con patrones de distribución espacial. Med Clin (Barc). 2004;122:68-72.        [ Links ]

29. Arizo V, Botella Rocamora P, Abellán C, Martínez-Beneito MA, García Blasco MJ, Amador A, et al. Aproximación a las enfermedades raras en la Comunidad Valenciana, 1999-2003. Valencia: Conselleria de Sanitat, Generalitat Valenciana; 2005.        [ Links ]

30. Botella Rocamora P, Zurriaga O, Posada de la Paz M, Martínez-Beneito M, Prieto Bel E, Robustillo Rodela A, et al. Atlas nacional provincial de enfermedades raras. 1999-2003. Valencia: Repier; 2006.        [ Links ]

31. EPER 2006 [citado 20 Jul 2007]. Disponible en: http://eper.ec.europa.eu/eper/        [ Links ]

32. Thomas RW. Spatial epidemiology. Londres: Pion; 1990.        [ Links ]

33. Sáez M, Saurina C. Estadística y epidemiología espacial. Girona: Edicions a Petició; 2007.        [ Links ]

34. Elliott P, Wartenberg D. Spatial epidemiology: current approaches and future challenges. Environm Health Perspect. 2004;112:998-1006.        [ Links ]

35. Ghosh M, Rao JMK. Small area estimation: an appraisal. Statistical Science. 1994;9:55-76.        [ Links ]

36. Cuzick J, Elliott P. Small-area studies: purpose and methods. En: Elliott P, Cuzick J, English D, Stern R, editors. Geographical and environmental epidemiology: methods for small area studies. Oxford: Oxford University Press; 1992. p. 14-21.        [ Links ]

37. Quinn M. Introduction. En: Office for National Statistics (ONS). Cancer Atlas of the United Kingdom and Ireland 1991-2000. Londres: ONS; 2005. p. 1-5 [citado 20 Jul 2007]. Disponible en: http://www.statistics.gov.uk/statbase/Product.asp?vlnk=14059        [ Links ]

38. Wilkinson P, Thakrar B, Shaddick G, Stevenson S, Patternden S, Landon M, et al. Cancer incidence and mortality around the Pan Britannica Industries pesticide factory, Waltham Abbey. Occup Environ Med. 1997;54:101-7.        [ Links ]

39. Lawson AB, Browne WJ, Vidal-Rodeiro CL. Disease mapping with WinBUGS and MLwiN. Chichester: John Wiley & Sons; 2003.        [ Links ]

40. Clayton D, Kaldor J. Empirical Bayes estimates for age-standardized relative risks. Biometrics. 1987;43:671-81.        [ Links ]

41. Böhning D. Computer-assisted analysis of mixtures and applications: meta-analysis, disease mapping and others. Boca Raton: Chapman and Hall; 1999.        [ Links ]

42. Lawson AB, Biggeri AB, Böhning D, Lesaffre E, Viel JF, Clark A, et al. Disease mapping models: an empirical evaluation. Stat Med. 2000;19:2217-41.        [ Links ]

43. Besag J, York J, Mollie A. Bayesian image restoration, with two applications in spatial statistics. Ann Institute Statistical Mathematics. 1991;43:1-59.        [ Links ]

44. Mollie A. Bayesian mapping of disease. En: Gilks WR, Richardson S, Spiegelhalter DJ, editors. Markov chain Monte Carlo in practice. Nueva York: Chapman & Hall; 1996. p. 359-79.        [ Links ]

45. Richardson S, Thomson A, Best NG, Elliot P. Interpreting posterior relative risk estimates in disease-mapping studies. Environ Health Persp. 2004;112:1016-25.        [ Links ]

46. Domínguez-Berjón MF, Borrell C, Cano-Serral G, Esnaola S, Nolasco A, Pasarín MI, et al. Construcción de un índice de privación a partir de datos censales en grandes ciudades (Proyecto MEDEA). Barcelona: Proyecto MEDEA (mimeo); 2007.        [ Links ]

47. Besag J. Spatial interaction and the statistical analysis of lattice systems (with discussion). J Royal Statistical Society, Series B. 1974;36:192-236.        [ Links ]

48. Clayton DG, Bernadinelli L, Montomoli C. Spatial correlation in ecological analysis. Int J Epidemiol. 1993;22:1193-202.        [ Links ]

49. Gelman A. Prior distributions for variance parameters in hierarchical models. Bayesian Analysis. 2005;1:1-19.        [ Links ]

50. Gilks WR, Richardson S, Spiegelhalter DJ, editors. Markov chain Monte Carlo in practice. Nueva York: Chapman & Hall; 1996.        [ Links ]

51. Elliot P, Best NG. Geographical patterns of disease. En: Armitage P, Colton T, editors. Encyclopedia of biostatistics. 2nd ed. Chichester: John Wiley & Sons; 2005 [citado 20 Jul 2007]. Disponible en: http://eu.wiley.com/legacy/wileychi/eob/articles.html        [ Links ]

52. Geman S, Geman D. Stochastic relaxation, Gibbs distributions and the Bayesian restoration of images. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1984;12:609-28.        [ Links ]

53. Tanner MA, Wong WH. The calculation of posterior distributions by data augmentation. JASA. 1987;82:528-49.        [ Links ]

54. Gelfand AE, Smith AFM. Sampling-based approaches to calculating marginal densities. JASA. 1990;85:398-409.        [ Links ]

55. Spiegelhalter D, Thomas A, Best N, Lunn D. WinBUGS User Manual. Version 1.4. Enero 2003 [citado 20 Jul 2007]. Disponible en: http://www.mrc-bsu.cam.ac.uk/bugs        [ Links ]

56. Brooks SP, Gelman A. Alternative's methods for monitoring convergence of iterative simulations. J Computational Graphical Statistics. 1998;7:434-55.        [ Links ]

57. Carlin BP, Gelman A, Neal RM. Markov chain Monte Carlo in practice: a roundtable discussion. The American Statistician. 1998;52:93-100.        [ Links ]

58. Gelman A, Carlin JB, Stern HS, Rubin DB. Bayesian data analysis. Londres: Chapman & Hall/CRC; 2004.        [ Links ]

59. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2007, ISBN 3-900051-07-0 [citado 20 Jul 2007]. Disponible en: http://www.R-project.org        [ Links ]

60. Gelman A, Sturtz S, Ligges U. The R2WinBUGS Package [citado 20 Jul 2007]. Disponible en: http://www.stat.columbia.edu/~gelman/bugsR/        [ Links ]

61. Levin-Koh NJ, Bivand R. The maptools package [citado 20 Jul 2007]. Disponible en: http://cran.r-project.org/        [ Links ]

62. Neuwirth E. The RColorBrewer package [citado 20 Jul 2007]. Disponible en: http://cran.r-project.org/        [ Links ]

63. Grup de Recerca en Estadística, Economia Aplicada i Salut (GRECS). Universitat de Girona [citado 1 Sep 2007]. Disponible en: http://www.3.udg.edu/fcee/economia/english/grecs.htm

64. López AD, Hollinshaw EN, Piha T. A descriptive model of cigarette epidemic in developed countries. Tobacco Control. 1994;3:242-7.        [ Links ]

65. Regidor E, Gutiérrez-Fisac JL, Calle E, Navarro P, Domínguez V. Trends in cigarrete smoking in Spain by social class. Prev Med. 2001;33:241-8.        [ Links ]

66. Borrell C, Rué M, Pasarín MI, Rohlfs I, Ferrando J, Fernández E. Trends in social class inequalities in health status, health-related behaviors, and health services utilization in a Southern European urban area (1983-1994). Prev Med. 2000;31: 691-701.        [ Links ]

67. Giskes K, Kunst AE, Benach J, Borrell C, Costa G, Dahl E, et al. Trends in smoking behaviour between 1985 and 2000 in nine European countries by education. J Epidemiol Comm Health. 2005;59:395-401.        [ Links ]

68. Bolumar-Montrull F, Vioque-López J, Rebagliato-Ruso M. Estudios descriptivos ecológicos. Estudios transversales. En: Piedrola-Gil G, editor. Medicina preventiva y salud pública. 10ª ed. Barcelona: Masson, 2002. p. 87-100.        [ Links ]

69. Wakefield J, Shaddick G. Health-exposure modelling and the ecological fallacy. Biostatistics. 2006;7:438-55.        [ Links ]

70. Portnov BA, Dubnov J, Barchana M. On ecological fallacy, assessment errors stemming from misguided variable selection, and the effect of aggregation on the outcome of epidemiological study. J Exposure Sci Environ Epidemiol. 2007; 17:106-21.        [ Links ]

71. Briggs D, Fecht D, De Hoogh K. Census data issues for epidemiology and health risk assessment: experiences from the Small Area Health Statistics Unit. JRSS, Series A. 2007; 170:355-78.        [ Links ]

Dirección para correspondencia:
M. Antònia Barceló.
Grup de Recerca en Estadística, Economia Aplicada i Salut (GRECS) y CIBERESP.
Universitat de Girona. Campus de Montilivi. ]]> 17071 Girona. España.
Correo electrónico: antonia.barcelo@udg.edu

Recibido: 5 de septiembre de 2007.
Aceptado: 31 de enero de 2008.