Chisq test
X2modcomp(largeModel, smallModel, betaH = 0, details = 0, ...)# S3 method for default
X2modcomp(largeModel, smallModel, betaH = 0, details = 0, ...)
An lmer model
An lmer model or a restriction matrix
A number or a vector of the beta of the hypothesis,
e.g. L beta=L betaH. If smallModel is a model object then betaH=0.
If larger than 0 some timing details are printed.
Additional arguments, currently not used.
Ulrich Halekoh uhalekoh@health.sdu.dk, Søren Højsgaard sorenh@math.aau.dk
(fm0 <- lmer(Reaction ~ (Days|Subject), sleepstudy)) (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) (fm2 <- lmer(Reaction ~ Days + I(Days^2) + (Days|Subject), sleepstudy))
X2modcomp(fm1, "Days") X2modcomp(fm1, ~.-Days) L1 <- cbind(0, 1)
X2modcomp(fm1, fm0) anova(fm1, fm0)
TBW