空间回归滞后的标准误差

Question

我正在使用空间滞后和durbin回归模型，我想估计影响的标准误差。关于如何做到这一点的任何想法？

使用durbin模型

的下面的可重复示例

# data  
  data(oldcol) 

# neighbours lists
  lw <- nb2listw(COL.nb, style="W")

# regression
  fit_durb <- lagsarlm(CRIME ~ INC + HOVAL, data=COL.OLD, type="Durbin",
                       listw=lw, method="eigen", 
                       zero.policy=T, na.action="na.omit")
# power traces
  W <- as(lw, "CsparseMatrix")
  trMC <- trW(W, type="MC", listw = lw)

# Impacts
  imp <- summary(impacts(fit_durb, tr=trMC, R=1000), zstats=TRUE, short=TRUE)

Answer 1

您应该能够使用存储在imp对象中的MC样本来获取标准错误，例如：

test1<-lapply(imp$sres, function(x){apply(x, 2, mean)})
test2<-lapply(imp$sres, function(x){apply(x, 2, sd)}
test1$direct/test2$direct

给出与imp

返回的z值相同的值