data(bminz)
fit = vgam(BMI ~ s(age, df=c(4,2)), fam=lms.bcn(zero=1), data=bminz)
fit@post$cdf[1:5]
cdf(fit)[1:5] # Same
fit@y[1:5]
fitted(fit)[1:5,]
cdf(fit, data.frame(age=c(31.5,39), BMI=c(28.4,24)))
Run the code above in your browser using DataLab