A comparison of conditional autoregressive models used in Bayesian disease mapping

Spatiotemporal patterns and socioeconomic determinants of pulmonary tuberculosis in Dongguan city, China, during 2011-2020: an ecological study

OBJECTIVE

Pulmonary tuberculosis (PTB) is a critical challenge worldwide, particularly in China. This study aimed to explore the spatiotemporal transmission patterns and socioeconomic factors of PTB in Dongguan city, China.

METHODS/DESIGN

An ecological study based on the reported new PTB cases between 2011 and 2020 was conducted in Dongguan city, China. The spatiotemporal analysis methods were used to explore the long-term trend, spatiotemporal transmission pattern and socioeconomic factors of PTB.

MAIN OUTCOME MEASURES

The number of new PTB cases.

PARTICIPANTS

We collected 35 756 new PTB cases, including 23 572 males and 12 184 females.

RESULTS

The seasonal-trend decomposition indicated a significant downward trend for PTB with a significant peak in 2017 and 2018, and local spatial autocorrelation showed more and more high-high clusters in the central and north-central towns with high incidence. The multivariate spatial time series analysis revealed that the endemic component had a leading role in driving PTB transmission, with a high total effect value being 189.40 (95% CI: 171.65-207.15). A Bayesian spatiotemporal model revealed that PTB incidence is positively associated with the agricultural population ratio (relative risk (RR) =1.074), gender ratio (RR=1.104) and the number of beds in medical institutions (RR=1.028).

CONCLUSIONS

These findings revealed potential spatiotemporal variability and spatial aggregation of PTB, so targeted preventive strategies should be made in different towns based on spatiotemporal transmission patterns and risk factors.

Improving the spatial and temporal resolution of burden of disease measures with Bayesian models

This paper contributes to the field by addressing the critical issue of enhancing the spatial and temporal resolution of health data. Although Bayesian methods are frequently employed to address this challenge in various disciplines, the application of Bayesian spatio-temporal models to burden of disease (BOD) studies remains limited. Our novelty lies in the exploration of two existing Bayesian models that we show to be applicable to a wide range of BOD data, including mortality and prevalence, thereby providing evidence to support the adoption of Bayesian modeling in full BOD studies in the future. We illustrate the benefits of Bayesian modeling with an Australian case study involving asthma and coronary heart disease. Our results showcase the effectiveness of Bayesian approaches in increasing the number of small areas for which results are available and improving the reliability and stability of the results compared to using data directly from surveys or administrative sources.

Spatiotemporal Bayesian modeling of the risk of congenital syphilis in São Paulo, SP, Brazil

The aim of this study is to analyze the spatiotemporal risk of congenital syphilis (CS) in high-prevalence areas in the city of São Paulo, SP, Brazil, and to evaluate its relationship with socioeconomic, demographic, and environmental variables. An ecological study was conducted based on secondary CS data with spatiotemporal components collected from 310 areas between 2010 and 2016. The data were modeled in a Bayesian context using the integrated nested Laplace approximation (INLA) method. Risk maps showed an increasing CS trend over time and highlighted the areas that presented the highest and lowest risk in each year. The model showed that the factors positively associated with a higher risk of CS were the Gini index and the proportion of women aged 18-24 years without education or with incomplete primary education, while the factors negatively associated were the proportion of women of childbearing age and the mean per capita income.

Geographical distribution of antimicrobial exposure among very preterm and very low birth weight infants: A nationwide database study in Japan

OBJECTIVES

To examine spatial effects in neonatal care, we conducted a retrospective cohort study to investigate the geographical distribution of antimicrobial exposure among very preterm and very low birth weight infants in Japan.

STUDY DESIGN

We utilized a nationwide claims database in Japan to extract prescriptions of injectable antimicrobials for 41,423 very preterm and very low birth weight infants admitted within the first two days of life from April 2010 to March 2021. We identified frequently prescribed antimicrobials, revealed early neonatal exposure and neonatal exposure to each antimicrobial agent by 47 prefectures in Japan, and evaluated their spatial autocorrelation using global and local Moran's I statistics. We then scrutinized regional disparities in antimicrobial drug prescriptions.

RESULTS

The top 10 antimicrobials prescribed to very preterm and very low birth weight infants in Japan were ampicillin, amikacin, gentamicin, cefotaxime, fluconazole, ampicillin combination, micafungin, cefmetazole, cefazolin, and vancomycin. We identified northern cold spots for fluconazole exposure and southern hot spots for ampicillin, amikacin, gentamicin, and cefmetazole exposure. Geographical heterogeneity in the selection of antibacterial and antimycotic agents was observed.

CONCLUSION

Our study revealed the geographical distribution of antimicrobial exposure among very preterm and very low birth weight infants in Japan, thus disclosing its spatial effects. Further research addressing the spatial effects of neonatal care is needed to understand how drug exposure affects the outcomes of preterm infants.

Causal decomposition maps: An exploratory tool for designing area-level interventions aimed at reducing health disparities

Methods for decomposition analyses have been developed to partition between-group differences into explained and unexplained portions. In this paper, we introduce the concept of causal decomposition maps, which allow researchers to test the effect of area-level interventions on disease maps before implementation. These maps quantify the impact of interventions that aim to reduce differences in health outcomes between groups and illustrate how the disease map might change under different interventions. We adapt a new causal decomposition analysis method for the disease mapping context. Through the specification of a Bayesian hierarchical outcome model, we obtain counterfactual small area estimates of age-adjusted rates and reliable estimates of decomposition quantities. We present two formulations of the outcome model, with the second allowing for spatial interference of the intervention. Our method is utilized to determine whether the addition of gyms in different sets of rural ZIP codes could reduce any of the rural-urban difference in age-adjusted colorectal cancer incidence rates in Iowa ZIP codes.

Estimation of the relationship between meteorological factors and measles using spatiotemporal Bayesian model in Shandong Province, China

BACKGROUND

Measles-containing vaccine (MCV) has been effective in controlling the spread of measles. Some countries have declared measles elimination. But recently years, the number of cases worldwide has increased, posing a challenge to the global goal of measles eradication. This study estimated the relationship between meteorological factors and measles using spatiotemporal Bayesian model, aiming to provide scientific evidence for public health policy to eliminate measles.

METHODS

Descriptive statistical analysis was performed on monthly data of measles and meteorological variables in 136 counties of Shandong Province from 2009 to 2017. Spatiotemporal Bayesian model was used to estimate the effects of meteorological factors on measles, and to evaluate measles risk areas at county level. Case population was divided into multiple subgroups according to gender, age and occupation. The effects of meteorological factors on measles in subgroups were compared.

RESULTS

Specific meteorological conditions increased the risk of measles, including lower relative humidity, temperature, and atmospheric pressure; higher wind velocity, sunshine duration, and diurnal temperature variation. Taking lowest value (Q1) as reference, RR (95%CI) for higher temperatures (Q2-Q4) were 0.79 (0.69-0.91), 0.54 (0.44-0.65), and 0.48 (0.38-0.61), respectively; RR (95%CI) for higher relative humidity (Q2-Q4) were 0.76 (0.66-0.88), 0.56 (0.47-0.67), and 0.49 (0.38-0.63), respectively; RR (95%CI) for higher wind velocity (Q2-Q4) were 1.43 (1.25-1.64), 1.85 (1.57-2.18), 2.00 (1.59-2.52), respectively. 22 medium-to-high risk counties were identified, mainly in northwestern, southwestern and central Shandong Province. The trend was basically same in the effects of meteorological factors on measles in subgroups, but the magnitude of the effects was different.

CONCLUSIONS

Meteorological factors have an important impact on measles. It is crucial to integrate these factors into public health policies for measles prevention and control in China.

A Systematic Review of Areal Units and Adjacency Used in Bayesian Spatial and Spatio-Temporal Conditional Autoregressive Models in Health Research

