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