⁹The spautolm function was originally written to replicate results in Waller, Gotway (2004, p. 378), and weights have now been added to the errorsarlm to match (relative differences in coefficient estimates are 1e-8, because of different underlying numerical methods).