Using Stan to model geostatistical count data with distance matrices

Switching to spatial regression modeling, I here show how to simulate discrete observations over a latent Gaussian Process spatial autocorrelation function. I then demonstrate how to use Stan to estimate parameters of the model using a spatial distance matrix as input.