誰か良い知恵をお貸し下さい。

遺伝的アルゴリズムを利用してラッパー法による
特徴選択をしたいと考えています。
例えば、5個の説明変数の組み合わせで2値の判別分析モデルを
構築するとして、下記のように2値の5個の遺伝子で染色体を表現します。

0 1 1 0 1

上記であれば、説明変数2、3、5の組み合わせでモデルを構築します。
分類の正答率やRMSEを評価関数として遺伝的アルゴリズムは
コーディングできそうなのですが、0と1で表現される列だけ
抜き出してデータフレーム(?)を作成して、説明変数全体の
サブセットを表現する方法が分かりません。

最初はsubset関数とcomb関数を使おうと思いましたが上手くいきませんでした。

どなた様かうまい解決方法をご存じないでしょうか？