powered by
pneumo
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