data = read.table('clipboard')
names(data)

y=datos$V1
x1=datos$V2
x2=datos$V3
x3=datos$V4
x4=datos$V5
x5=x1*x2
x6=x1*x3
x7=x1^2
x8=x2^2
x9=x3^2
x10=x4^2

m.general=lm(y~x1+x2+x3+x3+x4+x5+x6+x7+x8+x9+x10) summary(m.general) best=regsubsets(y~x1+x2+x3+x3+x4+x5+x6+x7+x8+x9+x10,datos,nb‌​est=3) summary(best)
m.r1=lm(y~x5,datos)
summary(m.r1)
AIC(m.r1)
press(m.r1)
x=data$V6
y=data$V1
m.r1Cp=leaps(x,y,method=c("Cp","adjr2","rsq"))