v<-list()

 for(i in seq(from=0,to=25, by=2)){
 +                    result<-NULL
 +                            XyMatrixtoTest<-na.omit(cbind(X_Data[2],X_Data[4],X_Data[6],X_Data[12],X_Data[17],diff( lagpad(Y_Data[,5],k=i)),2))
 +                               result<-bestglm(Xy=XyMatrixtoTest, family = gaussian, IC = "BIC", t = "default",  CVArgs = "default", qLevel = 0.99, TopModels = T,  method = "exhaustive", intercept = FALSE, weights = NULL, nvmax = "default", RequireFullEnumerationQ = FALSE)
 +                             XyMatrixtoTest<-matrix(ncol=ncol(XyMatrixtoTest), rep(NA, prod(dim(XyMatrixtoTest))))
 + v[[i+1]]<-result
+ }