x2__modcomp {pbkrtest} | R Documentation |
Chisq test
Description
Chisq test
Usage
X2modcomp(largeModel, smallModel, betaH = 0, details = 0, ...)
## Default S3 method:
X2modcomp(largeModel, smallModel, betaH = 0, details = 0, ...)
Arguments
largeModel |
An |
smallModel |
An |
betaH |
A number or a vector of the beta of the hypothesis,
e.g. L beta=L betaH. If |
details |
If larger than 0 some timing details are printed. |
... |
Additional arguments, currently not used. |
Details
TBW
Author(s)
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))
Test for no effect of Days in fm1, i.e. test fm0 under fm1
X2modcomp(fm1, "Days") X2modcomp(fm1, ~.-Days) L1 <- cbind(0, 1)
X2modcomp(fm1, L1) ## FIXME
X2modcomp(fm1, fm0) anova(fm1, fm0)
[Package pbkrtest version 0.5.4 Index]