问题标签 [spdep]

我对使用“spdep”包的 SAR 模型的执行时间有疑问。

我将两个相同维度的不同数据集传递给同一个函数，它需要非常不同的时间（几秒与几小时）。

我把我的代码写给你，如果你有任何想法，请告诉我。

谢谢基亚拉

我想使用mgcvR 中的包来运行模型

其中 spdf 是 SpatialPolygonDataFrame。

我在理解如何创建nb对象时遇到问题，因为数据df具有重复geoid值（每个月一个）以及我运行的时间

我得到错误

谢谢

我正在尝试安装包 NormqPCR，我必须更新两个包，igraph 和 spdep，但直到现在我还不能，我不知道我必须做什么。这是代码。

有人知道可行的解决方案吗？谢谢！！

正如标题所说，我正在尝试使用该spdep软件包为 527k 点数据集计算 Local Moran，并根据距离创建邻域。我正在做的广义过程如下：

但是我遇到了一个问题，我无法使用 来创建社区dnearneigh，因为数据集太大，社区由 200-1000 个点组成。

我尝试了这里描述的解决方案，我得到了一个数据框，第一行带有 ID，第二行带有包含相邻点 ID 的列表（即）：

但是，我不知道如何创建nb所需的对象nb2listw，并且四处挖掘对我没有太大帮助。

有没有办法将此数据框转换为nb对象？如果我能够这样做，权重矩阵会像邻居一样难以创建吗？有没有另一种方法来计算这个卷的数据集的本地莫兰？

在@GavinSimpson 的 supperb 博客之后，我正在尝试适应 2015 年波兰地方政府选举结果。 https://www.fromthebottomoftheheap.net/2017/10/19/first-steps-with-mrf-smooths/ 我用 6 位标识符的 shp 数据加入 xls（可能有前导 0 s）。我将其保留为文本变量。编辑，我简化了标识符，现在使用从 1 到 nrow 的序列来简化我的问题。

这是错误消息：

选举$idx 是一个因素。我用它来给 boroughs_nb 命名，以绝对确定我有相同数量的级别。我究竟做错了什么？

编辑：满足错误消息中提到的条件：

我的数据集中有 9,150 个多边形。我试图运行一个空间自回归模型（SAR）spdep来测试我的结果变量的空间依赖性。运行模型后，我想检查直接/间接影响，但遇到了一个错误，似乎与权重矩阵中邻居的长度不等于n.

我尝试运行与 SLX 模型（空间滞后 X）完全相同的方程，并且impacts()工作正常，即使我的集合中有一些没有邻居的多边形。我用谷歌搜索并查看了spdep文档，但找不到有关如何解决此错误的线索。

我正在尝试从伦敦下层超级输出区 (LSOA) 的多边形列表中构建邻居列表。我从伦敦数据存储区获得了我的 shapefile，通过此处下载：https ://data.london.gov.uk/dataset/statistical-gis-boundary-files-london

当我阅读 2004 shapefile ( LSOA_2004_London_Low_Resolution.shp) 时，它运行良好：

然而，当我阅读 2011 年 shapefile ( LSOA_2011_London_gen_MHW.shp) 时，似乎认为大多数 LSOA 没有邻居。

当我使用具有相同代码的 2011 版本时，我得到以下信息：

有谁知道我做错了什么？显然应该有邻居的 LSOA 显示为 0 个链接。我不知道为什么...！

通过以下可重现的示例更容易解释我要做什么：

这段代码生成以下图

我想做的是用与边缘相同的颜色为节点（或顶点）着色。可能吗？

先感谢您。

下面的会话信息（用于可重复性）。

我对 R 相当陌生，我正在尝试光栅化 spdep 的 localG 函数的输出。这段代码：

产生错误unable to find an inherited method for function 'rasterize' for signature '"localG", "RasterLayer"'

我尝试将 localG 类更改为 data.frame，但它创建了一个仍然不会光栅化的 1 列矩阵。

总结一下：我应该怎么做才能获得 localG 输出的栅格？

提前致谢！

我正在尝试使用该lagsarlm命令估计空间 SAR 模型，但出现以下错误：

lagsarlm(f1, data = df, spatialList, tol.solve = 1e-30) 中的错误：滞后因变量中的 NA 另外：警告消息：在 lag.listw(listw, y, zero.policy = zero.policy) 中：滞后值中的 NA

显然，这是因为我有岛屿（与任何其他观察无关的观察）。如果我在不创建岛的情况下运行以下代码（通过删除 W[1,1:50] <- 0），则代码可以正常工作。

这只是一个玩具示例。鉴于我的真实矩阵有几个岛，关于如何进行估计的任何想法？

非常感谢