pneumo
data frame has 8 rows and 4 columns.
Exposure time is explanatory, and there are 3 ordinal response variables.data(pneumo)
# Fit the proportional odds model, p.179, in McCullagh and Nelder (1989)
pneumo = transform(pneumo, let = log(exposure.time))
vglm(cbind(normal, mild, severe) ~ let, propodds, pneumo)
Run the code above in your browser using DataLab