Phân tích Bayes không đối xứng trong R

Tôi đang tìm kiếm một hướng dẫn tốt về phân cụm dữ liệu trong Rviệc sử dụng quy trình dirichlet phân cấp (HDP) (một trong những phương pháp Bayesian không tham số gần đây và phổ biến).

Có DPpackage(IMHO, toàn diện nhất trong tất cả những cái có sẵn) trong Rphân tích Bayes không đối xứng. Nhưng tôi không thể hiểu các ví dụ được cung cấp trongR News hoặc trong tài liệu tham khảo gói đủ tốt để mã HDP.

Bất kỳ trợ giúp hoặc con trỏ được đánh giá cao.

Một triển khai C ++ của HDP cho mô hình chủ đề có sẵn ở đây (vui lòng xem phía dưới để biết mã C ++)

Dưới đây là một số nguồn tài nguyên trực tuyến tôi thấy thú vị mà không đi sâu vào chi tiết (và tôi không phải là chuyên gia của chủ đề này):

• Các quy trình Dirichlet phân cấp , bởi Teh et al. (2005)
• Quy trình Dirichlet Một hướng dẫn nhẹ nhàng , của El-Arini (2008)
• Bayesian Nonparametrics , bởi Rosasco (2010)
• Phương pháp Bayes không tham số , bởi Ghahramani (2005)

Các tài liệu tham khảo dứt khoát dường như là

N. Hjort, C. Holmes, P. Müller và S. Walker, biên tập viên. Bayesian Nonparametrics . Số 28 trong sê-ri Cambridge về Toán học thống kê và xác suất. Nhà xuất bản Đại học Cambridge, 2010.

Về R, dường như có một số gói khác đáng để khám phá nếu DPpackage không phù hợp với nhu cầu của bạn, ví dụ: dpmixsim , BHC hoặc mbsc được tìm thấy trên Rseek.org .

Hai liên kết này cung cấp một số ví dụ mã R (và C) khi thực hiện hỗn hợp DP thông thường:

• http://ice.uchicago.edu/2008_presentations/Rossi/d mật_estimation_with_DP_priors.ppt
• http://www.duke.edu/~neelo003/r/DP02.r

Tìm thấy một tài liệu tham khảo khác. Chương 15 có mã winbugs DP:

• http://www.ics.uci.edu/~wjohnson/BIDA/BIDABook.html