comodex: Model comparison

In pbkrtest: Parametric Bootstrap, Kenward-Roger and Satterthwaite Based Methods for Test in Mixed Models

Model comparison

Description

Wrapper for functions KRmodcomp, SATmodcomp, PBmodcomp, X2modcomp

Usage

comodex(
  largeModel,
  smallModel,
  test = "x2",
  control = list(),
  details = 0,
  ...
)

## S3 method for class 'lmerMod'
comodex(
  largeModel,
  smallModel,
  test = "x2",
  control = list(),
  details = 0,
  ...
)

## Default S3 method:
comodex(
  largeModel,
  smallModel,
  test = "x2",
  control = list(),
  details = 0,
  ...
)

Arguments

largeModel	A model object
smallModel	A model object, a formula or a restriction matrix
test	A character string
control	A list controlling the model comparions.
details	should details be printed
...	Additional arguments to be passed on to other methods

Author(s)

Søren Højsgaard

Examples

(lmm0 <- lmer(Reaction ~ (Days|Subject), sleepstudy))
(lmm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
(lmm2 <- lmer(Reaction ~ Days + I(Days^2) + (Days|Subject), sleepstudy))

lm1 <- lm(dist ~ speed + I(speed^2), data=cars)
lm0 <- lm(dist ~ speed, data=cars)

comodex(lmm2, lmm1, test="x2")
comodex(lmm2, lmm1, test="kr")
comodex(lmm2, lmm1, test="sat")
comodex(lmm2, lmm1, test="PB", control=list(nsim=50, cl=1))
comodex(lmm2, .~. - I(Days^2))

comodex(lm1, lm0)
comodex(lm1, lm0, test="pb", control=list(nsim=50, cl=1))

pbkrtest documentation built on June 8, 2025, 11:20 a.m.