Advancements in Bayesian spatial and spatio-temporal modelling have been observed in recent years. Despite this, there are unresolved issues about the choice of appropriate spatial unit and adjacency matrix in disease mapping. There is limited systematic review evidence on this topic. This review aimed to address these problems. We searched seven databases to find published articles on this topic. A modified quality assessment tool was used to assess the quality of studies. A total of 52 studies were included, of which 26 (50.0%) were on infectious diseases, 10 (19.2%) on chronic diseases, 8 (15.5%) on maternal and child health, and 8 (15.5%) on other health-related outcomes. Only 6 studies reported the reasons for using the specified spatial unit, 8 (15.3%) studies conducted sensitivity analysis for prior selection, and 39 (75%) of the studies used Queen contiguity adjacency. This review highlights existing variation and limitations in the specification of Bayesian spatial and spatio-temporal models used in health research. We found that majority of the studies failed to report the rationale for the choice of spatial units, perform sensitivity analyses on the priors, or evaluate the choice of neighbourhood adjacency, all of which can potentially affect findings in their studies.

A practical illustration of spatial smoothing methods for disconnected regions with INLA: spatial survey on overweight and obesity in Malaysia

BACKGROUND

National prevalence could mask subnational heterogeneity in disease occurrence, and disease mapping is an important tool to illustrate the spatial pattern of disease. However, there is limited information on techniques for the specification of conditional autoregressive models in disease mapping involving disconnected regions. This study explores available techniques for producing district-level prevalence estimates for disconnected regions, using as an example childhood overweight in Malaysia, which consists of the Peninsular and Borneo regions separated by the South China Sea. We used data from Malaysia National Health and Morbidity Survey conducted in 2015. We adopted Bayesian hierarchical modelling using the integrated nested Laplace approximation (INLA) program in R-software to model the spatial distribution of overweight among 6301 children aged 5-17 years across 144 districts located in two disconnected regions. We illustrate different types of spatial models for prevalence mapping across disconnected regions, taking into account the survey design and adjusting for district-level demographic and socioeconomic covariates.

RESULTS

The spatial model with split random effects and a common intercept has the lowest Deviance and Watanabe Information Criteria. There was evidence of a spatial pattern in the prevalence of childhood overweight across districts. An increasing trend in smoothed prevalence of overweight was observed when moving from the east to the west of the Peninsular and Borneo regions. The proportion of Bumiputera ethnicity in the district had a significant negative association with childhood overweight: the higher the proportion of Bumiputera ethnicity in the district, the lower the prevalence of childhood overweight.

CONCLUSION

This study illustrates different available techniques for mapping prevalence across districts in disconnected regions using survey data. These techniques can be utilized to produce reliable subnational estimates for any areas that comprise of disconnected regions. Through the example, we learned that the best-fit model was the one that considered the separate variations of the individual regions. We discovered that the occurrence of childhood overweight in Malaysia followed a spatial pattern with an east-west gradient trend, and we identified districts with high prevalence of overweight. This information could help policy makers in making informed decisions for targeted public health interventions in high-risk areas.

Analysis of the injury severity of motor vehicle-pedestrian crashes at urban intersections using spatiotemporal logistic regression models

This paper conducted a comprehensive study on the injury severity of motor vehicle-pedestrian crashes at 489 urban intersections across a dense road network based on high-resolution accident data recorded by the police from 2010 to 2019 in Hong Kong. Given that accounting for the spatial and temporal correlations simultaneously among crash data can contribute to unbiased parameter estimations for exogenous variables and improved model performance, we developed spatiotemporal logistic regression models with various spatial formulations and temporal configurations. The results indicated that the model with the Leroux conditional autoregressive prior and random walk structure outperformed other alternatives in terms of goodness-of-fit and classification accuracy. According to the parameter estimates, pedestrian age, head injury, pedestrian location, pedestrian actions, driver maneuvers, vehicle type, first point of collision, and traffic congestion status significantly affected the severity of pedestrian injuries. On the basis of our analysis, a range of targeted countermeasures integrating safety education, traffic enforcement, road design, and intelligent traffic technologies were proposed to improve the safe mobility of pedestrians at urban intersections. The present study provides a rich and sound toolkit for safety analysts to deal with spatiotemporal correlations when modeling crashes aggregated at contiguous spatial units within multiple years.

The Potential of Surveillance Data for Dengue Risk Mapping: An Evaluation of Different Approaches in Cuba

To better guide dengue prevention and control efforts, the use of routinely collected data to develop risk maps is proposed. For this purpose, dengue experts identified indicators representative of entomological, epidemiological and demographic risks, hereafter called components, by using surveillance data aggregated at the level of Consejos Populares (CPs) in two municipalities of Cuba (Santiago de Cuba and Cienfuegos) in the period of 2010-2015. Two vulnerability models (one with equally weighted components and one with data-derived weights using Principal Component Analysis), and three incidence-based risk models were built to construct risk maps. The correlation between the two vulnerability models was high (tau > 0.89). The single-component and multicomponent incidence-based models were also highly correlated (tau ≥ 0.9). However, the agreement between the vulnerability- and the incidence-based risk maps was below 0.6 in the setting with a prolonged history of dengue transmission. This may suggest that an incidence-based approach does not fully reflect the complexity of vulnerability for future transmission. The small difference between single- and multicomponent incidence maps indicates that in a setting with a narrow availability of data, simpler models can be used. Nevertheless, the generalized linear mixed multicomponent model provides information of covariate-adjusted and spatially smoothed relative risks of disease transmission, which can be important for the prospective evaluation of an intervention strategy. In conclusion, caution is needed when interpreting risk maps, as the results vary depending on the importance given to the components involved in disease transmission. The multicomponent vulnerability mapping needs to be prospectively validated based on an intervention trial targeting high-risk areas.

A Systematic Review of Joint Spatial and Spatiotemporal Models in Health Research

With the advancement of spatial analysis approaches, methodological research addressing the technical and statistical issues related to joint spatial and spatiotemporal models has increased. Despite the benefits of spatial modelling of several interrelated outcomes simultaneously, there has been no published systematic review on this topic, specifically when such models would be useful. This systematic review therefore aimed at reviewing health research published using joint spatial and spatiotemporal models. A systematic search of published studies that applied joint spatial and spatiotemporal models was performed using six electronic databases without geographic restriction. A search with the developed search terms yielded 4077 studies, from which 43 studies were included for the systematic review, including 15 studies focused on infectious diseases and 11 on cancer. Most of the studies (81.40%) were performed based on the Bayesian framework. Different joint spatial and spatiotemporal models were applied based on the nature of the data, population size, the incidence of outcomes, and assumptions. This review found that when the outcome is rare or the population is small, joint spatial and spatiotemporal models provide better performance by borrowing strength from related health outcomes which have a higher prevalence. A framework for the design, analysis, and reporting of such studies is also needed.

Bayesian Spatial Modelling of HIV Prevalence in Jimma Zone, Ethiopia

Background: Although the human immunodeficiency virus (HIV) is spatially heterogeneous in Ethiopia, current regional estimates of HIV prevalence hide the epidemic’s heterogeneity. A thorough examination of the prevalence of HIV infection using district-level data could assist to develop HIV prevention strategies. The aims of this study were to examine the spatial clustering of HIV prevalence in Jimma Zone at district level and assess the effects of patient characteristics on the prevalence of HIV infection. Methods: The 8440 files of patients who underwent HIV testing in the 22 Districts of Jimma Zone between September 2018 and August 2019 were the source of data for this study. The global Moran’s index, Getis–Ord Gi* local statistic, and Bayesian hierarchical spatial modelling approach were applied to address the research objectives. Results: Positive spatial autocorrelation was observed in the districts and the local indicators of spatial analysis using the Getis–Ord statistic also identified three districts, namely Agaro, Gomma and Nono Benja, as hotspots, and two districts, namely Mancho and Omo Beyam, as coldspots with 95% and 90% confidence levels, respectively, for HIV prevalence. The results also showed eight patient-related characteristics that were considered in the study were associated with HIV prevalence in the study area. Furthermore, after accounting for these characteristics in the fitted model, there was no spatial clustering of HIV prevalence suggesting the patient characteristics had explained most of the heterogeneity in HIV prevalence in Jimma Zone for the study data. Conclusions: The identification of hotspot districts and the spatial dynamic of HIV infection in Jimma Zone at district level may allow health policymakers in the zone or Oromiya region or at national level to develop geographically specific strategies to prevent HIV transmission. Because clinic register data were used in the study, it is important to use caution when interpreting the results. The results are restricted to Jimma Zone districts and may not be generalizable to Ethiopia or the Oromiya region.

