## Fit model on oats data
data(oats.data)
## Fit simple model with only fixed effects.
LMM1 <- LMMsolve(fixed = yield ~ rep + gen,
data = oats.data)
## Obtain log-likelihood.
logLik(LMM1)
## Obtain log-likelihood without constant.
logLik(LMM1, includeConstant = FALSE)
Run the code above in your browser using DataLab