Solves a univariate mixed model of form y=X\beta+Mu+e
y=X\beta+Mu+e
solveRRBLUP(y, X, M)
y
a matrix with n rows and 1 column
X
a matrix with n rows and x columns
M
a matrix with n rows and m columns