A mixed model, mixed-effects model or mixed error-component model is a statistical model containing both fixed effects and random effects.^[1][2] These models are useful in a wide variety of disciplines in the physical, biological and social sciences. They are particularly useful in settings where repeated measurements are made on the same statistical units (longitudinal study), or where measurements are made on clusters of related statistical units.^[2] Mixed models are often preferred over traditional analysis of variance regression models because of their flexibility in dealing with missing values and uneven spacing of repeated measurements.^[3]The Mixed model analysis allows measurements to be explicitly modeled in a wider variety of correlation and variance-covariance structures.

This page will discuss mainly linear mixed-effects models (LMEM) rather than generalized linear mixed models or nonlinear mixed-effects models. ^[4]

Qualitative Description

Linear mixed models (LMMs) are statistical models that incorporate fixed and random effects to accurately represent non-independent data structures. LMM is an alternative to analysis of variance. Often, ANOVA assumes the independence of observations within each group, however, this assumption may not hold in non-independent data, such as multilevel/hierarchical, longitudinal, or correlated datasets.

Non-independent sets are ones in which the variability between outcomes is due to correlations within groups or between groups. Mixed models properly account for nest structures/hierarchical data structures where observations are influenced by their nested associations. For example, when studying education methods involving multiple schools, there are multiple levels of variables to consider. The individual level/lower level comprises individual students or teachers within the school. The observations obtained from this student/teacher is nested within their school. For example, Student A is a unit within the School A. The next higher level is the school. At the higher level, the school contains multiple individual students and teachers. The school level influences the observations obtained from the students and teachers. For Example, School A and School B are the higher levels each with its set of Student A and Student B respectively. This represents a hierarchical data scheme. A solution to modeling hierarchical data is using linear mixed models.

Representation of how data, related to education system, is non-independent and structured in nested/hierarchial levels.

LMMs allow us to understand the important effects between and within levels while incorporating the corrections for standard errors for non-independence embedded in the data structure.^[4][5]

History and current status

Representation of biased vs. unbiased data and the differences between LSR and LMM fitted estimates.

Ronald Fisher introduced random effects models to study the correlations of trait values between relatives.^[9] In the 1950s, Charles Roy Henderson provided best linear unbiased estimates of fixed effects and best linear unbiased predictions of random effects.^{[10][11][12][13]} Subsequently, mixed modeling has become a major area of statistical research, including work on computation of maximum likelihood estimates, non-linear mixed effects models, missing data in mixed effects models, and Bayesian estimation of mixed effects models. Mixed models are applied in many disciplines where multiple correlated measurements are made on each unit of interest. They are prominently used in research involving human and animal subjects in fields ranging from genetics to marketing, and have also been used in baseball ^[14] and industrial statistics.^[15] The mixed linear model association has improved the prevention of false positive associations. Populations are deeply interconnected and the relatedness structure of population dynamics is extremely difficult to model without the use of mixed models. Linear mixed models may not, however, be the only solution. LMM's have a constant-residual variance assumption that is sometimes violated when accounting or deeply associated continuous and binary traits. ^[16]

Definition

In matrix notation a linear mixed model can be represented as

${\displaystyle {\boldsymbol {y}}=X{\boldsymbol {\beta }}+Z{\boldsymbol {u}}+{\boldsymbol {\epsilon }}}$

where

• ${\displaystyle {\boldsymbol {y}}}$ is a known vector of observations, with mean ${\displaystyle E({\boldsymbol {y}})=X{\boldsymbol {\beta }}}$;
• ${\displaystyle {\boldsymbol {\beta }}}$ is an unknown vector of fixed effects;
• ${\displaystyle {\boldsymbol {u}}}$ is an unknown vector of random effects, with mean ${\displaystyle E({\boldsymbol {u}})={\boldsymbol {0}}}$ and variance–covariance matrix ${\displaystyle \operatorname {var} ({\boldsymbol {u}})=G}$;
• ${\displaystyle {\boldsymbol {\epsilon }}}$ is an unknown vector of random errors, with mean ${\displaystyle E({\boldsymbol {\epsilon }})={\boldsymbol {0}}}$ and variance ${\displaystyle \operatorname {var} ({\boldsymbol {\epsilon }})=R}$;
• ${\displaystyle X}$ is the known design matrix for the fixed effects relating the observations ${\displaystyle {\boldsymbol {y}}}$ to ${\displaystyle {\boldsymbol {\beta }}}$, respectively
• ${\displaystyle Z}$ is the known design matrix for the random effects relating the observations ${\displaystyle {\boldsymbol {y}}}$ to ${\displaystyle {\boldsymbol {u}}}$, respectively.

