CoChat
Sign Up
Cited by 2,094Statistical Methods and Bayesian InferenceAgriculture, Soil, Plant ScienceStatistical Methods and Applications

The DOI matches its Crossref record; 92 of the 92 references it cites were found in a registry or an index.

A brief introduction to mixed effects modelling and multi-model inference in ecology

Open access
Harrison, Xavier A., Donaldson, Lynda, Correa-Cano, Maria Eugenia, Evans, Julian, Fisher, David N., Goodwin, Cecily E.D., Robinson, Beth S., Hodgson, David J., Inger, Richard·2018·PeerJ
DOI 10.7717/peerj.4794PMC PMC5970551
Open this paper in CoChat — free account

Abstract

The use of linear mixed effects models (LMMs) is increasingly common in the analysis of biological data. Whilst LMMs offer a flexible approach to modelling a broad range of data types, ecological data are often complex and require complex model structures, and the fitting and interpretation of such models is not always straightforward. The ability to achieve robust biological inference requires that practitioners know how and when to apply these tools. Here, we provide a general overview of current methods for the application of LMMs to biological data, and highlight the typical pitfalls that can be encountered in the statistical modelling process. We tackle several issues regarding methods of model selection, with particular reference to the use of information theory and multi-model inference in ecology. We offer practical solutions and direct the reader to key references that provide further technical detail for those seeking a deeper understanding. This overview should serve as a widely accessible code of best practice for applying LMMs to complex biological problems and model structures, and in doing so improve the robustness of conclusions drawn from studies investigating ecological and evolutionary questions.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.7717/peerj.4794

Cites

92 of 92 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Differences between Random Intercept vs Random Slope Models. (A) A random-intercepts model where the outcome variable y is a function of predictor x, with a random intercept for group ID (coloured lines). Because all groups have been constrained to have a common slope, their regression lines are parallel. Solid lines are the regression lines fitted to the data. Dashed lines trace the regression lines back to the y intercept. Point colour corresponds to group ID of the data point. The black line represents the global mean value of the distribution of random effects. (B) A random intercepts and random slopes model, where both intercepts and slopes are permitted to vary by group. Random slope models give the model far more flexibility to fit the data, but require a lot more data to obtain accurate estimates of separate slopes for each group.
Figure 1 Differences between Random Intercept vs Random Slope Models. (A) A random-intercepts model where the outcome variable y is a function of predictor x, with a random intercept for group ID (coloured lines). Because all groups have been constrained to have a common slope, their regression lines are parallel. Solid lines are the regression lines fitted to the data. Dashed lines trace the regression lines back to the y intercept. Point colour corresponds to group ID of the data point. The black line represents the global mean value of the distribution of random effects. (B) A random intercepts and random slopes model, where both intercepts and slopes are permitted to vary by group. Random slope models give the model far more flexibility to fit the data, but require a lot more data to obtain accurate estimates of separate slopes for each group.CC BY · PubMed Central
Figure 2 The effect of collinearity on model parameter estimates. We simulated 10,000 iterations of a model y ∼ x1 + x2, where x1 had a positive effect on y (βx1 = 1, vertical dashed line). x2 is collinear with x1 with either a moderate (r = 0.5). (A) or strong correlation (r = 0.9). (B) With moderate collinearity, estimation of βx1 is precise, but certainty of the sign of βx2 is low. When collinearity is strong, estimation of βx1 is far less precise, with 14% of simulations estimating a negative coefficient for the effect of x1. For more elaborate versions of these simulations, see Freckleton (2011).
Figure 2 The effect of collinearity on model parameter estimates. We simulated 10,000 iterations of a model y ∼ x1 + x2, where x1 had a positive effect on y (βx1 = 1, vertical dashed line). x2 is collinear with x1 with either a moderate (r = 0.5). (A) or strong correlation (r = 0.9). (B) With moderate collinearity, estimation of βx1 is precise, but certainty of the sign of βx2 is low. When collinearity is strong, estimation of βx1 is far less precise, with 14% of simulations estimating a negative coefficient for the effect of x1. For more elaborate versions of these simulations, see Freckleton (2011).CC BY · PubMed Central
Figure 3 Using Simulation to Assess Model Fit for GLMMs. (A) Histogram of the proportion of zeroes in 10,000 datasets simulated from a Poisson GLMM. Vertical red line shows the proportion of zeroes in our real dataset. There is no strong evidence of zero-inflation for these data. (B) Histogram of the sum of squared Pearson residuals for 1,000 parametric bootstraps where the Poisson GLMM has been re-fitted to the data at each step. Vertical red line shows the test statistic for the original model, which lies well outside the simulated frequency distribution. The ratio of the real statistic to the simulated data can be used to calculate a mean dispersion statistic and 95% confidence intervals, which for these data is mean 3.16, (95% CI [2.77–3.59]). Simulating from models provides a simple yet powerful set of tools for assessing model fit and robustness.
Figure 3 Using Simulation to Assess Model Fit for GLMMs. (A) Histogram of the proportion of zeroes in 10,000 datasets simulated from a Poisson GLMM. Vertical red line shows the proportion of zeroes in our real dataset. There is no strong evidence of zero-inflation for these data. (B) Histogram of the sum of squared Pearson residuals for 1,000 parametric bootstraps where the Poisson GLMM has been re-fitted to the data at each step. Vertical red line shows the test statistic for the original model, which lies well outside the simulated frequency distribution. The ratio of the real statistic to the simulated data can be used to calculate a mean dispersion statistic and 95% confidence intervals, which for these data is mean 3.16, (95% CI [2.77–3.59]). Simulating from models provides a simple yet powerful set of tools for assessing model fit and robustness.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Research works printed but not linked: 5

  1. Cox, 1989, The Analysis of Binary Data
  2. Crawley, 2013, The R Book
  3. Noble, 2017, Planned missing data design: stronger inferences increased research efficiency and improved animal welfare in ecology and evolution, bioRxiv, 247064, 10.1101/247064 · DOI
  4. Peters, 1991, A Critique for Ecology
  5. Southwood, 2000, Ecological Methods

