I made a semivariogram in R using package gstat,
variogram() function. I want to see if there is spatial autocorrelation in the residuals of my model (species abundance as a function of habitat, across sites spaced a few km to 900 km apart, using a glmm).
My units are in km, and so my interpretation is that the range is just over 100 km until spatial autocorrelation is no longer an "issue". I am wondering if someone can explain why the nugget seems so high? Does this mean that even at similar locations, there is still a relatively high degree of difference? Or, does this wavy variogram mean I should adjust my number of lags and lag distance until I get a more typical shape?
In order to investigate a bit further, I also used the function
variog() in package geoR, and used
breaks=seq(0,100,10), to try and look just at the closer distances (using same points and same model residuals). This one indicates that the closest points are more different, which also doesn't make sense. Maybe this indicates there is no spatial autocorrelation, and that my model accounts for this already.
I found this excellent source, "Geostats without tears", and on page 51 has some good advice on fitting variograms. By this advice, my first one seems to have the correct range. So this goes back to the first question - how do I interpret this?