如何在矩阵内搜索少量特征

Question

我正在尝试做一些简单的组合，我有一个X矩阵由2000个变量组成，YI希望每次从X中选择4个变量并进行回归并仅保留那些R值高于R的集合0.9。我找到了另一个基于3个矩阵的解决方案，想要的是一个矩阵，我试图修改代码，但有些我没有把它弄好

这就是我试过的

tic
X1 = rand(100,2000);
Y  = rand(100,1);
in = ones(100,1);
R = zeros(size(:,1),1);
for c = 1:size(:,1)
    X    = [in X1(:,i)];
    yhat = X*(X\Y);
    ybar = mean(Y);
    ssr  = norm(yhat - ybar)^2;  
    sst  = norm(Y - ybar)^2; 
    R(c) = ssr/sst;
end
toc