We develop a two-stage spatial model for clustering disease risk.

First stage produces potential cluster structures, second stage fits these within a CAR model.

Model outperforms alternatives in a simulation study.

Applied to COPD mortalities in English local authorities and identifies north–south divide.