Spatial random variable with 9,738 coordinates causes R to crash about hmsc

[gtikhonov]

You have absolutely too many spatial locations in your model for vanilla Gaussian Process (GP) -priored latent factors to work well. Try formulations of HmscRandomLevel(...) that rely on approximate GP structures. For instance,

rL.spatial = HmscRandomLevel(sData = xycoords, sMethod ="NNGP")

See the help page of HmscRandomLevel(...) for more details or this paper (https://doi.org/10.1002/ecy.2929) for statistical presentation of implemented methods.

[devilevie]

Thanks, that's really helpful. I'll try the new code and have a look at the paper, hopefully I can get it working.