Statistical Applications in Geographical Health Studies

Aquesta tesi està formada per dues parts relacionades amb l'estudi de la salut d'una regió geogràfica dividida en un conjunt de zones (àrees petites). La primera part es basa en un estudi amb informació de salut agregada per cadascuna de les àrees que formen la regió d'estudi. En concret, es tracta d'una aplicació de mapes de salut (disease mapping), que utilitza mètodes Bayesians empírics per generar un Atles de mortalitat en àrees petites de Catalunya en el període 1984-1998. La segona part utilitza una nova perspectiva basada en la integració de les dades agregades i individuals de salut per cadascuna de les zones que formen la regió d'estudi, mitjançant equacions d'estimació (estimating equations). Aquesta nova perspectiva és una extensió de la regressió geogràfica.L'elaboració de la primera part d'aquesta tesi està justificada per diferents raons. En primer lloc, els atles de salut i en general els mapes d'indicadors de salut, ens han mostrat la seva gran utilitat per identificar les localitzacions geogràfiques de les malalties, formular hipòtesis sobre les causes de la malaltia i monitoritzar intervencions en salut pública. En segon lloc, els atles de mortalitat en àrees petites presenten la distribució del risc relatiu per les causes de mortalitat més importants utilitzant mapes amb un alt nivell de resolució geogràfica.El primer objectiu d'aquesta tesi va ser construir un atles de mortalitat en 289 àrees petites (municipis o municipis agregats) de la Comunitat Autònoma de Catalunya i 66 àrees bàsiques de salut de la ciutat de Barcelona (l'àrea petita analitzada amb una major població) per al període 1984-1998. Per obtenir els indicadors de salut en àrees petites s'han utilitzat mètodes Bayesians. Aquests mapes presenten, en un format de doble pàgina, els riscs relatius ajustats per edat, les àrees significatives d'alt i baix risc, el risc relatiu de la ciutat de Barcelona respecte a Catalunya i internament respecte a Barcelona, el risc relatiu per grups d'edat (0-64 i 65) i addicionalment l'evolució temporal del risc relatiu en cada àrea resumida en un únic mapa. En concret, per estudiar l'evolució del risc relatiu de mortalitat s'inclou: 1) l'evolució del risc relatiu en el període d'estudi de cada àrea comparada amb la tendència global de Catalunya i 2) l'evolució absoluta del risc relatiu a cada àrea. Segons el nostre coneixement, aquesta és la primera vegada que aquests dos tipus d'informació es combinen en un únic mapa. A més, aquest és el primer Atles que presenta informació sobre la distribució geogràfica de zones que formen àrees petites de gran població, com ciutats d'un país, i inclou l'esperança de vida obtinguda amb mètodes Bayesians empírics.La segona part d'aquesta tesi és útil per estudis epidemiològics on s'inclouen variables d'exposició i confusió que poden tenir diferents fonts de variabilitat (variabilitat dins les poblacions i entre les poblacions). Específicament, els anàlisis individuals que valoren la relació entre la malaltia i l'exposició dins d'una població són útils quan l'exposició presenta variabilitat dins la població. Quan aquesta variabilitat és limitada, la força dels anàlisis individuals es debilita. En aquesta situació, un anàlisis de dades agregades de la malaltia entre poblacions, amb una mostra de dades individuals d'exposició, pot ser eficaç en l'estimació de l'efecte d'exposició si aquest presenta gran variabilitat entre poblacions. No obstant, encara que es pugui conèixer quina de les dues variacions domina en la variable d'exposició, es poden considerar conjuntament variables d'exposició i/o confusió amb diferents tipus de variació. El segon objectiu d'aquesta tesi va ser considerar una nova perspectiva, combinació dels anàlisis de dades individuals i agregades, basat en equacions d'estimació (perspectiva population-based estimating equation (PBEE)). En funció de la variabilitat que domina en la exposició, la anàlisis proposada pren força de la perspectiva basada en dades individuals i agrades de salut, per estimar els efectes d'exposició. Es van realitzar estudis de simulació en diferents escenaris per a mostrar el poder de la perspectiva proposada en l'estimació dels efectes d'exposició d'interès.Finalment, esperem que els mètodes i els diferents aspectes utilitzats en aquesta tesi puguin ser d'utilitat per a aquells investigadors que vulguin millorar l'estudi de la salut a l'espai i temps.Esta tesis esta formada por dos partes relacionadas con el estudio de la salud en una región geográfica dividida en un conjunto de zonas (áreas pequeñas). La primera parte considera un estudio con información de salud agregada para cada una de las áreas que forman la región analizada. En concreto, se trata de una aplicación de mapas de salud (disease mapping), consistente en el uso de métodos Bayesianos empíricos para generar un Atlas de mortalidad en áreas pequeñas de Cataluña en el periodo 1984-1998. La segunda parte considera un nuevo enfoque que realiza una integración de los datos agregados e individuales de salud para cada una de las zonas que forman la región en estudio, mediante ecuaciones de estimación (estimating equations). Se considera que este nuevo enfoque es una extensión de la regresión geográfica. La elaboración de la primera parte de esta tesis esta justificada por diferentes razones. Primero, los atlas de salud y en general los mapas de indicadores de salud, han mostrado su gran utilidad para identificar localizaciones geográficas de las enfermedades, formular hipótesis sobre las causas de la enfermedad y monitorizar intervenciones en salud pública. En segundo lugar, los atlas de mortalidad en áreas pequeñas presentan la distribución del riesgo relativo para las causas de mortalidad más importantes usando mapas con un alto nivel de resolución geográfica. El primer objetivo de esta tesis fue construir un atlas de mortalidad en 289 áreas pequeñas (municipios o municipios agregados) de la Comunidad Autónoma de Cataluña y 66 áreas básicas de salud de la ciudad de Barcelona (el área pequeña analizada con mayor población) para el periodo 1984-1998. Para obtener los indicadores de salud en las áreas pequeñas se han aplicado métodos Bayesianos. Estos mapas presentan, en un formato de página doble, los riesgos relativos ajustados por edad, las áreas significativas de alto y bajo riesgo, el riesgo relativo de la ciudad de Barcelona con respecto a Cataluña e internamente con respecto a Barcelona, el riesgo relativo por grupos de edad (0-64 y 65) y adicionalmente la evolución temporal del riesgo relativo en cada área resumida en un único mapa. En concreto, para estudiar la evolución del riesgo relativo de mortalidad se incluye: 1) la evolución del riesgo relativo en el periodo de estudio de cada área comparada con la tendencia global de Cataluña y 2) la evolución absoluta del riesgo relativo en cada área. Según nuestro conocimiento, esta es la primera vez que ambos tipos de información se combinan en un único mapa. Además, este es el primer Atlas que presenta información sobre la distribución geográfica de zonas que forman áreas pequeñas de gran población, como ciudades de un país, e incluye la esperanza de vida obtenida mediante métodos Bayesianos empíricos. La segunda parte de esta tesis es útil en estudios epidemiológicos donde se incluyen variables de exposición y confusión que pueden tener diferentes fuentes de variabilidad (variabilidad dentro de las poblaciones y entre poblaciones). Específicamente, los análisis individuales que valoran la relación entre enfermedad y exposición dentro de una población son útiles cuando la exposición presenta variabilidad dentro de la población. Cuando dicha variabilidad es limitada el poder de los análisis individuales se reduce. En esta situación, un análisis de datos agregados de enfermedad entre poblaciones, con una muestra de datos individuales de exposición, puede ser eficaz en la estimación del efecto de exposición si este presenta gran variabilidad entre poblaciones. No obstante, aunque se pueda conocer cual de las dos variaciones domina en la variable de exposición, se pueden considerar conjuntamente variables de exposición y/o confusión con diferentes tipos de variación. El segundo objetivo de esta tesis fue considerar un nuevo enfoque, combinación de los análisis de datos individuales y agregados, basado en ecuaciones de estimación (enfoque population-based estimating equation (PBEE)). Dependiendo de la variabilidad que domina en dicha exposición, el análisis propuesto toma fuerza de los enfoques basados en datos individuales y agregados de salud, para estimar los efectos de exposición. Estudios de simulación bajo diferentes escenarios fueron realizados para mostrar el poder del enfoque propuesto en la estimación de los efectos de exposición de interés.Finalmente, esperamos que los métodos y diferentes aspectos empleados en esta tesis puedan ser de utilidad para aquellos investigadores que quieran mejorar el estudio de la salud en el espacio y en el tiempo.This thesis consists of two related parts based on the study of health in a geographical region divided in a set of zones (small areas). The first part considers studies based on health information aggregated for each area into which the region under study has been divided. Specifically, it is a disease mapping application, based on generation of an Atlas of mortality in small areas of Catalonia over the period 1984-1998, using empirical Bayes methods. The second part considers an innovative approach, based on an integration of aggregated and individual health data in each of the zones of the region under study, using an estimating equation approach. Specifically, we consider this new approach as an extension of geographical regression. The elaboration of the first part of this thesis is justified for different reasons. First, health atlases and the mapping of health indicators in general, has demonstrated its great utility in identifying geographical localizations of health problems, in formulation of hypotheses about disease causes, and in monitoring public health interventions. Second, most atlases of mortality at the small area level present patterns of relative mortality risk for the most important causes of death using maps with a high level of geographical resolution. The first goal of this thesis was to construct a mortality Atlas involving a decomposition of the Autonomous Community of Catalonia into 289 small areas (municipalities or aggregates thereof) and 66 primary health areas of Barcelona city (being a small area but with a large population) for the period 1984-1998. For Catalonia as a whole, these maps presented, using a double-page format, the age adjusted relative risk, significantly high and low relative risk areas, relative risk in Barcelona City with respect to Catalonia and internally with respect to Barcelona, relative risk by age group (0-64 and 65) and additionally the relative risk evolution over time in each area summarized in an single map, using spatial and temporal information modeled through Bayesian methods. Specifically, the atlas uses a strategy to include both: 1) relative risk evolution throughout the study period of each area compared to the average trend for all Catalonia and 2) the absolute relative risk evolution of each area. To our knowledge, this is the first time that both types of information have been combined in a single map. In addition, this is the first Atlas that presents information about geographical patterns in zones within small areas having a large population such as the cities of a country and includes life expectancy obtained with an empirical Bayes approach.The second part of this thesis can be useful in epidemiological studies where we include exposure and confounding variables that may have different sources of within and between-population variability. Specifically, analyses of individual disease-exposure data within a population are useful when exposure of interest varies sufficiently within the population. When the within-population variance of exposure is limited power of the individual-data analysis within a population is reduced. In such situations, aggregated-data analyses of disease data across populations, with a sample of individual exposure data from populations, can be powerful in estimating the exposure effect if between-population variation of exposure is large. However, although we may have knowledge of which variations dominate in each variable, exposure and/or confounding variables with different types of variation can be considered jointly. The second goal of this thesis was to consider a new analytical framework that is a combination of the individual- and aggregated-data analyses, based on an estimating equation approach ("population-based estimating equation" (PBEE) approach). The proposed analysis utilizes strengths from individual and aggregated health data approaches in the estimation of the exposure effect of interest, depending on which of the exposure variations (within- vs. between-population) dominates. Simulation studies under different scenarios were performed to show the strengths of the proposed approach in the estimation of the exposure effects of interest.Finally, we hope that some of the methods and topics employed may be of use to researchers who want to improve the study of health in space and time

