有一个100*6的矩阵X，X=（x1 x2 x3 x4 x5 x6),每个xi有100个观测，然后要对X的列的所有情况和y做回归，就是第一列有可能选有可能不选，第二列有可能选有可能不选，共需要做2^6次方个回归，请问各位大神这个R代码怎么写哇
关于X和Y：
miu<-c(0,0,0,0,0,0) #Xi~N(0,I6)
sigma<-diag(c(1,1,1,1,1,1))#I6
library(MASS)
x<-mvrnorm(100,miu,sigma)#生成100个服从miu和sigma的正态分布随机数
u<-rnorm(100,0,4)#随机误差项εi~N(0,4)
beta<-matrix(c(1,0.8,0.6,0,0,0),6,1)#生成β列向量
y<-x%*%beta+u#响应变量y=xβ+ε
以及请问各位：为什么我用这个com<-combn(c(x[,1],x[,2],x[,3],x[,4],x[,5],x[,6]),p2)在6列里选两列会出来2万种以及如何实现矩阵中在6列里选出两列