randef function

Compute posterior estimates of random effect