[R-sig-ME] some questions about "lmer" in GLMM
David Duffy
David.Duffy at qimr.edu.au
Wed May 9 01:32:23 CEST 2012
On Tue, 8 May 2012, Maddie Wang wrote:
Hi.
> I have a question about lmer for GLMM that I cannot find an answer
> on-line for quite a bit time. Suppose the response y is from a binomial
> distribution. The GLMM can be written as
> logit(pi)=x*\beta+z1b1+z2b2+z3b3, where pi is the expectation of
> binomial proportion, \beta is the fixed coefficients, b1,b2 and b3 are
> three random effects, where b1 is from a multivariate normal MN(0, \tau1
> * I_(3X3)), b2 is from a multivariate normal MN(0, \tau2 * I_(3X3)), and
> b3 is from a multivariate normal MN(0, \tau3 * I_(4X4)).
I'm having trouble understanding your setup. Why is b3 4X4 - are these
crossed indicator variables with 3,3,4 levels? And if so, what design and
how much data? Are the taus prespecified somehow?
Cheers, David Duffy.
More information about the R-sig-mixed-models
mailing list