# Which model to choose - lmer vs. glmer

I am trying to model y, a continuous variable that only takes positive values with fixed and random effects.

This is my first approach, using lmer():

log(y) ~ x1 + x2 + x3 + (1|plot) + (1|id)


The fit is not perfect but the residuals look roughly OK I guess.

In my second approach I am using glmer():

y ~ x1 + x2 + x3 + (1|plot) + (1|id)


with family = Gamma(link = "inverse")

The fit gets a bit better. But I am unsure if using the much more complex glmer is justified or if I should stick to lmer? On what grounds should I decide with which model to go?

Edit Simulated Residuals for the GLMM: Simulated Residuals for the LMM: • I think we need to know more about the problem. What is it you are modelling? Jun 26 '19 at 16:06
• I am modelling an indicator for plant growth, which is calculated as a ratio of two measurements. Jun 26 '19 at 16:11
• Ratios can be tricky and sometimes are better with the numerator and denominator modeled separately. Following up on the comment from @DemetriPananos, please say more about the data that contribute to calculating that ratio. Note that both models above show much lower fitted than observed values at the higher values, and the cloud of indistinguishable points seems correspondingly to have a slope much lower than 1. A standard plot of residuals against fitted values would be more informative. A normal distribution of errors doesn't help much if the basic fit isn't good.
– EdM
Jun 26 '19 at 16:39
• @EdM I have added the plot of residuals against fitted values. Jun 27 '19 at 5:18
• First of all, what is the goal of your modeling ? Predictive or Explanatory ? Do you want your coefficient estimates to make sense or just having the better prediction ? When you have that, you will have to select the criteria: either some error criteria for predictive or criteria like AIC, BIC, p-value, etc. for explanatory. After that, you can answer your question properly, you choose the model that minimize your chosen criteria. Jul 2 '19 at 0:38

However, you need to be careful with their use if you have missing at random missing data in your outcome variable y; for more on this check here.