# 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