Verteilungen auf dem Simplex mit korrelierten Komponenten

8

Ich suche nach einer Art Verteilung über den Simplex, in der Komponenten auf ordinale Weise korreliert werden. Das heißt, wenn p=(p1,...,pJ) aus unserer Verteilung auf dem Simplex gezogen wird, ich möchte pi positiv mit seinen Nachbarn korreliert werden und , sag. Ein Vanille-Dirichlet kann diese Anforderung eindeutig nicht erfüllen. Eine Option, nehme ich an, ist eine Mischung aus Dirichlet-Verteilungen; Wenn beispielsweise , könnte man D ( 1 , 1 , 0) nehmenpi+1pi1J=4D(1,1,0,0)+D(0,1,1,0)+D(0,0,1,1) oder ähnliches, um eine Korrelation zu induzieren, aber ich frage mich, ob es etwas Natürlicheres gibt. Eine andere Möglichkeit, nehme ich an, ist eine Verteilung auf{1,2,...,J} , sagen Sief(j|η) , setzen Sie eine Verteilung aufη takepj=f(j|η) . So könnte ich zum BeispielηBeta(α,β) und .f(j|η)=(Jj)ηj(1η)Jj

Auf jeden Fall möchte ich, dass alles, was ich am Ende habe, so gut wie möglich nachvollziehbar ist. Die Mischung aus Dirichlets ist ansprechend, weil ich eine nette bedingte Konjugation für mich in Gang bringen könnte, aber es ist nicht klar, wie man die Dinge einrichtet. Diese Frage spricht über die logistische Normalverteilung, aber ich weiß nicht viel darüber; ist es für die Bayes'sche Folgerung nachvollziehbar?

Natürlich sind die Komponenten eines Dirichlets bereits negativ korreliert, und die Frage nach einer "positiven Korrelation" ist wahrscheinlich nicht vollständig kohärent, da wenn es groß ist, von Natur aus den größten Teil der Masse einnimmt und somit die Wahrscheinlichkeit von erzwingt seine Nachbarn sollen klein sein. Vielleicht meine ich, dass p_i positiv mit p_ {i + 1} / \ sum_ {j \ ne i} p_j korreliertpipipi+1/jipj . Hoffentlich reicht die Frage aus, damit die Leute wissen, was ich will und mir helfen können.

Kerl
quelle

Antworten:

7

Eine Möglichkeit, ein zufälliges ohne die durch die negativen Kovarianzen der Dirichlet-Verteilung auferlegten Einschränkungen auf dem Simplex leben zu lassen, besteht darin, , für , wobei die Matrix den Rang . Durch Hinzufügen der Einschränkung kann jede dimensionale Normalverteilung zugewiesen werden.θ=(θ1,,θk)ϕi=j=1kcijlogθji=1,,k1(k1)×kC=(cij)k1i=1kθi=1k1ϕ=(ϕ1,,ϕk1) .

Die Bayes'sche Folgerung lässt sich innerhalb dieser reichen Klasse von Verteilungen nachvollziehen, die Aitchison in einer Reihe von Arbeiten vorgestellt und untersucht hat

Zeitschrift der Royal Statistical Society, ,B44 , 139-177 (1982),

Zeitschrift der Royal Statistical Society, ,B47 , 136-146 (1985);

und in seinem Buch

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

Zen
quelle