Tendencia temporal de la incapacidad laboral permanente en España (1992-2010)

Fundamentos: La incapacidad permanente (IP) es un problema de salud laboral en el que intervienen factores laborales, sanitarios o legales, los cuales han variado en España a lo largo de los últimos años. El objetivo del presente estudio es describir la tendencia temporal de la IP en España desde 1992 a 2010 por gravedad de la incapacidad y diferentes variables sociodemográficas y de afiliación a la seguridad social. Métodos: A partir de los datos de la Seguridad Social sobre afiliación y altas de IP en España se describió la tendencia temporal de la incidencia anual de la IP en el periodo comprendido entre 1992 y 2010. Para valorar las tendencias se estimó el porcentaje de cambio anual (PCA) y el intervalo de confianza al 95% mediante un modelo de regresión binomial negativa. Resultados: La tendencia de la incidencia de IP se subdividió en tres periodos. El primero (1992 a 1997) presentó una incidencia aproximadamente constante de 75 casos por 10.000 trabajadores (PCA1= -0,1; IC 95%:-1,6 a 1,4). Seguidamente las tasas descendieron durante un período de tres años (PCA2:-12,1%; IC 95%:-13,9 a -10,4), situándose alrededor de 50 nuevos casos por 10.000 trabajadores en el año 2000. A partir de este año la incidencia de IP se mantuvo aproximadamente constante hasta el año 2010 (PCA3=0,4; IC 95%:-0,2 a 1,0). Conclusiones: A pesar del incremento de trabajadores afiliados a la seguridad social en el periodo estudiado, la incidencia de IP se mantiene bastante estable, con un descenso entre 1997 y 2000 que se podría atribuir a las medidas de evaluación y seguimiento de las IP adoptados en 1995

