Multivariate Generalized Linear Mixed Models for Ranking Sports Teams

Maximum likelihood estimates are obtained via an EM algorithm with either a first-order or a fully exponential Laplace approximation as documented by Broatch and Karl (2018) , Karl, Yang, and Lohr (2014) , and by Karl (2012) . Karl and Zimmerman use this package to illustrate how the home field effect estimator from a mixed model can be biased under nonrandom scheduling.