Likelihood functions of spatial autoregressive models with normal but heteroskedastic disturbances have been already derived [Anselin (1988, ch.6)]. But there is no implementation for maximum likelihood estimation of these likelihood functions in general (heteroskedastic disturbances) cases....

The likelihood functions for spatial autoregressive models with normal but heteroskedastic disturbances have been derived [Anselin (1988, ch.6)], but there is no implementation of maximum likelihood estimation for these likelihood functions in general cases with heteroskedastic disturbances....