Estimation

The joint density of ${\displaystyle {\boldsymbol {y}}}$ and ${\displaystyle {\boldsymbol {u}}}$ can be written as: ${\displaystyle f({\boldsymbol {y}},{\boldsymbol {u}})=f({\boldsymbol {y}}|{\boldsymbol {u}})\,f({\boldsymbol {u}})}$. Assuming normality, ${\displaystyle {\boldsymbol {u}}\sim {\mathcal {N}}({\boldsymbol {0}},G)}$, ${\displaystyle {\boldsymbol {\epsilon }}\sim {\mathcal {N}}({\boldsymbol {0}},R)}$ and ${\displaystyle \mathrm {Cov} ({\boldsymbol {u}},{\boldsymbol {\epsilon }})={\boldsymbol {0}}}$, and maximizing the joint density over ${\displaystyle {\boldsymbol {\beta }}}$ and ${\displaystyle {\boldsymbol {u}}}$, gives Henderson's "mixed model equations" (MME) for linear mixed models:^[10][12][17]

${\displaystyle {\begin{pmatrix}X'R^{-1}X&X'R^{-1}Z\\Z'R^{-1}X&Z'R^{-1}Z+G^{-1}\end{pmatrix}}{\begin{pmatrix}{\hat {\boldsymbol {\beta }}}\\{\hat {\boldsymbol {u}}}\end{pmatrix}}={\begin{pmatrix}X'R^{-1}{\boldsymbol {y}}\\Z'R^{-1}{\boldsymbol {y}}\end{pmatrix}}}$

The solutions to the MME, ${\displaystyle \textstyle {\hat {\boldsymbol {\beta }}}}$ and ${\displaystyle \textstyle {\hat {\boldsymbol {u}}}}$ are best linear unbiased estimates and predictors for ${\displaystyle {\boldsymbol {\beta }}}$ and ${\displaystyle {\boldsymbol {u}}}$, respectively. This is a consequence of the Gauss–Markov theorem when the conditional variance of the outcome is not scalable to the identity matrix. When the conditional variance is known, then the inverse variance weighted least squares estimate is best linear unbiased estimates. However, the conditional variance is rarely, if ever, known. So it is desirable to jointly estimate the variance and weighted parameter estimates when solving MMEs.

One method used to fit such mixed models is that of the expectation–maximization algorithm (EM) where the variance components are treated as unobserved nuisance parameters in the joint likelihood.^[18] Currently, this is the method implemented in statistical software such as Python (statsmodels package) and SAS (proc mixed), and as initial step only in R's nlme package lme(). The solution to the mixed model equations is a maximum likelihood estimate when the distribution of the errors is normal.^[19][20]

Fixed, mixed, and random effects influence linear regression models.

There are several other methods to fit mixed models, including using an MEM initially, and then Newton-Raphson (used by R package nlme^[21]'s lme()), penalized least squares to get a profiled log likelihood only depending on the (low-dimensional) variance-covariance parameters of ${\displaystyle {\boldsymbol {u}}}$, i.e., its cov matrix ${\displaystyle {\boldsymbol {G}}}$, and then modern direct optimization for that reduced objective function (used by R's lme4^[22] package lmer() and the Julia package MixedModels.jl) and direct optimization of the likelihood (used by e.g. R's glmmTMB). Notably, while the canonical form proposed by Henderson is useful for theory, many popular software packages use a different formulation for numerical computation in order to take advantage of sparse matrix methods (e.g. lme4 and MixedModels.jl).

See also

• Nonlinear mixed-effects model
• Fixed effects model
• Generalized linear mixed model
• Linear regression
• Mixed-design analysis of variance
• Multilevel model
• Random effects model
• Repeated measures design
• Empirical Bayes method

References

Further reading

• Gałecki, Andrzej; Burzykowski, Tomasz (2013). Linear Mixed-Effects Models Using R: A Step-by-Step Approach. New York: Springer. ISBN 978-1-4614-3900-4.
• Milliken, G. A.; Johnson, D. E. (1992). Analysis of Messy Data: Vol. I. Designed Experiments. New York: Chapman & Hall.
• West, B. T.; Welch, K. B.; Galecki, A. T. (2007). Linear Mixed Models: A Practical Guide Using Statistical Software. New York: Chapman & Hall/CRC.