r - R中二进制数据的空间自相关测试

Question

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

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

我的数据是这样的：

 gridnumber  species  
     1       1     
     2       0
     3       0
     4       1
     ……      

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