If the ordinary least squares (OLS) diagnostics discussed in the previous chapter indicate the existence of spatial lag or spatial error dependence, the researcher will wish to model the type of dependence indicated by these diagnostics. If the OLS diagnostics indicate the presence of a diffusion process, the researcher will wish to estimate a spatial lag model via maximum likelihood (ML) estimation or an instrumental variables specification incorporating instruments for the spatially lagged dependent variable. Alternatively, if the OLS diagnostics indicate the existence of spatial error dependence, the researcher may choose to estimate a more fully specified OLS model to model the spatial dependence or may choose to employ a ML or generalized method of moments (GMM) approach incorporating the spatial dependence in the errors.

The spatial dependence diagnosed via the diagnostics discussed in Chapter 5 may alternatively be produced by spatial heterogeneity in the effects of covariates. If this is the only source of spatial dependence, modeling this heterogeneity will be sufficient to capture the spatial dependence. As a consequence, any specification search should also consider the possibility of spatial heterogeneity, which is the focus of Chapter 7. This chapter will first, however, examine alternative approaches for modeling spatial dependence if spatial heterogeneity is not present.

This chapter begins by examining ML estimation of spatial lag models that derives from Ord (1975). Next, I explore alternative instrumental variables and GMM estimators for spatial lag dependence. Next, I turn to approaches for estimating spatial error models. I conclude by considering areas of concern in the estimation of spatial models. These include estimators for large sample sizes and diagnostics for continued spatial dependence.

MAXIMUM LIKELIHOOD SPATIAL LAG ESTIMATION

The mixed regressive, spatial autoregressive model, or spatial lag model, extends the pure spatial autoregressive model considered in Section 3.2 to include also the set of covariates and associated parameters:

y = ρWy+Xβ+ε

where X is again an N by K matrix of observations on the covariates, β is a K by 1 vector of parameters, and the remaining notation is as discussed in Section 3.2.