Fit a spatio-temporal generalised linear mixed model to data, with a common temporal main effect and separate spatial surfaces with individual variances.
Useful links