Calculation of the GAMs
stepGAM<-function(X,Y,option="gam") {