Suggested bug fix

[plewis]

In pglmm-utils.R, line 622 is

H <- Y - X %*% B

You need to add the following line before line 622

Y <- as.matrix(Y)

otherwise a "non-conformable array" error is issued due to the fact that Y is not a column vector.