We introduce a Bayesian hierarchical regression model that extends the traditional least-squares regression model used to estimate gravity or spatial interaction relations involving origin-destination flows. Spatial interaction models attempt to explain variation in flows from n origin regions...