将spdep
用于R(空间计量经济学的工具集合),我想自己建立一个回归矩阵。矩阵是两个矩阵的组合,一个是权重矩阵,另一个是协变矩阵。
构建我需要的矩阵的命令叫做create_WX()
。我似乎错误地使用它,因为它创建了一个只有NA
的列。
library(spdep)
data(oldcol)
# Build the weights matrix
lw <- nb2listw(COL.nb, style="W")
# Build the covariates matrix
X <- COL.OLD[, c("INC", "HOVAL")]
# Combine
WX <- create_WX(X, lw)
head(WX)
输出
.INC .HOVAL
[1,] 46.67233 NA
[2,] 32.81675 NA
[3,] 35.45417 NA
[4,] 35.40625 NA
[5,] 35.45850 NA
[6,] 45.36475 NA
应该没有一个NA
。我想知道为什么他们在那里,我怎么能避免这个问题?遗憾的是,create_WX()
的文档已嵌入errorsarlm()
并缩减为
create_WX
会产生空间滞后的RHS变量,并且会暴露出来 用于模型拟合函数。