Abstract:There is currently much interest in conducting spatial analyses of health outcomes at the small-area scale. This requires sophisticated statistical techniques, usually involving Bayesian models, to smooth the underlying risk estimates because the data are typically sparse. However, questions have been raised about the performance of these models for recovering the "true" risk surface, about the influence of the prior structure specified, and about the amount of smoothing of the risks that is actually performed. We describe a comprehensive simulation study designed to address these questions. Our results show that Bayesian disease-mapping models are essentially conservative, with high specificity even in situations with very sparse data but low sensitivity if the raised-risk areas have only a moderate (less than 2-fold) excess or are not based on substantial expected counts (> 50 per area). Semiparametric spatial mixture models typically produce less smoothing than their conditional autoregressive counterpart when there is sufficient information in the data (moderate-size expected count and/or high true excess risk). Sensitivity may be improved by exploiting the whole posterior distribution to try to detect true raised-risk areas rather than just reporting and mapping the mean posterior relative risk. For the widely used conditional autoregressive model, we show that a decision rule based on computing the probability that the relative risk is above 1 with a cutoff between 70 and 80% gives a specific rule with reasonable sensitivity for a range of scenarios having moderate expected counts (approximately 20) and excess risks (approximately 1.5- to 2-fold). Larger (3-fold) excess risks are detected almost certainly using this rule, even when based on small expected counts, although the mean of the posterior distribution is typically smoothed to about half the true value.

Nonparametric geostatistical risk mapping

Nonparametric conditional risk mapping under heteroscedasticity

Bayesian statistics in spatial epidemiology]

Likelihood-based Estimation of Spatial Intensity and Variation in Disease Risk from Locations Observed with Error

Spatial hierarchical modeling of threshold exceedances using rate mixtures

New Bayesian method for estimation of Value at Risk and Conditional Value at Risk

Nonparametric quantile regression for spatio-temporal processes

Interpreting posterior relative risk estimates in disease-mapping studies

Bayesian Geostatistical Modeling for Discrete-Valued Processes

Spatial extremes and stochastic geometry for Gaussian-based peaks-over-threshold processes

Bayesian nonparametric models for spatially indexed data of mixed type

A Study on the Method for Determining Statistical Risk Boundaries for Spatial Data with Application to Traffic Accident Data

Estimating Systemic Risk within Financial Networks: A Two-Step Nonparametric Method

Bayesian nonparametric nonhomogeneous Poisson process with applications to USGS earthquake data

Patterns in Spatio-Temporal Extremes

A new Bayesian method for estimation of value at risk and conditional value at risk

Interval-Valued Kriging Models for Geostatistical Mapping with Imprecise Inputs

Efficient global estimation of conditional-value-at-risk through stochastic kriging and extreme value theory

Bayesian inference for high-dimensional nonstationary Gaussian processes

Big problems in spatio-temporal disease mapping: methods and software

Scalable Bayesian modeling for smoothing disease risks in large spatial data sets