Classical LMM using r() around random effects, e.g., y ~ x * r(z). Support for various coding schemes. Sums-of-squares type I/II/III. Support for REML/ML.