pneumo <- transform(pneumo, let = log(exposure.time))
fit <- vglm(cbind(normal, mild, severe) ~ let, acat, data = pneumo)
coef(fit) # 8-vector
Influence(fit) # 8 x 4
all(abs(colSums(Influence(fit))) < 1e-6) # TRUE
Run the code above in your browser using DataLab