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

A new two-way map between time domain and numerical magnitudes or values domain (v-dom) provides a new solution to heteroscedasticity. Since sorted logs of squared fitted residuals are monotonic in the v-dom, we obtain a parsimonious fit there. Two theorems prove consistency, asymptotic...