Lasso pada Model Regresi Binomial Negatif

8

Apakah ada yang bisa saya lakukan LASSO dengan Regresi Binomial Negatif pada R? Saya melakukan regresi binomial negatif pada dataset saya karena data terlalu tersebar untuk memaksakan regresi poisson. Sementara itu, saya juga menghadapi beberapa masalah multikolinieritas. Saya sudah mencoba menggunakan glmnetdengan family = poisson, namun data tersebut tidak pas dengan baik (untuk kedua alpha = 0 dan alpha = 1) ... Sejujurnya aku tidak tahu apa yang harus dilakukan untuk menganalisis kekacauan besar ini data: /

Terima kasih

EDIT: di sini adalah tabel varians-kovarians dari kecocokan binomial negatif

       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 Jin-Dominique
sumber
1
Saya pikir Anda akan ingin mengambil yang ini ke Stack Overflow .
Nick Stauner
1
Sejujurnya, saya tidak yakin apakah pertanyaan ini bahkan akan on-topik pada SO; Anda mungkin ingin bertanya pada r-help listserv .
gung - Reinstate Monica
1
Pertanyaan ini tampaknya di luar topik karena menyangkut apakah analisis tertentu dapat dijalankan di R.
gung - Reinstate Monica
4
itu juga akan dipalu pada SO karena itu hanya "bagaimana saya bisa?" pertanyaan, daripada pertanyaan pemrograman tertentu ... Bisakah Anda memberikan sedikit lebih banyak konteks? Saya akan tergoda untuk melakukan quasi-Poisson fit (yaitu, pas model sebagai laso Poisson, misalnya dengan glmnetpaket, kemudian membuat penyesuaian post hoc untuk kesalahan standar dari parameter berdasarkan perkiraan sisa penyimpangan ...)
Ben Bolker
1
Sudahkah Anda mencoba model quasipoisson? Berdasarkan uraian yang sangat kabur, saya pikir masalah substantif Anda mungkin ada hubungannya dengan masalah singularitas. Tunjukkan pada kami apa yang vcov(fit)memberi, fitsebagai objek glm Anda.
AdamO

Jawaban:

5

LASSO dan metode hukuman lainnya untuk binomial negatif dan nol-inflatif negatif binom disediakan oleh mpathpaket dalam R, seperti yang telah dicatat pada halaman Cross Validated yang lebih baru . Namun, satu jawaban pada halaman itu menunjukkan beberapa kesulitan dalam menggunakan mpath. Sebuah publikasi terbaru menggambarkan aplikasi dari mpathpaket; sebuah sketsa dalam paket R mereproduksi analisis data publikasi yang.

EdM
sumber