# Horseshoe Priors for Spatial Small Area Estimation: Regular Variation, Tail Robustness, and Deep Learning

arXiv: https://arxiv.org/abs/2606.30659 · HTML: https://mlsanity.com/paper/2606.30659/

## Key facts

- No citing papers in the map yet (data through 14 Oct 2026).
- Builds on 5 of its 37 references that are in the map.
- First posted on arXiv 16 Jun 2026.

## Builds on (most cited first, 5 of 5)

- [Conformalized Quantile Regression](https://mlsanity.com/paper/1905.03222/) (2019)
- [Lasso Meets Horseshoe : A Survey](https://mlsanity.com/paper/1706.10179/) (2017)
- [Deep Learning for Spatio-Temporal Modeling: Dynamic Traffic Flows and High Frequency Trading](https://mlsanity.com/paper/1705.09851/) (2017)
- [An intuitive Bayesian spatial model for disease mapping that accounts for scaling](https://mlsanity.com/paper/1601.01180/) (2016)
- [Default Bayesian analysis with global-local shrinkage priors](https://mlsanity.com/paper/1510.03516/) (2015)

## Abstract (arXiv)

Small area estimation borrows strength across domains to repair the poor precision of direct survey estimators. Two philosophies dominate the area-level literature. The first, descending from Ghosh and Rao (1994), borrows strength through structured Gaussian smoothing: an intrinsic conditional autoregression or its BYM2 reparameterization pools each area towards its neighbours. The second borrows strength globally but acts locally through a heavy-tailed global-local prior on area effects, of which the horseshoe of Carvalho et al. (2010) is the canonical instance; Tang et al. (2018) first brought this idea to small area estimation. We study the horseshoe Fay-Herriot model with known unequal sampling variances and make four contributions. First, a tail-robustness theorem: through a heteroscedastic Tweedie identity the posterior mean leaves strongly signalled areas essentially unshrunk, bounding the influence of an outlying direct estimate, unlike Gaussian random-effect models. Second, standardizing by the known design variances transfers the minimax contraction and credible-set theory of the homoscedastic sequence model to the heteroscedastic Fay-Herriot problem; the posterior contracts at the nearly-black minimax rate, with a matching lower bound. Third, we give an account of when structured smoothing and global-local shrinkage each win. Fourth, an O(m) Gibbs sampler, simulations, and an analysis of the Scottish lip cancer data confirm the account: on strongly spatial data the smoother predicts held-out districts best, yet the horseshoe flags exceptional districts that smoothing suppresses. Throughout we argue, following the regular-variation theory of Bhadra et al. (2016), that these properties make the horseshoe a sound default prior for area effects: it borrows strength aggressively yet lets genuinely exceptional areas speak for themselves, with no tuning and no neighbourhood graph.