Tendencia temporal de la incapacidad laboral permanente en España (1992-2010)

Fundamentos: La incapacidad permanente (IP) es un problema de salud laboral en el que intervienen factores laborales, sanitarios o legales, los cuales han variado en España a lo largo de los últimos años. El objetivo del presente estudio es describir la tendencia temporal de la IP en España desde 1992 a 2010 por gravedad de la incapacidad y diferentes variables sociodemográficas y de afiliación a la seguridad social. Métodos: A partir de los datos de la Seguridad Social sobre afiliación y altas de IP en España se describió la tendencia temporal de la incidencia anual de la IP en el periodo comprendido entre 1992 y 2010. Para valorar las tendencias se estimó el porcentaje de cambio anual (PCA) y el intervalo de confianza al 95% mediante un modelo de regresión binomial negativa. Resultados: La tendencia de la incidencia de IP se subdividió en tres periodos. El primero (1992 a 1997) presentó una incidencia aproximadamente constante de 75 casos por 10.000 trabajadores (PCA1= -0,1; IC 95%:-1,6 a 1,4). Seguidamente las tasas descendieron durante un período de tres años (PCA2:-12,1%; IC 95%:-13,9 a -10,4), situándose alrededor de 50 nuevos casos por 10.000 trabajadores en el año 2000. A partir de este año la incidencia de IP se mantuvo aproximadamente constante hasta el año 2010 (PCA3=0,4; IC 95%:-0,2 a 1,0). Conclusiones: A pesar del incremento de trabajadores afiliados a la seguridad social en el periodo estudiado, la incidencia de IP se mantiene bastante estable, con un descenso entre 1997 y 2000 que se podría atribuir a las medidas de evaluación y seguimiento de las IP adoptados en 1995

