Bay phân cấp cho R hoặc Python

Question

Các mô hình Bay phân cấp thường được sử dụng trong Marketing, Khoa học Chính trị và Toán Kinh tế. Tuy nhiên, gói duy nhất tôi biết là bayesm, thực sự là bạn đồng hành với một cuốn sách (Thống kê Bayesian và Tiếp thị, bởi Rossi, et al.) Tôi có thiếu gì đó không? Có một gói phần mềm cho R hay Python thực hiện công việc ngoài kia, và/hoặc một ví dụ đã làm việc trong ngôn ngữ liên quan không?

Answer 1

Dưới đây là bốn cuốn sách về mô hình phân cấp và phân tích Bayes được viết bằng mã R trong suốt sách.

Hierarchical Modeling và phân tích cho không gian dữ liệu (chuyên khảo về thống kê và Xác suất Ứng dụng) (Bìa cứng) http://www.amazon.com/gp/product/158488410X

phân tích dữ liệu Sử dụng Regression và đa cấp/Mô hình phân cấp (Paperback) http://www.amazon.com/Analysis-Regression-Multilevel-Hierarchical-Models/dp/052168689X/ref=pd_sim_b_1

Bayesian Tính với R (Sử dụng R) (Bìa mềm) http://www.amazon.com/Bayesian-Computation-R-Use/dp/0387922970/ref=pd_bxgy_b_img_c

Lập mô hình phân cấp cho khoa học môi trường: Phương pháp thống kê và áp dụng ions (Oxford Biology) (Bìa mềm) (Tôi giả định rằng mã này có mã R vì cả hai tác giả đều sử dụng R rộng rãi)

Tôi biết một số sách python có khả năng phân tích đa biến (ví dụ như Trí tuệ tập thể) nhưng tôi chưa nhìn thấy bất kỳ điều gì thực sự nghiên cứu kỹ lưỡng về mô hình phân tầng hoặc phân cấp.

Answer 2

Có gói OpenBUGS và R helper. Kiểm tra trang web Gelman cho cuốn sách của ông, trong đó có hầu hết các liên kết có liên quan:

• http://www.stat.columbia.edu/~gelman/software/
• Example of computation in R and Bugs

Về phía Python, tôi chỉ biết PyMC:

• http://code.google.com/p/pymc/
• An example statistical model

EDIT: Đã thêm liên kết vào phụ lục thích hợp từ sách của Gelman, có sẵn trực tuyến, ví dụ sử dụng R và BUGS.

Answer 3

Có một vài mô hình phân cấp trong số MCMCpack cho R, kiến ​​thức của tôi là mẫu lấy mẫu nhanh nhất cho nhiều loại mô hình phổ biến. (Tôi đã viết mô hình [mục phản hồi phân cấp] [2] trong đó.)

[RJAGS] [3] tên của nó nghe như thế nào. Mã hóa mô hình .bug, cung cấp dữ liệu trong R và gọi Jags từ R.

Answer 4

Gói lme4, ước tính mô hình phân cấp sử dụng phương pháp thường xuyên, có hàm gọi là mcmcsamp cho phép bạn lấy mẫu từ hậu sau phân phối mô hình bằng MCMC. Điều này hiện chỉ hoạt động cho các mô hình tuyến tính, khá tiếc.

Answer 5

trong python, hãy thử PyMC.Có một ví dụ về mô hình đa cấp với nó ở đây: http://groups.google.com/group/pymc/browse_thread/thread/c6ce37a80edf7f85/1bfd9138c8db891d

Answer 6

Tôi áp dụng mô hình Bayes thứ bậc trong R kết hợp với JAGS (Linux) hoặc đôi khi WinBUGS (Windows hoặc Wine). Hãy xem cuốn sách của Andrew Gelman, như đã nói ở trên.