Splus にコードがありますが、それを R に変換する必要がありますが、これは大したことではありません。ただし、私は両方のソフトウェアに非常に慣れていません。これは私が苦労しているコードです:
a <- regsubsets(x, y, wt = wt, method = "adjr2", nbest=nbest,
names = dimnames(x)[[2]], int=F)
x は独立変数の行列で、最初の length(keep1) 列は常に BMA に保持される変数に対応します (ベイジアン モデルの平均化 -- これは重要ではありません。基本的に x は行列です)。y も行列です。
当初、これは「リープ」関数を使用することでしたが、行列が大きすぎたため、R から regsubsets を使用するように依頼されました。現在、「adjr2」は「あいまいまたは認識されないメソッド名」であると言われています。これは、leaps関数で認識されているメソッドであると確信していますか? leaps 関数ではなく regsubset 関数を使用しているため、使用する別のコマンドはありますか? 助けてください！