Exposure to ergonomic risk by occupation in spanish working: a data exploitation of MatEmESp

Fundamentos: La exposición a riesgos ergonómicos es la principal causa de daños de origen laboral. El objetivo es describir la prevalencia de exposición a carga física por ocupación en población laboral española y su relación con las condiciones de empleo y características sociodemográficas de los trabajadores. Métodos: Se utilizó la información contenida en la Matriz Empleo-Exposición Española (MatEmESp) con datos entre 1997 y 2005. Se describieron las prevalencias de exposición a riesgos ergonómicos, identificándose las ocupaciones con mayores prevalencias y se analizó la relación con las condiciones de empleo y las características sociodemográficas por ocupación mediante coeficientes de correlación de Spearman y diagramas de dispersión. Resultados: Los movimientos repetidos fueron el riesgo ergonómico declarado con mayor frecuencia (prevalencias en torno al 60%). Las mayores prevalencias de exposición a posturas forzadas, manipulación de cargas, movimientos repetidos y trabajo sedentario se dieron en "albañiles" (96%), en "peones de construcción" (89%), en "trabajadores de artes gráficas" (95%) y "auxiliares administrativos" (98%), respectivamente. Encontramos una fuerte relación (p50%). Conclusiones: En general, la prevalencia de exposición a riesgos ergonómicos se mantiene muy elevada en el periodo analizado. Las ocupaciones menos cualificadas presentaron habitualmente mayores prevalencias de exposición

