Re: [R] Specify multiple nested random effects in lme with heteroskedastic variance across group

2018-02-21 Thread Bert Gunter
You should post this on the r-sig-mixed-models list. You are more likely to get a helpful response there. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic

[R] Specify multiple nested random effects in lme with heteroskedastic variance across group

2018-02-20 Thread suraj keshri
I want to fit a random effects model with two separate nested random effects. I can easily do this using the `lmer` package in R. Here's how: model<-lmer(y ~ 1 + x + (1 | oid/gid) + (1 | did/gid), data=data) Here, I'm fitting a random intercept for `oid` nested within `gid` and `did` nested