pneumo <- transform(pneumo, let = log(exposure.time))
(fit <- vglm(cbind(normal, mild, severe) ~ let, propodds, pneumo))
fit@y # Sample proportions (not recommended)
depvar(fit) # Better than using fit@y
weights(fit, type = "prior") # Number of observations
Run the code above in your browser using DataLab