Chapter 11 Linear Mixed-Effects Models
- Mixed-effects model
y=Xβ+Zu+e
- the elements of β are considered to be non-random and are called fixed effects.
- the elements of u are random variables and are called random effects.
- We assume that E(e)=0,Var(e)=R,E(u)=0,Var(u)=G,Cov(u,e)=0
- E(y)=Xβ, Var(y)=ZGZ′+R≡Σ, E(y∣u)=Xβ+Zu
When there are m random factors, we can partition Z and u as Z=[Z1,…,Zm] and u=[u′1,…,u′m]′ so that Zu=∑mj=1Zjuj.
- We often assume that all random effects are mutually independent and random effects associated with jth random factor have variance σ2j. Then, Var(y)=ZGZ′+R=∑Rj=1σ2jZjZ′j+σ2eI.
- The unknown variance parameters σ2j,σ2e are called variance components.
Experimental Design Terminology
- Experiment: An investigation in which the investigator applies some treatments to experimental units and then observes the effect of the treatments on the experimental units by measuring one or more response variables.
- Treatment: a condition or set of conditions applied to experimental units in an experiment.
- Experimental Unit: the physical entity to which a treatment is randomly assigned and independently applied.
- Response Variable: a characteristic of an experimental unit that is measured after treatment and analyzed to assess the effects of treatments on experimental units.
- Observational Unit: the unit on which a response variable is measured.
- Completely Randomized Design (CRD) – experimental design in which, for given number of experiment units per treatment, all possible assignments of treatments to experimental units are equally likely.
- Block – a group of experimental units that, prior to treatment, are expected to be more like one another (with respect to one or more response variables) than experimental units in general.
- Randomized Complete Block Design (RCBD) – experimental design in which separate and completely randomized treatment assignments are made for each of multiple blocks in such a way that all treatments have at least one experimental unit in each block.
Whenever an experiment involves multiple observations per experimental unit, it is important to include a random effect for each experimental unit.
Without a random effect for each experimental unit, a one-to-one correspondence between observations and experimental units is assumed.