我想测试物种的空间自相关存在/不存在记录。我的数据在GIS中包含> 130,000个网格，并且有大约700种物种的存在记录。
我已经读过正常Moran的$ I $无法处理这类数据，但是包spdep中的连接计数方法可以做到这一点。但是，我是R的新手，我仍然无法理解joincount.mc或joincount.test的帮助中的信息和代码。
我的数据是这样的：
gridnumber species
1 1
2 0
3 0
4 1
……
我知道如何将shp.file读入R，我知道我必须计算数据的权重，但spdep的以下步骤超出了我的能力范围。