Age-space-time CAR models in Bayesian disease mapping

Using mortality to predict incidence for rare and lethal cancers in very small areas

Incidence and mortality figures are needed to get a comprehensive overview of cancer burden. In many countries, cancer mortality figures are routinely recorded by statistical offices, whereas incidence depends on regional cancer registries. However, due to the complexity of updating cancer registries, incidence numbers become available 3 or 4 years later than mortality figures. It is, therefore, necessary to develop reliable procedures to predict cancer incidence at least until the period when mortality data are available. Most of the methods proposed in the literature are designed to predict total cancer (except nonmelanoma skin cancer) or major cancer sites. However, less frequent lethal cancers, such as brain cancer, are generally excluded from predictions because the scarce number of cases makes it difficult to use univariate models. Our proposal comes to fill this gap and consists of modeling jointly incidence and mortality data using spatio-temporal models with spatial and age shared components. This approach allows for predicting lethal cancers improving the performance of individual models when data are scarce by taking advantage of the high correlation between incidence and mortality. A fully Bayesian approach based on integrated nested Laplace approximations is considered for model fitting and inference. A validation process is also conducted to assess the performance of alternative models. We use the new proposals to predict brain cancer incidence rates by gender and age groups in the health units of Navarre and Basque Country (Spain) during the period 2005-2008.

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.

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.

Beyond standardized mortality ratios; some uses of smoothed age-specific mortality rates on small areas studies

Background

Most epidemiological risk indicators strongly depend on the age composition of populations, which makes the direct comparison of raw (unstandardized) indicators misleading because of the different age structures of the spatial units of study. Age-standardized rates (ASR) are a common solution for overcoming this confusing effect. The main drawback of ASRs is that they depend on age-specific rates which, when working with small areas, are often based on very few, or no, observed cases for most age groups. A similar effect occurs with life expectancy at birth and many more epidemiological indicators, which makes standardized mortality ratios (SMR) the omnipresent risk indicator for small areas epidemiologic studies.

Methods

To deal with this issue, a multivariate smoothing model, the M-model, is proposed in order to fit the age-specific probabilities of death (PoDs) for each spatial unit, which assumes dependence between closer age groups and spatial units. This age–space dependence structure enables information to be transferred between neighboring consecutive age groups and neighboring areas, at the same time, providing more reliable age-specific PoDs estimates.

Results

Three case studies are presented to illustrate the wide range of applications that smoothed age specific PoDs have in practice . The first case study shows the application of the model to a geographical study of lung cancer mortality in women. This study illustrates the convenience of considering age–space interactions in geographical studies and to explore the different spatial risk patterns shown by the different age groups. Second, the model is also applied to the study of ischaemic heart disease mortality in women in two cities at the census tract level. Smoothed age-standardized rates are derived and compared for the census tracts of both cities, illustrating some advantages of this mortality indicator over traditional SMRs. In the latest case study, the model is applied to estimate smoothed life expectancy (LE), which is the most widely used synthetic indicator for characterizing overall mortality differences when (not so small) spatial units are considered.

Conclusion

Our age–space model is an appropriate and flexible proposal that provides more reliable estimates of the probabilities of death, which allow the calculation of enhanced epidemiological indicators (smoothed ASR, smoothed LE), thus providing alternatives to traditional SMR-based studies of small areas.

A systematic review of Bayesian spatial-temporal models on cancer incidence and mortality

OBJECTIVES

This study aimed to review the types and applications of fully Bayesian (FB) spatial-temporal models and covariates used to study cancer incidence and mortality.

METHODS

This systematic review searched articles published within Medline, Embase, Web-of-Science and Google Scholar between 2014 and 2018.

RESULTS

A total of 38 studies were included in our study. All studies applied Bayesian spatial-temporal models to explore spatial patterns over time, and over half assessed the association with risk factors. Studies used different modelling approaches and prior distributions for spatial, temporal and spatial-temporal interaction effects depending on the nature of data, outcomes and applications. The most common Bayesian spatial-temporal model was a generalized linear mixed model. These models adjusted for covariates at the patient, area or temporal level, and through standardization.