Monitoring social determinants of health inequalities: the impact of unemployment among vulnerable

Surveillance of social determinants of health inequalities is an essential but still underdeveloped issue in public health. Existing research has identified unemployment as an important social determinant of health inequalities. This cross-sectional study investigates the impact of unemployment on mental health outcomes among vulnerable groups, using the 2006 Catalonian Health Survey (N=8591). The authors estimate the prevalence ratios and differences (excess of prevalence) for poor mental health in the unemployed and employed, with 95 percent confidence intervals. After taking into account the interactions among social mechanisms of inequality and related factors, the authors identified seven vulnerable groups to monitor. Primary findings indicate that unemployment has a greater adverse effect on the mental health of male manual workers, single mothers, main-earner women, and manual workers without unemployment benefits for both sexes. Findings support the need to devote more research to the surveillance of unemployment as a social determinant of health inequalities, to identify additional unemployment indicators, and to consider how various social mechanisms of inequality interact with each other to produce health inequalities among vulnerable groups.This work was supported by the CIBER Epidemiologia y Salud Pública (CIBERESP), Spai

The value of using administrative data in public health research: the continuous working life Sample

La utilización de datos administrativos es práctica habitual en la investigación en salud pública. El objetivo de este trabajo es describir la Muestra continua de vidas laborales (MCVL) y su utilidad en la investigación en salud. La MCVL está construida sobre registros que generan todos los contactos con la Seguridad Social (laborales, incapacidades, etc.), más datos fiscales (percepciones dinerarias, ingresos, etc.) y del padrón (nivel educativo, país de nacimiento, etc.), sin que sea posible la identificación individual. La MCVL se inició en 2004 con el 4% (1,1 millón de personas) del total de población cotizante o perceptora de la Seguridad Social. Los datos de las personas que forman parte de la MCVL son actualizados anualmente, y se sustituyen las bajas producidas. Este diseño continuo permite generar una cohorte con información de la vida laboral y económica, y valorar su asociación con la incapacidad laboral. El supuesto de conexión con la historia clínica posibilitaría el análisis de otros resultados de salud

Migration & self-rated health: comparing occupational cohorts in California & Spain (MICASA and ITSAL Projects)

International migration is a growing global phenomenon. The magnitude of the global population living outside their countries of origin substantiates the value of considering potential public health issues and their population-wide burden. As migration has yet to be generally accepted as an exposure in and of itself, and encompasses a wide range of experiences and health effects, a measure of overall health is well suited for this research. This study compares self-rated health between two independently collected occupational cohorts as part of the ITSAL study in Spain and the MICASA study in Mendota, California, USA. We observed greater gender balance in the MICASA sample than the ITSAL sample, where there was a substantial male majority. Mexican-born workers in the MICASA sample tended to be older, less educated, and more likely to work in agriculture than their Moroccan-born counterparts in the ITSAL study. We also observed a higher prevalence of poor self-rated health in the MICASA sample compared to the ITSAL sample. Differences may be due to sampling and data collection issues, cultural issues and the subjectivity of self-rated health as an outcome, as well as actual health differences. Further research is needed to determine common and distinct migration-related public health issues.Funding for the MICASA study was provided by the National Institute for Occupational Safety and Health (2U500H007550 and RO1OH009293) and The California Endowment. Funding for the ITSAL project was provided by the Fondo de Investigación Sanitaria, Instituto de Salud Carlos III in 2005 (PI050497, PI052202, PI052334), 2006 (PI061701), and 2007 (PI0790470)

Migration & self-rated health: comparing occupational cohorts in California & Spain (MICASA and ITSAL Projects)

