---
title: 'SpatialPlus: Unbiased Estimation in Spatial Models'
url: https://www.emergentmind.com/topics/spatialplus
type: topic
---

# SpatialPlus: Unbiased Estimation in Spatial Models

SpatialPlus (“Spatial+”) comprises a set of statistical methodologies devised to mitigate bias arising from spatial confounding in spatial regression models. In these models, collinearity between spatially indexed covariates and spatial random effects often introduces significant bias in fixed-effect estimates, particularly when both are smooth functions of space. Spatial+ achieves bias reduction by actively decorrelating covariates from the spatial structure, thereby isolating fixed-effect inference from the influence of spatial smoothing. The core concept originated with frequentist two-stage adjustment using thin-plate splines, has been extended to Bayesian hierarchical modeling with smoothness parameter coupling, and has inspired spectral domain and areal data generalizations.

## 1. Spatial Confounding and Smoothing-Induced Bias

Spatial confounding occurs when observed covariates used in regression share similar spatial patterns with random effects or smoothers included to account for spatial residuals. In the canonical model,
\[
Y_i = X_i^\top\beta + f(s_i) + \varepsilon_i,\quad \varepsilon_i \sim N(0,\sigma^2),
\]
$X$ and $f$ may be approximately collinear, allowing the smoother $f$ to absorb variation attributable to $X$. As articulated by Rice (1986), penalized smoothing (e.g., with thin-plate splines) introduces a *smoothing-induced bias* in $\hat\beta$, which can decay at the same or a slower rate than the estimator's standard error, so that smoothing-induced bias can dominate in finite samples [2009.09420].

Generalized spatial regression formulations—including those using partial splines, GMRF/ICAR random effects, or multivariate areal models—are equally susceptible. The identifiability of $\beta$ is compromised by shared large-scale and/or high-frequency spatial features between covariates and random effects, a phenomenon not remedied merely by increased sample size.

## 2. The Spatial+ Methodology: Frequentist Formulation

The Spatial+ approach, as introduced by Dupont, Wood, and Augustin, involves a two-stage residualization procedure that decouples the covariate from the spatial pattern captured by the smoother. The prototypical procedure is:
1. For each covariate $X$, fit the spatial model $X = f_X(s) + \varepsilon_X$ using thin-plate splines and extract the residuals $\tilde X = X - \widehat f_X(s)$.
2. Fit the primary regression model using $\tilde X$ in place of $X$:
\[
Y_i = \tilde X_i^\top\beta + f^+(s_i) + \varepsilon_i
\]
Bias analysis shows that the leading-order smoothing-induced bias in $\hat\beta$ disappears to first order, resulting in an estimator that is asymptotically unbiased:
\[
\hat\beta^+ - \beta = o(n^{-1/2}) + O\left((\lambda \lambda_x)^{1/2}\right), \quad \text{with } \Var(\hat\beta^+) \sim \frac{\sigma^2}{n\sigma_x^2}
\]
When tuning parameters $\lambda$ and $\lambda_x$ are chosen optimally, the bias decays faster than sampling variability, unlike in the unadjusted model [2009.09420].

This residualization ensures that, under the natural spline basis, $\tilde X$ is (approximately) orthogonal to $f^+$, preventing the smoother from "stealing" covariate effects.

## 3. Extensions: Bayesian Spatial+, Joint Smoothing Priors, and Uncertainty Quantification

A major limitation of frequentist Spatial+ is the lack of propagation of uncertainty from covariate residualization to the primary response model. The Bayesian Spatial+ framework jointly models the covariate and response:
- **Stage 1** (Covariate): $x(\mathbf{s}_i) = \beta_0^x + f^x(\mathbf{s}_i) + \varepsilon^x_i$
- **Stage 2** (Response): $y(\mathbf{s}_i) = \beta_0 + \beta\,\varepsilon^x_i + f(\mathbf{s}_i) + \varepsilon_i$
where $\varepsilon^x_i = x(\mathbf{s}_i) - \beta_0^x - f^x(\mathbf{s}_i)$.

All model parameters—including smoothers' coefficients, smoothing penalties, and regression coefs—are sampled simultaneously in a hierarchical setting (typically via HMC/NUTS in RStan). This enables *full posterior* propagation of uncertainty, meaning that interval estimates for $\beta$ reflect all sources of variability, in contrast to the plug-in two-stage approach [2309.05496].

