logLik.lme {nlme} | R Documentation |
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.
## S3 method for class 'lme': logLik(object, REML, ...)
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 the method of estimation used, that is
TRUE if and only object was fitted with method =
"REML" (the default for these fitting functions) . |
... |
some methods for this generic require additional arguments. None are used in this method. |
the (restricted) log-likelihood of the model represented by
object
evaluated at the estimated coefficients.
Jose Pinheiro and Douglas Bates
Harville, D.A. (1974) “Bayesian Inference for Variance Components Using Only Error Contrasts”, Biometrika, 61, 383–385.
Pinheiro, J.C., and Bates, D.M. (2000) "Mixed-Effects Models in S and S-PLUS", Springer.
lme
,gls
,
logLik.corStruct
,
logLik.glsStruct
,
logLik.lmeStruct
,
logLik.lmList
,
logLik.reStruct
,
logLik.varFunc
,
fm1 <- lme(distance ~ Sex * age, Orthodont, random = ~ age, method = "ML") logLik(fm1) logLik(fm1, REML = TRUE)