> library(foreign)
> ds <- read.dta("c:/arthritis.dta")
> attach(ds)
> library(MASS)
> model1 <- polr(factor(y4) ~ age + trt)
> summary(model1)
Re-fitting to get Hessian
Call:
polr(formula = factor(y4) ~ age + trt)
Coefficients:
Value Std. Error t value
age 0.02048 0.009833 2.083
trt -0.60791 0.214223 -2.838
Intercepts:
Value Std. Error t value
1|2 -1.2119 0.5317 -2.2793
2|3 0.5251 0.5249 1.0003
3|4 2.1494 0.5381 3.9944
4|5 4.1365 0.6122 6.7571
Residual Deviance: 814.1095
AIC: 826.1095
(10 observations deleted due to missingness)