International migration is a growing global phenomenon. The magnitude of the global population living outside their countries of origin substantiates the value of considering potential public health issues and their population-wide burden. As migration has yet to be generally accepted as an exposure in and of itself, and encompasses a wide range of experiences and health effects, a measure of overall health is well suited for this research. This study compares self-rated health between two independently collected occupational cohorts as part of the ITSAL study in Spain and the MICASA study in Mendota, California, USA. We observed greater gender balance in the MICASA sample than the ITSAL sample, where there was a substantial male majority. Mexican-born workers in the MICASA sample tended to be older, less educated, and more likely to work in agriculture than their Moroccan-born counterparts in the ITSAL study. We also observed a higher prevalence of poor self-rated health in the MICASA sample compared to the ITSAL sample. Differences may be due to sampling and data collection issues, cultural issues and the subjectivity of self-rated health as an outcome, as well as actual health differences. Further research is needed to determine common and distinct migration-related public health issues.Funding for the MICASA study was provided by the National Institute for Occupational Safety and Health (2U500H007550 and RO1OH009293) and The California Endowment. Funding for the ITSAL project was provided by the Fondo de Investigación Sanitaria, Instituto de Salud Carlos III in 2005 (PI050497, PI052202, PI052334), 2006 (PI061701), and 2007 (PI0790470)

Estimación de la carga de enfermedad por incapacidad laboral permanente en España durante el período 2009-2012

Fundamentos: La incapacidad laboral permanente (IP) supone una importante carga social y económica. El objetivo de este estudio fue estimar la carga de enfermedad debida a la IP en España. Métodos: A partir de la Muestra Continua de Vidas Laborales (MCVL) se imputó como diagnóstico de la IP (2009-2012) el de la incapacidad temporal (IT) iniciada en 2009 en los individuos que cumplieron dos criterios: larga duración de la IT y corta duración del período entre la IT y la IP. Para los casos con diagnóstico imputado se calculó la carga de enfermedad por IP a partir de los años potenciales de vida laboral perdidos (APVLP), definidos como la diferencia entre la edad de jubilación (65 años) y la edad de inicio de la IP. Resultados: Según la MCVL, 163.135 (13,6%) personas sufrieron de IT en 2009, de las cuáles a 4.738 (0,39%) se le reconoció una IP entre 2009-2012. Se pudo establecer el diagnóstico de IP en 3.073 (64,9%) casos. Las causas más frecuentes de IP fueron los trastornos musculo-esqueléticos (27,4%) y mentales (11,6%) que, asimismo, produjeron mayor cantidad de APVLP por IP: 11.086 (26,5%) y 7.052 (16,9%), respectivamente. Los trastornos mentales destacaron por producir IP a edades más jóvenes (mediana=49 años). Conclusiones: Los trastornos mentales y los musculo-esqueléticos representaron la mayor carga de enfermedad en individuos en situación de IP

Estimación de la carga de enfermedad por incapacidad laboral permanente en España durante el período 2009-2012

Fundamentos: La incapacidad laboral permanente (IP) supone una importante carga social y económica. El objetivo de este estudio fue estimar la carga de enfermedad debida a la IP en España. Métodos: A partir de la Muestra Continua de Vidas Laborales (MCVL) se imputó como diagnóstico de la IP (2009-2012) el de la incapacidad temporal (IT) iniciada en 2009 en los individuos que cumplieron dos criterios: larga duración de la IT y corta duración del período entre la IT y la IP. Para los casos con diagnóstico imputado se calculó la carga de enfermedad por IP a partir de los años potenciales de vida laboral perdidos (APVLP), definidos como la diferencia entre la edad de jubilación (65 años) y la edad de inicio de la IP. Resultados: Según la MCVL, 163.135 (13,6%) personas sufrieron de IT en 2009, de las cuáles a 4.738 (0,39%) se le reconoció una IP entre 2009-2012. Se pudo establecer el diagnóstico de IP en 3.073 (64,9%) casos. Las causas más frecuentes de IP fueron los trastornos musculo-esqueléticos (27,4%) y mentales (11,6%) que, asimismo, produjeron mayor cantidad de APVLP por IP: 11.086 (26,5%) y 7.052 (16,9%), respectivamente. Los trastornos mentales destacaron por producir IP a edades más jóvenes (mediana=49 años). Conclusiones: Los trastornos mentales y los musculo-esqueléticos representaron la mayor carga de enfermedad en individuos en situación de IP