REML=FALSE
, returns the log-likelihood value of the linear
mixed-effects model represented by
object
evaluated at the
estimated coefficients; else, the restricted log-likelihood evaluated
at the estimated coefficients is returned.
logLik(object, REML)
lme
, representing
a fitted linear mixed-effects model.
TRUE
the restricted
log-likelihood is returned, else, if
FALSE
, the log-likelihood
is returned. Defaults to
FALSE
.
object
evaluated at the estimated coefficients.
Harville, D.A. (1974) "Bayesian Inference for Variance Components Using Only Error Contrasts", Biometrika, 61, 383-385.
fm1 <- lme(distance ~ Sex * age, Orthodont, random = ~ age, method = "ML") logLik(fm1) logLik(fm1, REML = TRUE)