getDiag function

Extract posterior summaries of random effects