Cited by

Sorted by citations

See all 2,094

Related

See all

Metadata

Authors
Harrison, X. A., Donaldson, L., Correa-Cano, M. E., Evans, J., Fisher, D. N., Goodwin, C. E. D., Robinson, B. S., Hodgson, D. J., & Inger, R.
Affiliations
  • Zoological Society of London (GB)
  • University of Exeter (GB)
  • Environment and Sustainability Institute, University of Exeter, Penryn, UK
  • University of Ottawa (CA)
Corresponding author
Harrison, Xavier A.
ORCID
Venue
PeerJ
Publisher
PeerJ
Published
2018-05-23
Volume
6
Pages
e4794
ISSN
2167-8359 (ISSN-L)
DOI
10.7717/peerj.4794
PMC
PMC5970551
OpenAlex
W2783564415
Cited by
2,094
References
91
Citations by year
2026: 227 · 2025: 228 · 2024: 271 · 2023: 246 · 2022: 327 · 2021: 383 · 2020: 261 · 2019: 116 · 2018: 32
FWCI
164.31
Type
journal article
Language
en
Open access
Licence
Locations
  • PeerJ (publishedVersion, cc-by)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Aberdeen University Research Archive (Aberdeen University) (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
Funders
  • Institute of Zoology Research Fellowship
  • NERC studentship (NE/H02249X/1, NE/L501669/1)
  • University of Exeter and the Animal and Plant Health as part of ‘Wildlife Research Co-Operative’
  • CONACYT (The Mexican National Council for Science and Technology)
  • SEP (The Mexican Ministry of Education)
  • Forestry Commission
Topics
  • Physical Sciences › Mathematics › Statistics and Probability › Statistical Methods and Bayesian Inference
  • Life Sciences › Agricultural and Biological Sciences › Agronomy and Crop Science › Agriculture, Soil, Plant Science
  • Physical Sciences › Mathematics › Statistics and Probability › Statistical Methods and Applications
Keywords
linear mixed models, multi-model inference, model selection, information theory, ecological data analysis, statistical pitfalls

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review