負の二項回帰モデルのなげなわ

とにかく、Rで負の二項回帰を使用してLASSOを実行できますか？データが分散しすぎてポアソン回帰を適用できないため、データセットで負の二項回帰を実行しています。その間、多重共線性の問題にも直面しています。私はすでにで使用しようとしglmnetましたfamily = poissonが、データはあまりうまく適合していません（アルファ= 0とアルファ= 1の両方）...この大きなデータの混乱を分析するために何をすべきか正直にわかりません：/

ありがとうございました

編集：これは負の二項適合の分散共分散表です

       8.392729e+18  1.239178e+06  -3.624090e+05  1.896258e+17  -3.702521e+17
       1.239178e+06  1.119052e-04   5.201989e-06 -1.877590e+05  -2.558095e+05
      -3.624090e+05  5.201989e-06   5.179343e-06 -8.021543e+04  -1.436381e+05
       1.896258e+17 -1.877590e+05  -8.021543e+04  2.193290e+17   6.413947e+16
      -3.702521e+17 -2.558095e+05  -1.436381e+05  6.413947e+16   2.142183e+17

r regression generalized-linear-model

— ジンドミニク
ソース

これをStack Overflowに持っていきたいと思います。

— Nick Stauner

正直なところ、この質問がSOについても話題になるかどうかはわかりません。r-help listservで質問することもできます。

— ガン-モニカの回復

それは、特定の分析はR.で実行できるかどうかについてですので、この質問は、オフトピックのように見える

— GUNG -復活モニカ

また、「どうすればよいのか」というだけの理由で、SOに打ち込まれます。特定のプログラミングの質問ではなく質問...もう少しコンテキストを教えてもらえますか？疑似ポアソンフィット（つまり、glmnetパッケージを使用してモデルをポアソンラッソとしてフィットさせ、推定された残差逸脱度に基づいてパラメーターの標準誤差に事後調整を行う）に誘惑されます...）

— Ben Bolker、2013

次に、準ポアソンモデルを試しましたか？非常にあいまいな説明に基づいて、私はあなたの実質的な問題はおそらく特異点の問題に関係していると思います。何vcov(fit)を与えるか、fitあなたのglmオブジェクトであるかを示してください。

— AdamO 2013

LASSOおよび負の二項式とゼロインフレートされた負の二項式に対するその他のペナルティ付きの方法はmpath、R のパッケージで提供されています。ただし、そのページの回答の1つは、の使用が困難であることを示していmpathます。最近の出版物はのアプリケーションを示しmpathたパッケージを、ビネット Rパッケージには、その出版物のデータ解析を再生します。

— EdM
ソース