Plot lmer model with Effects package
Dear Marte, On Thu, 25 Oct 2012 16:11:07 +0200
Marte Lilleeng <mlilleeng at gmail.com> wrote:
Hi everyone! I have a simple model that i would like to plot with 95% CIs. It is like follows: m1<-lmer(Richness~Grazing+I(Grazing^2)+(1|Plot),family=poisson) By using the effects package I get two plots, one for the linear term and one for the squared term. Q1: Can I get all in one? I.e. with one line for the whole model?
It's better to specify the model as Richness~poly(Grazing, 2)+(1|Plot) so that effect() understands that the quadratic is a single term in the model. Alternatively, I believe that Effect() would work as you want for either form.
Q2: Can I also visualize the random effects?
I'm not sure what that would mean. I suppose that you could look at the distribution of the BLUPs. I hope this helps, John ------------------------------------------------ John Fox Sen. William McMaster Prof. of Social Statistics Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/
I would be very happy for your answers! -- Mvh Marte S. Lilleeng
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.