Role of social determinants of health in differential respiratory exposure and health outcomes among children

BACKGROUND

Attributes defining the Social Determinants of Health (SDoH) are associated with disproportionate exposures to environmental hazards and differential health outcomes among communities. The dynamics between SDoH, disproportionate environmental exposures, and differential health outcomes are often specific to micro-geographic areas.

METHODS

This study focused on children less than 20 years of age who lived in Douglas County, Nebraska, during 2016-2019. To assess the role of SDoH in differential exposures, we evaluated the association between SDoH metrics and criteria pollutant concentrations and the association between SDoH and pediatric asthma exacerbations to quantify the role of SDoH in differential pediatric asthma outcomes. The Bayesian Poisson regression model with spatial random effects was used to evaluate associations.

RESULTS

We identified significant positive associations between the annual mean concentration of criteria pollutants (carbon monoxide, particulate matter_2.5, nitrogen dioxide, sulfur dioxide) with race (Non-Hispanic Black and Hispanic/Latino), financial stability, and literacy. Additionally, there were significant positive associations between higher rates of pediatric asthma emergency department visits and neighborhoods with more Non-Hispanic Black children, children without health insurance coverage, and households without access to a vehicle.

CONCLUSIONS

Non-Hispanic Black and Hispanic/Latino children living in Douglas County, NE experience disproportionately higher exposure to criteria pollutant concentrations. Additionally, higher rates of asthma exacerbations among Non-Hispanic Black children could be due to reduced access to respiratory care that is potentially the result of financial instability and vehicle access. These results could inform city planners and health care providers to mitigate respiratory risks among these higher at-risk populations.

Revisiting Gaussian Markov random fields and Bayesian disease mapping

We revisit several conditionally formulated Gaussian Markov random fields, known as the intrinsic conditional autoregressive model, the proper conditional autoregressive model, and the Leroux et al. conditional autoregressive model, as well as convolution models such as the well known Besag, York and Mollie model, its (adaptive) re-parameterization, and its scaled alternatives, for their roles of modelling underlying spatial risks in Bayesian disease mapping. Analytic and simulation studies, with graphic visualizations, and disease mapping case studies, present insights and critique on these models for their nature and capacities in characterizing spatial dependencies, local influences, and spatial covariance and correlation functions, and in facilitating stabilized and efficient posterior risk prediction and inference. It is illustrated that these models are Gaussian (Markov) random fields of different spatial dependence, local influence, and (covariance) correlation functions and can play different and complementary roles in Bayesian disease mapping applications.

Deconvolution algorithms for inference of the cell-type composition of the spatial transcriptome

The spatial transcriptome has enabled researchers to resolve transcriptome expression profiles while preserving information about cell location to better understand the complex biological processes that occur in organisms. Due to technical limitations, the current high-throughput spatial transcriptome sequencing methods (known as next-generation sequencing with spatial barcoding methods or spot-based methods) cannot achieve single-cell resolution. A single measurement site, called a spot, in these technologies frequently contains multiple cells of various types. Computational tools for determining the cellular composition of a spot have emerged as a way to break through these limitations. These tools are known as deconvolution tools. Recently, a couple of deconvolution tools based on different strategies have been developed and have shown promise in different aspects. The resulting single-cell resolution expression profiles and/or single-cell composition of spots will significantly affect downstream data mining; thus, it is crucial to choose a suitable deconvolution tool. In this review, we present a list of currently available tools for spatial transcriptome deconvolution, categorize them based on the strategies they employ, and explain their advantages and limitations in detail in order to guide the selection of these tools in future studies.

A joint Bayesian spatiotemporal risk prediction model of COVID-19 incidence, IC admission, and death with application to Sweden

The three closely related COVID-19 outcomes of incidence, intensive care (IC) admission and death, are commonly modelled separately leading to biased estimation of the parameters and relatively poor forecasts. This paper presents a joint spatiotemporal model of the three outcomes based on weekly data that is used for risk prediction and identification of hotspots. The paper applies a pure spatiotemporal model consisting of structured and unstructured spatial and temporal effects and their interaction capturing the effects of the unobserved covariates. The pure spatiotemporal model limits the data requirements to the three outcomes and the population at risk per spatiotemporal unit. The empirical study for the 21 Swedish regions for the period 1 January 2020-4 May 2021 confirms that the joint model predictions outperform the separate model predictions. The fifteen-week-ahead spatiotemporal forecasts (5 May-11 August 2021) show a significant decline in the relative risk of COVID-19 incidence, IC admission, death and number of hotspots.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00168-022-01191-1.

A unified Gaussian copula methodology for spatial regression analysis

Spatially referenced data arise in many fields, including imaging, ecology, public health, and marketing. Although principled smoothing or interpolation is paramount for many practitioners, regression, too, can be an important (or even the only or most important) goal of a spatial analysis. When doing spatial regression it is crucial to accommodate spatial variation in the response variable that cannot be explained by the spatially patterned explanatory variables included in the model. Failure to model both sources of spatial dependence-regression and extra-regression, if you will-can lead to erroneous inference for the regression coefficients. In this article I highlight an under-appreciated spatial regression model, namely, the spatial Gaussian copula regression model (SGCRM), and describe said model's advantages. Then I develop an intuitive, unified, and computationally efficient approach to inference for the SGCRM. I demonstrate the efficacy of the proposed methodology by way of an extensive simulation study along with analyses of a well-known dataset from disease mapping.

A Bayesian modelling framework to quantify multiple sources of spatial variation for disease mapping

Spatial connectivity is an important consideration when modelling infectious disease data across a geographical region. Connectivity can arise for many reasons, including shared characteristics between regions and human or vector movement. Bayesian hierarchical models include structured random effects to account for spatial connectivity. However, conventional approaches require the spatial structure to be fully defined prior to model fitting. By applying penalized smoothing splines to coordinates, we create two-dimensional smooth surfaces describing the spatial structure of the data while making minimal assumptions about the structure. The result is a non-stationary surface which is setting specific. These surfaces can be incorporated into a hierarchical modelling framework and interpreted similarly to traditional random effects. Through simulation studies, we show that the splines can be applied to any symmetric continuous connectivity measure, including measures of human movement, and that the models can be extended to explore multiple sources of spatial structure in the data. Using Bayesian inference and simulation, the relative contribution of each spatial structure can be computed and used to generate hypotheses about the drivers of disease. These models were found to perform at least as well as existing modelling frameworks, while allowing for future extensions and multiple sources of spatial connectivity.

Spatially informed cell-type deconvolution for spatial transcriptomics

Many spatially resolved transcriptomic technologies do not have single-cell resolution but measure the average gene expression for each spot from a mixture of cells of potentially heterogeneous cell types. Here, we introduce a deconvolution method, conditional autoregressive-based deconvolution (CARD), that combines cell-type-specific expression information from single-cell RNA sequencing (scRNA-seq) with correlation in cell-type composition across tissue locations. Modeling spatial correlation allows us to borrow the cell-type composition information across locations, improving accuracy of deconvolution even with a mismatched scRNA-seq reference. CARD can also impute cell-type compositions and gene expression levels at unmeasured tissue locations to enable the construction of a refined spatial tissue map with a resolution arbitrarily higher than that measured in the original study and can perform deconvolution without an scRNA-seq reference. Applications to four datasets, including a pancreatic cancer dataset, identified multiple cell types and molecular markers with distinct spatial localization that define the progression, heterogeneity and compartmentalization of pancreatic cancer.

