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....