We develop Bayesian model selection via fractional Bayes factors to simultaneously assess spatial dependence and select regressors in Gaussian hierarchical models with intrinsic conditional autoregressive (ICAR) random effects. Selection of covariates structure is difficult, as confounding creates a tension between fixed Researchers have commonly performed separately for effects models. Simulta...