Пространственный автокорреляционный тест для двоичных данных в R

Question

Я хочу проверить записи присутствия / отсутствия вида для пространственной автокорреляции. Мои данные содержат>130 000 сеток в ГИС и около 700 записей присутствия видов.

Я читал, что обычный $ I $ Морана не может иметь дело с такого рода данными, но что метод счетчика соединений в пакете spdep может сделать это. Тем не менее, я новичок в R и я до сих пор не могу понять информацию и код в справке для joincount.mc или же joincount.test,

Мои данные таковы:

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

Я умею читать shp.file в R и я знаю, что должен рассчитать вес моих данных, но следующие шаги с spdep вне моей способности.