##load data
data(mesa.data.res)
##Extract pre-computed cross-validated predictions
pred.cv <- mesa.data.res$pred.cv
##compute long term averages of predictions and observations
pred.lta <- compute.ltaCV(pred.cv)
##we can now compare observed and predicted averages at each site
plot(pred.lta[,"obs"],pred.lta[,"pred"])
abline(0,1)
Run the code above in your browser using DataLab