Ho fatto un semivariogramma in R usando il pacchetto gstat, variogram()funzione. Voglio vedere se c'è un'autocorrelazione spaziale nei residui del mio modello (l'abbondanza di specie in funzione dell'habitat, attraverso siti distanziati a pochi km a 900 km di distanza, usando un glmm).
Le mie unità sono in km, e quindi la mia interpretazione è che la distanza è di poco più di 100 km fino a quando l'autocorrelazione spaziale non è più un "problema". Mi chiedo se qualcuno può spiegare perché la pepita sembra così alta? Questo significa che anche in luoghi simili, c'è ancora un grado relativamente elevato di differenza? Oppure, questo variogramma ondulato significa che dovrei regolare il mio numero di ritardi e la distanza di ritardo fino a ottenere una forma più tipica?
Per approfondire un po ', ho anche usato la funzione variog()nel pacchetto geoR, e ho usato breaks=seq(0,100,10), per cercare di guardare solo le distanze più vicine (usando gli stessi punti e gli stessi residui del modello). Questo indica che i punti più vicini sono più diversi, il che non ha senso. Forse questo indica che non c'è autocorrelazione spaziale e che il mio modello già tiene conto di ciò.
Ho trovato questa eccellente fonte, "Geostati senza lacrime" , e a pagina 51 ha alcuni buoni consigli sull'adattamento di vari diagrammi . Con questo consiglio, il mio primo sembra avere la gamma corretta. Quindi questo risale alla prima domanda: come posso interpretarlo?