options(contrasts = c("contr.treatment", "contr.poly"))
data(soup)
m1 <- clm(SURENESS ~ PROD, scale = ~PROD, data = soup,
link = "logistic")
m2 <- update(m1, link = "probit")
m3 <- update(m1, link = "cloglog")
m4 <- update(m1, link = "loglog")
anova(m1, update(m1, scale = ~.-PROD))
mT1 <- update(m1, threshold = "symmetric")
## Fit model from polr example:
data(housing, package = "MASS")
fm1 <- clm(Sat ~ Infl + Type + Cont, weights = Freq, data = housing)
update(fm1, location = ~ . - Type, scale =~ Cont)
Run the code above in your browser using DataLab