MCMC methods for putative pollution source problems in environmental epidemiology

Responses to COVID-19 with probabilistic programming

The COVID-19 pandemic left its unique mark on the twenty-first century as one of the most significant disasters in history, triggering governments all over the world to respond with a wide range of interventions. However, these restrictions come with a substantial price tag. It is crucial for governments to form anti-virus strategies that balance the trade-off between protecting public health and minimizing the economic cost. This work proposes a probabilistic programming method to quantify the efficiency of major initial non-pharmaceutical interventions. We present a generative simulation model that accounts for the economic and human capital cost of adopting such strategies, and provide an end-to-end pipeline to simulate the virus spread and the incurred loss of various policy combinations. By investigating the national response in 10 countries covering four continents, we found that social distancing coupled with contact tracing is the most successful policy, reducing the virus transmission rate by 96% along with a 98% reduction in economic and human capital loss. Together with experimental results, we open-sourced a framework to test the efficacy of each policy combination.

A latent model to detect multiple clusters of varying sizes

This article develops a latent model and likelihood-based inference to detect temporal clustering of events. The model mimics typical processes generating the observed data. We apply model selection techniques to determine the number of clusters, and develop likelihood inference and a Monte Carlo expectation-maximization algorithm to estimate model parameters, detect clusters, and identify cluster locations. Our method differs from the classical scan statistic in that we can simultaneously detect multiple clusters of varying sizes. We illustrate the methodology with two real data applications and evaluate its efficiency through simulation studies. For the typical data-generating process, our methodology is more efficient than a competing procedure that relies on least squares.

An autoregressive point source model for spatial processes

We suggest a parametric modeling approach for nonstationary spatial processes driven by point sources. Baseline near-stationarity, which may be reasonable in the absence of a point source, is modeled using a conditional autoregressive (CAR) Markov random field. Variability due to the point source is captured by our proposed autoregressive point source (ARPS) model. Inference proceeds according to the Bayesian hierarchical paradigm, and is implemented using Markov chain Monte Carlo (MCMC) methods. The parametric approach allows a formal test of effectiveness of the point source. Application is made to a real dataset on electric potential measurements in a field containing a metal pole and the finding is that our approach captures the pole's impact on small-scale variability of the electric potential process.

Bayesian statistics in medicine: a 25 year review

This review examines the state of Bayesian thinking as Statistics in Medicine was launched in 1982, reflecting particularly on its applicability and uses in medical research. It then looks at each subsequent five-year epoch, with a focus on papers appearing in Statistics in Medicine, putting these in the context of major developments in Bayesian thinking and computation with reference to important books, landmark meetings and seminal papers. It charts the growth of Bayesian statistics as it is applied to medicine and makes predictions for the future. From sparse beginnings, where Bayesian statistics was barely mentioned, Bayesian statistics has now permeated all the major areas of medical statistics, including clinical trials, epidemiology, meta-analyses and evidence synthesis, spatial modelling, longitudinal modelling, survival modelling, molecular genetics and decision-making in respect of new technologies.

Disease cluster detection: a critique and a Bayesian proposal

This paper reviews issues in the analysis of non-focussed clustering, and proposes a novel approach to cluster modelling that can be used in a surveillance context. The novel approach involves the use of local likelihood models for the analysis of clustering in small area health data. Local likelihood is used when interdependence between data events at locations is modelled directly, as opposed to the modelling of a hidden process of cluster centres. This approach allows the use of conventional posterior sampling. It also allows a less parameterized approach to the form of clusters detected. The idea of a spatially dependent lasso which provides the local maxima for the aggregation of locations is considered as an approximation. The methods are applied to a well known data set and compared with Satscan, and a conditional logistic Bayesian model.

L’analyse de cluster en épidémiologie géographique : utilisation de plusieurs méthodes statistiques et comparaison de leurs résultats

BACKGROUND

The increasing interest in environmental epidemiology has been followed by the development of many statistical tests for detecting disease clustering near a point source. The objectives of this study were to compare several tests to detect disease clustering, among which modelisation using Markov Chain Monte Carlo methods.

METHODS

We compared six statistical methods for detecting disease clustering of bladder cancer around an industrial centre of Isère (France) for the period 1983-1997: Stone's test, score test, and two log-linear modelisations (with and without corrections for extra-Poisson variations) using two ways of parameters estimation (maximum likelihood and Markov Chain Monte Carlo methods).

RESULTS

The results of the Stone test and the score test are not in favour of a higher risk of bladder cancer around the considered point source. The conclusions brought by the log linear modelisations are the same, but the results obtained using the Markov Chain Monte Carlo Method are very dependant of prior distributions determined for the different parameters.

CONCLUSION

Markov Chain Monte Carlo methods, which allow taking into account complex geographical effects, seem well adapted to cluster analysis in geographical epidemiology. However, they remain difficult to implement.

A hierarchical model for spatially clustered disease rates

Maps of regional disease rates are potentially useful tools in examining spatial patterns of disease and for identifying clusters. Bayes and empirical Bayes approaches to this problem have proven useful in smoothing crude maps of disease rates. In recent years, models including both spatially correlated random effects and spatially unstructured random effects have been very popular. The spatially correlated random effects have been proposed in an attempt to capture a general clustering in the data. As an alternative, we propose replacing the spatially structured random effect with fixed clustering effects associated with particular areas. A reversible jump Markov chain Monte Carlo (RJMCMC) algorithm for posterior inference is described. We illustrate the model using the well-known New York leukaemia data.

Spatial statistical methods in health

The study of the geographical distribution of disease incidence and its relationship to potential risk factors (referred to here as "geographical epidemiology") has provided, and continues to provide, rich ground for the application and development of statistical methods and models. In recent years increasingly powerful and versatile statistical tools have been developed in this application area. This paper discusses the general classes of problem in geographical epidemiology and reviews the key statistical methods now being employed in each of the application areas identified. The paper does not attempt to exhaustively cover all possible methods and models, but extensive references are provided to further details and to additional approaches. The overall aim is to provide a picture of the "current state of the art" in the use of spatial statistical methods in epidemiological and public health research. Following the review of methods, the main software environments which are available to implement such methods are discussed. The paper concludes with some brief general reflections on the epidemiological and public health implications of the use of spatial statistical methods in health and on associated benefits and problems.

Cluster modelling of disease incidence via RJMCMC methods: a comparative evaluation. Reversible jump Markov chain Monte Carlo

The spatial modelling of small area health data has, for some time, included spatial autocorrelation as a random effect. This effect is non-specific and global and does not address the location of clusters of disease (a specific task). This paper addresses the need for specific and non-specific random effects within spatial epidemiology. In addition, individual frailty is also considered important and a computational algorithm based on reversible jump Markov chain Monte Carlo (RJMCMC) methods is described.

Bayesian detection and modeling of spatial disease clustering

Many current statistical methods for disease clustering studies are based on a hypothesis testing paradigm. These methods typically do not produce useful estimates of disease rates or cluster risks. In this paper, we develop a Bayesian procedure for drawing inferences about specific models for spatial clustering. The proposed methodology incorporates ideas from image analysis, from Bayesian model averaging, and from model selection. With our approach, we obtain estimates for disease rates and allow for greater flexibility in both the type of clusters and the number of clusters that may be considered. We illustrate the proposed procedure through simulation studies and an analysis of the well-known New York leukemia data.