[ ] ~need to remove dependency on BlockDiagonals for compatibility with MixedModels 4.0~ looks like this isn't necessary and I can't put off the dependency reduction for later
[x] I need to document rewriting the Z indicator matrices as a single orthogonal matrix
This doesn't quite work yet