1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
\name{anova.clm}
%%\alias{anova}
\alias{anova.clm}
\title{ANODE Tables and Likelihood ratio test of cumulative link models}
\description{
Type I, II, and III analysis of deviance (ANODE) tables for
cumulative link models and
comparison of cumulative link models with likelihood ratio tests.
Models may differ by terms in location, scale and nominal
formulae, in link, threshold function.
}
\usage{
\method{anova}{clm}(object, ..., type = c("I", "II", "III", "1", "2", "3"))
}
\arguments{
\item{object}{a \code{\link{clm}} object.
}
\item{\dots}{optionally one or more additional \code{\link{clm}} objects.
}
\item{type}{the type of hypothesis test if \code{anova} is called with a
single model; ignored if more than one model is passed to the method.
}
}
\details{
The ANODE table returned when \code{anova} is called with a single model apply only to
terms in \code{formula}, that is, terms in \code{nominal} and \code{scale} are
ignored.
}
\value{
An analysis of deviance table based on Wald chi-square test if called with a
single model and a comparison of
models with likelihood ratio tests if called with more than one model.
}
\author{Rune Haubo B Christensen}
\seealso{
\code{\link[ordinal]{clm}}
}
\examples{
## Analysis of deviance tables with Wald chi-square tests:
fm <- clm(rating ~ temp * contact, scale=~contact, data=wine)
anova(fm, type="I")
anova(fm, type="II")
anova(fm, type="III")
options(contrasts = c("contr.treatment", "contr.poly"))
m1 <- clm2(SURENESS ~ PROD, scale = ~PROD, data = soup,
link = "logistic")
## anova
anova(m1, update(m1, scale = ~.-PROD))
mN1 <- clm2(SURENESS ~ 1, nominal = ~PROD, data = soup,
link = "logistic")
anova(m1, mN1)
anova(m1, update(m1, scale = ~.-PROD), mN1)
## Fit model from polr example:
if(require(MASS)) {
fm1 <- clm2(Sat ~ Infl + Type + Cont, weights = Freq, data = housing)
anova(fm1, update(fm1, scale =~ Cont))
}
}
\keyword{models}
|
