The SPATIALREG Procedure

Spatial Error Models

The spatial error model (SEM) accounts for spatial dependence in the error terms rather than in the dependent variable. Let y Subscript i denote the observation associated with the spatial unit bold s Subscript i for i equals 1 comma 2 comma ellipsis comma n. For these spatial units, let bold upper W be an n times n spatial weights matrix. Further, let bold x Subscript i be a p times 1 vector that denotes values of p regressors recorded at unit bold s Subscript i.

The SEM model can be described in vector form by using the following two-stage formulation (LeSage and Pace 2009),

where bold y equals left-parenthesis y 1 comma y 2 comma ellipsis comma y Subscript n Baseline right-parenthesis prime and bold-italic epsilon equals left-parenthesis epsilon 1 comma epsilon 2 comma ellipsis comma epsilon Subscript n Baseline right-parenthesis prime, with epsilon Subscript i Baseline tilde Overscript normal i normal i normal d Endscripts upper N left-parenthesis 0 comma sigma squared right-parenthesis. Moreover, bold upper X is an n times p matrix where each row consists of bold x prime Subscript i. In addition, bold-italic beta is a p times 1 parameter vector.

The log-likelihood function for the SEM model is

where bold upper B equals bold upper I Subscript n Baseline minus lamda bold upper W.

For the SEM model, the gradients are

In addition to the interaction effects among error terms, you might also want to include exogenous interaction effects in the model. In such cases, you need to consider a spatial Durbin error model (SDEM).

Last updated: August 11, 2020