问题标签 [spdep]

我是 R 编程的新手。我目前正在研究一个可以改进 spdep 包中的 mstree() 的研究项目。我有一个关于如何使 Polygon()、Polygons()、SpatialPolygons() 和 SpationPolygonsDataFrame() 工作的问题。以下是我尝试过的代码：

...和错误！我上面的代码有什么问题？我在互联网上搜索和搜索以找出使用 Polygon()、Polygons()、SpatialPolygons()、SpatialPolygonsDataFrame() 和 poly2nb() 的正确方法，但我只是不知道如何使这些函数工作！我已经阅读了 R 手册，但它并没有太多帮助......不过，一件有趣的事情是，当我查看 Polygons 对象时，Polygons 插槽是空的；同样，当我查看 Sr 时，SpatialPolygons 插槽和 Polygons 插槽是空的。谢谢，

我有 4000 条人工林的数量记录。我需要计算整个种植园的 Moran's I。我使用ape库，因为据说 spdep速度较慢。我的代码是这样的：

当我运行代码时，会出现类似溢出的错误。

使用 ff 库

更多错误信息：

• 我怎样才能得到整个种植园的计算？

• 我应该使用sdep而不是猿库吗？

• 并行库如何解决这个问题？

提前感谢胡安

我正在尝试使用 mat2listw 函数在 R 中创建一个权重对象。我有一个在 Excel 中创建并读入 R 的非常大的空间权重矩阵（大约 22,000x22,000），我现在正在尝试实现：

我收到以下错误：

这里出了什么问题？我当前的矩阵都是 0 和 1，没有缺失值和负元素。我错过了什么？

我会很感激任何建议。在此先感谢您的帮助！

我有一个点列表，我想使用 Moran's I 并通过将感兴趣区域除以 4 x 4 样方来检查自相关。

现在，我在 Google 上找到的每个示例（例如http://www.ats.ucla.edu/stat/r/faq/morans_i.htm）都使用某种测量值作为 Moran's I 函数的第一个输入，无论哪个使用了库（我查看了 ape 和 spdep 包）。

但是，我所拥有的只是我想要检查相关性的点本身。

问题是，尽管这听起来很有趣（或悲伤），但我不知道我在这里做什么。我不是一个（空间）统计专家，我只想知道使用 Moran's I 的点集合是分散的、聚集的还是随机的。

我的方法正确吗？如果不是，我在哪里做错了什么？

谢谢

这是我到目前为止所拥有的：

我想测试一个物种的存在/不存在记录的空间自相关性。我的数据在 GIS 中包含超过 130,000 个网格和大约 700 个物种的存在记录。

我读过普通的 Moran 的 $I$ 不能处理这种数据，但是 package 中的 join count 方法spdep可以做到。但是，我是新手R，我仍然无法理解joincount.mcor帮助中的信息和代码joincount.test。

我的数据是这样的：

我知道如何读shp.file入R，我知道我必须计算我的数据的权重，但以下步骤spdep超出了我的能力。

我在具有 64 位操作系统和 7.8 GB RAM 的 PC 上使用 spdep 包的 errorsarlm() 函数。我的管理员分配的配额是 12 GB。

我的 errorsarlm() 的输入是：

1) 一个非对称权重矩阵，从 .gwt 转换而来（通过 nb2listw 函数；1.7 mb）。对于 85684 个区域（观察次数），它是 ak=4 最近邻矩阵：

2) 一个 246.3 mb 的 CSV 数据集，我在 errorsarlm() 模型中使用了 80 个变量。每个变量有 85684 个观测值。

当我基于 80 个变量运行简单线性回归 (lm) 时，我没有任何问题。但是，当我使用默认参数运行 errorsalm 模型时，我会立即收到以下消息：

traceback() 告诉我：

使用矩阵 nb 对象的 sum(card()) 进行矩阵密度检查，该对象是从带有 read.gwt2nb 的 .gwt 转换而来的，它告诉我：

由于错误消息的立即出现，内存使用量不会增加。是否有某种我已经超过的 spdep 允许的最大矩阵大小？还是有不同的解释？

当我在我的errorsarlm（）中包含method="LU"和时，模型运行良好并给出结果。method "MC"

希望有人可以帮助我理解这一点。

谢谢，

戴安娜

我正在尝试spdep在 OS X 10.9 Mavericks 上安装 R 包（），R 3.0.1（sessionInfo在底部）。该软件包是从http://cran.r-project.org/web/packages/spdep/index.html（源文件spdep_0.5-71.tar.gz）下载的，我尝试使用R CMD INSTALL.

问题是终端产生此错误：

我尝试重新安装 Xcode 和我已经拥有的命令行工具。而且我还在终端中安装了gfortran-4.8.2-darwin13.tar.bz2doing sudo tar fvxz gfortran-4.8.2-darwin13.tar.bz2 。但它仍然给我同样的错误。

有什么解决办法吗？

谢谢！

我有一个问题spdep()。从函数产生的非缺失距离矩阵开始

然后我们变成具有线性逆的权重

现在

我检查以确保没有丢失重量：

既然没有，请继续：

导致以下错误：

这是我第一次spdep在 R 中使用这个包，所以我希望有人可以帮助我做一些事情。

让我们考虑以下几点。

本质上，我想通过连接彼此最近的点来构建路径。例如，如果我们想要长度为 3 的路径并且我们从点 1 开始，那么我们应该接下来到点 11，然后最后到点 10。这是我手动获得的其他路径。

我正在尝试使用spdep包中的函数获取相同的信息。这是我得到的：

结果图几乎是我想要的，除了一些不应该连接的点被连接起来。例如，点 11 和 12 不应该连接，因为我们不能直接从点 11 到 12，反之亦然。我可能没有n2正确创建。我如何解决它？以及如何查看类似于我之前在矩阵中的邻居列表？

我想spdep::lagsarlm在 R 中估计一个模型（空间自回归）。我的观察结果（n=447）是多边形，每个代表柏林的一个行政区域。

然而，问题在于这些地区的居民数量差异很大（在 500 到 32000 之间）。因此，我想用它的居民数量来衡量每个观测值。这很容易，因为lm它接受可选参数weights=......

我怎样才能做类似的事情spdep::lagsarlm？有解决方法吗？