Abstract
Summary: A two dimensional autoregressive moving average spatial model is used to analyse spatial interaction. Maximum likelihood estimates of the unknown parameters are derived as the solution of a system of nonlinear equations, and are shown to be best asymptotic normal. One important computational procedure is discussed. The argument is extended to the general regression model with autoregressive moving average residuals. Explicit computational formulae are given.