1
|
Gause EL, Schumacher AE, Ellyson AM, Withers SD, Mayer JD, Rowhani-Rahbar A. An introduction to bayesian spatial smoothing methods for disease mapping: modeling county firearm suicide mortality rates. Am J Epidemiol 2024; 193:1002-1009. [PMID: 38375682 DOI: 10.1093/aje/kwae005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Revised: 11/20/2023] [Accepted: 02/13/2024] [Indexed: 02/21/2024] Open
Abstract
This article introduces bayesian spatial smoothing models for disease mapping-a specific application of small area estimation where the full universe of data is known-to a wider audience of public health professionals using firearm suicide as a motivating example. Besag, York, and Mollié (BYM) Poisson spatial and space-time smoothing models were fitted to firearm suicide counts for the years 2014-2018. County raw death rates in 2018 ranged from 0 to 24.81 deaths per 10 000 people. However, the highest mortality rate was highly unstable, based on only 2 deaths in a population of approximately 800, and 80.5% of contiguous US counties experienced fewer than 10 firearm suicide deaths and were thus suppressed. Spatially smoothed county firearm suicide mortality estimates ranged from 0.06 to 4.05 deaths per 10 000 people and could be reported for all counties. The space-time smoothing model produced similar estimates with narrower credible intervals as it allowed counties to gain precision from adjacent neighbors and their own counts in adjacent years. bayesian spatial smoothing methods are a useful tool for evaluating spatial health disparities in small geographies where small numbers can result in highly variable rate estimates, and new estimation techniques in R software have made fitting these models more accessible to researchers.
Collapse
|
2
|
Wang F, Duan C, Li Y, Huang H, Shia BC. Spatiotemporal varying coefficient model for respiratory disease mapping in Taiwan. Biostatistics 2023; 25:40-56. [PMID: 36484310 DOI: 10.1093/biostatistics/kxac046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 11/13/2022] [Accepted: 11/16/2022] [Indexed: 12/17/2023] Open
Abstract
Respiratory diseases have been global public health problems for a long time. In recent years, air pollutants as important risk factors have drawn lots of attention. In this study, we investigate the influence of $\pm2.5$ (particulate matters in diameter less than 2.5 ${\rm{\mu }} m$) on hospital visit rates for respiratory diseases in Taiwan. To reveal the spatiotemporal pattern of data, we propose a Bayesian disease mapping model with spatially varying coefficients and a parametric temporal trend. Model fitting is conducted using the integrated nested Laplace approximation, which is a widely applied technique for large-scale data sets due to its high computational efficiency. The finite sample performance of the proposed method is studied through a series of simulations. As demonstrated by simulations, the proposed model can improve both the parameter estimation performance and the prediction performance. We apply the proposed model on the respiratory disease data in 328 third-level administrative regions in Taiwan and find significant associations between hospital visit rates and $\pm2.5$.
Collapse
Affiliation(s)
- Feifei Wang
- Center for Applied Statistics, School of Statistics, Renmin University of China, Beijing, 100872, China
| | - Congyuan Duan
- School of Mathematics, Sun Yat-Sen University, Guangdong, 510275, China
| | - Yang Li
- Center for Applied Statistics, School of Statistics, and Statistical Consulting Center, Renmin University of China, Beijing, 100872, China
| | - Hui Huang
- School of Mathematics, Sun Yat-Sen University, Guangdong, 510275, China
| | - Ben-Chang Shia
- AI Development Center of Taiwan Institute of Artificial Intelligence, Fu Jen Catholic University, New Taipei City, 24205, Taiwan
| |
Collapse
|
3
|
Chen CC, Lo GJ, Chan TC. Spatial Analysis on Supply and Demand of Adult Surgical Masks in Taipei Metropolitan Areas in the Early Phase of the COVID-19 Pandemic. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2022; 19:6704. [PMID: 35682289 PMCID: PMC9179980 DOI: 10.3390/ijerph19116704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Revised: 05/26/2022] [Accepted: 05/27/2022] [Indexed: 12/10/2022]
Abstract
This study aimed to assess the gap between the supply and demand of adult surgical masks under limited resources. Owing to the implementation of the real-name mask rationing system, the historical inventory data of aggregated mask consumption in a pharmacy during the early period of the COVID-19 outbreak (April and May 2020) in Taiwan were analyzed for supply-side analysis. We applied the Voronoi diagram and areal interpolation methods to delineate the average supply of customer counts from a pharmacy to a village (administrative level). On the other hand, the expected number of demand counts was estimated from the population data. The relative risk (RR) of supply, which is the average number of adults served per day divided by the expected number in a village, was modeled under a Bayesian hierarchical framework, including Poisson, negative binomial, Poisson spatial, and negative binomial spatial models. We observed that the number of pharmacies in a village is associated with an increasing supply, whereas the median annual per capita income of the village has an inverse relationship. Regarding land use percentages, percentages of the residential and the mixed areas in a village are negatively associated, while the school area percentage is positively associated with the supply in the Poisson spatial model. The corresponding uncertainty measurement: villages where the probability exceeds the risk of undersupply, that is, Pr (RR < 1), were also identified. The findings of the study may help health authorities to evaluate the spatial allocation of anti-epidemic resources, such as masks and rapid test kits, in small areas while identifying priority areas with the suspicion of undersupply in the beginning stages of outbreaks.
Collapse
Affiliation(s)
- Chien-Chou Chen
- Department of Statistics and Information Science, Fu Jen Catholic University, New Taipei City 242, Taiwan; (C.-C.C.); (G.-J.L.)
| | - Guo-Jun Lo
- Department of Statistics and Information Science, Fu Jen Catholic University, New Taipei City 242, Taiwan; (C.-C.C.); (G.-J.L.)
| | - Ta-Chien Chan
- Research Center for Humanities and Social Sciences, Academia Sinica, Taipei 115, Taiwan
- Institute of Public Health, School of Medicine, National Yang Ming Chiao Tung University, Taipei 112, Taiwan
| |
Collapse
|
4
|
Makarenko C, Pedro AS, Paiva NS, Souza-Santos R, Medronho RDA, Gibson G. Identificação de áreas de risco e fatores associados à epidemia de sarampo de 2019 no Estado de São Paulo, Brasil. CAD SAUDE PUBLICA 2022; 38:e00039222. [DOI: 10.1590/0102-311xpt039222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Accepted: 08/04/2022] [Indexed: 11/27/2022] Open
Abstract
O objetivo foi analisar a ocorrência de clusters e fatores associados ao ressurgimento de casos de sarampo da maior epidemia do período pós-eliminação, ocorrida no Estado de São Paulo, Brasil, em 2019. Fatores sociossanitários e assistenciais foram analisados por modelos de Poisson inflacionado de zero (ZIP) e ZIP com efeito espacial estruturado e não estruturado. A estatística de varredura SCAN foi usada para analisar a ocorrência de clusters de casos. Foram identificados clusters de casos de alto risco em municípios que compõem a região intermediária de São Paulo. No modelo ZIP, foram observadas como fatores de risco no nível municipal as variáveis chefes de domicílio menores de 18 anos (RR ajustado = 1,39; ICr95%: 1,27-1,53), desigualdade na distribuição de renda (RR ajustado = 36,67; ICr95%: 26,36-51,15), desocupação em maiores de 18 anos (RR ajustado = 1,10; ICr95%: 1,08-1,12) e iluminação pública inexistente (RR ajustado = 1,05; ICr95%: 1,04-1,05). Nos modelos ZIP com efeito espacial estruturado e não estruturado, foram identificados como fatores de risco os indicadores chefes de domicílio menores de 18 anos (RR ajustado = 1,36; ICr95%: 1,04-1,90) e desigualdade na distribuição dos rendimentos do trabalho (RR ajustado = 3,12; ICr95%: 1,02-9,48). Em ambos os modelos, a cobertura de agentes de saúde se apresentou como fator de proteção. Os achados reforçam a importância de intensificar as ações de vigilância de sarampo articuladas à Estratégia Saúde da Família, especialmente em áreas de maior vulnerabilidade social, para garantir coberturas vacinais equânimes e satisfatórias e reduzir o risco de reemergência da doença.
Collapse
Affiliation(s)
| | | | | | | | | | - Gerusa Gibson
- Universidade Federal do Rio de Janeiro, Brazil; Fundação Oswaldo Cruz, Brazil
| |
Collapse
|
5
|
Costa C, Santana P. Trends of amenable deaths due to healthcare within the European Union countries. Exploring the association with the economic crisis and education. SSM Popul Health 2021; 16:100982. [PMID: 34926783 PMCID: PMC8648806 DOI: 10.1016/j.ssmph.2021.100982] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 11/23/2021] [Accepted: 11/25/2021] [Indexed: 02/03/2023] Open
Abstract
The study of premature deaths from causes that are generally preventable given the current availability of healthcare - called amenable deaths due to healthcare - provides information on the quality of services. However, they are not only impacted by healthcare characteristics: other factors are also likely to influence. Therefore, identifying the association between amenable deaths due to healthcare and health determinants, such as education, might be the key to preventing these deaths in the future. Still unclear however, is how this works and how amenable deaths due to healthcare are distributed and evolve within the European Union (EU) below the national level. We therefore studied the geographical and temporal patterns of amenable deaths due to healthcare in the 259 EU regions from 1999 to 2016, including the 2007-2008 financial crisis and the post-2008 economic downturn, and identified whether any association with education exists. A cross-sectional ecological study was carried out. Using a hierarchical Bayesian model, we estimated the average smoothed Standardized Mortality Ratios (sSMR). A regression model was also applied to measure the relative risks (RR) at 95% credible intervals for cause-specific mortality association with education. Results show that amenable deaths due to healthcare decreased globally. Nevertheless, the decrease is not the same across all regions, and inequalities within countries do persist, with lower mortality ratios seen in regions from Central European countries and higher mortality ratios in regions from Eastern European countries. Also, the evolution trend reveals that after the financial crisis, the number of these deaths increased in regions across almost all EU countries. Moreover, educational disparities in mortality emerged, and a statistical association was found between amenable deaths due to healthcare and early exit from education and training. These results confirm that identifying and understanding the background of regional differences may lead to a better understanding of the amenable deaths due to healthcare and allow for the application of more effective policies.
Collapse
Affiliation(s)
- Claudia Costa
- Centre of Studies in Geography and Spatial Planning (CEGOT), Department of Geography and Tourism, University of Coimbra, Portugal
| | - Paula Santana
- Centre of Studies in Geography and Spatial Planning (CEGOT), Department of Geography and Tourism, University of Coimbra, Portugal
| |
Collapse
|
6
|
Carroll R, Prentice CR. Community vulnerability and mobility: What matters most in spatio-temporal modeling of the COVID-19 pandemic? Soc Sci Med 2021; 287:114395. [PMID: 34530217 PMCID: PMC8434688 DOI: 10.1016/j.socscimed.2021.114395] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 08/12/2021] [Accepted: 09/09/2021] [Indexed: 12/02/2022]
Abstract
Community vulnerability is widely viewed as an important aspect to consider when modeling disease. Although COVID-19 does disproportionately impact vulnerable populations, human behavior as measured by community mobility is equally influential in understanding disease spread. In this research, we seek to understand which of four composite measures perform best in explaining disease spread and mortality, and we explore the extent to which mobility account for variance in the outcomes of interest. We compare two community mobility measures, three composite measures of community vulnerability, and one composite measure that combines vulnerability and human behavior to assess their relative feasibility in modeling the US COVID-19 pandemic. Extensions – via temporally dependent fixed effect coefficients – of the commonly used Bayesian spatio-temporal Poisson disease mapping models are implemented and compared in terms of goodness of fit as well as estimate precision and viability. A comparison of goodness of fit measures nearly unanimously suggests the human behavior-based models are superior. The duration at residence mobility measure indicates two unique and seemingly inverse relationships between mobility and the COVID-19 pandemic: the findings indicate decreased COVID-19 presence with decreased mobility early in the pandemic and increased COVID-19 presence with decreased mobility later in the pandemic. The early indication is likely influenced by a large presence of state-issued stay at home orders and self-quarantine, while the later indication likely emerges as a consequence of holiday gatherings in a country under limited restrictions. This study implements innovative statistical methods and furnishes results that challenge the generally accepted notion that vulnerability and deprivation are key to understanding disparities in health outcomes. We show that human behavior is equally, if not more important to understanding disease spread. We encourage researchers to build upon the work we start here and continue to explore how other behaviors influence the spread of COVID-19.
Collapse
Affiliation(s)
- Rachel Carroll
- Department of Mathematics and Statistics, University of North Carolina Wilmington, 601 S College Rd., Wilmington, NC, USA.
| | - Christopher R Prentice
- Department of Public and International Affairs, University of North Carolina Wilmington, 601 S College Rd., Wilmington, NC, USA
| |
Collapse
|
7
|
Stipancic J, Racine EB, Labbe A, Saunier N, Miranda-Moreno L. Massive GNSS data for road safety analysis: Comparing crash models for several Canadian cities and data sources. ACCIDENT; ANALYSIS AND PREVENTION 2021; 159:106232. [PMID: 34186470 DOI: 10.1016/j.aap.2021.106232] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Revised: 05/26/2021] [Accepted: 05/31/2021] [Indexed: 06/13/2023]
Abstract
Mobile sensors are a useful data source with applications in several transportation fields. Though cost of collection, transmission, and storage has limited studies on driving data and safety, this can be overcome through usage-based insurance (UBI). In UBI programs, drivers are monitored, and their premiums are adjusted based on driver-level surrogate safety measures (SSMs) related to exposure and driving style. Contextual link-level SSMs (volume, speed, or density) could further improve discount calibration. This study quantifies relationships between contextual SSMs and crashes and includes the validation of previous results (correlations between SSMs and crashes and statistical models estimated using smartphone-collected data from Quebec City) and the comparison of three Canadian cities (using UBI data from Quebec City, Montreal, and Ottawa). Extracted SSMs were compared to large volumes of historical crash frequency data using Spearman's Rank Correlation Coefficient and then implemented into spatial Bayesian crash models. Results from the UBI data generally matched those from the previous study, with observed correlations mirroring previous results in direction (braking, congestion, and speed variation are positively associated with crash frequency while mean speed is negatively associated) while correlation strength was slightly higher. Furthermore, these results were consistent between cities. For the crash modelling, repeatability of previous results in Quebec City was moderately good for the UBI data. Importantly for large-scale implementation, models estimated using UBI data were largely consistent between cities. This work provides an important contribution to the existing literature, clearly demonstrating how contextual safety measures could be applied to benefit UBI practices.
Collapse
Affiliation(s)
- Joshua Stipancic
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec H3T 2A7, Canada.
| | - Etienne B Racine
- Intactlab - Data Science, Intact Insurance, Suite 100, 2020 Boulevard Robert-Bourassa, Montréal, Québec H3T 2A7, Canada.
| | - Aurélie Labbe
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec H3T 2A7, Canada; Interuniversity Research Centre on Enterprise Networks, Logistics and Transportation (CIRRELT) Pavillon André Aisenstadt, Room 3520 2920 Chemin de la Tour Université de Montréal, Montréal, Quebec H3T 1J4, Canada.
| | - Nicolas Saunier
- Department of Civil, Geological and Mining Engineering, Polytechnique Montréal, C.P. 6079, succ. Centre-Ville, Montréal, Québec H3C 3A7, Canada; Interuniversity Research Centre on Enterprise Networks, Logistics and Transportation (CIRRELT) Pavillon André Aisenstadt, Room 3520 2920 Chemin de la Tour Université de Montréal, Montréal, Quebec H3T 1J4, Canada.
| | - Luis Miranda-Moreno
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 268, Macdonald Engineering Building, 817 Sherbrooke Street West, Montréal, Québec H3A 0C3, Canada; Interuniversity Research Centre on Enterprise Networks, Logistics and Transportation (CIRRELT) Pavillon André Aisenstadt, Room 3520 2920 Chemin de la Tour Université de Montréal, Montréal, Quebec H3T 1J4, Canada.
| |
Collapse
|
8
|
Using spatial and temporal modeling to visualize the effects of U.S. state issued stay at home orders on COVID-19. Sci Rep 2021; 11:13939. [PMID: 34230582 PMCID: PMC8260658 DOI: 10.1038/s41598-021-93433-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Accepted: 06/23/2021] [Indexed: 01/03/2023] Open
Abstract
Coronavirus disease 2019 dominated and augmented many aspects of life beginning in early 2020. Related research and data generation developed alongside its spread. We developed a Bayesian spatio-temporal Poisson disease mapping model for estimating real-time characteristics of the coronavirus disease in the United States. We also created several dashboards for visualization of the statistical model for fellow researchers and simpler spatial and temporal representations of the disease for consumption by analysts and data scientists in the policymaking community in our region. Findings suggest that the risk of confirmed cases is higher for health regions under partial stay at home orders and lower in health regions under full stay at home orders, when compared to before stay at home orders were declared. These results confirm the benefit of state-issued stay at home orders as well as suggest compliance to the directives towards the older population for adhering to social distancing guidelines.
Collapse
|
9
|
Addressing challenges in routine health data reporting in Burkina Faso through Bayesian spatiotemporal prediction of weekly clinical malaria incidence. Sci Rep 2020; 10:16568. [PMID: 33024162 PMCID: PMC7538437 DOI: 10.1038/s41598-020-73601-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Accepted: 09/07/2020] [Indexed: 11/15/2022] Open
Abstract
Sub-Saharan African (SSA) countries’ health systems are often vulnerable to unplanned situations that can hinder their effectiveness in terms of data completeness and disease control. For instance, in Burkina Faso following a workers' strike, comprehensive data on several diseases were unavailable for a long period in 2019. Weather, seasonal-malaria-chemoprevention (SMC), free healthcare, and other contextual data, which are purported to influence malarial disease, provide opportunities to fit models to describe the clinical malaria data and predict the disease spread. Bayesian spatiotemporal modeling was applied to weekly malaria surveillance data from Burkina Faso (2011–2018) while considering the effects of weather, health programs and contextual factors. Then, a prediction was used to deal with weekly missing data for the entire year of 2019, and SMC and free healthcare effects were quantified. Our proposed model accurately predicted weekly clinical malaria incidence (correlation coefficient, r = 0.90). The distribution of clinical malaria incidence was heterogeneous across the country. Overall, national predicted clinical malaria incidence in 2019 (605 per 1000 [95% CrI: 360–990]) increased by 24.7% compared with the year 2015. SMC and the interaction between free healthcare and health facility attendance were associated with a reduction in clinical malaria incidence. Our modeling approach could be a useful tool for strengthening health systems’ resilience by addressing data completeness and could support SSA countries in developing appropriate targets and indicators to facilitate the subnational control effort.
Collapse
|
10
|
Kanankege KST, Alvarez J, Zhang L, Perez AM. An Introductory Framework for Choosing Spatiotemporal Analytical Tools in Population-Level Eco-Epidemiological Research. Front Vet Sci 2020; 7:339. [PMID: 32733923 PMCID: PMC7358365 DOI: 10.3389/fvets.2020.00339] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2020] [Accepted: 05/15/2020] [Indexed: 12/04/2022] Open
Abstract
Spatiotemporal visualization and analytical tools (SATs) are increasingly being applied to risk-based surveillance/monitoring of adverse health events affecting humans, animals, and ecosystems. Different disciplines use diverse SATs to address similar research questions. The juxtaposition of these diverse techniques provides a list of options for researchers who are new to population-level spatial eco-epidemiology. Here, we are conducting a narrative review to provide an overview of the multiple available SATs, and introducing a framework for choosing among them when addressing common research questions across disciplines. The framework is comprised of three stages: (a) pre-hypothesis testing stage, in which hypotheses regarding the spatial dependence of events are generated; (b) primary hypothesis testing stage, in which the existence of spatial dependence and patterns are tested; and (c) secondary-hypothesis testing and spatial modeling stage, in which predictions and inferences were made based on the identified spatial dependences and associated covariates. In this step-wise process, six key research questions are formulated, and the answers to those questions should lead researchers to select one or more methods from four broad categories of SATs: (T1) visualization and descriptive analysis; (T2) spatial/spatiotemporal dependence and pattern recognition; (T3) spatial smoothing and interpolation; and (T4) geographic correlation studies (i.e., spatial modeling and regression). The SATs described here include both those used for decades and also other relatively new tools. Through this framework review, we intend to facilitate the choice among available SATs and promote their interdisciplinary use to support improving human, animal, and ecosystem health.
Collapse
Affiliation(s)
- Kaushi S. T. Kanankege
- Department of Veterinary Population Medicine, College of Veterinary Medicine, University of Minnesota, Saint Paul, MN, United States
| | - Julio Alvarez
- Departamento de Sanidad Animal, Centro de Vigilancia Sanitaria Veterinaria (VISAVET), Facultad de Veterinaria, Universidad Complutense, Madrid, Spain
| | - Lin Zhang
- Division of Biostatistics, School of Public Health, University of Minnesota, Minneapolis, MN, United States
| | - Andres M. Perez
- Department of Veterinary Population Medicine, College of Veterinary Medicine, University of Minnesota, Saint Paul, MN, United States
| |
Collapse
|
11
|
Otiende VA, Achia TN, Mwambi HG. Bayesian hierarchical modeling of joint spatiotemporal risk patterns for Human Immunodeficiency Virus (HIV) and Tuberculosis (TB) in Kenya. PLoS One 2020; 15:e0234456. [PMID: 32614847 PMCID: PMC7332062 DOI: 10.1371/journal.pone.0234456] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 05/27/2020] [Indexed: 11/25/2022] Open
Abstract
The simultaneous spatiotemporal modeling of multiple related diseases strengthens inferences by borrowing information between related diseases. Numerous research contributions to spatiotemporal modeling approaches exhibit their strengths differently with increasing complexity. However, contributions that combine spatiotemporal approaches to modeling of multiple diseases simultaneously are not so common. We present a full Bayesian hierarchical spatio-temporal approach to the joint modeling of Human Immunodeficiency Virus and Tuberculosis incidences in Kenya. Using case notification data for the period 2012–2017, we estimated the model parameters and determined the joint spatial patterns and temporal variations. Our model included specific and shared spatial and temporal effects. The specific random effects allowed for departures from the shared patterns for the different diseases. The space-time interaction term characterized the underlying spatial patterns with every temporal fluctuation. We assumed the shared random effects to be the structured effects and the disease-specific random effects to be unstructured effects. We detected the spatial similarity in the distribution of Tuberculosis and Human Immunodeficiency Virus in approximately 29 counties around the western, central and southern regions of Kenya. The distribution of the shared relative risks had minimal difference with the Human Immunodeficiency Virus disease-specific relative risk whereas that of Tuberculosis presented many more counties as high-risk areas. The flexibility and informative outputs of Bayesian Hierarchical Models enabled us to identify the similarities and differences in the distribution of the relative risks associated with each disease. Estimating the Human Immunodeficiency Virus and Tuberculosis shared relative risks provide additional insights towards collaborative monitoring of the diseases and control efforts.
Collapse
Affiliation(s)
- Verrah A. Otiende
- Department of Mathematical Sciences, Pan African University Institute of Basic Sciences Technology and Innovation, Nairobi, Kenya
- * E-mail: ,
| | - Thomas N. Achia
- School of Mathematics, Statistics & Computer Science, University of KwaZulu Natal, Pietermaritzburg, South Africa
| | - Henry G. Mwambi
- School of Mathematics, Statistics & Computer Science, University of KwaZulu Natal, Pietermaritzburg, South Africa
| |
Collapse
|
12
|
Rouamba T, Samadoulougou S, Tinto H, Alegana VA, Kirakoya-Samadoulougou F. Severe-malaria infection and its outcomes among pregnant women in Burkina Faso health-districts: Hierarchical Bayesian space-time models applied to routinely-collected data from 2013 to 2018. Spat Spatiotemporal Epidemiol 2020; 33:100333. [PMID: 32370941 PMCID: PMC7613547 DOI: 10.1016/j.sste.2020.100333] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Revised: 09/15/2019] [Accepted: 12/27/2019] [Indexed: 11/12/2022]
Abstract
Fine-scale hotspots detection is crucial for optimum delivery of essential health-services for reducing severe malaria in pregnancy (MiP) and death cases in Burkina Faso. This study used hierarchical-Bayesian Spatio-temporal modeling to explore space-time patterns and pinpoint health-districts with an exceedance probability of severe MiP incidence and fatality rate. Study also assessed effect of health-district service delivery (readiness) on severe-MiP outcomes. Severe-MiP fatality rate declined considerably while its incidence rate remained unchanged between January-2013 and December-2018. Severe-MiP cases persisted throughout the year with peaks between August and November. These peaks increased 2.5-fold the fatality rate. Furthermore, severe-MiP fatality was higher in health-districts classified as low-readiness (IRR = 2.469, 95%CrI: 1.632–3.738). However, the fatality rate decreased significantly with proper coverage with three doses for intermittent-preventive-treatment with sulphadoxine-pyrimethamine. Severe-MiP burden was heterogeneous spatially and temporally. The study suggested that health-programs should increase health-districts readiness and optimize resource allocation in high burden areas and months.
Collapse
Affiliation(s)
- Toussaint Rouamba
- Clinical Research Unit of Nanoro, Institut de Recherche en Sciences de la Santé, 42, Avenue Kumda-Yonre, Centre National de la Recherche Scientifique et Technologique, 11 BP 218 Ouaga CMS 11, Ouagadougou, Burkina Faso; Centre d'Epidémiologie, Biostatistique et Recherche Clinique, Ecole de Santé Publique, Université Libre de Bruxelles (ULB), Route de Lennik, 808 B-1070, Bruxelles, Belgique.
| | - Sekou Samadoulougou
- Evaluation Platform on Obesity Prevention, Quebec Heart and Lung Institute Research Center, Quebec City, QC G1V 4G5, Canada; Centre for Research on Planning and Development (CRAD), Laval University, Quebec, G1V 0A6, Canada.
| | - Halidou Tinto
- Clinical Research Unit of Nanoro, Institut de Recherche en Sciences de la Santé, 42, Avenue Kumda-Yonre, Centre National de la Recherche Scientifique et Technologique, 11 BP 218 Ouaga CMS 11, Ouagadougou, Burkina Faso
| | - Victor A Alegana
- Population Health Theme, Kenya Medical Research Institute-Wellcome Trust Research Programme, Nairobi, Kenya; Geography and Environmental Science, University of Southampton, Southampton SO17 1BJ, UK
| | - Fati Kirakoya-Samadoulougou
- Centre d'Epidémiologie, Biostatistique et Recherche Clinique, Ecole de Santé Publique, Université Libre de Bruxelles (ULB), Route de Lennik, 808 B-1070, Bruxelles, Belgique.
| |
Collapse
|
13
|
Stipancic J, Miranda-Moreno L, Strauss J, Labbe A. Pedestrian safety at signalized intersections: Modelling spatial effects of exposure, geometry and signalization on a large urban network. ACCIDENT; ANALYSIS AND PREVENTION 2020; 134:105265. [PMID: 31704639 DOI: 10.1016/j.aap.2019.105265] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/31/2015] [Revised: 08/09/2019] [Accepted: 08/09/2019] [Indexed: 06/10/2023]
Abstract
Intersections represent the most dangerous sites in the road network for pedestrians: not only is modal separation often impossible, but elements of geometry, traffic control, and built environment further exacerbate crash risk. Evaluating the safety impact of intersection features requires methods to quantify relationships between different factors and pedestrian injuries. The purpose of this paper is to model the effects of exposure, geometry, and signalization on pedestrian injuries at urban signalized intersections using a Full Bayes spatial Poisson Log-Normal model that accounts for unobserved heterogeneity and spatial correlation. Using the Integrated Nested Laplace Approximation (INLA) technique, this work leverages a rich database of geometric and signalization variables for 1864 intersections in Montreal, Quebec. To collect exposure data, short-term pedestrian and vehicle counts were extrapolated to AADT using developed expansion factors. Results of the model confirmed the positive relationship between pedestrian and vehicle volumes and pedestrian injuries. Curb extensions, raised medians, and exclusive left turn lanes were all found to reduce pedestrian injuries, while the total number of lanes and the number of commercial entrances were found to increase them. Pedestrian priority phases reduced injuries while the green straight arrow increased injuries. Lastly, the posterior expected number of crashes was used to identify hotspots. The proposed ranking criteria identified many intersections close to the city centre where the expected number of crashes is highest and intersections along arterials with lower pedestrian volumes where individual pedestrian risk is elevated. Understanding the effects of intersection geometry and pedestrian signalization will aid in ensuring the safety of pedestrians at signalized intersections.
Collapse
Affiliation(s)
- Joshua Stipancic
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec, H3T 2A7, Canada.
| | - Luis Miranda-Moreno
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 268, Macdonald, Engineering Building, 817 Sherbrooke Street West, Montréal, Québec, H3A 0C3, Canada.
| | - Jillian Strauss
- Polytechnique Montréal, C.P. 6079, succ. Centre-Ville, Montréal, Québec, H3C 3A7, Canada.
| | - Aurélie Labbe
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec, H3T 2A7, Canada.
| |
Collapse
|
14
|
Zhang Y, Wang X, Li Y, Ma J. Spatiotemporal Analysis of Influenza in China, 2005-2018. Sci Rep 2019; 9:19650. [PMID: 31873144 PMCID: PMC6928232 DOI: 10.1038/s41598-019-56104-8] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 12/04/2019] [Indexed: 12/14/2022] Open
Abstract
Influenza is a major cause of morbidity and mortality worldwide, as well as in China. Knowledge of the spatial and temporal characteristics of influenza is important in evaluating and developing disease control programs. This study aims to describe an accurate spatiotemporal pattern of influenza at the prefecture level and explore the risk factors associated with influenza incidence risk in mainland China from 2005 to 2018. The incidence data of influenza were obtained from the Chinese Notifiable Infectious Disease Reporting System (CNIDRS). The Besag York Mollié (BYM) model was extended to include temporal and space-time interaction terms. The parameters for this extended Bayesian spatiotemporal model were estimated through integrated nested Laplace approximations (INLA) using the package R-INLA in R. A total of 702,226 influenza cases were reported in mainland China in CNIDRS from 2005–2018. The yearly reported incidence rate of influenza increased 15.6 times over the study period, from 3.51 in 2005 to 55.09 in 2008 per 100,000 populations. The temporal term in the spatiotemporal model showed that much of the increase occurred during the last 3 years of the study period. The risk factor analysis showed that the decreased number of influenza vaccines for sale, the new update of the influenza surveillance protocol, the increase in the rate of influenza A (H1N1)pdm09 among all processed specimens from influenza-like illness (ILI) patients, and the increase in the latitude and longitude of geographic location were associated with an increase in the influenza incidence risk. After the adjusting for fixed covariate effects and time random effects, the map of the spatial structured term shows that high-risk areas clustered in the central part of China and the lowest-risk areas in the east and west. Large space-time variations in influenza have been found since 2009. In conclusion, an increasing trend of influenza was observed from 2005 to 2018. The insufficient flu vaccine supplements, the newly emerging influenza A (H1N1)pdm09 and expansion of influenza surveillance efforts might be the major causes of the dramatic changes in outbreak and spatio-temporal epidemic patterns. Clusters of prefectures with high relative risks of influenza were identified in the central part of China. Future research with more risk factors at both national and local levels is necessary to explain the changing spatiotemporal patterns of influenza in China.
Collapse
Affiliation(s)
- Yewu Zhang
- Center for Public Health Surveillance and Information Service, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Xiaofeng Wang
- Center for Public Health Surveillance and Information Service, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Yanfei Li
- Center for Public Health Surveillance and Information Service, Chinese Center for Disease Control and Prevention, Beijing, China
| | - Jiaqi Ma
- Center for Public Health Surveillance and Information Service, Chinese Center for Disease Control and Prevention, Beijing, China.
| |
Collapse
|
15
|
Otiende V, Achia T, Mwambi H. Bayesian modeling of spatiotemporal patterns of TB-HIV co-infection risk in Kenya. BMC Infect Dis 2019; 19:902. [PMID: 31660883 PMCID: PMC6819548 DOI: 10.1186/s12879-019-4540-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 10/09/2019] [Indexed: 02/01/2023] Open
Abstract
Background Tuberculosis (TB) and Human Immunodeficiency Virus (HIV) diseases are globally acknowledged as a public health challenge that exhibits adverse bidirectional relations due to the co-epidemic overlap. To understand the co-infection burden we used the case notification data to generate spatiotemporal maps that described the distribution and exposure hypotheses for further epidemiologic investigations in areas with unusual case notification levels. Methods We analyzed the TB and TB-HIV case notification data from the Kenya national TB control program aggregated for forty-seven counties over a seven-year period (2012–2018). Using spatiotemporal poisson regression models within the Integrated Nested Laplace Approach (INLA) paradygm, we modeled the risk of TB-HIV co-infection. Six competing models with varying space-time formulations were compared to determine the best fit model. We then assessed the geographic patterns and temporal trends of coinfection risk by mapping the posterior marginal from the best fit model. Results Of the total 608,312 TB case notifications, 194,129 were HIV co-infected. The proportion of TB-HIV co-infection was higher in females (39.7%) than in males (27.0%). A significant share of the co-infection was among adults aged 35 to 44 years (46.7%) and 45 to 54 years (42.1%). Based on the Bayesian Defiance Information (DIC) and the effective number of parameters (pD) comparisons, the spatiotemporal model allowing space-time interaction was the best in explaining the geographical variations in TB-HIV coinfection. The model results suggested that the risk of TB-HIV coinfection was influenced by infrastructure index (Relative risk (RR) = 5.75, Credible Interval (Cr.I) = (1.65, 19.89)) and gender ratio (RR = 5.81e−04, Cr. I = (1.06e−04, 3.18e−03). The lowest and highest temporal relative risks were in the years 2016 at 0.9 and 2012 at 1.07 respectively. The spatial pattern presented an increased co-infection risk in a number of counties. For the spatiotemporal interaction, only a few counties had a relative risk greater than 1 that varied in different years. Conclusions We identified elevated risk areas for TB/HIV co-infection and fluctuating temporal trends which could be because of improved TB case detection or surveillance bias caused by spatial heterogeneity in the co-infection dynamics. Focused interventions and continuous TB-HIV surveillance will ensure adequate resource allocation and significant reduction of HIV burden amongst TB patients.
Collapse
Affiliation(s)
- Verrah Otiende
- Department of Mathematical Sciences, Pan African University Institute of Basic Sciences Technology and Innovation, Nairobi, Kenya.
| | - Thomas Achia
- School of Mathematics, Statistics & Computer Science, University of KwaZulu-Natal, Pietermaritzburg, South Africa
| | - Henry Mwambi
- School of Mathematics, Statistics & Computer Science, University of KwaZulu-Natal, Pietermaritzburg, South Africa
| |
Collapse
|
16
|
Stipancic J, Miranda-Moreno L, Saunier N, Labbe A. Network screening for large urban road networks: Using GPS data and surrogate measures to model crash frequency and severity. ACCIDENT; ANALYSIS AND PREVENTION 2019; 125:290-301. [PMID: 30818096 DOI: 10.1016/j.aap.2019.02.016] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Revised: 02/07/2019] [Accepted: 02/13/2019] [Indexed: 06/09/2023]
Abstract
Crash frequency and injury severity are independent dimensions defining crash risk in road safety management and network screening. Traditional screening techniques model crashes using regression and historical crash data, making them intrinsically reactive. In response, surrogate measures of safety have become a popular proactive alternative. The purpose of this paper is to develop models for crash frequency and severity incorporating GPS-derived surrogate safety measures (SSMs) as predictive variables. SSMs based on vehicle manoeuvres and traffic flow were extracted from data collected in Quebec City. The mixed multivariate outcome is estimated using two models; a Full Bayes Spatial Negative Binomial model for crash frequency estimated using the Integrated Nested Laplace Approximation approach and a fractional Multinomial Logit model for crash severity. Model outcomes are combined to generate posterior expected crash frequency at each severity level and rank sites based on crash cost. The crash frequency model was accurate at the network scale, with the majority of proposed SSMs statistically significant at 95% confidence and the direction of their effect generally consistent with previous research. In the crash severity model, fewer variables were significant, yet the direction of the effect of all significant variables was again consistent with previous research. Correlations between rankings predicted by the mixed multivariate model and by the crash data were adequate for intersections (0.46) but were poorer for links (0.25). The ability to prioritize sites based on GPS data and SSMs rather than historical crash data represents a substantial contribution to the field of road safety.
Collapse
Affiliation(s)
- Joshua Stipancic
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 391, Macdonald Engineering Building, 817 Sherbrooke Street West, Montréal, Québec, H3A 0C3, Canada.
| | - Luis Miranda-Moreno
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 268, Macdonald Engineering Building, 817 Sherbrooke Street West, Montréal, Québec, H3A 0C3, Canada.
| | - Nicolas Saunier
- Department of Civil, Geological and Mining Engineering, Polytechnique Montréal, C.P. 6079, succ. Centre-Ville, Montréal, Québec, H3C 3A7, Canada.
| | - Aurélie Labbe
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec, H3T 2A7, Canada.
| |
Collapse
|
17
|
Corpas-Burgos F, Botella-Rocamora P, Martinez-Beneito MA. On the convenience of heteroscedasticity in highly multivariate disease mapping. TEST-SPAIN 2019. [DOI: 10.1007/s11749-019-00628-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
18
|
Stipancic J, Miranda-Moreno L, Saunier N, Labbe A. Surrogate safety and network screening: Modelling crash frequency using GPS travel data and latent Gaussian Spatial Models. ACCIDENT; ANALYSIS AND PREVENTION 2018; 120:174-187. [PMID: 30142497 DOI: 10.1016/j.aap.2018.07.013] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Revised: 07/03/2018] [Accepted: 07/05/2018] [Indexed: 06/08/2023]
Abstract
Improving road safety requires accurate network screening methods to identify and prioritize sites in order to maximize the effectiveness of implemented countermeasures. In screening, hotspots are commonly identified using statistical models and ranking criteria derived from observed crash data. However, collision databases are subject to errors, omissions, and underreporting. More importantly, crash-based methods are reactive and require years of crash data. With the arrival of new technologies including Global Positioning System (GPS) trajectory data, proactive surrogate safety methods have gained popularity as an alternative approach for screening. GPS-enabled smartphones can collect reliable and spatio-temporally rich driving data from regular drivers using an inexpensive, simple, and user-friendly tool. However, few studies to date have analyzed large volumes of smartphone GPS data and considered surrogate-safety modelling techniques for network screening. The purpose of this paper is to propose a surrogate safety screening approach based on smartphone GPS data and a Full Bayesian modelling framework. After processing crash data and GPS data collected in Quebec City, Canada, several surrogate safety measures (SSMs), including vehicle manoeuvres (hard braking) and measures of traffic flow (congestion, average speed, and speed variation), were extracted. Then, spatial crash frequency models incorporating the extracted SSMs were proposed and validated. A Latent Gaussian Spatial Model was estimated using the Integrated Nested Laplace Approximation (INLA) technique. While the INLA Negative Binomial models outperformed alternative models, incorporating spatial correlations provided the greatest improvement in model fit. Relationships between SSMs and crash frequency established in previous studies were generally supported by the modelling results. For example, hard braking, congestion, and speed variation were all positively linked to crash counts at the intersection level. Network screening based on SSMs presents a substantial contribution to the field of road safety and works towards the elimination of crash data in evaluation and monitoring.
Collapse
Affiliation(s)
- Joshua Stipancic
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 391, Macdonald Engineering Building, 817 Sherbrooke Street West, Montréal, Québec, H3A 0C3, Canada.
| | - Luis Miranda-Moreno
- Department of Civil Engineering and Applied Mechanics, McGill University, Room 268, Macdonald Engineering Building, 817 Sherbrooke Street West, Montréal, Québec, H3A 0C3, Canada.
| | - Nicolas Saunier
- Department of Civil, Geological and Mining Engineering, Polytechnique Montréal, C.P. 6079, Succ. Centre-Ville, Montréal, Québec, H3C 3A7, Canada.
| | - Aurelie Labbe
- Department of Decision Sciences, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montréal, Québec, H3T 2A7, Canada.
| |
Collapse
|
19
|
Bakka H, Rue H, Fuglstad G, Riebler A, Bolin D, Illian J, Krainski E, Simpson D, Lindgren F. Spatial modeling with R‐INLA: A review. ACTA ACUST UNITED AC 2018. [DOI: 10.1002/wics.1443] [Citation(s) in RCA: 121] [Impact Index Per Article: 20.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Haakon Bakka
- CEMSE Division King Abdullah University of Science and Technology Thuwal Saudi Arabia
| | - Håvard Rue
- CEMSE Division King Abdullah University of Science and Technology Thuwal Saudi Arabia
| | - Geir‐Arne Fuglstad
- Department of Mathematical Sciences Norwegian University of Science and Technology Trondheim Norway
| | - Andrea Riebler
- Department of Mathematical Sciences Norwegian University of Science and Technology Trondheim Norway
| | - David Bolin
- Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Gothenburg Sweden
| | - Janine Illian
- CREEM, School of Mathematics and Statistics University of St Andrews St. Andrews UK
| | - Elias Krainski
- Departamento de Estatística Universidade Federal do Paraná Paraná Brazil
| | - Daniel Simpson
- Department of Statistical Sciences University of Toronto Toronto Canada
| | - Finn Lindgren
- School of Mathematics University of Edinburgh Edinburgh UK
| |
Collapse
|
20
|
Carroll R, Zhao S. Gaining relevance from the random: Interpreting observed spatial heterogeneity. Spat Spatiotemporal Epidemiol 2018; 25:11-17. [PMID: 29751888 DOI: 10.1016/j.sste.2018.01.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Revised: 12/26/2017] [Accepted: 01/10/2018] [Indexed: 10/18/2022]
Abstract
In Bayesian disease mapping, spatial random effects are used to account for confounding in the data so that reasonable estimates for the fixed effects can be obtained. Typically, the spatial random effects are mapped and qualitative comments are made related to an increase or decrease in risk for certain areas. The approach outlined here illustrates how a quantitative secondary assessment can be applied to make more useful and applicable inference related to these spatial random effects. We are able to recover important but unmeasured or unincluded risk factors via a secondary model fit. Results from the secondary model fit can determine association between spatial region-level risk factors and the estimated spatial random effects. We believe this work presents a useful, quantitative technique highlighting the importance and applicability of spatial random effects as well as illustrates how these methods lead to more interpretable conclusions.
Collapse
Affiliation(s)
- Rachel Carroll
- Biostatistics and Computational Biology Branch, National Institute of Environmental Health Sciences, 111 TW Alexander Dr., Research Triangle Park, NC 27709 USA.
| | - Shanshan Zhao
- Biostatistics and Computational Biology Branch, National Institute of Environmental Health Sciences, 111 TW Alexander Dr., Research Triangle Park, NC 27709 USA
| |
Collapse
|
21
|
Riebler A, Sørbye SH, Simpson D, Rue H. An intuitive Bayesian spatial model for disease mapping that accounts for scaling. Stat Methods Med Res 2018; 25:1145-65. [PMID: 27566770 DOI: 10.1177/0962280216660421] [Citation(s) in RCA: 174] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In recent years, disease mapping studies have become a routine application within geographical epidemiology and are typically analysed within a Bayesian hierarchical model formulation. A variety of model formulations for the latent level have been proposed but all come with inherent issues. In the classical BYM (Besag, York and Mollié) model, the spatially structured component cannot be seen independently from the unstructured component. This makes prior definitions for the hyperparameters of the two random effects challenging. There are alternative model formulations that address this confounding; however, the issue on how to choose interpretable hyperpriors is still unsolved. Here, we discuss a recently proposed parameterisation of the BYM model that leads to improved parameter control as the hyperparameters can be seen independently from each other. Furthermore, the need for a scaled spatial component is addressed, which facilitates assignment of interpretable hyperpriors and make these transferable between spatial applications with different graph structures. The hyperparameters themselves are used to define flexible extensions of simple base models. Consequently, penalised complexity priors for these parameters can be derived based on the information-theoretic distance from the flexible model to the base model, giving priors with clear interpretation. We provide implementation details for the new model formulation which preserve sparsity properties, and we investigate systematically the model performance and compare it to existing parameterisations. Through a simulation study, we show that the new model performs well, both showing good learning abilities and good shrinkage behaviour. In terms of model choice criteria, the proposed model performs at least equally well as existing parameterisations, but only the new formulation offers parameters that are interpretable and hyperpriors that have a clear meaning.
Collapse
Affiliation(s)
- Andrea Riebler
- Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim, Norway
| | - Sigrunn H Sørbye
- Department of Mathematics and Statistics, UiT The Arctic University of Norway, Tromsø, Norway
| | - Daniel Simpson
- Department of Mathematical Sciences, University of Bath, Bath, UK
| | - Håvard Rue
- Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
22
|
Goicoa T, Adin A, Etxeberria J, Militino AF, Ugarte MD. Flexible Bayesian P-splines for smoothing age-specific spatio-temporal mortality patterns. Stat Methods Med Res 2017; 28:384-403. [PMID: 28847210 DOI: 10.1177/0962280217726802] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
In this paper age-space-time models based on one and two-dimensional P-splines with B-spline bases are proposed for smoothing mortality rates, where both fixed relative scale and scale invariant two-dimensional penalties are examined. Model fitting and inference are carried out using integrated nested Laplace approximations, a recent Bayesian technique that speeds up computations compared to McMC methods. The models will be illustrated with Spanish breast cancer mortality data during the period 1985-2010, where a general decline in breast cancer mortality has been observed in Spanish provinces in the last decades. The results reveal that mortality rates for the oldest age groups do not decrease in all provinces.
Collapse
Affiliation(s)
- T Goicoa
- 1 Department of Statistics and Operations Research, Public University of Navarre, Spain.,2 Institute for Advanced Materials (InaMat), Public University of Navarre, Spain.,3 Research Network on Health Services in Chronic Diseases (REDISSEC), Spain
| | - A Adin
- 1 Department of Statistics and Operations Research, Public University of Navarre, Spain.,2 Institute for Advanced Materials (InaMat), Public University of Navarre, Spain
| | - J Etxeberria
- 1 Department of Statistics and Operations Research, Public University of Navarre, Spain.,2 Institute for Advanced Materials (InaMat), Public University of Navarre, Spain.,4 Consortium for Biomedical Research in Epidemiology and Public Health (CIBERESP), Spain
| | - A F Militino
- 1 Department of Statistics and Operations Research, Public University of Navarre, Spain.,2 Institute for Advanced Materials (InaMat), Public University of Navarre, Spain
| | - M D Ugarte
- 1 Department of Statistics and Operations Research, Public University of Navarre, Spain.,2 Institute for Advanced Materials (InaMat), Public University of Navarre, Spain
| |
Collapse
|
23
|
Panczak R, Moser A, Held L, Jones PA, Rühli FJ, Staub K. A tall order: Small area mapping and modelling of adult height among Swiss male conscripts. ECONOMICS AND HUMAN BIOLOGY 2017; 26:61-69. [PMID: 28284175 DOI: 10.1016/j.ehb.2017.01.005] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2016] [Revised: 01/10/2017] [Accepted: 01/11/2017] [Indexed: 06/06/2023]
Abstract
Adult height reflects an individual's socio-economic background and offers insights into the well-being of populations. Height is linked to various health outcomes such as morbidity and mortality and has consequences on the societal level. The aim of this study was to describe small-area variation of height and associated factors among young men in Switzerland. Data from 175,916 conscripts (aged between 18.50 and 20.50 years) was collected between 2005 and 2011, which represented approximately 90% of the corresponding birth cohorts. These were analysed using Gaussian hierarchical models in a Bayesian framework to investigate the spatial pattern of mean height across postcodes. The models varied both in random effects and degree of adjustment (professional status, area-based socioeconomic position, and language region). We found a strong spatial structure for mean height across postcodes. The range of height differences between mean postcode level estimates was 3.40cm according to the best fitting model, with the shorter conscripts coming from the Italian and French speaking parts of Switzerland. There were positive socioeconomic gradients in height at both individual and area-based levels. Spatial patterns for height persisted after adjustment for individual factors, but not when language region was included. Socio-economic position and cultural/natural boundaries such as language borders and mountain passes are shaping patterns of height for Swiss conscripts. Small area mapping of height contributes to the understanding of its cofactors.
Collapse
Affiliation(s)
- Radoslaw Panczak
- Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland; Institute of Social and Preventive Medicine, University of Bern, Finkenhubelweg 11, CH-3012 Bern, Switzerland
| | - André Moser
- Institute of Social and Preventive Medicine, University of Bern, Finkenhubelweg 11, CH-3012 Bern, Switzerland; Department of Geriatrics, Inselspital, Bern University Hospital, University of Bern, CH-3012 Bern, Switzerland
| | - Leonhard Held
- Epidemiology, Biostatistics and Prevention Institute, University of Zurich, Hirschengraben 84, CH-8001 Zurich, Switzerland
| | - Philip A Jones
- Department of Geography, Swansea University, Wallace Building, Singleton Park, Swansea SA2 8PP, UK
| | - Frank J Rühli
- Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Kaspar Staub
- Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland.
| |
Collapse
|
24
|
Carroll R, Lawson AB, Faes C, Kirby RS, Aregay M, Watjou K. Spatio-temporal Bayesian model selection for disease mapping. ENVIRONMETRICS 2016; 27:466-478. [PMID: 28070156 PMCID: PMC5217709 DOI: 10.1002/env.2410] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Spatio-temporal analysis of small area health data often involves choosing a fixed set of predictors prior to the final model fit. In this paper, we propose a spatio-temporal approach of Bayesian model selection to implement model selection for certain areas of the study region as well as certain years in the study time line. Here, we examine the usefulness of this approach by way of a large-scale simulation study accompanied by a case study. Our results suggest that a special case of the model selection methods, a mixture model allowing a weight parameter to indicate if the appropriate linear predictor is spatial, spatio-temporal, or a mixture of the two, offers the best option to fitting these spatio-temporal models. In addition, the case study illustrates the effectiveness of this mixture model within the model selection setting by easily accommodating lifestyle, socio-economic, and physical environmental variables to select a predominantly spatio-temporal linear predictor.
Collapse
Affiliation(s)
- R Carroll
- Department of Public Health, Medical University of South Carolina
- Corresponding author, Dr. R Carroll, Department of Public Health, Medical University of South Carolina, 135 Cannon St, Charleston, SC 29425, USA,
| | - AB Lawson
- Department of Public Health, Medical University of South Carolina
| | - C Faes
- Interuniversity Institute for Statistics and Statistical Bioinformatics, Hasselt University
| | - RS Kirby
- Department of Community and Family Health, University of South Florida
| | - M Aregay
- Department of Public Health, Medical University of South Carolina
| | - K Watjou
- Interuniversity Institute for Statistics and Statistical Bioinformatics, Hasselt University
| |
Collapse
|
25
|
Tahden M, Manitz J, Baumgardt K, Fell G, Kneib T, Hegasy G. Epidemiological and Ecological Characterization of the EHEC O104:H4 Outbreak in Hamburg, Germany, 2011. PLoS One 2016; 11:e0164508. [PMID: 27723830 PMCID: PMC5056673 DOI: 10.1371/journal.pone.0164508] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2016] [Accepted: 09/25/2016] [Indexed: 11/19/2022] Open
Abstract
In 2011, a large outbreak of entero-hemorrhagic E. coli (EHEC) and hemolytic uremic syndrome (HUS) occurred in Germany. The City of Hamburg was the first focus of the epidemic and had the highest incidences among all 16 Federal States of Germany. In this article, we present epidemiological characteristics of the Hamburg notification data. Evaluating the epicurves retrospectively, we found that the first epidemiological signal of the outbreak, which was in form of a HUS case cluster, was received by local health authorities when already 99 EHEC and 48 HUS patients had experienced their first symptoms. However, only two EHEC and seven HUS patients had been notified. Middle-aged women had the highest risk for contracting the infection in Hamburg. Furthermore, we studied timeliness of case notification in the course of the outbreak. To analyze the spatial distribution of EHEC/HUS incidences in 100 districts of Hamburg, we mapped cases' residential addresses using geographic information software. We then conducted an ecological study in order to find a statistical model identifying associations between local socio-economic factors and EHEC/HUS incidences in the epidemic. We employed a Bayesian Poisson model with covariates characterizing the Hamburg districts as well as incorporating structured and unstructured spatial effects. The Deviance Information Criterion was used for stepwise variable selection. We applied different modeling approaches by using primary data, transformed data, and preselected subsets of transformed data in order to identify socio-economic factors characterizing districts where EHEC/HUS outbreak cases had their residence.
Collapse
Affiliation(s)
- Maike Tahden
- Biometry and Data Management, Leibniz Institute for Prevention Research and Epidemiology—BIPS, Bremen, Germany
- Department of Psychology and Cluster of Excellence “Hearing4all”, Carl von Ossietzky University Oldenburg, Oldenburg, Germany
| | - Juliane Manitz
- Department for Statistics and Econometrics, University of Goettingen, Goettingen, Germany
| | - Klaus Baumgardt
- Division for Environmental Monitoring, Institut fuer Hygiene und Umwelt, Hamburg, Germany
| | - Gerhard Fell
- Centre for Infectious Diseases Epidemiology, Institut fuer Hygiene und Umwelt, Hamburg, Germany
| | - Thomas Kneib
- Department for Statistics and Econometrics, University of Goettingen, Goettingen, Germany
| | - Guido Hegasy
- Centre for Infectious Diseases Epidemiology, Institut fuer Hygiene und Umwelt, Hamburg, Germany
- * E-mail:
| |
Collapse
|
26
|
Incorporating High-Dimensional Exposure Modelling into Studies of Air Pollution and Health. STATISTICS IN BIOSCIENCES 2016; 9:559-581. [PMID: 29225714 PMCID: PMC5711999 DOI: 10.1007/s12561-016-9150-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2015] [Accepted: 05/20/2016] [Indexed: 11/11/2022]
Abstract
Performing studies on the risks of environmental hazards on human health requires accurate estimates of exposures that might be experienced by the populations at risk. Often there will be missing data and in many epidemiological studies, the locations and times of exposure measurements and health data do not match. To a large extent this will be due to the health and exposure data having arisen from completely different data sources and not as the result of a carefully designed study, leading to problems of both ‘change of support’ and ‘misaligned data’. In such cases, a direct comparison of the exposure and health outcome is often not possible without an underlying model to align the two in the spatial and temporal domains. The Bayesian approach provides the natural framework for such models; however, the large amounts of data that can arise from environmental networks means that inference using Markov Chain Monte Carlo might not be computationally feasible in this setting. Here we adapt the integrated nested Laplace approximation to implement spatio–temporal exposure models. We also propose methods for the integration of large-scale exposure models and health analyses. It is important that any model structure allows the correct propagation of uncertainty from the predictions of the exposure model through to the estimates of risk and associated confidence intervals. The methods are demonstrated using a case study of the levels of black smoke in the UK, measured over several decades, and respiratory mortality.
Collapse
|
27
|
Hanandita W, Tampubolon G. Geography and social distribution of malaria in Indonesian Papua: a cross-sectional study. Int J Health Geogr 2016; 15:13. [PMID: 27072128 PMCID: PMC4830039 DOI: 10.1186/s12942-016-0043-y] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2015] [Accepted: 03/30/2016] [Indexed: 11/10/2022] Open
Abstract
Background Despite being one of the world’s most affected regions, only little is known about the social and spatial distributions of malaria in Indonesian Papua. Existing studies tend to be descriptive in nature; their inferences are prone to confounding and selection biases. At the same time, there remains limited malaria-cartographic activity in the region. Analysing a subset (N = 22,643) of the National Basic Health Research 2007 dataset (N = 987,205), this paper aims to quantify the district-specific risk of malaria in Papua and to understand how socio-demographic/economic factors measured at individual and district levels are associated with individual’s probability of contracting the disease. Methods We adopt a Bayesian hierarchical logistic regression model that accommodates not only the nesting of individuals within the island’s 27 administrative units but also the spatial autocorrelation among these locations. Both individual and contextual characteristics are included as predictors in the model; a normal conditional autoregressive prior and an exchangeable one are assigned to the random effects. Robustness is then assessed through sensitivity analyses using alternative hyperpriors. Results We find that rural Papuans as well as those who live in poor, densely forested, lowland districts are at a higher risk of infection than their counterparts. We also find age and gender differentials in malaria prevalence, if only to a small degree. Nine districts are estimated to have higher-than-expected malaria risks; the extent of spatial variation on the island remains notable even after accounting for socio-demographic/economic risk factors. Conclusions Although we show that malaria is geography-dependent in Indonesian Papua, it is also a disease of poverty. This means that malaria eradication requires not only biological (proximal) interventions but also social (distal) ones.
Collapse
Affiliation(s)
- Wulung Hanandita
- Cathie Marsh Institute for Social Research (CMIST), University Manchester, Oxford Road, Manchester, M13 9PL, UK.
| | - Gindo Tampubolon
- Cathie Marsh Institute for Social Research (CMIST), University Manchester, Oxford Road, Manchester, M13 9PL, UK
| |
Collapse
|
28
|
Panczak R, Held L, Moser A, Jones PA, Rühli FJ, Staub K. Finding big shots: small-area mapping and spatial modelling of obesity among Swiss male conscripts. BMC OBESITY 2016; 3:10. [PMID: 26918194 PMCID: PMC4758017 DOI: 10.1186/s40608-016-0092-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/18/2015] [Accepted: 02/10/2016] [Indexed: 12/03/2022]
Abstract
BACKGROUND In Switzerland, as in other developed countries, the prevalence of overweight and obesity has increased substantially since the early 1990s. Most of the analyses so far have been based on sporadic surveys or self-reported data and did not offer potential for small-area analyses. The goal of this study was to investigate spatial variation and determinants of obesity among young Swiss men using recent conscription data. METHODS A complete, anonymized dataset of conscription records for the 2010-2012 period were provided by Swiss Armed Forces. We used a series of Bayesian hierarchical logistic regression models to investigate the spatial pattern of obesity across 3,187 postcodes, varying them by type of random effects (spatially unstructured and structured), level of adjustment by individual (age and professional status) and area-based [urbanicity and index of socio-economic position (SEP)] characteristics. RESULTS The analysed dataset consisted of 100,919 conscripts, out of which 5,892 (5.8 %) were obese. Crude obesity prevalence increased with age among conscripts of lower individual and area-based SEP and varied greatly over postcodes. Best model's estimates of adjusted odds ratios of obesity on postcode level ranged from 0.61 to 1.93 and showed a strong spatial pattern of obesity risk across the country. Odds ratios above 1 concentrated in central and north Switzerland. Smaller pockets of elevated obesity risk also emerged around cities of Geneva, Fribourg and Lausanne. Lower estimates were observed in North-East and East as well as south of the Alps. Importantly, small regional outliers were observed and patterning did not follow administrative boundaries. Similarly as with crude obesity prevalence, the best fitting model confirmed increasing risk of obesity with age and among conscripts of lower professional status. The risk decreased with higher area-based SEP and, to a lesser degree - in rural areas. CONCLUSION In Switzerland, there is a substantial spatial variation in obesity risk among young Swiss men. Small-area estimates of obesity risk derived from conscripts records contribute to its understanding and could be used to design further studies and interventions.
Collapse
Affiliation(s)
- Radoslaw Panczak
- />Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
- />Institute of Social and Preventive Medicine, University of Bern, Finkenhubelweg 11, CH-3012 Bern, Switzerland
| | - Leonhard Held
- />Epidemiology, Biostatistics and Prevention Institute, University of Zurich, Hirschengraben 84, CH-8001 Zurich, Switzerland
| | - André Moser
- />Institute of Social and Preventive Medicine, University of Bern, Finkenhubelweg 11, CH-3012 Bern, Switzerland
| | - Philip A. Jones
- />Department of Geography, Swansea University, Wallace Building, Singleton Park, Swansea, SA2 8PP UK
| | - Frank J. Rühli
- />Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Kaspar Staub
- />Institute of Evolutionary Medicine, University of Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| |
Collapse
|
29
|
A latent process model for forecasting multiple time series in environmental public health surveillance. Stat Med 2016; 35:3085-100. [DOI: 10.1002/sim.6904] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2015] [Revised: 11/26/2015] [Accepted: 01/21/2016] [Indexed: 01/19/2023]
|
30
|
Costa JV, Silveira LVDA, Donalísio MR. Análise espacial de dados de contagem com excesso de zeros aplicado ao estudo da incidência de dengue em Campinas, São Paulo, Brasil. CAD SAUDE PUBLICA 2016; 32:e00036915. [DOI: 10.1590/0102-311x00036915] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2015] [Accepted: 04/06/2016] [Indexed: 11/21/2022] Open
Abstract
Resumo: A incidência de dengue ocorre predominantemente em áreas urbanas das cidades. Identificar o padrão de distribuição espacial da doença no nível local contribui na formulação de estratégias de controle e prevenção da doença. A análise espacial de dados de contagem para pequenas áreas comumente viola as suposições dos modelos tradicionais de Poisson, devido à quantidade excessiva de zeros. Neste estudo, comparou-se o desempenho de quatro modelos de contagem utilizados no mapeamento de doenças: Poisson, Binomial negativa, Poisson inflacionado de zeros e Binomial negativa inflacionado de zeros. Os métodos foram comparados em um estudo de simulação. Os modelos analisados no estudo de simulação foram aplicados em um estudo ecológico espacial, aos dados de dengue agregados por setores censitários, do Município de Campinas, São Paulo, Brasil, em 2007. A análise espacial foi conduzida por modelos hierárquicos bayesianos. O modelo de Poisson inflacionado de zeros apresentou melhor desempenho para estimar o risco relativo de incidência de dengue nos setores censitários.
Collapse
|
31
|
Abstract
BACKGROUND This study quantifies the spatiotemporal risk of pedestrian and bicyclist injury in New York City at the census tract level over a recent 10-year period, identifies areas of increased risk, and evaluates the role of socioeconomic and traffic-related variables in injury risk. METHODS Crash data on 140,835 pedestrian and bicyclist injuries in 1908 census tracts from 2001 to 2010 were obtained from the New York City Department of Transportation. We analyzed injury counts within census tracts with Bayesian hierarchical spatial models using integrated nested Laplace approximations. The model included variables for social fragmentation, median household income, and average vehicle speed and traffic density, as well as a spatially unstructured random effect term, a spatially structured conditional autoregression term, a first-order random walk-correlated time variable, and an interaction term for time and place. Incidence density ratios, credible intervals, and probability exceedances were calculated and mapped. RESULTS The yearly rate of crashes involving injuries to "pedestrians" (including bicyclists) decreased 16.2% over the study period, from 23.7 per 10,000 population to 16.2 per 10,000. The temporal term in the spatiotemporal model indicated that much of the decrease over the study period occurred during the first 4 years of the study period. Despite an overall decrease, the model identified census tracts that were at persistently high risk of pedestrian injury throughout the study period, as well as areas that experienced sporadic annual increases in risk. Aggregate social, economic, and traffic-related measures were associated with pedestrian injury risk at the ecologic level. Every 1-unit increase in a standardized social fragmentation index was associated with a 19% increase in pedestrian injury risk (incidence density ratio = 1.19 [95% credible interval = 1.16 - 1.23]), and every 1 standardized unit increase in traffic density was associated with a 20% increase in pedestrian injury risk (1.20 [1.15 - 1.26]). Each 10-mile-per-hour increase in average traffic speed in a census tract was associated with a 24% decrease in pedestrian injury risk (0.76 [0.69 - 0.83]). CONCLUSIONS The risk of a pedestrian or bicyclist being struck by a motor vehicle in New York City decreased from 2001 to 2004 and held fairly steady thereafter. Some census tracts in the city did not benefit from overall reductions or experienced sporadic years of increased risk compared with the city as a whole. Injury risk at the census tract level was associated with social, economic, and traffic-related factors.
Collapse
|
32
|
Sparks C. An examination of disparities in cancer incidence in Texas using Bayesian random coefficient models. PeerJ 2015; 3:e1283. [PMID: 26421245 PMCID: PMC4586809 DOI: 10.7717/peerj.1283] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Accepted: 09/09/2015] [Indexed: 01/05/2023] Open
Abstract
Disparities in cancer risk exist between ethnic groups in the United States. These disparities often result from differential access to healthcare, differences in socioeconomic status and differential exposure to carcinogens. This study uses cancer incidence data from the population based Texas Cancer Registry to investigate the disparities in digestive and respiratory cancers from 2000 to 2008. A Bayesian hierarchical regression approach is used. All models are fit using the INLA method of Bayesian model estimation. Specifically, a spatially varying coefficient model of the disparity between Hispanic and Non-Hispanic incidence is used. Results suggest that a spatio-temporal heterogeneity model best accounts for the observed Hispanic disparity in cancer risk. Overall, there is a significant disadvantage for the Hispanic population of Texas with respect to both of these cancers, and this disparity varies significantly over space. The greatest disparities between Hispanics and Non-Hispanics in digestive and respiratory cancers occur in eastern Texas, with patterns emerging as early as 2000 and continuing until 2008.
Collapse
Affiliation(s)
- Corey Sparks
- Department of Demography, The University of Texas at San Antonio , San Antonio, TX , USA
| |
Collapse
|
33
|
Rastaghi S, Jafari-Koshki T, Mahaki B. Application of Bayesian Multilevel Space-Time Models to Study Relative Risk of Esophageal Cancer in Iran 2005-2007 at a County Level. Asian Pac J Cancer Prev 2015; 16:5787-92. [PMID: 26320452 DOI: 10.7314/apjcp.2015.16.14.5787] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Reported age standardized incidence rates for esophageal cancer in Iran are 0.88 and 6.15 for females and males, at fifth and the eighth ranks, respectively, of cancers overall. The present study aimed to map relative risk using more realistic and less problematic methods than common estimators. MATERIALS AND METHODS In this ecological investigation, the studied population consisted of all esophageal cancer patients in Iran from 2005 to 2007. The Bayesian multilevel space-time model with three levels of county, province, and time was used to measure the relative risk of esophageal cancer. Analyses were conducted using R package INLA. RESULTS The total number of registered patients was 7,160. According to the results, the three-level model with adjustment for risk factors of physical activity and smoking had the best fit among all models .The overall temporal trend was significantly increasing. At county level, Ahar, Marand, Salmas, Bojnoord, Saghez, Sarakhs, Shahroud and Torbatejam had the highest relative risks. Physical activity was found to have significant direct association with risk of developing esophageal cancer. CONCLUSIONS Given to great variation across geographical areas, many different factors affect the incidence of esophageal cancer. Conducting further studies at the individual level in areas with high incidence could provide more detailed information on risk factors of esophageal cancer.
Collapse
Affiliation(s)
- Sedigheh Rastaghi
- Department of Biostatistics, School of Public Health, Isfahan University of Medical Sciences, Isfahan, Iran E-mail :
| | | | | |
Collapse
|
34
|
Kang SY, McGree J, Baade P, Mengersen K. A Case Study for Modelling Cancer Incidence Using Bayesian Spatio-Temporal Models. AUST NZ J STAT 2015. [DOI: 10.1111/anzs.12127] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Su Yun Kang
- Mathematical Sciences School; Queensland University of Technology; GPO Box 2434 Brisbane QLD 4001 Australia
- CRC for Spatial Information; 204 Lygon Street Carlton Vic. 3053 Australia
| | - James McGree
- Mathematical Sciences School; Queensland University of Technology; GPO Box 2434 Brisbane QLD 4001 Australia
- CRC for Spatial Information; 204 Lygon Street Carlton Vic. 3053 Australia
| | - Peter Baade
- Viertel Centre for Research in Cancer Control; Cancer Council Queensland; Gregory Terrace Fortitude Valley Australia
- School of Public Health; Queensland University of Technology; GPO Box 2434 Brisbane QLD 4001 Australia
- Griffith Health Institute; Griffith University; Brisbane QLD 4001 Australia
| | - Kerrie Mengersen
- Mathematical Sciences School; Queensland University of Technology; GPO Box 2434 Brisbane QLD 4001 Australia
- CRC for Spatial Information; 204 Lygon Street Carlton Vic. 3053 Australia
| |
Collapse
|
35
|
Carroll R, Lawson AB, Faes C, Kirby RS, Aregay M, Watjou K. Comparing INLA and OpenBUGS for hierarchical Poisson modeling in disease mapping. Spat Spatiotemporal Epidemiol 2015; 14-15:45-54. [PMID: 26530822 DOI: 10.1016/j.sste.2015.08.001] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/25/2015] [Revised: 08/05/2015] [Accepted: 08/06/2015] [Indexed: 11/18/2022]
Abstract
The recently developed R package INLA (Integrated Nested Laplace Approximation) is becoming a more widely used package for Bayesian inference. The INLA software has been promoted as a fast alternative to MCMC for disease mapping applications. Here, we compare the INLA package to the MCMC approach by way of the BRugs package in R, which calls OpenBUGS. We focus on the Poisson data model commonly used for disease mapping. Ultimately, INLA is a computationally efficient way of implementing Bayesian methods and returns nearly identical estimates for fixed parameters in comparison to OpenBUGS, but falls short in recovering the true estimates for the random effects, their precisions, and model goodness of fit measures under the default settings. We assumed default settings for ground truth parameters, and through altering these default settings in our simulation study, we were able to recover estimates comparable to those produced in OpenBUGS under the same assumptions.
Collapse
Affiliation(s)
- R Carroll
- Department of Public Health, Medical University of South Carolina, 135 Cannon St, Charleston, SC 29425, USA.
| | - A B Lawson
- Department of Public Health, Medical University of South Carolina, 135 Cannon St, Charleston, SC 29425, USA
| | - C Faes
- Interuniversity Institute for Statistics and Statistical Bioinformatics, Hasselt University, Agoralaan 1, 3590 Diepenbeek, Belgium
| | - R S Kirby
- Department of Community and Family Health, University of Southern Florida, 13201 Bruce B Downs Blvd MDC 56, Tampa, FL 33612, USA
| | - M Aregay
- Department of Public Health, Medical University of South Carolina, 135 Cannon St, Charleston, SC 29425, USA
| | - K Watjou
- Interuniversity Institute for Statistics and Statistical Bioinformatics, Hasselt University, Agoralaan 1, 3590 Diepenbeek, Belgium
| |
Collapse
|
36
|
Iftimi A, Montes F, Santiyán AM, Martínez-Ruiz F. Space-time airborne disease mapping applied to detect specific behaviour of varicella in Valencia, Spain. Spat Spatiotemporal Epidemiol 2015; 14-15:33-44. [PMID: 26530821 DOI: 10.1016/j.sste.2015.07.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/08/2015] [Revised: 05/14/2015] [Accepted: 07/09/2015] [Indexed: 10/23/2022]
Abstract
Airborne diseases are one of humanity's most feared sicknesses and have regularly caused concern among specialists. Varicella is an airborne disease which usually affects children before the age of 10. Because of its nature, varicella gives rise to interesting spatial, temporal and spatio-temporal patterns. This paper studies spatio-temporal exploratory analysis tools to detect specific behaviour of varicella in the city of Valencia, Spain, from 2008 to 2013. These methods have shown a significant association between the spatial and the temporal component, confirmed by the space-time models applied to the data. High relative risk of varicella is observed in economically disadvantaged regions, areas less involved in vaccination programmes.
Collapse
|
37
|
Estimation of district-level under-5 mortality in Zambia using birth history data, 1980-2010. Spat Spatiotemporal Epidemiol 2014; 11:89-107. [PMID: 25457599 DOI: 10.1016/j.sste.2014.09.002] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/21/2014] [Revised: 09/16/2014] [Accepted: 09/18/2014] [Indexed: 11/24/2022]
Abstract
Birth history data-the primary source of data on under-5 mortality in developing countries-are infrequently used for subnational estimates due to concerns over small sample sizes. In this study we consider different methods for analyzing birth history data in combination with various small area models. We construct a simulation environment to assess the performance of different combinations of birth history methods and small area models in terms of bias, efficiency, and coverage. We find that performance is highly dependent on the birth history method applied and how temporal trends are accounted for. We estimated trends in district-level under-5 mortality in Zambia from 1980 to 2010 using the best-performing model. We find that under-5 mortality is highly variable within Zambia: there was a 1.8-fold difference between the lowest and highest levels in 2010, and declines over the period 1980 to 2010 ranged from less than 5% to more than 50%.
Collapse
|
38
|
Assareh H, Ou L, Chen J, Hillman K, Flabouris A, Hollis SJ. Geographic variation of failure-to-rescue in public acute hospitals in New South Wales, Australia. PLoS One 2014; 9:e109807. [PMID: 25310260 PMCID: PMC4195695 DOI: 10.1371/journal.pone.0109807] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2013] [Accepted: 09/14/2014] [Indexed: 12/21/2022] Open
Abstract
Despite the wide acceptance of Failure-to-Rescue (FTR) as a patient safety indicator (defined as the deaths among surgical patients with treatable complications), no study has explored the geographic variation of FTR in a large health jurisdiction. Our study aimed to explore the spatiotemporal variations of FTR rates across New South Wales (NSW), Australia. We conducted a population-based study using all admitted surgical patients in public acute hospitals during 2002-2009 in NSW, Australia. We developed a spatiotemporal Poisson model using Integrated Nested Laplace Approximation (INLA) methods in a Bayesian framework to obtain area-specific adjusted relative risk. Local Government Area (LGA) was chosen as the areal unit. LGA-aggregated covariates included age, gender, socio-economic and remoteness index scores, distance between patient residential postcode and the treating hospital, and a quadratic time trend. We studied 4,285,494 elective surgical admissions in 82 acute public hospitals over eight years in NSW. Around 14% of patients who developed at least one of the six FTR-related complications (58,590) died during hospitalization. Of 153 LGAs, patients who lived in 31 LGAs, accommodating 48% of NSW patients at risk, were exposed to an excessive adjusted FTR risk (10% to 50%) compared to the state-average. They were mostly located in state's centre and western Sydney. Thirty LGAs with a lower adjusted FTR risk (10% to 30%), accommodating 8% of patients at risk, were mostly found in the southern parts of NSW and Sydney east and south. There were significant spatiotemporal variations of FTR rates across NSW over an eight-year span. Areas identified with significantly high and low FTR risks provide potential opportunities for policy-makers, clinicians and researchers to learn from the success or failure of adopting the best care for surgical patients and build a self-learning organisation and health system.
Collapse
Affiliation(s)
- Hassan Assareh
- Simpson Centre for Health Services Research, Australian Institute of Health Innovation & South Western Sydney Clinical School, University of New South Wales, Sydney, New South Wales, Australia
- Epidemiology, Western Sydney Local Health District, Sydney, New South Wales, Australia
| | - Lixin Ou
- Simpson Centre for Health Services Research, Australian Institute of Health Innovation & South Western Sydney Clinical School, University of New South Wales, Sydney, New South Wales, Australia
| | - Jack Chen
- Simpson Centre for Health Services Research, Australian Institute of Health Innovation & South Western Sydney Clinical School, University of New South Wales, Sydney, New South Wales, Australia
| | - Kenneth Hillman
- Simpson Centre for Health Services Research, Australian Institute of Health Innovation & South Western Sydney Clinical School, University of New South Wales, Sydney, New South Wales, Australia
| | - Arthas Flabouris
- Intensive Care Unit, Royal Adelaide Hospital, Adelaide, South Australia, Australia
| | - Stephanie J. Hollis
- Simpson Centre for Health Services Research, Australian Institute of Health Innovation & South Western Sydney Clinical School, University of New South Wales, Sydney, New South Wales, Australia
| |
Collapse
|
39
|
Affiliation(s)
- Andrew B. Lawson
- Department of Public Health Sciences Medical University of South Carolina Charleston SC USA
| |
Collapse
|
40
|
Ugarte MD, Adin A, Goicoa T, Militino AF. On fitting spatio-temporal disease mapping models using approximate Bayesian inference. Stat Methods Med Res 2014; 23:507-30. [PMID: 24713158 DOI: 10.1177/0962280214527528] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Spatio-temporal disease mapping comprises a wide range of models used to describe the distribution of a disease in space and its evolution in time. These models have been commonly formulated within a hierarchical Bayesian framework with two main approaches: an empirical Bayes (EB) and a fully Bayes (FB) approach. The EB approach provides point estimates of the parameters relying on the well-known penalized quasi-likelihood (PQL) technique. The FB approach provides the posterior distribution of the target parameters. These marginal distributions are not usually available in closed form and common estimation procedures are based on Markov chain Monte Carlo (MCMC) methods. However, the spatio-temporal models used in disease mapping are often very complex and MCMC methods may lead to large Monte Carlo errors and a huge computation time if the dimension of the data at hand is large. To circumvent these potential inconveniences, a new technique called integrated nested Laplace approximations (INLA), based on nested Laplace approximations, has been proposed for Bayesian inference in latent Gaussian models. In this paper, we show how to fit different spatio-temporal models for disease mapping with INLA using the Leroux CAR prior for the spatial component, and we compare it with PQL via a simulation study. The spatio-temporal distribution of male brain cancer mortality in Spain during the period 1986-2010 is also analysed.
Collapse
Affiliation(s)
| | - Aritz Adin
- Department of Statistics and O. R., Public University of Navarre, Spain
| | - Tomas Goicoa
- Department of Statistics and O. R., Public University of Navarre, Spain Research Network on Health Services in Chronic Diseases (REDISSEC), Spain
| | | |
Collapse
|
41
|
Wang XF, Li Y. Bayesian inferences for beta semiparametric-mixed models to analyze longitudinal neuroimaging data. Biom J 2014; 56:662-77. [DOI: 10.1002/bimj.201300003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2013] [Revised: 11/19/2013] [Accepted: 12/06/2013] [Indexed: 11/07/2022]
Affiliation(s)
- Xiao-Feng Wang
- Department of Quantitative Health Sciences/Biostatistics Section; Cleveland Clinic Lerner Research Institute; Cleveland OH 44195 USA
| | - Yingxing Li
- The Wang Yanan Institute for Studies in Economics; Xiamen University; Xiamen China
| |
Collapse
|
42
|
Trends in socioeconomic inequalities in ischemic heart disease mortality in small areas of nine Spanish cities from 1996 to 2007 using smoothed ANOVA. J Urban Health 2014; 91:46-61. [PMID: 23564269 PMCID: PMC3907633 DOI: 10.1007/s11524-013-9799-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
The aim of this study was to analyze the evolution of socioeconomic inequalities in mortality due to ischemic heart diseases (IHD) in the census tracts of nine Spanish cities between the periods 1996-2001 and 2002-2007. Among women, there are socioeconomic inequalities in IHD mortality in the first period which tended to remain stable or even increase in the second period in most of the cities. Among men, in general, no socioeconomic inequalities have been detected for this cause in either of the periods. These results highlight the importance of intra-urban inequalities in mortality due to IHD and their evolution over time.
Collapse
|
43
|
Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol 2013; 7:39-55. [DOI: 10.1016/j.sste.2013.07.003] [Citation(s) in RCA: 116] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
44
|
Sturtz S, Ickstadt K. Comparison of Bayesian methods for flexible modeling of spatial risk surfaces in disease mapping. Biom J 2013; 56:5-22. [PMID: 24130078 DOI: 10.1002/bimj.201200176] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2012] [Revised: 07/16/2013] [Accepted: 07/25/2013] [Indexed: 11/06/2022]
Abstract
Bayesian hierarchical models usually model the risk surface on the same arbitrary geographical units for all data sources. Poisson/gamma random field models overcome this restriction as the underlying risk surface can be specified independently to the resolution of the data. Moreover, covariates may be considered as either excess or relative risk factors. We compare the performance of the Poisson/gamma random field model to the Markov random field (MRF)-based ecologic regression model and the Bayesian Detection of Clusters and Discontinuities (BDCD) model, in both a simulation study and a real data example. We find the BDCD model to have advantages in situations dominated by abruptly changing risk while the Poisson/gamma random field model convinces by its flexibility in the estimation of random field structures and by its flexibility incorporating covariates. The MRF-based ecologic regression model is inferior. WinBUGS code for Poisson/gamma random field models is provided.
Collapse
Affiliation(s)
- Sibylle Sturtz
- Department of Medical Biometry, Institute for Quality and Efficiency in Health Care, Im Mediapark 8, 50670, Köln, Germany; Faculty of Statistics, TU Dortmund University, 44221, Dortmund, Germany
| | | |
Collapse
|
45
|
Kang SY, McGree J, Mengersen K. The impact of spatial scales and spatial smoothing on the outcome of bayesian spatial model. PLoS One 2013; 8:e75957. [PMID: 24146799 PMCID: PMC3795684 DOI: 10.1371/journal.pone.0075957] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2013] [Accepted: 08/19/2013] [Indexed: 12/02/2022] Open
Abstract
Discretization of a geographical region is quite common in spatial analysis. There have been few studies into the impact of different geographical scales on the outcome of spatial models for different spatial patterns. This study aims to investigate the impact of spatial scales and spatial smoothing on the outcomes of modelling spatial point-based data. Given a spatial point-based dataset (such as occurrence of a disease), we study the geographical variation of residual disease risk using regular grid cells. The individual disease risk is modelled using a logistic model with the inclusion of spatially unstructured and/or spatially structured random effects. Three spatial smoothness priors for the spatially structured component are employed in modelling, namely an intrinsic Gaussian Markov random field, a second-order random walk on a lattice, and a Gaussian field with Matérn correlation function. We investigate how changes in grid cell size affect model outcomes under different spatial structures and different smoothness priors for the spatial component. A realistic example (the Humberside data) is analyzed and a simulation study is described. Bayesian computation is carried out using an integrated nested Laplace approximation. The results suggest that the performance and predictive capacity of the spatial models improve as the grid cell size decreases for certain spatial structures. It also appears that different spatial smoothness priors should be applied for different patterns of point data.
Collapse
Affiliation(s)
- Su Yun Kang
- Mathematical Sciences School, Queensland University of Technology, Brisbane, Queensland, Australia
- Cooperative Research Centre for Spatial Information, Melbourne, Victoria, Australia
| | - James McGree
- Mathematical Sciences School, Queensland University of Technology, Brisbane, Queensland, Australia
- Cooperative Research Centre for Spatial Information, Melbourne, Victoria, Australia
| | - Kerrie Mengersen
- Mathematical Sciences School, Queensland University of Technology, Brisbane, Queensland, Australia
- Cooperative Research Centre for Spatial Information, Melbourne, Victoria, Australia
| |
Collapse
|
46
|
Bilancia M, Demarinis G. Bayesian scanning of spatial disease rates with integrated nested Laplace approximation (INLA). STAT METHOD APPL-GER 2013. [DOI: 10.1007/s10260-013-0241-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
47
|
Alegana VA, Atkinson PM, Wright JA, Kamwi R, Uusiku P, Katokele S, Snow RW, Noor AM. Estimation of malaria incidence in northern Namibia in 2009 using Bayesian conditional-autoregressive spatial-temporal models. Spat Spatiotemporal Epidemiol 2013; 7:25-36. [PMID: 24238079 PMCID: PMC3839406 DOI: 10.1016/j.sste.2013.09.001] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/23/2013] [Revised: 08/05/2013] [Accepted: 09/05/2013] [Indexed: 10/29/2022]
Abstract
As malaria transmission declines, it becomes increasingly important to monitor changes in malaria incidence rather than prevalence. Here, a spatio-temporal model was used to identify constituencies with high malaria incidence to guide malaria control. Malaria cases were assembled across all age groups along with several environmental covariates. A Bayesian conditional-autoregressive model was used to model the spatial and temporal variation of incidence after adjusting for test positivity rates and health facility utilisation. Of the 144,744 malaria cases recorded in Namibia in 2009, 134,851 were suspected and 9893 were parasitologically confirmed. The mean annual incidence based on the Bayesian model predictions was 13 cases per 1000 population with the highest incidence predicted for constituencies bordering Angola and Zambia. The smoothed maps of incidence highlight trends in disease incidence. For Namibia, the 2009 maps provide a baseline for monitoring the targets of pre-elimination.
Collapse
Affiliation(s)
- Victor A Alegana
- Malaria Public Health Department, KEMRI-Wellcome Trust-University of Oxford Collaborative Programme, P.O. Box 43640, 00100 GPO Nairobi, Kenya; Centre for Geographical Health Research, Geography and Environment, University of Southampton, Highfield, Southampton SO17 1BJ, UK.
| | | | | | | | | | | | | | | |
Collapse
|
48
|
Maheu-Giroux M, Castro MC. Impact of community-based larviciding on the prevalence of malaria infection in Dar es Salaam, Tanzania. PLoS One 2013; 8:e71638. [PMID: 23977099 PMCID: PMC3743749 DOI: 10.1371/journal.pone.0071638] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2013] [Accepted: 07/01/2013] [Indexed: 12/18/2022] Open
Abstract
Background The use of larval source management is not prioritized by contemporary malaria control programs in sub-Saharan Africa despite historical success. Larviciding, in particular, could be effective in urban areas where transmission is focal and accessibility to Anopheles breeding habitats is generally easier than in rural settings. The objective of this study is to assess the effectiveness of a community-based microbial larviciding intervention to reduce the prevalence of malaria infection in Dar es Salaam, United Republic of Tanzania. Methods and Findings Larviciding was implemented in 3 out of 15 targeted wards of Dar es Salaam in 2006 after two years of baseline data collection. This intervention was subsequently scaled up to 9 wards a year later, and to all 15 targeted wards in 2008. Continuous randomized cluster sampling of malaria prevalence and socio-demographic characteristics was carried out during 6 survey rounds (2004–2008), which included both cross-sectional and longitudinal data (N = 64,537). Bayesian random effects logistic regression models were used to quantify the effect of the intervention on malaria prevalence at the individual level. Effect size estimates suggest a significant protective effect of the larviciding intervention. After adjustment for confounders, the odds of individuals living in areas treated with larviciding being infected with malaria were 21% lower (Odds Ratio = 0.79; 95% Credible Intervals: 0.66–0.93) than those who lived in areas not treated. The larviciding intervention was most effective during dry seasons and had synergistic effects with other protective measures such as use of insecticide-treated bed nets and house proofing (i.e., complete ceiling or window screens). Conclusion A large-scale community-based larviciding intervention significantly reduced the prevalence of malaria infection in urban Dar es Salaam.
Collapse
Affiliation(s)
- Mathieu Maheu-Giroux
- Department of Global Health and Population, Harvard School of Public Health, Boston, Massachusetts, United States of America
| | - Marcia C. Castro
- Department of Global Health and Population, Harvard School of Public Health, Boston, Massachusetts, United States of America
- * E-mail:
| |
Collapse
|
49
|
Blangiardo M, Cameletti M, Baio G, Rue H. Spatial and spatio-temporal models with R-INLA. Spat Spatiotemporal Epidemiol 2013; 4:33-49. [DOI: 10.1016/j.sste.2012.12.001] [Citation(s) in RCA: 158] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/11/2012] [Revised: 11/28/2012] [Accepted: 12/05/2012] [Indexed: 10/27/2022]
|
50
|
Abstract
Spatial statistics and spatial econometrics fill an impressive body of literature, clearly demonstrating the importance of the corresponding methods. Though the two fields are being referred to in the same breath, they differ substantially as they pursue quite different objectives and hence exhibit different interpretations. In this article we shed some light on the nature and relevance of those to demonstrate how the models can be employed with three data examples in a regional context. The discussion covers models for continuous responses and count data, real data examples and Monte Carlo simulations.
Collapse
|