Generalized Rank Dirichlet distributions

Abstract

We study a new parametric family of distributions on the ordered simplex ∇d−1={y∈Rd:y1≥⋯≥yd≥0,∑k=1dyk=1}, which we call Generalized Rank Dirichlet (GRD) distributions. Their density is proportional to ∏k=1dykak−1 for a parameter a=(a1,…,ad)∈Rd satisfying ak+ak+1+⋯+ad>0 for k=2,…,d. The density is similar to the Dirichlet distribution, but is defined on ∇d−1, leading to different properties. In particular, certain components ak can be negative. Random variables Y=(Y1,…,Yd) with GRD distributions have previously been used to model capital distribution in financial markets and more generally can be used to model ranked order statistics of weight vectors. We obtain for any dimension dexplicit expressions for moments of order M∈N for the Yk’s and moments of all orders for the log gaps Zk=logYk−1−logYk when a1+⋯+ad=−M. Additionally, we propose an algorithm to exactly simulate random variates in this case. In the general case a1+⋯+ad∈R we obtain series representations for these quantities and provide an approximate simulation algorithm.