# Fit a GAM using vgam() and compare it with the Waitakere Ranges one
fit.h <- vgam(agaaus ~ s(altitude, df = 2), binomialff, data = hunua)
## Not run:
# plot(fit.h, se = TRUE, lcol = "orange", scol = "orange",
# llwd = 2, slwd = 2, main = "Orange is Hunua, Blue is Waitakere") ## End(Not run)
head(predict(fit.h, hunua, type = "response"))
fit.w <- vgam(agaaus ~ s(altitude, df = 2), binomialff, data = waitakere)
## Not run:
# plot(fit.w, se = TRUE, lcol = "blue", scol = "blue", add = TRUE) ## End(Not run)
head(predict(fit.w, hunua, type = "response")) # Same as above?
Run the code above in your browser using DataLab