# NOT RUN {
Mental <- read.table("http://www.stat.ufl.edu/~aa/glm/data/Mental.dat",
header = TRUE) # Make take a while to load in
Mental$impair <- ordered(Mental$impair)
pfit3 <- vglm(impair ~ ses + life, data = Mental,
cumulative(link = "probit", reverse = FALSE, parallel = TRUE))
coef(pfit3, matrix = TRUE)
ordsup(pfit3) # The 'ses' variable is binary
# Fit a crude LM
fit7 <- vglm(as.numeric(impair) ~ ses + life, uninormal, data = Mental)
coef(fit7, matrix = TRUE) # 'sd' is estimated by MLE
ordsup(fit7)
ordsup(fit7, all.vars = TRUE) # Some output may not be meaningful
ordsup(fit7, confint = TRUE, method = "profile")
# }
Run the code above in your browser using DataLab