The article investigates the finite sample properties of estimators for spatial autoregressive models where the disturbance terms may follow a spatial autoregressive process. In particular we investigate the finite sample behavior of the feasible generalized spatial two-stage least squares (FGS2SLS) estimator introduced by Kelejian and Prucha (1998), the maximum likelihood (ML) estimator, as well as that of several other estimators. We find that the FGS2SLS estimator is virtually as efficient as the ML estimator. This is important because the ML estimator is computationally burdensome, and may even be forbidding in large samples, while the FGS2SLS estimator remains computationally feasible in large samples.