Closed orlando-sabogal closed 1 year ago
No, the approach is not recommended at all. Either use sphet or a Bayesian approach giving the marginal posterior distribution. I'll drop the confusing documentation. tary
is $y - \rho W y$ and similarly for tarX
in the spatial error model case. Note that tary
etc. only occur in spdep in documentation for localmoran.exact()
and localmoran.sad()
; were you using out of date package versions?
Thanks for your fast response. I'll go with sphet. I am using the latest version of the package.
Please also consider installing the updated development version from github rather than the otdated version on CRAN.
This is probably an issue more related to the documentation. I asked this question on stack overflow but I think it is helpful to bring it here.
I am fitting a Spatial Error Model using the errorsarlm() function in the spdep library. The Breusch-Pagan test for spatial models, calculated using the bptest.sarlm() function, suggest the presence of heteroskedasticity.
A natural next step would be to get the robust standard error estimates and update the p-values. In the documentation of the bptest.sarlm() function says the following:
and the following code (as reference) is presented:
where error.col is the spatial error model estimated.
Now, I can easily adapt the code to my problem and get the robust standard errors. Nevertheless, I was wondering: