Implements Expectation/Conditional Maximization Either (ECME), rapidly converging algorithms, and Bayesian inference for linear mixed models following Schafer (1998), "Some Improved Procedures for Linear Mixed Models", Department of Statistics, The Pennsylvania State University.