library(MASS)
library(leaps)
library(corrplot)
library(glmnet)

reg.data <- read.table("tp3_a18_reg_app.txt")

reg.fit <-regsubsets(y ~ .,data=reg.data,method="exhaustive",nvmax=50, really.big = T)
subsets <- summary(reg.fit)

reg.mask <- subsets$which[24,2:51]

reg.train.X <- reg.data[,reg.mask]
reg.train.y <- reg.data[,y]