coalminers <- transform(coalminers, Age = (age - 42) / 5)
fit <- vglm(cbind(nBnW,nBW,BnW,BW) ~ Age, binom2.or(zero = NULL),
data = coalminers)
is.zero(fit)
is.zero(coef(fit, matrix = TRUE))
Run the code above in your browser using DataLab