R code used to produce the graph shown below.

library(mgcv) 
library(isdals) 
data(fev) 
model <- gam(FEV ~ Gender + Smoke + s(Age) + s(Ht), data=fev) 
summary(model)