This paper focuses on panel data models combining spatial dependence with a nested (hierarchical) structure. We use a generalized moments estimator to estimate the spatial autoregressive parameter and the variance components of the disturbance process. A spatial counterpart of the Cochrane‐Orcutt transformation leads to a feasible generalized least squares procedure to estimate the regression parameters. Monte Carlo simulations show that our estimators perform well in terms of root mean square error compared to the maximum likelihood estimator. The approach is applied to English house price data for districts nested within counties.