# Load data
data(met)
# Response
Y = Obs[,6]
# ID, covariate (Year) and spatial parameters (Block-Row-Col)
DTA = Obs[,1:5]
# Fit model
fit = gmm(y=Y,gen=Gen,dta=DTA)
# Summary
cor(fit$hat,fit$obs)
# Plot
plot(fit$hat,fit$obs,col=as.numeric(Obs$Year)+1,pch=20)
legend('topleft',legend=2013:2015,col=2:4,pch=20)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.