gremlin (version 0.1.0.1)

anova.gremlin: anova() for gremlin objects

Description

REML Likelihood Ratio Tests for gremlin models using anova()

Usage

# S3 method for gremlin
anova(object, ..., model.names = NULL)

Arguments

object

An object of class ‘gremlin’.

…

Additional objects of class ‘gremlin’.

model.names

Optional character vector with model names to be used in the anova table

Value

A data.frame containing the nested comparison of model objects via a REML likelihood ratio test.

Examples

# NOT RUN {
mod11 <- gremlinR(WWG11 ~ sex - 1,
random = ~ calf,
data = Mrode11,
Gstart = matrix(0.1), Rstart = matrix(0.4),
maxit = 50, v = 2, algit = "EM")
logLik(mod11)
# }