added posterior predictive checks

This allows us to check model fit with Bayesian p-values and to
calculate the sums of square residuals. Right now it is only for
overall cover estimates; could be broken up by regression like in Adler
et al. 2012.