Spatial patterning of occupational stress and its related factors in Iranian critical care nurses using a hierarchical Bayesian technique

BACKGROUND: In each community, health problems’ patterns and geographical changes are of prime importance to determine high and low-risk areas. OBJECTIVES: This study aimed to investigate the Spatial patterning of occupational stress and its related factors in Iranian critical care nurses using a hierarchical Bayesian technique METHODS: The current research was a cross-sectional descriptive-analytical study. The data includes the number of critical care unit nurses who show a high stress level based on a questionnaire. We used variables such as age, gender, collaboration status, working time, marital status, clinical experience, education, supervisor support, stress score, and working on holiday days for this study. The survey participants had to be at least 18 years old, a registered nurse, and working in the intensive care unit (ICU). OpenBUGS version 3.2.3 was used to implement the Bayesian hierarchical Poisson model and find partial patterning of occupational stress and its related factors. RESULTS: The final sample size was 17414 nurses. The overall prevalence of occupational stress in ICU nurses was estimated at 70%. The lowest and highest prevalence was 65.8% in the North Khorasan province and 75.2% in Golestan province. Occupational stress had a statistically significant association with collaboration status, but with demographic variables, shift work, supportive supervisor, and working on holidays had no statistically significant association. CONCLUSIONS: According to the findings, it is necessary to eliminate or reduce job stress and increase efficiency in Iranian nurses, encourage teamwork and collaboration as an essential element of a healthy workplace environment.

Bayesian disease mapping: Past, present, and future

On the occasion of the Spatial Statistics' 10th Anniversary, I reflect on the past and present of Bayesian disease mapping and look into its future. I focus on some key developments of models, and on recent evolution of multivariate and adaptive Gaussian Markov random fields and their impact and importance in disease mapping. I reflect on Bayesian disease mapping as a subject of spatial statistics that has advanced to date, and continues to grow, in scope and complexity alongside increasing needs of analytic tools for contemporary health science research, such as spatial epidemiology, population and public health, and medicine. I illustrate (potential) utility and impact of some of the disease mapping models and methods for analysing and monitoring communicable disease such as the COVID-19 infection risks during an ongoing pandemic.

Geographical variation of COVID-19 vaccination coverage, ethnic diversity and population composition in Flanders

The vaccination coverage in Flanders is high, but some regions show lower vaccination willingness as compared to the overall vaccination coverage. Beginning November of 2021, the vaccination rate in Flanders was above 93% in age groups above 45 years, and around 85% in the age groups 12 to 44 years. Apart from Flanders as a whole, focus here is on the health sector Maasland, which has a slightly lower vaccination rate, especially in the age groups 12 to 44 years. In the Maasland region, located on the eastern border of Flanders, there are between 1% and 10% less vaccinated individuals than expected according to the vaccination rate in the whole of Flanders, with lowest vaccination rates in the south of the Maasland region. We study the impact of ethnic diversity in the population, population composition with respect to the ethnicity of individuals (in the sense of how the local population composition differs from the Flemish average), and socio-economic status on the vaccination rate at the level of the statistical sector, apart from the effect of age. We explain the statistical methods to investigate geographical differences and illustrate how one can deal with incomplete information in vaccination registries. Ethnic diversity in a region is associated with lower vaccination rates, as is a lower regional socio-economic status. The composition of the population in Maasland is associated with a 35% reduction in the odds to get vaccinated as compared to the overall Flemish population.

Accessibility to rabies centers and human rabies post-exposure prophylaxis rates in Cambodia: A Bayesian spatio-temporal analysis to identify optimal locations for future centers

Rabies is endemic in Cambodia. For exposed humans, post-exposure prophylaxis (PEP) is very effective in preventing this otherwise fatal disease. The Institut Pasteur du Cambodge (IPC) in Phnom Penh was the primary distributor of PEP in Cambodia until 2018. Since then, and to increase distribution of PEP, two new centers have been opened by IPC in the provinces of Battambang and Kampong Cham. Data on bitten patients, who sometimes bring the head of the biting animal for rabies analyses, have been recorded by IPC since 2000. However, human cases are not routinely recorded in Cambodia, making it difficult to establish a human burden of disease and generate a risk map of dog bites to inform the selection of future PEP center locations in high-risk areas. Our aim was to assess the impact of accessibility to rabies centers on the yearly rate of PEP patients in the population and generate a risk map to identify the locations where new centers would be the most beneficial to the Cambodian population. To accomplish this, we used spatio-temporal Bayesian regression models with the number of PEP patients as the outcome. The primary exposure variable considered was travel time to the nearest IPC center. Secondary exposure variables consisted of travel time to a provincial capital and urban proportion of the population. Between 2000 and 2016, a total of 293,955 PEP patient records were identified. Our results showed a significant negative association between travel time to IPC and the rate of PEP patients: an increase in one hour travel time from the living location to IPC PEP centers leads to a reduction in PEP rate of 70% to 80%. Five provinces were identified as the most efficient locations for future centers to maximize PEP accessibility: Banteay Meanchey, Siem Reap, Takeo, Kampot and Svay Rieng. Adding a PEP center in every provincial capital would increase the proportion of Cambodians living within 60 minutes of a PEP center from 26.6% to 64.9%, and living within 120 minutes from 52.8% to 93.3%, which could save hundreds of lives annually.

Rabies is a fatal viral disease that affects the nervous system. It is endemic in many countries in Africa and Asia where free roaming dogs form a reservoir. Transmission to humans occurs most often through a dog bite. However, post-exposure prophylaxis (PEP), if administered before symptom onset, is highly effective at preventing the disease. In Cambodia, a few number of centers offer PEP, with the Institut Pasteur du Cambodge in Phnom Penh being the main one. These few locations lead to limited accessibility for rural areas distant from Phnom Penh and underestimations of the dog-bite burden and PEP needs. Through statistical modelling, we measured the impact of accessibility on the number of PEP patients and predicted the impact of opening new centers in other locations. We found that travel time was significantly associated with the rate of PEP patients. IPC opened new rabies centers in Battambang and Kampong Cham provinces in 2018 and 2019, respectively, and we identified four provinces where future openings would be the most beneficial: Banteay Meanchey, Siem Reap, Takeo, Kampot and Svay Rieng. This study is part of a broader drive to eradicate rabies in Cambodia by 2030 through increased PEP infrastructure and control measures in the dog population.

Spatial-temporal generalized additive model for modeling COVID-19 mortality risk in Toronto, Canada

This article presents a spatial-temporal generalized additive model for modeling geo-referenced COVID-19 mortality data in Toronto, Canada. A range of factors and spatial-temporal terms are incorporated into the model. The non-linear and interactive effects of the neighborhood-level factors, i.e., population density and average of income, are modeled as a two-dimensional spline smoother. The change of spatial pattern over time is modeled as a three-dimensional tensor product smoother. By fitting this model, the space-time effect can uncover the underlying spatial-temporal pattern that is not explained by the covariates. The performance of the modeling method based on the individual data is also compared to the modeling methods based on the aggregated data in terms of in-sample and out-of-sample predictive checking. The results suggest that the individual-level based analysis provided a better overall model fit and higher predictive accuracy for detecting epidemic peaks in this application as compared to the analysis based on the aggregated data.

Evaluation of spatial Bayesian Empirical Likelihood models in analysis of small area data

