coalminers <- transform(coalminers, Age = (age - 42) / 5)
fit <- vglm(cbind(nBnW, nBW, BnW, BW) ~ sm.bs(Age),
binom2.or(zero = NULL), data = coalminers)
par(mfrow = c(1, 3))
plot(fit, type = "vgam", se = TRUE, ylim = c(-3, 2), las = 1)
plot(fit, type = "vgam", se = TRUE, which.cf = 1:2,
lcol = "blue", scol = "orange", ylim = c(-3, 2))
plot(fit, type = "vgam", se = TRUE, which.cf = 1:2,
lcol = "blue", scol = "orange", overlay = TRUE)
Run the code above in your browser using DataLab