In various applications with large spatial regions, the relationship between response variable and covariates is expected to exhibit complex patterns. We propose a spatially clustered varying coefficient model, where regression coefficients are allowed vary smoothly within each cluster but change abruptly across boundaries of adjacent clusters, we develop unified approach for simultaneous estim...