If we have a vector of probabilities $\textbf{p} = (p_1,...,p_n)$, with $p_i \in [0,1]$ and $\sum p_i = 1$, is it possible to model $\textbf{p}$ as a joint random variable with some multivariate distribution?
Each $p_i$ itself could perhaps be modeled in a variety of ways (although the lower and upper bounds complicate things a bit), but how could we possibly deal with the fact that they must sum to $1$ if we want a joint distribution?