Deviance of an LMMsolve object — deviance.LMMsolve • LMMsolver

Obtain the deviance of a model fitted using LMMsolve.

Usage

# S3 method for LMMsolve
deviance(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 in lme4 and SAS. In asreml the constant is omitted.
...: some methods for this generic require additional arguments. None are used in this method.

Value

The deviance of the fitted model.

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 deviance.
logLik(LMM1)
#> [1] -34.95557

## Obtain deviance. without constant.
logLik(LMM1, includeConstant = FALSE)
#> [1] 7.315605