
The pneumo
data frame has 8 rows and 4 columns.
Exposure time is explanatory, and there are 3 ordinal response variables.
data(pneumo)
This data frame contains the following columns:
a numeric vector, in years
a numeric vector, counts
a numeric vector, counts
a numeric vector, counts
These were collected from coalface workers. In the original data set, the two most severe categories were combined.
McCullagh, P. and Nelder, J. A. (1989) Generalized Linear Models, 2nd ed. London: Chapman & Hall.
# NOT RUN {
# 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