
I have a two level-model with a random intercept and two random slope components. I can fit this with lme and then do plot(ranef(model), layout=c(3,1,1)), and somewhere in the code, it recognises that the x-axis scales are different and does the plot accordingly, and as expected.

With lmer, things seem more complicated. I can do:

rr1 <- ranef(model, postVar=TRUE)
dotplot(rr1,scales = list(x = list(relation = 'free'), layout=c(3,1,1)))[["groupingfactor"]]

However it seems like the explicit specification of relation=free then causes layout to be ignored, generally leading to a 2x2 panel plot, which causes too much bunching in the groups to be able to read them, and generally not using the space on the page in an optimal manner.

It is possible to go for
dotplot(rr1 , layout=c(3,1,1)))[["groupingfactor"]]
and get the layout right but the x-axis scales are then equal for each panel, which will force the random slopes to be plotted on the same scale as the random intercepts :-(

Has anyone else also had these problems, and perhaps found a solution?

regards

Mike Dunbar




-- 
This message (and any attachments) is for the recipient ...{{dropped:10}}


