# Rank-1 stereotype model of Anderson (1984)
n = nrow(pneumo)
pneumo = transform(pneumo, let=log(exposure.time), x1=runif(n), x2=runif(n))
fit = rrvglm(cbind(normal,mild,severe) ~ let + x1 + x2, multinomial, pneumo)
coef(fit, mat=TRUE)
Coef(fit)
print(Coef(fit), digits=3)
Run the code above in your browser using DataLab