Log-Likelihood of an lme Object

DESCRIPTION:

If 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.

USAGE:

logLik(object, REML) 

REQUIRED ARGUMENTS:

object
an object inheriting from class lme, representing a fitted linear mixed-effects model.

OPTIONAL ARGUMENTS:

REML
an optional logical value. If TRUE the restricted log-likelihood is returned, else, if FALSE, the log-likelihood is returned. Defaults to FALSE.

VALUE:

the (restricted) log-likelihood of the linear mixed-effects model represented by object evaluated at the estimated coefficients.

REFERENCES:

Harville, D.A. (1974) "Bayesian Inference for Variance Components Using Only Error Contrasts", Biometrika, 61, 383-385.

SEE ALSO:

EXAMPLES:

fm1 <- lme(distance ~ Sex * age, Orthodont, random = ~ age, method = "ML") 
logLik(fm1) 
logLik(fm1, REML = TRUE)