
Last chance! 50% off unlimited learning
Sale ends in
data(bminz)
There are some outliers in the data set.
A variable gender
would be useful, and may be added later.
with(bminz, plot(age, BMI, col="blue"))
fit = vgam(BMI ~ s(age, df=c(2,4,2)), fam=lms.yjn, bminz, tr=TRUE)
qtplot(fit, pcol="blue", tcol="brown", lcol="brown")
Run the code above in your browser using DataLab