Regression modelCausal inferenceQuasi-experimental / causal inferenceModel

Spatial Coarsened Exact Matching (Spatial CEM)

Also known as: Spatial CEM, Geographic CEM, Spatial exact matching, CEM with spatial covariates

OriginatorIacus, King & Porro (CEM foundation, 2012); extended to spatial contexts by applied spatial econometriciansYear2012 (CEM foundation); spatial extension in applied literature 2015-presentSources2Related methods6

Spatial Coarsened Exact Matching applies the Coarsened Exact Matching framework to study designs involving geographic units — neighbourhoods, census tracts, municipalities, or grid cells. Covariates are coarsened into discrete bins and units are matched exactly on those bins, with spatial attributes (location, adjacency, geographic characteristics) incorporated as matching dimensions to control for spatial confounding.

Key highlights

  • Guarantees exact covariate balance within matched strata, eliminating model-dependent extrapolation across geographic or covariate space.
  • Transparently controls for spatial confounding by including geographic covariates as mandatory matching dimensions.
  • Non-parametric: does not require specifying a propensity score model, reducing sensitivity to functional-form assumptions.
  • The pruning of unmatched strata makes the scope of inference explicit — researchers know exactly which portion of the spatial sample supports the estimate.
  • Straightforward to implement in R (MatchIt, CEM packages) and easily combined with spatial data tools (sf, sp, GeoDa).

Intuition

This section is available to Pro members. Upgrade to Pro

How it works

This section is available to Pro members. Upgrade to Pro

When to use it

Use Spatial CEM when your units are geographic areas or locations, treatment assignment varies spatially, and you suspect that location-specific characteristics (urban density, proximity to infrastructure, regional policy regimes) confound the treatment-outcome relationship. It is particularly well-suited when the number of spatial units is moderate (roughly 50 to several thousand), covariate distributions overlap sufficiently between treated and control areas, and you want transparent, assumption-light balance rather than model-based propensity scores. Do not use it when geographic units are too few to form matched strata, when spatial spillovers from treated to control units are pervasive (SUTVA violations), or when exact stratum matching would prune so many units that the remaining sample is unrepresentative.

Strengths & limitations

Strengths
  • Guarantees exact covariate balance within matched strata, eliminating model-dependent extrapolation across geographic or covariate space.
  • Transparently controls for spatial confounding by including geographic covariates as mandatory matching dimensions.
  • Non-parametric: does not require specifying a propensity score model, reducing sensitivity to functional-form assumptions.
  • The pruning of unmatched strata makes the scope of inference explicit — researchers know exactly which portion of the spatial sample supports the estimate.
  • Straightforward to implement in R (MatchIt, CEM packages) and easily combined with spatial data tools (sf, sp, GeoDa).
Limitations
  • Coarsening choices (bin boundaries) are subjective and can materially affect which units are matched and what effect is estimated.
  • When many spatial covariates are included or bins are fine-grained, strata become sparse, pruning large fractions of the sample and reducing external validity.
  • Spatial CEM controls for observed spatial confounders but remains vulnerable to unobserved spatial heterogeneity not captured by included covariates.
  • Does not directly model spatial spillovers; if treatment in one area affects outcomes in adjacent control areas, the ATT estimate is contaminated.
  • Matching on spatial bins ignores residual within-stratum geographic variation, which may still be confounding.

Common pitfalls

This section is available to Pro members. Upgrade to Pro

Applications

This section is available to Pro members. Upgrade to Pro

Frequently asked

How does Spatial CEM differ from standard CEM?

Standard CEM matches on any set of observed covariates. Spatial CEM explicitly includes geographic covariates — such as region, urban-rural classification, distance deciles, or spatial cluster indicators — as mandatory matching dimensions. This ensures matched pairs are both observationally and geographically comparable, controlling for location-based confounding that standard CEM would ignore if spatial variables were not included.

How do I choose coarsening bins for spatial covariates?

Bins should reflect substantively meaningful geographic strata set before looking at outcomes — for example, administrative regions, urban-suburban-rural classification, or distance quartiles from a landmark. Avoid data-driven bin selection that maximises matched sample size, as this introduces specification searching. Sensitivity analyses across plausible coarsening schemes are recommended to assess robustness.

What if many units are pruned during matching?

Large pruning means the treated and control covariate distributions overlap only in a narrow region of the geographic-covariate space. The resulting ATT is valid for that overlapping subpopulation but may not generalise to all treated areas. Report the fraction pruned and characterise how retained and pruned treated units differ, so readers understand the scope of inference.

Should I use spatially clustered standard errors?

Yes, if residuals are spatially autocorrelated within matched strata. After matching, test for residual spatial autocorrelation (e.g., Moran's I on residuals). If significant, cluster standard errors at the spatial block or geographic unit level, or use spatial HAC estimators. Ignoring spatial autocorrelation leads to anti-conservative inference.

Can Spatial CEM handle spatial spillovers between treated and control units?

No — like all standard matching estimators, Spatial CEM assumes the Stable Unit Treatment Value Assumption (SUTVA): a unit's outcome depends only on its own treatment status. If treated areas generate externalities in neighbouring control areas (e.g., displacement effects, knowledge spillovers), the control group outcomes are contaminated. In such cases, consider spatial difference-in-differences with spillover-robust estimators, or define geographic buffers that exclude potentially contaminated control units.

Sources

  1. 1.
    Iacus, S. M., King, G., & Porro, G. (2012). Causal Inference without Balance Checking: Coarsened Exact Matching. Political Analysis, 20(1), 1-24.
  2. 2.
    Anselin, L., & Rey, S. J. (Eds.) (2014). Modern Spatial Econometrics in Practice: A Guide to GeoDa, GeoDaSpace and PySAL. GeoDa Press.
    ISBN 978-0986342103

You have read it. What now?

Cite this page

ScholarGate. (2026, June 3). Spatial Coarsened Exact Matching. ScholarGate. https://scholargate.app/causal-inference/spatial-coarsened-exact-matching

Spatial Coarsened Exact Matching (Spatial CEM) | ScholarGate