Lasso về mô hình hồi quy nhị thức âm

Có cách nào để tôi có thể thực hiện LASSO với Hồi quy nhị thức âm trên R không? Tôi đang thực hiện hồi quy nhị thức âm trên tập dữ liệu của mình vì dữ liệu quá phân tán để áp đặt hồi quy poisson. Trong khi đó, tôi cũng đang phải đối mặt với một số vấn đề về đa cộng đồng. Tôi đã thử sử dụng glmnetvới family = poisson, nhưng dữ liệu không phù hợp lắm (cho cả alpha = 0 và alpha = 1) ... Thật lòng tôi không biết phải làm gì để phân tích mớ dữ liệu lớn này: /

cảm ơn bạn

EDIT: đây là bảng phương sai hiệp phương sai của nhị thức âm

       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
nguồn

Tôi nghĩ bạn sẽ muốn đưa cái này vào Stack Overflow .

— Nick Stauner

Thành thật mà nói, tôi không chắc câu hỏi này thậm chí sẽ có chủ đề về SO hay không; bạn có thể muốn hỏi về listerv r-help .

— gung - Phục hồi Monica

Câu hỏi này dường như không có chủ đề vì đó là về việc liệu một phân tích cụ thể có thể được chạy trong R.

— gung - Tái lập Monica

nó cũng sẽ bị tấn công bởi vì nó chỉ là "làm thế nào tôi có thể?" câu hỏi, thay vì một câu hỏi lập trình cụ thể ... Bạn có thể đưa ra một chút bối cảnh không? Tôi sẽ bị cám dỗ để làm một sự phù hợp quasi-Poisson (tức là, phù hợp với những mô hình như một Poisson Lasso, ví dụ như với các glmnetgói phần mềm, sau đó đưa ra một bài hoc điều chỉnh các sai số chuẩn của các thông số dựa trên lệch còn lại ước tính ...)

— Ben Bolker

Bạn đã thử một mô hình quasipoisson sau đó? Dựa trên mô tả rất mơ hồ, tôi nghĩ rằng vấn đề thực sự của bạn có thể liên quan đến một vấn đề đơn lẻ. Cho chúng tôi thấy những gì vcov(fit)cho, fitlà đối tượng glm của bạn.

— AdamO

LASSO và các phương pháp bị phạt khác đối với nhị thức âm tính và nhị phân âm được thổi phồng bằng mpathgói được cung cấp trong R, như đã được ghi chú trên trang Xác thực chéo gần đây . Tuy nhiên, một câu trả lời trên trang đó cho thấy một số khó khăn trong việc sử dụng mpath. Một ấn phẩm gần đây minh họa một ứng dụng của mpathgói; một họa tiết trong gói R tái tạo phân tích dữ liệu của ấn phẩm đó.

— EdM
nguồn