Function to obtain restricted log-likelihood and the first derivatives of the log-likelihood, given values for the penalty parameters
mLogLik(object, theta)
A data.frame with logL and the first derivatives of log-likelihood
an object of class LMMsolve
a matrix with values of precision parameters theta.