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 to n destination regions resulting in a sample of N = n 2 observations that reflect an n by n flow matrix converted to a vector. Explanatory variables typically include origin and destination characteristics as well as distance between each region and all other regions. Our extension introduces latent spatial effects parameters structured to follow a spatial autoregressive process. Individual effects parameters are included in the model to reflect latent or unobservable influences at work that are unique to each region treated as an origin and destination. That is, we estimate 2n individual effects parameters using the sample of N = n 2 observations. We illustrate the method using a sample of commodity flows between 18 Spanish regions during the 2002 period.