[R] Gibbs sampling

[Emma Tan]

Hi,

I am trying to use the glmm() function in the
GLMMGibbs package to fit a crossed random effects
model.  I want to include two random effects and I'm
not sure how to go about doing this.  Do I include two
Ra() functions for the two random effects I wish to
include?  I did try this, adding both random
components to the glmm() function but only seemed to
get an estimate of one random error component and not
both.

Thanks for any help you can give!

Emma