Spatial Coarsened Exact Matching (Spatial CEM)
Spatial Coarsened Exact Matching Estimator · Also known as: Spatial CEM, Geographic CEM, Spatial exact matching, CEM with spatial covariates
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.
Read the full method
Sign in with a free account to read this section.
Method map
The neighbourhood of related methods — select a node to explore.
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
- 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).
- 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.
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
- Iacus, S. M., King, G., & Porro, G. (2012). Causal Inference without Balance Checking: Coarsened Exact Matching. Political Analysis, 20(1), 1-24. DOI: 10.1093/pan/mpr013 ↗
- Anselin, L., & Rey, S. J. (Eds.) (2014). Modern Spatial Econometrics in Practice: A Guide to GeoDa, GeoDaSpace and PySAL. GeoDa Press. ISBN: 978-0986342103
How to cite this page
ScholarGate. (2026, June 3). Spatial Coarsened Exact Matching Estimator. ScholarGate. https://scholargate.app/en/causal-inference/spatial-coarsened-exact-matching
Which method?
Set this method beside its closest kin and read them side by side — the library lays the books on the table; the choice is yours.
- Coarsened Exact MatchingCausal inference↔ compare
- Difference-in-DifferencesEconometrics↔ compare
- Propensity Score MatchingResearch Statistics↔ compare
- Spatial Doubly Robust EstimationCausal inference↔ compare
- Spatial Propensity Score MatchingCausal inference↔ compare
- Spatial Regression Discontinuity DesignCausal inference↔ compare