Phân phối siêu nhân cho các tham số (ma trận tỷ lệ và mức độ tự do) của một điều ước trước một ma trận hiệp phương sai nghịch đảo


9

Tôi đang ước tính một số ma trận hiệp phương sai nghịch đảo của một tập hợp các phép đo trên các nhóm dân cư khác nhau bằng cách sử dụng một wishart trước trong jags / rjags / R.

Thay vì chỉ định một ma trận tỷ lệ và mức độ tự do trên ma trận hiệp phương sai trước (phân phối wishart), tôi muốn sử dụng một siêu nhân trên ma trận tỷ lệ và mức độ tự do, để chúng có thể được ước tính từ sự thay đổi giữa các quần thể.

Tôi chưa tìm thấy nhiều tài liệu về các siêu nhân cho ma trận quy mô và mức độ tự do. Hầu hết các tài liệu dường như dừng việc phân cấp theo sự lựa chọn của hiệp phương sai / hiệp phương sai và / hoặc tập trung vào việc ước tính một ma trận hiệp phương sai duy nhất thay vì một số ma trận hiệp phương sai trên các quần thể khác nhau.

Bất kỳ đề xuất nào về cách thực hiện điều này - các phân phối siêu nhân được đề xuất sử dụng cho ma trận tỷ lệ và mức độ tự do của phân phối wishart là gì? Có một số tài liệu về điều này mà tôi đang thiếu?

Câu trả lời:


5

DPpackage của R cho phép một hệ thống phân cấp đi xa như bạn đang đề xuất trên ma trận tỷ lệ trong DPd mật độ chức năng. Bạn có thể xem những gì họ làm trong hướng dẫn của họ hoặc trong họa tiết liên quan để có được một số ý tưởng. Đặt là ma trận hiệp phương sai. Nó đặt và trong đó nghịch đảo-Wishart với mức độ tự do và có nghĩa là trong đó là thứ nguyên của dữ liệu. Điều này thoạt nhìn hơi ngược với tôi lúc đầu, nhưng nếu bạn chơi với mật độ bạn có thể thấy nó là liên hợp. Mật độ Wishart không có vẻ hứa hẹn cho bất cứ điều gì phân tíchΣΣIW(ν1,Ψ1)Ψ1IW(ν2,Ψ2)IW(ν,Ψ)νΨ1νp1pν . Bạn luôn có thể đặt bất cứ thứ gì lên và sử dụng một bước của Metropolis-Hastings.ν

EDIT : Tôi chỉ nhận thấy bạn đang sử dụng jags. Có một cơ hội tốt tôi nghĩ rằng nó sẽ khó chịu nếu bạn cố gắng đặt bất kỳ ưu tiên nào vào , mặc dù nghịch đảo-Wishart là liên hợp. Việc triển khai BUGS có thể không ổn định về những gì họ cho phép cho các bản phân phối đa biến của họ, vì vậy nó có thể không biết cách thực hiện cập nhật liên hợp. Tôi không biết chắc chắn mặc dù.Ψ1


Thật không may, tôi nghĩ rằng bạn đúng về jags. Tôi gặp lỗi "Không thể tìm thấy bộ lấy mẫu thích hợp" khi tôi đặt trước cho Phi1 đó là một nghịch đảo () được áp dụng cho một người lùn () rv
user4733
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.