Obtain the Restricted Maximum Log-Likelihood of a model fitted using LMMsolve.
Usage
# S3 method for LMMsolve
logLik(object, includeConstant = TRUE, ...)
Arguments
- object
an object of class LMMsolve
- includeConstant
Should the constant in the restricted log-likelihood be included. Default is
TRUE
, as for example inlme4
and SAS. Inasreml
the constant is omitted.- ...
some methods for this generic require additional arguments. None are used in this method.
Examples
## Fit model on john.alpha data from agridat package.
data(john.alpha, package = "agridat")
## Fit simple model with only fixed effects.
LMM1 <- LMMsolve(fixed = yield ~ rep + gen,
data = john.alpha)
## Obtain log-likelihood.
logLik(LMM1)
#> [1] -34.95557
## Obtain log-likelihood without constant.
logLik(LMM1, includeConstant = FALSE)
#> [1] 7.315605