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