(→개요) |
|||
3번째 줄: | 3번째 줄: | ||
* "Ordered Logistic Or Probit Regression" | * "Ordered Logistic Or Probit Regression" | ||
<source lang='r' | <source lang='r' notebook> | ||
library(MASS) | library(MASS) | ||
options(contrasts = c("contr.treatment", "contr.poly")) | options(contrasts = c("contr.treatment", "contr.poly")) | ||
house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = housing) | house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = housing) | ||
house.plr | house.plr | ||
</source> | |||
<source lang='r' notebook> | |||
summary(house.plr, digits = 3) | summary(house.plr, digits = 3) | ||
</source> | |||
<source lang='r' notebook> | |||
summary(update(house.plr, method = "probit", Hess = TRUE), digits = 3) | summary(update(house.plr, method = "probit", Hess = TRUE), digits = 3) | ||
</source> | |||
<source lang='r' notebook> | |||
summary(update(house.plr, method = "loglog", Hess = TRUE), digits = 3) | summary(update(house.plr, method = "loglog", Hess = TRUE), digits = 3) | ||
</source> | |||
<source lang='r' notebook> | |||
summary(update(house.plr, method = "cloglog", Hess = TRUE), digits = 3) | summary(update(house.plr, method = "cloglog", Hess = TRUE), digits = 3) | ||
</source> | |||
<source lang='r' notebook> | |||
predict(house.plr, housing, type = "p") | predict(house.plr, housing, type = "p") | ||
addterm(house.plr, ~.^2, test = "Chisq") | addterm(house.plr, ~.^2, test = "Chisq") | ||
house.plr2 <- stepAIC(house.plr, ~.^2) | house.plr2 <- stepAIC(house.plr, ~.^2) | ||
house.plr2$anova | house.plr2$anova | ||
</source> | |||
<source lang='r' notebook> | |||
anova(house.plr, house.plr2) | anova(house.plr, house.plr2) | ||
</source> | |||
<source lang='r' notebook> | |||
house.plr <- update(house.plr, Hess=TRUE) | house.plr <- update(house.plr, Hess=TRUE) | ||
pr <- profile(house.plr) | pr <- profile(house.plr) |
2021년 5월 7일 (금) 19:25 판
1 개요
- R polr()
- "Ordered Logistic Or Probit Regression"
R
Reload
Copy
library(MASS)
options(contrasts = c("contr.treatment", "contr.poly"))
house.plr <- polr(Sat ~ Infl + Type + Cont, weights = Freq, data = housing)
house.plr
Loading
Copy
summary(house.plr, digits = 3)
Loading
Copy
summary(update(house.plr, method = "probit", Hess = TRUE), digits = 3)
Loading
Copy
summary(update(house.plr, method = "loglog", Hess = TRUE), digits = 3)
Loading
Copy
summary(update(house.plr, method = "cloglog", Hess = TRUE), digits = 3)
Loading
Copy
predict(house.plr, housing, type = "p")
addterm(house.plr, ~.^2, test = "Chisq")
house.plr2 <- stepAIC(house.plr, ~.^2)
house.plr2$anova
Loading
Copy
anova(house.plr, house.plr2)
Loading
Copy
house.plr <- update(house.plr, Hess=TRUE)
pr <- profile(house.plr)
confint(pr)
plot(pr)
pairs(pr)
Loading
2 같이 보기
3 참고
편집자 Jmnote Jmnote bot
로그인하시면 댓글을 쓸 수 있습니다.