"R polr()"의 두 판 사이의 차이

3번째 줄: 3번째 줄:
* "Ordered Logistic Or Probit Regression"
* "Ordered Logistic Or Probit Regression"


<source lang='r' run>
<source lang='r' notebook>
options(echo = TRUE)
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)
## slightly worse fit from
</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)
## although it is not really appropriate, can fit
</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 참고