Large-scale modelling and forecasting of ambulance calls in northern Sweden using spatio-temporal log-Gaussian Cox processes

Abstract: Although ambulance call data typically come in the form of spatio-temporal point patterns, point process-based modelling approaches presented in the literature are scarce. In this paper, we study a unique set of Swedish spatio-temporal ambulance call data, which consist of the spatial (GPS) locations of the calls (within the four northernmost regions of Sweden) and the associated days of occurrence of the calls (January 1, 2014, to December 31, 2018). Motivated by the nature of the data, we here employ log-Gaussian Cox processes (LGCPs) for the spatio-temporal modelling and forecasting of the calls. To this end, we propose a K-means clustering based bandwidth selection method for the kernel estimation of the spatial component of the separable spatio-temporal intensity function. The temporal component of the intensity function is modelled using Poisson regression, using different calendar covariates, and the spatio-temporal random field component of the random intensity of the LGCP is fitted using the Metropolis-adjusted Langevin algorithm. Spatial hot-spots have been found in the south-eastern part of the study region, where most people in the region live and our fitted model/forecasts manage to capture this behavior quite well. Also, there is a significant association between the expected number of calls and the day-of-the-week and the season-of-the-year. A non-parametric second-order analysis indicates that LGCPs seem to be reasonable models for the data. Finally, we find that the fitted forecasts generate simulated future spatial event patterns that quite well resemble the actual future data.