A novel feature of Bayesian Spatial+ is the *joint prior on smoothing parameters*, ensuring the response surface $f$ is smoother than the covariate surface $f^x$:
\[
\tau^2 \sim \mathrm{Weibull}(0.5, \zeta),\quad \tau_x^2 = \frac{\sigma_x^2}{\sigma^2}\tau^2 + \xi, \quad \xi > 0
\]
This restricts the response smoother from capturing high-frequency spatial features that are not present in the covariate residuals, directly mitigating spectral confounding and preventing overfitting in the second stage.

## 4. Spectral and Areal Data Generalizations

Areal models, especially in the context of multivariate counts, are prone to spatial confounding due to similar issues of collinearity between covariates and spatial random effects. The *simplified spatial+* approach as in Urdangarin et al. [2308.11260] operates via spectral decomposition:
- Decompose the areal adjacency/precision matrix $Q$ into eigenvectors: $Q=U\Delta U'$.
- Covariate $X$ is split into large-scale ($X_L=U_LU_L'X$) and short-scale ($X_S=U_SU_S'X$) components, where $U_L$ holds $k$ eigenvectors corresponding to the smoothest spatial direction.
- The regression model uses only $X_S$, thus orthogonalizing the fixed effect against the spatial random effect $\Theta$.

This method is particularly computationally efficient (implemented via INLA) and is adaptable to multivariate, multi-covariate settings, albeit with the need to select the cutoff $k$.

## 5. Empirical Performance and Simulation Studies

Simulation studies consistently demonstrate that:
- Standard spatial regression (uncorrected) yields substantial bias in $\hat\beta$ under moderate to strong confounding.
- Both frequentist and Bayesian Spatial+ methods substantially reduce bias; the Bayesian variant achieves near-nominal coverage for $\beta$ and better MSE, especially under moderate or large sample sizes [2309.05496, 2009.09420, 2308.11260].
- When confounding arises from small-scale spatial structure, all methods suffer nonvanishing bias, though Bayesian Spatial+ maintains superior mean squared error properties.

The table below (Scenario 1 from [2308.11260]) illustrates the gain:

| Model              | True $\beta_1$ | Mean Est. | Bias    | 95% Coverage |
|--------------------|:--------------:|:---------:|:-------:|:------------:|
| M–Spatial          | -0.15          | -0.2996   | -0.1496 | 35.7%        |
| M–SpatPlus(10)     | -0.15          | -0.1863   | -0.0363 | 95.3%        |

In real-world applications (e.g., German precipitation and Indian district-level crime counts), Spatial+ adjustments yield enlarged effect estimates and wider, more reliable credible intervals, preventing overconfident inference and improving interpretability [2309.05496, 2308.11260].

## 6. Implementation and Practical Considerations

Spatial+ is readily implementable with standard spline-based or areal spatial modeling software. For frequentist estimation, software such as `mgcv` enables specification of sequential residualization via thin-plate splines. For Bayesian inference, RStan with HMC/NUTS is used, with diagnostic checks on $\hat{R}$ and effective sample size [2309.05496]. INLA enables rapid inference for the simplified spatial+ model in areal and multivariate contexts.

The selection of the spectral cutoff $k$ in areal settings remains an open methodological question; practical heuristics suggest removing 7–20% of eigenvectors for smooth covariates, up to 40% for mixed spatial scales. Posterior standard errors can be inflated, and model selection metrics (e.g., DIC, WAIC) may inadequately distinguish among plausible values of $k$.

Spatial+ is directly extensible to non-Gaussian responses and to a wide range of spatial modeling frameworks, provided residualization is performed using the corresponding smoothing/quadratic penalty defined by the prior or spatial effect [2009.09420, 2308.11260].

## 7. Limitations and Extensions

Current limitations include the restriction of original Spatial+ to linear covariate effects, the need for separate residualization for each linear term, and incomplete theoretical guarantees for high-dimensional or nonlinear covariate scenarios. Automated or cross-validated selection of spectral cutoffs in areal implementations, as well as variance component inference and improved standard error calibration, are active areas of research [2308.11260]. Extensions to allow joint decorrelation of multiple covariates and to handle nonlinear/additive smooth effects remain open.

Spatial+ and its modern Bayesian generalizations provide a principled solution to the pervasive problem of spatial confounding, yielding unbiased or nearly-unbiased effect estimates, valid uncertainty quantification, and better fixed-effect inference in spatial regression and areal models [2009.09420, 2309.05496, 2308.11260].

Source: https://www.emergentmind.com/topics/spatialplus