Estimation of higher‐order spatial autoregressive cross‐section models with heteroscedastic disturbances