Distribusi pada simpleks dengan komponen berkorelasi

8

Saya mencari beberapa jenis distribusi melalui simpleks di mana komponen-komponennya berkorelasi secara ordinal. Yaitu, jika diambil dari distribusi kami pada simpleks, saya ingin berkorelasi positif dengan tetangganya dan , katakan . Vanilla Dirichlet jelas tidak dapat memenuhi persyaratan ini. Satu opsi yang saya kira adalah campuran dari distribusi Dirichlet; misalnya, ketika seseorang dapat mengambil atau sesuatu mirip dengan menginduksi korelasi, tetapi saya bertanya-tanya apakah ada sesuatu yang sedikit lebih alami. Pilihan lain yang saya kira adalah mengambil distribusip=(p1,...,pJ)pipi+1pi1J=4D(1,1,0,0)+D(0,1,1,0)+D(0,0,1,1){1,2,...,J} , katakanlah , beri distribusi pada ambil . Jadi saya dapat mengambil, misalnya, dan biarkan f (j | \ eta) = {J \ pilih j} \ eta ^ j (1 - \ eta) ^ {J - j} .f(j|η)ηpj=f(j|η)ηBeta(α,β)f(j|η)=(Jj)ηj(1η)Jj

Bagaimanapun, saya ingin apa pun yang saya bisa selesaikan mungkin. Campuran Dirichlet's menarik karena saya bisa mendapatkan konjugasi kondisional yang bagus untuk saya, tetapi tidak jelas bagaimana mengaturnya. Pertanyaan ini berbicara tentang distribusi normal logistik, tetapi saya tidak tahu banyak tentang itu; apakah bisa ditelusuri untuk kesimpulan Bayesian?

Tentu saja, komponen-komponen Dirichlet sudah berkorelasi negatif, dan meminta "korelasi positif" mungkin tidak sepenuhnya koheren karena jika pi besar maka secara alami, mengambil sebagian besar massa dan karenanya memaksa kemungkinan tetangganya menjadi kecil. Mungkin yang saya maksud adalah bahwa pi berkorelasi positif dengan pi+1/jipj . Semoga pertanyaan yang disampaikan sudah cukup bagi orang untuk mengetahui apa yang saya inginkan dan dapat membantu saya.

orang
sumber

Jawaban:

7

Salah satu cara untuk memiliki hidup di simpleks, tanpa batasan yang diberlakukan oleh kovarian negatif dari distribusi Dirichlet, adalah dengan mendefinisikan , untuk , di mana matriks memiliki peringkat . Menambahkan batasan , sebarang distribusi normal dimensi dapat ditugaskan ke .θ=(θ1,,θk)ϕi=j=1kcijlogθji=1,,k1(k1)×kC=(cij)k1i=1kθi=1k1ϕ=(ϕ1,,ϕk1)

Inferensi Bayesian dapat ditelusuri dalam kelas distribusi yang kaya ini yang diperkenalkan dan dipelajari oleh Aitchison dalam serangkaian makalah

Jurnal Masyarakat Statistik Kerajaan, , , 139-177 (1982),B44

Jurnal Masyarakat Statistik Kerajaan, , , 136-146 (1985);B47

dan di bukunya

The Statistical Analysis of Compositional Data . Chapman & Hall: London (1986).

Zen
sumber