我正在使用错误()函数在斯戴普库中拟合一个空间错误模型。利用bptest.sarlm()函数计算的空间模型的Breusch检验表明存在异方差.
自然的下一步将是获得稳健的标准错误估计和更新p-值。在bptest.sarlm()函数的文档中,说明如下:
技术上也可以通过使用sarlm对象的“lm.target”组件对标准误差估计进行异方差校正--在lmtest和三明治包中使用函数。
并提出了以下代码(作为参考):
lm.target <- lm(error.col$tary ~ error.col$tarX - 1)
if (require(lmtest) && require(sandwich)) {
print(coeftest(lm.target, vcov=vcovHC(lm.target, type="HC0"), df=Inf))} 其中error.col是估计的空间误差模型。
现在,我可以很容易地使代码适应我的问题,并获得健壮的标准错误。然而,我想知道:
发布于 2021-12-06 14:32:11
我在github上报道了这个问题,并且有一个Roger Bivand的答复:
不,根本不推荐这种方法。要么使用sphet,要么使用Bayesian方法给出边缘的后验分布。我会删除那些令人困惑的文档。在空间误差模型的情况下,$y - \rho $和tarX类似。请注意,localmoran.exact()和localmoran.sad()的文档中只在spdep中出现of等;您是否使用过期包版本?
https://stackoverflow.com/questions/70173881
复制相似问题