Стандартные ошибки воздействий в пространственной регрессии Лагсарла

Question

Я использую модели пространственной лаги и регрессии Дурбина, и я хотел бы оценить стандартные ошибки воздействий. Есть идеи, как это сделать?

Воспроизводимый пример ниже с использованием модели Дурбина

# 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

Вы должны иметь возможность использовать образцы MC, хранящиеся в объекте imp, чтобы получить стандартные ошибки, например:

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

дать те же значения Z, что и Imp