思うところあって（というか、非線形混合効果モデルの式がなかなか把握できないので）、線形代数のあたりを復習してみたりとか。
プログラミングのための線形代数はやっぱりいい本だなぁと思う。幾何的なイメージつかみやすいし、QR分解とかLU分解とかもわかりやすい
プログラミングのための線形代数以下、Rで。
LU分解
library(Matrix)
x <- Matrix(rnorm(9), 3, 3)
elu <- expand(lu(x))
elu$L
elu$U
QR分解
x <- matrix(rnorm(9), 3, 3)
y <- qr(x)
qr.Q(y)
qr.R(y)