CONCLUSIONS

Few studies (4) modelled patient-level clinical characteristics (11%), and the applications of an FB approach in the forecasting of spatial-temporally aligned cancer data were limited. This review highlighted the need for Bayesian spatial-temporal models to incorporate patient-level prognostic characteristics through the multi-level framework and forecast future cancer incidence and outcomes for cancer prevention and control strategies.

Advances in spatiotemporal models for non-communicable disease surveillance

Surveillance systems are commonly used to provide early warning detection or to assess an impact of an intervention/policy. Traditionally, the methodological and conceptual frameworks for surveillance have been designed for infectious diseases, but the rising burden of non-communicable diseases (NCDs) worldwide suggests a pressing need for surveillance strategies to detect unusual patterns in the data and to help unveil important risk factors in this setting. Surveillance methods need to be able to detect meaningful departures from expectation and exploit dependencies within such data to produce unbiased estimates of risk as well as future forecasts. This has led to the increasing development of a range of space-time methods specifically designed for NCD surveillance. We present an overview of recent advances in spatiotemporal disease surveillance for NCDs, using hierarchically specified models. This provides a coherent framework for modelling complex data structures, dealing with data sparsity, exploiting dependencies between data sources and propagating the inherent uncertainties present in both the data and the modelling process. We then focus on three commonly used models within the Bayesian Hierarchical Model (BHM) framework and, through a simulation study, we compare their performance. We also discuss some challenges faced by researchers when dealing with NCD surveillance, including how to account for false detection and the modifiable areal unit problem. Finally, we consider how to use and interpret the complex models, how model selection may vary depending on the intended user group and how best to communicate results to stakeholders and the general public.

A Bayesian spatiotemporal statistical analysis of out-of-hospital cardiac arrests

We propose a Bayesian spatiotemporal statistical model for predicting out-of-hospital cardiac arrests (OHCAs). Risk maps for Ticino, adjusted for demographic covariates, are built for explaining and forecasting the spatial distribution of OHCAs and their temporal dynamics. The occurrence intensity of the OHCA event in each area of interest, and the cardiac risk-based clustering of municipalities are efficiently estimated, through a statistical model that decomposes OHCA intensity into overall intensity, demographic fixed effects, spatially structured and unstructured random effects, time polynomial dependence, and spatiotemporal random effect. In the studied geography, time evolution and dependence on demographic features are robust over different categories of OHCAs, but with variability in their spatial and spatiotemporal structure. Two main OHCA incidence-based clusters of municipalities are identified.

Online relative risks/rates estimation in spatial and spatio-temporal disease mapping

BACKGROUND AND OBJECTIVE

Spatial and spatio-temporal analyses of count data are crucial in epidemiology and other fields to unveil spatial and spatio-temporal patterns of incidence and/or mortality risks. However, fitting spatial and spatio-temporal models is not easy for non-expert users. The objective of this paper is to present an interactive and user-friendly web application (named SSTCDapp) for the analysis of spatial and spatio-temporal mortality or incidence data. Although SSTCDapp is simple to use, the underlying statistical theory is well founded and all key issues such as model identifiability, model selection, and several spatial priors and hyperpriors for sensitivity analyses are properly addressed.

METHODS

The web application is designed to fit an extensive range of fairly complex spatio-temporal models to smooth the very often extremely variable standardized incidence/mortality risks or crude rates. The application is built with the R package shiny and relies on the well founded integrated nested Laplace approximation technique for model fitting and inference.

RESULTS

The use of the web application is shown through the analysis of Spanish spatio-temporal breast cancer data. Different possibilities for the analysis regarding the type of model, model selection criteria, and a range of graphical as well as numerical outputs are provided.

CONCLUSIONS

Unlike other software used in disease mapping, SSTCDapp facilitates the fit of complex statistical models to non-experts users without the need of installing any software in their own computers, since all the analyses and computations are made in a powerful remote server. In addition, a desktop version is also available to run the application locally in those cases in which data confidentiality is a serious issue.

Flexible Bayesian P-splines for smoothing age-specific spatio-temporal mortality patterns

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.