問題タブ [glmnet]

glmnet で次のコードを使用しています。

ただし、リッジ回帰で行われるように、せいぜい Lambda で係数を出力したいと思います。次の適合構造が表示されます。

しかし、最適な Lambda と対応する係数を取得できません。ご協力いただきありがとうございます。

この投稿 ( R の glmnet プロットで曲線にラベルを追加) の関数を変更して、次のようにプロットに凡例を追加しました。

ただし、プロットと凡例のテキスト ラベルが一致しません。変数「am」は明らかに間違った色になっています。エラーはどこにありますか? ご協力いただきありがとうございます。

glmnet を使用して、ペナルティ付きロジスティック回帰法を学習しようとしています。mtcars のサンプル データの車にオートマチック トランスミッションまたはマニュアルのどちらがあるかを予測しようとしています。私のコードはかなり単純だと思いますが、エラーが発生しているようです:

この最初のブロックは、mtcar を 80% のトレーニング セットと 20% のテスト セットに単純に分割します。

x データは応答のない行列形式であることがわかっているので、2 つのトレーニング セットを非応答行列 (train_x) と応答ベクトル (train_y) に分けます。

しかし、モデルをトレーニングしようとすると、

エラーが発生します：

何か不足していますか？

glmnet パッケージを使用して、並列化されたロジスティック リッジ回帰を実行したいと考えています。私のデータは大きな疎行列です (1,000 万の観測値と約 60,000 列)。

データのサブセット (観測と列サブセットの両方) に対して小さな試行を行ったところ、うまくいきました。次のコードは、私がやっていることと同等です:

ただし、マトリックス全体を入力すると、プロセスがクラッシュし、次のエラー メッセージが表示されます。

エラーの原因がわかりません。また、エラー メッセージが何を指しているのかわかりません。

8 コアの RStudio Linux サーバーで r を使用しています。

sessionInfo():

更新I：

エラー (機密性の問題) を生成するデータを共有できないため、示されているエラーではなく、メモリ オーバーフローを生成しようとした再現を試みたので、質問を再定式化します。

私が得たエラーメッセージはメモリに関連していますか、それとも何か他のものに関連していますか?

データセットのサイズを考えると、メモリ関連のエラーはオプションです。ただし、エラー メッセージは、ラムダ値内に複数の最小値があることに関連する内部の問題を示しています。メモリの問題ではない場合、どうすればよいですか?回避策はありますか?

glmnet (なげなわ用) とキャレット (k 分割交差検証用) を使用して、ロジスティック回帰モデルをデータに適合させようとしています。2 つの異なる構文を試しましたが、どちらもエラーをスローします。

データフレームをマトリックスに変換する必要がありますか?

応答変数は因数である必要がありますか、それとも 0/1 の整数である必要がありますか?

df_without データ フレームを含む .Rdata ファイルはこちらです。

sessionInfo()

私はデータセットに取り組んでおり、線形 LASSO/Ridge 回帰にglmnetを使用しています。

簡単にするために、私が使用しているモデルは次のとおりであると仮定しましょう。

クライアント向けのプレゼンテーションを準備していて、変数の T 統計とR 二乗値を示す必要があります。さらに、モデルの適合値に対して残差をプロットする必要もあります。

これを行う関数をゼロから作成する前に、これがライブラリで既にカバーされているかどうかを尋ねたかったのです。glmnet ビネットを確認しましたが、何も見つかりませんでした。

ご協力いただきありがとうございます！

glmnetMatLab のパッケージを使用しようとしています。glmnet は正常に動作しますが、cvglmnet(パッケージの別の機能である) 使用しようとすると、MatLab は「壊滅的な」エラーを返し、閉じる必要があります。呼び出している MEX ファイルと関係があります。以下にエラーメッセージ全体をコピーしました。どんな助けでも大歓迎です。