# Rank-1 stereotype model of Anderson (1984)
set.seed(111)
pneumo = transform(pneumo, let=log(exposure.time),
x1=runif(nrow(pneumo))) # x1 is some unrelated covariate
fit = rrvglm(cbind(normal, mild, severe) ~ let + x1,
multinomial, pneumo, Rank=1)
Coef(fit)
Run the code above in your browser using DataLab