Bayesian empirical likelihood (BEL) models are becoming increasingly popular as an attractive alternative to fully parametric models. However, they have only recently been applied to spatial data analysis for small area estimation. This study considers the development of spatial BEL models using two popular conditional autoregressive (CAR) priors, namely BYM and Leroux priors. The performance of the proposed models is compared with their parametric counterparts and with existing spatial BEL models using independent Gaussian priors and generalised Moran basis priors. The models are applied to two benchmark spatial datasets, simulation study and COVID-19 data. The results indicate promising opportunities for these models to capture new insights into spatial data. Specifically, the spatial BEL models outperform the parametric spatial models when the underlying distributional assumptions of data appear to be violated.

Identifying extreme COVID-19 mortality risks in English small areas: a disease cluster approach

The COVID-19 pandemic is having a huge impact worldwide and has highlighted the extent of health inequalities between countries but also in small areas within a country. Identifying areas with high mortality is important both of public health mitigation in COVID-19 outbreaks, and of longer term efforts to tackle social inequalities in health. In this paper we consider different statistical models and an extension of a recent method to analyze COVID-19 related mortality in English small areas during the first wave of the epidemic in the first half of 2020. We seek to identify hotspots, and where they are most geographically concentrated, taking account of observed area factors as well as spatial correlation and clustering in regression residuals, while also allowing for spatial discontinuities. Results show an excess of COVID-19 mortality cases in small areas surrounding London and in other small areas in North-East and and North-West of England. Models alleviating spatial confounding show ethnic isolation, air quality and area morbidity covariates having a significant and broadly similar impact on COVID-19 mortality, whereas nursing home location seems to be slightly less important.

Identifying extreme COVID-19 mortality risks in English small areas: a disease cluster approach

The COVID-19 pandemic is having a huge impact worldwide and has highlighted the extent of health inequalities between countries but also in small areas within a country. Identifying areas with high mortality is important both of public health mitigation in COVID-19 outbreaks, and of longer term efforts to tackle social inequalities in health. In this paper we consider different statistical models and an extension of a recent method to analyze COVID-19 related mortality in English small areas during the first wave of the epidemic in the first half of 2020. We seek to identify hotspots, and where they are most geographically concentrated, taking account of observed area factors as well as spatial correlation and clustering in regression residuals, while also allowing for spatial discontinuities. Results show an excess of COVID-19 mortality cases in small areas surrounding London and in other small areas in North-East and and North-West of England. Models alleviating spatial confounding show ethnic isolation, air quality and area morbidity covariates having a significant and broadly similar impact on COVID-19 mortality, whereas nursing home location seems to be slightly less important.

Spatial Statistical Models: An Overview under the Bayesian Approach

Spatial documentation is exponentially increasing given the availability of Big Data in the Internet of Things, enabled by device miniaturization and data storage capacity. Bayesian spatial statistics is a useful statistical tool to determine the dependence structure and hidden patterns in space through prior knowledge and data likelihood. However, this class of modeling is not yet well explored when compared to adopting classification and regression in machine-learning models, in which the assumption of the spatiotemporal independence of the data is often made, that is an inexistent or very weak dependence. Thus, this systematic review aims to address the main models presented in the literature over the past 20 years, identifying the gaps and research opportunities. Elements such as random fields, spatial domains, prior specification, the covariance function, and numerical approximations are discussed. This work explores the two subclasses of spatial smoothing: global and local.

Spatiotemporal Analysis of Overall Health in the United States Between 2010 and 2018

Background

Although many previous studies have documented spatial heterogeneity in health outcomes across the United States at different geographic scales, spatiotemporal analyses to understand overall health are scant.

Methodology

We used the County Health Rankings (CHR) data to analyze the three types of health outcomes, viz., overall health, length of life, and quality of life for 2010-2018 in the contiguous United States employing hierarchal Bayesian methods. Composite scores were created to proxy these outcomes utilizing predefined weights of several variables as recommended by CHR. Our methods assumed a convolution of spatially structured and unstructured errors to model the overall spatial error. Spatial effects were modeled using conditional autoregressive distribution.

Results

The substantial disparity in these health outcomes was evident, with counties having poorer health outcomes mostly concentrated in the southeastern United States. Models that incorporated county-level demographic and socioeconomic characteristics partially explained the observed spatial heterogeneity in health outcomes. Interestingly, there was no time effect in any of the outcomes suggesting a perpetuation of health disparity over the years.

Conclusions

County-specific health policy interventions that take into account the contextual factors might be beneficial in improving population health and breaking the perpetuation of health disparity.

From first to second wave: follow-up of the prospective COVID-19 cohort (KoCo19) in Munich (Germany)

BACKGROUND

In the 2nd year of the COVID-19 pandemic, knowledge about the dynamics of the infection in the general population is still limited. Such information is essential for health planners, as many of those infected show no or only mild symptoms and thus, escape the surveillance system. We therefore aimed to describe the course of the pandemic in the Munich general population living in private households from April 2020 to January 2021.

METHODS

The KoCo19 baseline study took place from April to June 2020 including 5313 participants (age 14 years and above). From November 2020 to January 2021, we could again measure SARS-CoV-2 antibody status in 4433 of the baseline participants (response 83%). Participants were offered a self-sampling kit to take a capillary blood sample (dry blood spot; DBS). Blood was analysed using the Elecsys® Anti-SARS-CoV-2 assay (Roche). Questionnaire information on socio-demographics and potential risk factors assessed at baseline was available for all participants. In addition, follow-up information on health-risk taking behaviour and number of personal contacts outside the household (N = 2768) as well as leisure time activities (N = 1263) were collected in summer 2020.

RESULTS

Weighted and adjusted (for specificity and sensitivity) SARS-CoV-2 sero-prevalence at follow-up was 3.6% (95% CI 2.9-4.3%) as compared to 1.8% (95% CI 1.3-3.4%) at baseline. 91% of those tested positive at baseline were also antibody-positive at follow-up. While sero-prevalence increased from early November 2020 to January 2021, no indication of geospatial clustering across the city of Munich was found, although cases clustered within households. Taking baseline result and time to follow-up into account, men and participants in the age group 20-34 years were at the highest risk of sero-positivity. In the sensitivity analyses, differences in health-risk taking behaviour, number of personal contacts and leisure time activities partly explained these differences.

CONCLUSION

The number of citizens in Munich with SARS-CoV-2 antibodies was still below 5% during the 2nd wave of the pandemic. Antibodies remained present in the majority of SARS-CoV-2 sero-positive baseline participants. Besides age and sex, potentially confounded by differences in behaviour, no major risk factors could be identified. Non-pharmaceutical public health measures are thus still important.

Bayesian spatiotemporal forecasting and mapping of COVID-19 risk with application to West Java Province, Indonesia

The coronavirus disease (COVID-19) has spread rapidly to multiple countries including Indonesia. Mapping its spatiotemporal pattern and forecasting (small area) outbreaks are crucial for containment and mitigation strategies. Hence, we introduce a parsimonious space-time model of new infections that yields accurate forecasts but only requires information regarding the number of incidences and population size per geographical unit and time period. Model parsimony is important because of limited knowledge regarding the causes of COVID-19 and the need for rapid action to control outbreaks. We outline the basics of Bayesian estimation, forecasting, and mapping, in particular for the identification of hotspots. The methodology is applied to county-level data of West Java Province, Indonesia.

Combining assumptions and graphical network into gene expression data analysis

Background

Analyzing gene expression data rigorously requires taking assumptions into consideration but also relies on using information about network relations that exist among genes. Combining these different elements cannot only improve statistical power, but also provide a better framework through which gene expression can be properly analyzed.

Material and methods

We propose a novel statistical model that combines assumptions and gene network information into the analysis. Assumptions are important since every test statistic is valid only when required assumptions hold. So, we propose hybrid p-values and show that, under the null hypothesis of primary interest, these p-values are uniformly distributed. These proposed hybrid p-values take assumptions into consideration. We incorporate gene network information into the analysis because neighboring genes share biological functions. This correlation factor is taken into account via similar prior probabilities for neighboring genes.

Results

