X: Covariables matrix # Y: Response vector

n=dim(X)[1] nvar=dim(X)[2]

stepGAM<-function(X,Y,option="gam") {
