種の存在/不在記録の空間的自己相関をテストしたいと考えています。私のデータには、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
は私の能力を超えています。