2012-04-16 32 views
9

Tôi tìm thấy một gói 'bivpois' cho R đánh giá một mô hình cho hai quá trình Poisson liên quan (ví dụ, số lượng mục tiêu của nhà và đội khách trong một trò chơi bóng đá). Tuy nhiên, gói này dường như không còn sử dụng được trong các phiên bản mới hơn của R.Biến đổi Poisson Bivariate trong R?

Có cách nào hợp lý để sửa đổi hàm glm() để thực hiện quy trình tương tự hay chạy gói cũ hơn này trên phiên bản R mới của tôi không? Tôi đã tìm thấy rất ít tài liệu về các loại quy trình này và đã tìm thấy rất ít về mặt triển khai dễ dàng trong các gói thống kê khác như STATA.

Mọi đề xuất sẽ được đánh giá cao.

+1

Nếu bạn mở để sử dụng OpenBUGS, có vẻ như tác giả của gói có thể bao gồm mã ở đây (và có thể thảo luận cách tiếp cận trên trang 291-293 của cuốn sách của mình) http://stat-athens.aueb.gr/~ jbn/winbugs_book/Ngay cả khi bạn không muốn sử dụng Bayesian bạn vẫn có thể kiểm tra cuốn sách bởi vì tôi nghĩ đôi khi anh ấy trình bày mã R thường xuyên, mặc dù tôi không chắc chắn 100% về điều đó. Ngoài ra, hãy xem xét gửi email cho anh ấy. –

+0

Google đã bật tính năng này http://www.ats.ucla.edu/stat/r/dae/poissonreg.htm – Seth

+0

@Seth, tôi không nghĩ rằng liên kết của bạn trả lời câu hỏi (OP muốn * bivariate * Hồi quy Poisson, không hồi quy Poisson trơn tru) –

Trả lời

6

Trong khi CRAN không lưu trữ nhị phân hiện tại của bivpois, bạn có thể tạo gói từ mã nguồn đã lưu trữ (xem http://cran.r-project.org/doc/manuals/R-exts.html#Checking-and-building-packages). Xây dựng bivpois 0,50-3,1 từ nguồn (có sẵn tại http://cran.r-project.org/src/contrib/Archive/bivpois/) hoạt động với tôi trên R 2.15.0 Windows x64. Nén nhị phân Windows được nén mà tôi đã tạo sẵn có tại đây: http://commondatastorage.googleapis.com/jthetzel-public/bivpois_0.50-3.1.zip.

+0

jthetzel, Cảm ơn thông tin. Tôi đã đi trước và tải về một phiên bản cũ của R và gói nén của tôi không hoạt động trên đó. Vì vậy, tôi đoán một câu hỏi followup là liệu có bất kỳ chức năng như 'predprob' (trong gói 'pscl'), hoặc nếu tôi phải tự làm dự đoán từ mô hình (họ không phải là thẳng về phía trước như tôi đã hy vọng). Cảm ơn một lần nữa. – Millsy11

+0

Bạn cũng có thể xem gói VGLM trên CRAN. –

+0

Xin lỗi, đó là gói VGAM! –

Các vấn đề liên quan