# NOT RUN {
fit <- vgam(cbind(nolnor, nolr, lnor, lr) ~ s(age, c(4, 3)),
binom2.or(exchan = TRUE, zero = NULL), data = chest.nz)
coef(fit, matrix = TRUE)
# }
# NOT RUN {
plot(fit, which.cf = 2, se = TRUE)
# }
Run the code above in your browser using DataLab