R USCounties 美国县的连续性矩阵

说明

该矩阵使用至少一个共享顶点或边的皇后准则给出了 3111 个美国县的邻近区域。

用法

data(USCounties)

格式

类 dsCMatrix 的 3111 \times 3111 稀疏对称矩阵，具有 9101 个非零条目。

例子

data(USCounties, package = "Matrix")
(n <- ncol(USCounties))
I <- .symDiagonal(n)

set.seed(1)
r <- 50L
rho <- 1 / runif(r, 0, 0.5)

system.time(MJ0 <- sapply(rho, function(mult)
    determinant(USCounties + mult * I, logarithm = TRUE)$modulus))

## Can be done faster by updating the Cholesky factor:

C1 <- Cholesky(USCounties, Imult = 2)
system.time(MJ1 <- sapply(rho, function(mult)
    determinant(update(C1, USCounties, mult), sqrt = FALSE)$modulus))
stopifnot(all.equal(MJ0, MJ1))

C2 <- Cholesky(USCounties, super = TRUE, Imult = 2)
system.time(MJ2 <- sapply(rho, function(mult)
    determinant(update(C2, USCounties, mult), sqrt = FALSE)$modulus))
stopifnot(all.equal(MJ0, MJ2))

来源

GAL 点阵文件‘usc_q.GAL”(2008 年检索自“http://sal.uiuc.edu/weights/zips/usc.zip’经 Luc Anselin 许可使用和分发)被读入R使用函数read.gal从包装中spdep.

使用包 spdep 和 spatialreg 中的函数 nb2listw 和 similar.listw 使用 row-standardized(然后对称)空间权重来增强邻居列表。使用 spatialreg 中的 as_dsTMatrix_listw 将生成的 listw 对象强制转换为 dsTMatrix 类，然后强制转换为 dsCMatrix 类。

参考

Ord, J. K. (1975). Estimation methods for models of spatial interaction. Journal of the American Statistical Association, 70(349), 120-126. doi:10.2307/2285387