# NOT RUN {
fit <- vgam(agaaus ~ s(altitude, df = 2), binomialff, waitakere)
head(predict(fit, waitakere, type = "response"))
# }
# NOT RUN {
plot(fit, se = TRUE, lcol = "orange", scol = "blue")
# }
Run the code above in your browser using DataLab