I am new with mixed effect models and I need your help please. I have plotted the below graph in ggplot: ggplot(tempEf,aes(TRTYEAR,CO2effect,group=Myc,col=Myc)) + facet_grid(~N) + geom_smooth(m...

I used to use the code below to calculate standardized coefficients of a lmer model. However, with the new version of lme the structure of the returned object has changed. How to adapt the functi...

I would like to use try() or tryCatch() or a function like this to detect if there is an error in my model called "fit1". If the model is fine, I want to use "fit1", otherwise I want to use "fit2" ...

I fitted a model in R with the lmer()-function from the lme4 package. I scaled the dependent variable: mod <- lmer(scale(Y) ~ X + (X | Z), data = df, REML = FALSE) I look at the fixed-ef...

I'm using lmer() in package lme4 to estimate mixed effects models. This works well, but now I want to run the estimation process for a fixed number of iterations, then resume the process by specify...

I am currently testing whether I should include certain random effects in my lmer model or not. I use the anova function for that. My procedure so far is to fit the model with a function call to lm...

I'm trying to understand the function lmer. I've found plenty of information about how to use the command, but not much about what it's actually doing (save for some cryptic comments here: http://w...

I'm trying to use lmerTest to have p-values for my fixed effects. I have 4 different random intercepts, 3 crossed and one nested : test.reml <- lmerTest::lmer(y ~ s1 + min + cot + min:cot + ge ...

Consider a simple example: > library(lme4) Loading required package: lattice Loading required package: Matrix > set.seed(1) > df <- data.frame(x = c(rnorm(7), NA), y = rep(c("A", "B"),...

I just updated lme4 to version 1.0-4 and when I run lmer() my mixed effects model, which was converging before, now prints this warning: Warning message: In (function (fn, par, lower = rep.int(-In...

I am trying to run diagnostic plots on an lmer model but keep hitting a wall. I'm not sure how much information I need to provide here, but here goes: The model is simple: best <- lmer(MSV_mm...

Is it possible to plot the random intercept or slope of a mixed model when it has more than one predictor? With one predictor I would do like this: #generate one response, two predictors and one ...