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