posterior predictive check for binomial glm with matrix response

posterior predictive check for binomial glm with matrix response #558