[プログラム・コード公開] コピペだけで実用的かつ実践的なLeast Absolute Shrinkage and Selection Operator (LASSO) (R言語)

【更新・値上げ中】好評につき再度 値上げしました。

仕事や研究において、変数選択が同時に可能な回帰分析を行うためにLeast Absolute Shrinkage and Selection Operator (LASSO) をする方もいらっしゃると思います。LASSOの実用的かつ実践的な方法はこちらに書きました。

しかし、LASSOのやり方はわかっても、実際にLASSOができるようになるわけではありません。ネットや本でLASSOのプログラミングを説明しているものはありますが、データの読み込み方とか結果の出し方とか、他にも調べてやらなくちゃいけないこと、多いんですよね・・・。手間と時間がかかります。

そこですぐにLASSOを実行可能なプログラムを作りました。下に示す形式のデータ(data.csv, data_prediction1.csv, data_prediction2.csv)さえ準備すれば、R言語でLASSOができます。

●必要なもの

・モデル構築用データのcsvファイル [data.csv]
下図のように、一番上が変数の名前、一番左がサンプルの名前です。一番左の変数が目的変数Yでありその右が説明変数Xです。サンプルの名前と変数の名前はすべて異なるものにしてください。

・予測用データ1のcsvファイル [data_prediction1.csv]
予測用データ1は目的変数Yの値が分かっている予測用データです。data.csvと説明変数の種類および個数を揃える必要があります。準備できない場合はdata.csvと同じもので名前をdata_prediction1.csvとしてください。下図のように、一番上が変数の名前、一番左がサンプルの名前です。一番左の変数が目的変数Yでありその右が説明変数Xです。サンプルの名前はすべて異なるものにしてください。

・予測用データ2のcsvファイル [data_prediction2.csv]
予測用データ2は目的変数Yの値が分かっていない予測用データです。data.csvと説明変数の種類および個数を揃える必要があります。準備できない場合はdata.csvの説明変数だけ取り出したもので名前をdata_prediction2.csvとしてください。下図のように、一番上が変数の名前、一番左がサンプルの名前です。サンプルの名前はすべて異なるものにしてください。

実行結果を下に示します。data_prediction2.csvの目的変数の予測値がPredictedY2.csvというファイルに保存され、標準回帰係数がRegressionCoefficients.csvというファイルに保存されます。標準回帰係数の値が0以外の変数が選択されたということになります。

R言語で必要なものは以下の通りです。
RStudio [こちらではVersion 0.99.879です]

もちろんLASSOからスタートしてさらにプログラミングを進めたいと考えている方にもぜひ利用していただければと思います。

R言語のプログラムは有料コンテンツとします。ただペットボトル一本分で、こちらに記載したLASSOをそのまますぐに実行できます。

こちらからプログラムのzipファイル自体はダウンロードできます。購入していただくと解凍のためのパスワードがありますのでそちらをご利用ください。