R polr()

1 개요

R polr()

• "Ordered Logistic Or Probit Regression"

library(MASS)
options(contrasts = c("contr.treatment", "contr.poly"))
house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = housing)
house.plr

summary(house.plr, digits = 3)

summary(update(house.plr, method = "probit", Hess = TRUE), digits = 3)

summary(update(house.plr, method = "loglog", Hess = TRUE), digits = 3)

summary(update(house.plr, method = "cloglog", Hess = TRUE), digits = 3)

predict(house.plr, housing, type = "p")

addterm(house.plr, ~.^2, test = "Chisq")

house.plr2 <- stepAIC(house.plr, ~.^2)

house.plr2$anova

anova(house.plr, house.plr2)

house.plr <- update(house.plr, Hess=TRUE)
pr <- profile(house.plr)
confint(pr)

plot(pr)

pairs(pr)

2 같이 보기

• R glm()
• R optim()
• R multinom()

3 참고

• https://www.rdocumentation.org/packages/MASS/versions/7.3-51.5/topics/polr