## Not run: with(bmi.nz, plot(age, BMI, col = "blue"))
# fit <- vgam(BMI ~ s(age, df = c(2, 4, 2)), lms.yjn, data = bmi.nz, trace = TRUE)
# qtplot(fit, pcol = "blue", tcol = "brown", lcol = "brown") ## End(Not run)
Run the code above in your browser using DataLab