With a series of simulations our approach is compared with other approaches. Area Under the ROC Curves (AUCs) are constructed to compare the different methodologies; the AUC based on our methodology is larger than others. For regression analysis, AUC from our proposed method contains AUCs of Spearman test and of Pearson test. In addition, true negative rates (TNRs) also known as specificities are higher with our approach than with the other approaches. For two group comparison analysis, for instance, with a sample size of n=10, specificity corresponding to our proposed methodology is 0.716146 and specificities for t-test and rank sum are 0.689223 and 0.69797, respectively. Our method that combines assumptions and network information into the analysis is shown to be more powerful.

Conclusions

These proposed procedures are introduced as a general class of methods that can incorporate procedure-selection, account for multiple-testing, and incorporate graphical network information into the analysis. We obtain very good performance in simulations, and in real data analysis.

Spatial and spatio-temporal methods for mapping malaria risk: a systematic review

Background

Approaches in malaria risk mapping continue to advance in scope with the advent of geostatistical techniques spanning both the spatial and temporal domains. A substantive review of the merits of the methods and covariates used to map malaria risk has not been undertaken. Therefore, this review aimed to systematically retrieve, summarise methods and examine covariates that have been used for mapping malaria risk in sub-Saharan Africa (SSA).

Methods

A systematic search of malaria risk mapping studies was conducted using PubMed, EBSCOhost, Web of Science and Scopus databases. The search was restricted to refereed studies published in English from January 1968 to April 2020. To ensure completeness, a manual search through the reference lists of selected studies was also undertaken. Two independent reviewers completed each of the review phases namely: identification of relevant studies based on the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines, data extraction and methodological quality assessment using a validated scoring criterion.

Results

One hundred and seven studies met the inclusion criteria. The median quality score across studies was 12/16 (range: 7–16). Approximately half (44%) of the studies employed variable selection techniques prior to mapping with rainfall and temperature selected in over 50% of the studies. Malaria incidence (47%) and prevalence (35%) were the most commonly mapped outcomes, with Bayesian geostatistical models often (31%) the preferred approach to risk mapping. Additionally, 29% of the studies employed various spatial clustering methods to explore the geographical variation of malaria patterns, with Kulldorf scan statistic being the most common. Model validation was specified in 53 (50%) studies, with partitioning data into training and validation sets being the common approach.

Conclusions

Our review highlights the methodological diversity prominent in malaria risk mapping across SSA. To ensure reproducibility and quality science, best practices and transparent approaches should be adopted when selecting the statistical framework and covariates for malaria risk mapping. Findings underscore the need to periodically assess methods and covariates used in malaria risk mapping; to accommodate changes in data availability, data quality and innovation in statistical methodology.

Association between environmental and climatic risk factors and the spatial distribution of cystic and alveolar echinococcosis in Kyrgyzstan

Background

Cystic and alveolar echinococcosis (CE and AE) are neglected tropical diseases caused by Echinococcus granulosus sensu lato and E. multilocularis, and are emerging zoonoses in Kyrgyzstan. In this country, the spatial distribution of CE and AE surgical incidence in 2014-2016 showed marked heterogeneity across communities, suggesting the presence of ecological determinants underlying CE and AE distributions.

Methodology/Principal findings

For this reason, in this study we assessed potential associations between community-level confirmed primary CE (no.=2359) or AE (no.=546) cases in 2014-2016 in Kyrgyzstan and environmental and climatic variables derived from satellite-remote sensing datasets using conditional autoregressive models. We also mapped CE and AE relative risk. The number of AE cases was negatively associated with 10-year lag mean annual temperature. Although this time lag should not be considered as an exact measurement but with associated uncertainty, it is consistent with the estimated 10–15-year latency following AE infection. No associations were detected for CE. We also identified several communities at risk for CE or AE where no disease cases were reported in the study period.

Conclusions/Significance

Our findings support the hypothesis that CE is linked to an anthropogenic cycle and is less affected by environmental risk factors compared to AE, which is believed to result from spillover from a wild life cycle. As CE was not affected by factors we investigated, hence control should not have a geographical focus. In contrast, AE risk areas identified in this study without reported AE cases should be targeted for active disease surveillance in humans. This active surveillance would confirm or exclude AE transmission which might not be reported with the present passive surveillance system. These areas should also be targeted for ecological investigations in the animal hosts.

Cystic and alveolar echinococcosis (CE and AE) are parasitic zoonoses that cause a substantial disease burden in Kyrgyzstan. The etiologic agents of these diseases are parasites in the genus Echinococcus. These parasites have complex life cycles which include mammalian definitive and intermediate hosts and a free-living egg stage in the environment. Consequently, environmental and climatic factors can affect the prevalence and geographical distribution of these diseases because such factors influence the parasites’ eggs survival and longevity, and can affect suitable habitats for the intermediate and definitive hosts. In this geographic correlation study, we assessed environmental and climatic determinants of the spatial distributions of CE and AE in Kyrgyzstan. We found that 10-year lag annual temperature plays an important role in AE distribution, whilst none of the variables assessed was found to significantly affect that of CE. Moreover, communities at risk where these diseases are potentially under- or misdiagnosed were identified. Our findings provide vital information for targeted, area-specific interventions in Kyrgyzstan, and add to the body of knowledge on the ecology of these neglected parasitic diseases that are emerging and reemerging in several regions in North America, Europe and Asia.

Community factors associated with local epidemic timing of respiratory syncytial virus: A spatiotemporal modeling study

Respiratory syncytial virus (RSV) causes a large burden of morbidity in young children and the elderly. Spatial variability in the timing of RSV epidemics provides an opportunity to probe the factors driving its transmission, including factors that influence epidemic seeding and growth rates. Using hospitalization data from Connecticut, New Jersey, and New York, we estimated epidemic timing at the ZIP code level using harmonic regression and then used a Bayesian meta-regression model to evaluate correlates of epidemic timing. Earlier epidemics were associated with larger household size and greater population density. Nearby localities had similar epidemic timing. Our results suggest that RSV epidemics grow faster in areas with more local contact opportunities, and that epidemic spread follows a spatial diffusion process based on geographic proximity. Our findings can inform the timing of delivery of RSV extended half-life prophylaxis and maternal vaccines and guide future studies on the transmission dynamics of RSV.

Modeling the Social and Spatial Proximity of Crime: Domestic and Sexual Violence Across Neighborhoods

OBJECTIVES

Our goal is to understand the social dynamics affecting domestic and sexual violence in urban areas by investigating the role of connections between area nodes, or communities. We use innovative methods adapted from spatial statistics to investigate the importance of social proximity measured based on connectedness pathways between area nodes. In doing so, we seek to extend the standard treatment in the neighborhoods and crime literature of areas like census blocks as independent analytical units or as interdependent primarily due to geographic proximity.

METHODS

In this paper, we develop techniques to incorporate two types of proximity, geographic proximity and commuting proximity in spatial generalized linear mixed models (SGLMM) in order to estimate domestic and sexual violence in Detroit, Michigan and Arlington County, Virginia. Analyses are based on three types of CAR models (the Besag, York, and Mollié (BYM), Leroux, and the sparse SGLMM models) and two types of SAR models (the spatial lag and spatial error models) to examine how results vary with different model assumptions. We use data from local and federal sources such as the Police Data Initiative and American Community Survey.

RESULTS

Analyses show that incorporating information on commuting ties, a non-spatially bounded form of social proximity, to spatial models contributes to better deviance information criteria (DIC) scores (a metric which explicitly accounts for model fit and complexity) in Arlington for sexual and domestic crime as well as overall crime. In Detroit, the fit is improved only for overall crime. The distinctions in model fit are less pronounced when using cross-validated mean absolute error (MAE) as a comparison criteria.

CONCLUSION

Overall, the results indicate variations across crime type, urban contexts, and modeling approaches. Nonetheless, in important contexts, commuting ties among neighborhoods are observed to greatly improve our understanding of urban crime. If such ties contribute to the transfer of norms, social support, resources, and behaviors between places, they may then transfer also the effects of crime prevention efforts.

