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
# This is new:
fit@family@infos()$muxypw # Proportion or count?
fit@family@infos()$Mux4vglmnet # vglmnet needs counts?
fit@family@infos()$roundmux # Product is integer-valued?
depvar(fit, muxypw = TRUE, roundmux = TRUE) # Counts
Run the code above in your browser using DataLab