# Perfect Gibbs Sampling of Order Constrained Non-IID Ordered Random Variates with Application to Bayesian Principal Components Analysis

Order statistics arising from m independent but not identically distributed random variables are typically constructed by arranging some X_1, X_2, …, X_m, with X_i having distribution function F_i(x), in increasing order denoted as X_(1)≤ X_(2)≤…≤ X_(m). In this case, X_(i) is not necessarily associated with F_i(x). Assuming one can simulate values from each distribution, one can generate such “non-iid" order statistics by simulating X_i from F_i, for i=1,2,…, m, and simply putting 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 F_i(x). This problem arises in Bayesian principal components analysis (BPCA) where the X_i are ordered eigenvalues that are a posteriori independent but not identically distributed. In this paper, 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.

• 3 publications
• 3 publications
research
05/16/2022

### Order statistics from exchangeable random variables are always sufficient

Let (X_1,…,X_n) be an exchangeable random vector with distribution funct...
research
07/05/2023

### On multivariate orderings of some general ordered random vectors

Ordered random vectors are frequently encountered in many problems. The ...
research
06/13/2018

### Optimal moment inequalities for order statistics from nonnegative random variables

We obtain the best possible upper bounds for the moments of a single ord...
research
05/10/2019

### Generating Random Samples from Non-Identical Truncated Order Statistics

We provide an efficient algorithm to generate random samples from the bo...
research
12/21/2018

### Efficient Calculation of the Joint Distribution of Order Statistics

We consider the problem of computing the joint distribution of order sta...
research
06/15/2021

### Diagonal sections of copulas, multivariate conditional hazard rates and distributions of order statistics for minimally stable lifetimes

As a motivating problem, we aim to study some special aspects of the mar...
research
03/11/2020

### Lifted samplers for partially ordered discrete state-spaces

A technique called lifting is employed in practice for avoiding that the...