Software application profile: the Rapid Inquiry Facility 4.0: an open access tool for environmental public health tracking

The Rapid Inquiry Facility 4.0 (RIF) is a new user-friendly and open-access tool, developed by the UK Small Area Health Statistics Unit (SAHSU), to facilitate environment public health tracking (EPHT) or surveillance (EPHS). The RIF is designed to help public health professionals and academics to rapidly perform exploratory investigations of health and environmental data at the small-area level (e.g. postcode or detailed census areas) in order to identify unusual signals, such as disease clusters and potential environmental hazards, whether localized (e.g. industrial site) or widespread (e.g. air and noise pollution). The RIF allows the use of advanced disease mapping methods, including Bayesian small-area smoothing and complex risk analysis functionalities, while accounting for confounders. The RIF could be particularly useful to monitor spatio-temporal trends in mortality and morbidity associated with cardiovascular diseases, cancers, diabetes and chronic lung diseases, or to conduct local or national studies on air pollution, flooding, low-magnetic fields or nuclear power plants.

District-level estimation of vaccination coverage: Discrete vs continuous spatial models

Health and development indicators (HDIs) such as vaccination coverage are regularly measured in many low- and middle-income countries using household surveys, often due to the unreliability or incompleteness of routine data collection systems. Recently, the development of model-based approaches for producing subnational estimates of HDIs using survey data, particularly cluster-level data, has been an active area of research. This is mostly driven by the increasing demand for estimates at certain administrative levels, for example, districts, at which many development goals are set and evaluated. In this study, we explore spatial modeling approaches for producing district-level estimates of vaccination coverage. Specifically, we compare discrete spatial smoothing models which directly model district-level data with continuous Gaussian process (GP) models that utilize geolocated cluster-level data. We adopt a fully Bayesian framework, implemented using the INLA and SPDE approaches. We compare the predictive performance of the models by analyzing vaccination coverage using data from two Demographic and Health Surveys (DHS), namely the 2014 Kenya DHS and the 2015-16 Malawi DHS. We find that the continuous GP models performed well, offering a credible alternative to traditional discrete spatial smoothing models. Our analysis also revealed that accounting for between-cluster variation in the continuous GP models did not have any real effect on the district-level estimates. Our results provide guidance to practitioners on the reliability of these model-based approaches for producing estimates of vaccination coverage and other HDIs.

Origin Sample Prediction and Spatial Modeling of Antimicrobial Resistance in Metagenomic Sequencing Data

