Spatio-temporal disease mapping comprises a wide range of models used to describe the distribution of a disease in space and its evolution in time. These models have been commonly formulated within a hierarchical Bayesian framework with two main approaches: an empirical Bayes (EB) and a fully Bayes (FB) approach. The EB approach provides point estimates of the parameters relying on the well-known penalized quasi-likelihood (PQL) technique. The FB approach provides the posterior distribution of the target parameters. These marginal distributions are not usually available in closed form and common estimation procedures are based on Markov chain Monte Carlo (MCMC) methods. However, the spatio-temporal models used in disease mapping are often very complex and MCMC methods may lead to large Monte Carlo errors and a huge computation time if the dimension of the data at hand is large. To circumvent these potential inconveniences, a new technique called integrated nested Laplace approximations (INLA), based on nested Laplace approximations, has been proposed for Bayesian inference in latent Gaussian models. In this paper, we show how to fit different spatio-temporal models for disease mapping with INLA using the Leroux CAR prior for the spatial component, and we compare it with PQL via a simulation study. The spatio-temporal distribution of male brain cancer mortality in Spain during the period 1986-2010 is also analysed.

The purpose of this article is to draw attention to the possible need for inclusion of interaction effects between regions and age groups in mapping studies. We propose a simple model for including such an interaction in order to develop a test for its significance. The assumption of an absence of such interaction effects is a helpful simplifying one. The measure of relative risk related to a particular region becomes easily and neatly summarized. Indeed, such a test seems warranted because it is anticipated that the simple model, which ignores such interaction, as is in common use, may at times be adequate. The test proposed is a score test and hence only requires fitting the simpler model. We illustrate our approaches using mortality data from British Columbia, Canada, over the 5-year period 1985-1989. For this data, the interaction effect between age groups and regions is quite large and significant.

SUMMARYAnalyzing the temporal evolution of the geographical distribution of mortality (or incidence) risks is an important area of research in disease mapping. It might help to better determining the risk factors involved in the studied disease, and to address important epidemiological questions about the stability of the estimated patterns of disease. Traditionally, risk smoothing is carried out using conditional autoregressive (CAR) models but very recently, penalized splines have also been considered in an Empirical Bayes (EB) spatial context to estimate large-scale spatial trends together with region random effects. In this paper, penalized splines for smoothing risks in both the spatial and the temporal dimensions will be applied. The mean squared error (MSE) of the log-risk predictor will be derived allowing for constructing confidence intervals for the risks. To illustrate the procedure mortality data due to brain cancer in continental Spain over the period 1996-2005 are analyzed.

Detecting change-points and trends are common tasks in the analysis of remote sensing data. Over the years, many different methods have been proposed for those purposes, including (modified) Mann–Kendall and Cox–Stuart tests for detecting trends; and Pettitt, Buishand range, Buishand U, standard normal homogeneity (Snh), Meanvar, structure change (Strucchange), breaks for additive season and trend (BFAST), and hierarchical divisive (E.divisive) for detecting change-points. In this paper, we describe a simulation study based on including different artificial, abrupt changes at different time-periods of image time series to assess the performances of such methods. The power of the test, type I error probability, and mean absolute error (MAE) were used as performance criteria, although MAE was only calculated for change-point detection methods. The study reveals that if the magnitude of change (or trend slope) is high, and/or the change does not occur in the first or last time-periods, the methods generally have a high power and a low MAE. However, in the presence of temporal autocorrelation, MAE raises, and the probability of introducing false positives increases noticeably. The modified versions of the Mann–Kendall method for autocorrelated data reduce/moderate its type I error probability, but this reduction comes with an important power diminution. In conclusion, taking a trade-off between the power of the test and type I error probability, we conclude that the original Mann–Kendall test is generally the preferable choice. Although Mann–Kendall is not able to identify the time-period of abrupt changes, it is more reliable than other methods when detecting the existence of such changes. Finally, we look for trend/change-points in land surface temperature (LST), day and night, via monthly MODIS images in Navarre, Spain, from January 2001 to December 2018.