Controlled Accuracy Gibbs Sampling of Order Constrained Non-IID Ordered Random Variates

Abstract

Order statistics arising from m independent but not identically distributed random variables are typically constructed by arranging some X1, X2, …, Xm, with Xi having distribution function Fi(x), in increasing order denoted as X(1) ≤ X(2) ≤ … ≤ X(m). In this case, X(i) is not necessarily associated with Fi(x). Assuming one can simulate values from each distribution, one can generate such "non-iid" order statistics by simulating Xi from Fi, for i=1,2,…, m, and arranging them in order. In this paper, we consider the problem of simulating ordered values X(1), X(2), …, X(m) such that the marginal distribution of X(i) is Fi(x). This problem arises in Bayesian principal components analysis (BPCA) where the Xi are ordered eigenvalues that are a posteriori independent but not identically distributed. We propose a novel coupling-from-the-past algorithm to "perfectly" (up to computable order of accuracy) simulate such order-constrained non-iid order statistics. We demonstrate the effectiveness of our approach for several examples, including the BPCA problem.

0

Discussion (0)

Sign in to join the discussion.

Loading comments…