The steady elaboration of the Metagenomic and Metadesign of Subways and Urban Biomes (MetaSUB) international consortium project raises important new questions about the origin, variation, and antimicrobial resistance of the collected samples. CAMDA (Critical Assessment of Massive Data Analysis, http://camda.info/) forum organizes annual challenges where different bioinformatics and statistical approaches are tested on samples collected around the world for bacterial classification and prediction of geographical origin. This work proposes a method which not only predicts the locations of unknown samples, but also estimates the relative risk of antimicrobial resistance through spatial modeling. We introduce a new component in the standard analysis as we apply a Bayesian spatial convolution model which accounts for spatial structure of the data as defined by the longitude and latitude of the samples and assess the relative risk of antimicrobial resistance taxa across regions which is relevant to public health. We can then use the estimated relative risk as a new measure for antimicrobial resistance. We also compare the performance of several machine learning methods, such as Gradient Boosting Machine, Random Forest, and Neural Network to predict the geographical origin of the mystery samples. All three methods show consistent results with some superiority of Random Forest classifier. In our future work we can consider a broader class of spatial models and incorporate covariates related to the environment and climate profiles of the samples to achieve more reliable estimation of the relative risk related to antimicrobial resistance.

Population-weighted exposure to air pollution and COVID-19 incidence in Germany

Many countries have enforced social distancing to stop the spread of COVID-19. Within countries, although the measures taken by governments are similar, the incidence rate varies among areas (e.g., counties, cities). One potential explanation is that people in some areas are more vulnerable to the coronavirus disease because of their worsened health conditions caused by long-term exposure to poor air quality. In this study, we investigate whether long-term exposure to air pollution increases the risk of COVID-19 infection in Germany. The results show that nitrogen dioxide (NO 2 ) is significantly associated with COVID-19 incidence, with a 1 μ g  m - 3 increase in long-term exposure to NO 2 increasing the COVID-19 incidence rate by 5.58% (95% credible interval [CI]: 3.35%, 7.86%). This result is consistent across various models. The analyses can be reproduced and updated routinely using public data sources and shared R code.

Spatial modeling of individual-level infectious disease transmission: Tuberculosis data in Manitoba, Canada

Geographically dependent individual level models (GD-ILMs) are a class of statistical models that can be used to study the spread of infectious disease through a population in discrete-time in which covariates can be measured both at individual and area levels. The typical ILMs to illustrate spatial data are based on the distance between susceptible and infectious individuals. A key feature of GD-ILMs is that they take into account the spatial location of the individuals in addition to the distance between susceptible and infectious individuals. As a motivation of this article, we consider tuberculosis (TB) data which is an infectious disease which can be transmitted through individuals. It is also known that certain areas/demographics/communities have higher prevalent of TB (see Section 4 for more details). It is also of interest of policy makers to identify those areas with higher infectivity rate of TB for possible preventions. Therefore, we need to analyze this data properly to address those concerns. In this article, the expectation conditional maximization algorithm is proposed for estimating the parameters of GD-ILMs to be able to predict the areas with the highest average infectivity rates of TB. We also evaluate the performance of our proposed approach through some simulations. Our simulation results indicate that the proposed method provides reliable estimates of parameters which confirms accuracy of the infectivity rates.

A hierarchical spatio-temporal model to analyze relative risk variations of COVID-19: a focus on Spain, Italy and Germany

The novel coronavirus disease (COVID-19) has spread rapidly across the world in a short period of time and with a heterogeneous pattern. Understanding the underlying temporal and spatial dynamics in the spread of COVID-19 can result in informed and timely public health policies. In this paper, we use a spatio-temporal stochastic model to explain the temporal and spatial variations in the daily number of new confirmed cases in Spain, Italy and Germany from late February 2020 to mid January 2021. Using a hierarchical Bayesian framework, we found that the temporal trends of the epidemic in the three countries rapidly reached their peaks and slowly started to decline at the beginning of April and then increased and reached their second maximum in the middle of November. However decline and increase of the temporal trend seems to show different patterns in Spain, Italy and Germany.

Methods of spatial cluster detection in rare childhood cancers: Benchmarking data and results from a simulation study on nephroblastoma

The potential existence of spatial clusters in childhood cancer incidence is a debated topic. Identification of rare disease clusters in general may help to better understand disease etiology and develop preventive strategies against such entities.

The incidence of newly diagnosed childhood malignancies under 15 years of age is 140/1,000,000. In this context, the subgroup of nephroblastoma represents an extremely rare entity with an annual incidence of 7/1,000,000. We evaluated widely used statistical approaches for spatial cluster detection in childhood cancer (Ref. Schündeln et al., 2021, Cancer Epidemiology). For the simulation study, random high risk clusters of 1 to 50 adjacent districts (NUTS-level 3, nomenclature des unités territoriales statistiques) were generated on the basis of the 402 German administrative districts. Each cluster was simulated with different relative risk levels (1 to 100). For each combination of cluster size and risk level 2000 iterations were performed. Simulated data was then analyzed by three local clustering tests: Besag-Newell method, spatial scan statistic and the Bayesian Besag-York-Mollié approach (fit by Integrated Nested Laplace Approximation). The performance characteristics of all three methods were systematically documented (sensitivity, specificity, positive/negative predictive values, exact- and minimum power, correct classification, positive/negative diagnostic likelihood and false positive/negative rate).

This data article links to a Mendeley online repository which includes the raw data of simulated high-risk clusters and simulated cases on the district level for an all-childhood-malignancy scenario as well as for cases of nephroblastoma. These data was used for the evaluation of the three cluster detection methods. The R code for simulation and analysis are available from GitHub.

The article also includes analyzed data summarizing the performance of the cluster detection tests in very rare disease entities, using the example of simulated nephroblastoma cases.

The raw data from the study can be used for benchmarking analyses applying different spatial statistical methods systematically and evaluating their performance characteristics comparatively. The analyzed data from the nephroblastoma example can be useful to interpret the performance of the three applied local cluster detection tests in the setting of extremely rare disease entities. As a practical application, data and R code can be used for performance analyses when planning to establish surveillance systems for rare disease entities.

Statistical methods for spatial cluster detection in childhood cancer incidence: A simulation study

BACKGROUND AND OBJECTIVE

The potential existence of spatial clusters in childhood cancer incidence is a debated topic. Identification of such clusters may help to better understand etiology and develop preventive strategies. We evaluated widely used statistical approaches to cluster detection in this context.

METHODS

Incidence of newly diagnosed childhood cancer (140/1,000,000 children under 15 years) and nephroblastoma (7/1,000,000) was simulated. Clusters of defined size (1-50) were randomly assembled on the district level in Germany. Each cluster was simulated with different relative risk levels (1-100). For each combination 2000 iterations were done. Simulated data was then analyzed by three local clustering tests: Besag-Newell method, spatial scan statistic and Bayesian Besag-York-Mollié with Integrated Nested Laplace Approximation approach. The operating characteristics (sensitivity, specificity, predictive values, power and correct classification) of all three methods were systematically described.

RESULTS

Performance varied considerably within and between methods, depending on the simulated setting. Sensitivity of all methods was positively associated with increasing size, incidence and RR of the high-risk area. Besag-York-Mollié showed highest specificity for minimally increased RR in most scenarios. The performance of all methods was lower in the nephroblastoma scenario compared with the scenario including all cancer cases.

CONCLUSION

This study illustrates the challenge to make reliable inferences on the existence of spatial clusters based on single statistical approaches in childhood cancer. Application of multiple methods, ideally with known operating characteristics, and a critical discussion of the joint evidence seems recommendable when aiming to identify high-risk clusters.

Spatial-temporal patterns of malaria incidence in Uganda using HMIS data from 2015 to 2019

BACKGROUND

As global progress to reduce malaria transmission continues, it is increasingly important to track changes in malaria incidence rather than prevalence. Risk estimates for Africa have largely underutilized available health management information systems (HMIS) data to monitor trends. This study uses national HMIS data, together with environmental and geographical data, to assess spatial-temporal patterns of malaria incidence at facility catchment level in Uganda, over a recent 5-year period.

METHODS

Data reported by 3446 health facilities in Uganda, between July 2015 and September 2019, was analysed. To assess the geographic accessibility of the health facilities network, AccessMod was employed to determine a three-hour cost-distance catchment around each facility. Using confirmed malaria cases and total catchment population by facility, an ecological Bayesian conditional autoregressive spatial-temporal Poisson model was fitted to generate monthly posterior incidence rate estimates, adjusted for caregiver education, rainfall, land surface temperature, night-time light (an indicator of urbanicity), and vegetation index.

RESULTS

An estimated 38.8 million (95% Credible Interval [CI]: 37.9-40.9) confirmed cases of malaria occurred over the period, with a national mean monthly incidence rate of 20.4 (95% CI: 19.9-21.5) cases per 1000, ranging from 8.9 (95% CI: 8.7-9.4) to 36.6 (95% CI: 35.7-38.5) across the study period. Strong seasonality was observed, with June-July experiencing highest peaks and February-March the lowest peaks. There was also considerable geographic heterogeneity in incidence, with health facility catchment relative risk during peak transmission months ranging from 0 to 50.5 (95% CI: 49.0-50.8) times higher than national average. Both districts and health facility catchments showed significant positive spatial autocorrelation; health facility catchments had global Moran's I = 0.3 (p < 0.001) and districts Moran's I = 0.4 (p < 0.001). Notably, significant clusters of high-risk health facility catchments were concentrated in Acholi, West Nile, Karamoja, and East Central - Busoga regions.

CONCLUSION

Findings showed clear countrywide spatial-temporal patterns with clustering of malaria risk across districts and health facility catchments within high risk regions, which can facilitate targeting of interventions to those areas at highest risk. Moreover, despite high and perennial transmission, seasonality for malaria incidence highlights the potential for optimal and timely implementation of targeted interventions.

Towards activity-based exposure measures in spatial analysis of pedestrian-motor vehicle crashes

BACKGROUND

Although numerous efforts have been devoted to exploring the effects of area-wide factors on the frequency of pedestrian crashes in neighborhoods over the past two decades, existing studies have largely failed to provide a full picture of the factors that contribute to the incidence of zonal pedestrian crashes, due to the unavailability of reliable exposure data and use of less sound analytical methods.

METHODS

Based on a crowdsourced dataset in Hong Kong, we first proposed a procedure to extract pedestrian trajectories from travel-diary survey data. We then aggregated these data to 209 neighborhoods and developed a Bayesian spatially varying coefficients model to investigate the spatially non-stationary relationships between the number of pedestrian-motor vehicle (PMV) crashes and related risk factors. To dissect the role of pedestrian exposure, the estimated coefficients of models with population, walking trips, walking time, and walking distance as the measure of pedestrian exposure were presented and compared.

RESULTS

Our results indicated substantial inconsistencies in the effects of several risk factors between the models of population and activity-based exposure measures. The model using walking trips as the measure of pedestrian exposure had the best goodness-of-fit. We also provided new insights that in addition to the unstructured variability, heterogeneity in the effects of explanatory variables on the frequency of PMV crashes could also arise from the spatially correlated effects. After adjusting for vehicle volume and pedestrian activity, road density, intersection density, bus stop density, and the number of parking lots were found to be positively associated with PMV crash frequency, whereas the percentage of motorways and median monthly income had negative associations with the risk of PMV crashes.

CONCLUSIONS

The use of population or population density as a surrogate for pedestrian exposure when modeling the frequency of zonal pedestrian crashes is expected to produce biased estimations and invalid inferences. Spatial heterogeneity should also not be negligible when modeling pedestrian crashes involving contiguous spatial units.

A conceptual framework to understand the role of built environment on traffic safety

INTRODUCTION

Many U.S. cities have adopted the Vision Zero strategy with the specific goal of eliminating traffic-related deaths and injuries. To achieve this ambitious goal, safety professionals have increasingly called for the development of a safe systems approach to traffic safety. This approach calls for examining the macrolevel risk factors that may lead road users to engage in errors that result in crashes. This study explores the relationship between built environment variables and crash frequency, paying specific attention to the environmental mediating factors, such as traffic exposure, traffic conflicts, and network-level speed characteristics.

METHODS

Three years (2011-2013) of crash data from Mecklenburg County, North Carolina, were used to model crash frequency on surface streets as a function of built environment variables at the census block group level. Separate models were developed for total and KAB crashes (i.e., crashes resulting in fatalities (K), incapacitating injuries (A), or non-incapacitating injuries (B)) using the conditional autoregressive modeling approach to account for unobserved heterogeneity and spatial autocorrelation present in data.

RESULTS

Built environment variables that are found to have positive associations with both total and KAB crash frequencies include population, vehicle miles traveled, big box stores, intersections, and bus stops. On the other hand, the number of total and KAB crashes tend to be lower in census block groups with a higher proportion of two-lane roads and a higher proportion of roads with posted speed limits of 35 mph or less.

CONCLUSIONS

This study demonstrates the plausible mechanism of how the built environment influences traffic safety. The variables found to be significant are all policy-relevant variables that can be manipulated to improve traffic safety. Practical Applications: The study findings will shape transportation planning and policy level decisions in designing the built environment for safer travels.