powered by
pneumo
data(pneumo)
cumulative
# 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, data = pneumo)
Run the code above in your browser using DataLab