logLik.lme {nlme} | R Documentation

## 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

    ## S3 method for class 'lme':
    logLik(object, REML, ...)

### Arguments

`object` | an object inheriting from class `lme`, representing a fitted linear mixed-effects model.
`REML` | an optional logical value. If `TRUE` the restricted log-likelihood is returned, else, if `FALSE`, the log-likelihood is returned. Defaults to `FALSE`.
`...` | some methods for this generic require additional arguments. None are used in this method.

### Value

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

### Author(s)

Jose Pinheiro and Douglas Bates

### References

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

### See Also

`lme`,`gls`

### Examples

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

[Package Contents]