Hello,

I am working with a GLM where my random effects are factor variables f1 and
f2. I am new to MCMCglmm, so i'm not sure that I have the following syntax
correct:

mod1 <- MCMCglmm(outcome ~ X1 + X2, random=~f1 + f2, family="categorical",
...)

When I look at mod1$VCV I only see three columns. One for f1, one for f2 and
"units". This makes me think that f1 and f2 aren't being treated as factors?
is(f1) does show it's a factor, and trying ~as.factor(f1) doesn't seem to
help (Syntax error).

How can I 1) confirm that it's treating f1 and f2 as factors (with many
levels) and 2) get estimates of the random effects out for those factors?

Thanks,
-Robin

	[[alternative HTML version deleted]]


