Article summary
This paper develops and provides software for a class of spatial hierarchical (multilevel) models for binary outcomes, motivated by situations where outcome dependence arises both from diffusion among lower-level units and from spatially structured unobserved group effects. The authors propose two model variants—one that treats unobserved group effects as spatially independent and another that treats them as spatially dependent—and use a Bayesian estimation approach. They evaluate performance with Monte Carlo simulations and illustrate the approach by analyzing causes of civil rights protests in the United States in the 1960s.