[R-sig-Geo] Mean Effects SAR Model

Larry Layne ljlayne at unm.edu
Tue Apr 4 00:36:48 CEST 2006


--On Monday, April 03, 2006 Roger Bivand wrote:

> OK. I can replicate this - there is still a bug in handling intercept-only
> models in lagsarlm() with method="SparseM"; it was partly fixed in January
> 2005. The equivalent model for errorsarlm() works:
>
> and the residuals are the same as for the lag model, as you can see. I'll
> try to fix the problem tomorrow.
>
> Thanks for a helpful report.

Thanks very much for your help. Using residuals(lagsarlm(Y ~ 
1,method="eigen")) computes what I am looking for. Not sure why I was being 
so dense with respect to specifying an intercept-only model. One of those 
days.

Larry




More information about the R-sig-Geo mailing list