Conditional Rank Statistics

• Conditional Rank Statistics (CRS) are inference methods that transform data via conditional ranks, allowing analysis based on order and group structure.
• They enable nonparametric estimation, conditional regression, and robust variable selection through techniques like ranked set sampling and maximally selected rank statistics.
• CRS methods provide solid asymptotic theory and computational efficiency, making them valuable in economic, biomedical, and high-dimensional survival analyses.

Conditional Rank Statistics (CRS) refer generically to inferential or modeling approaches that utilize order, ranking, or permutation structure within or across strata defined by covariates or sample design, typically by conditioning either on ranks or covariates. CRS formalizes distributional and inferential properties (e.g., of statistics, estimators, or predictions) that arise when the usual parameters of interest are replaced by or connected to their conditional ranks, or when rank-based procedures are adapted to handle conditioning on observed group structure or covariate values. Applications span nonparametric estimation of distribution functions from ranked samples, conditional rank regression models, random forest split criteria, and modern high-dimensional coverage criteria indexed by rank.

1. Definitions and Conceptual Foundations

Several canonical constructions exemplify Conditional Rank Statistics:

• Conditional Sample Ranks: For a continuous outcome $Y$ and covariates $X$, the conditional rank is defined as $U = F_{Y|X}(Y|X)$, where $F_{Y|X}(y|x)$ denotes the conditional cumulative distribution function (CDF) of $Y$ given $X$ (Chernozhukov et al., 2024). The mapping $Y \mapsto U$ yields $U \sim \mathrm{Uniform}(0,1)$ conditionally on $X$.
• Ranked Set Sampling (RSS) and Judgment Post-Stratification (JPS): Observations $(X_i, R_i)$ are recorded, with $X$0 assigned as the (possibly error-prone) rank of $X$1 within a hypothetical sample of size $X$2. Conditionally, $X$3 is distributed as the $X$4th order statistic in a sample from an unknown $X$5 (Duembgen et al., 2013).
• Conditional Linear Rank Statistics in Random Forests: For structure discovery in survival forests, linear or maximally selected rank statistics are calculated conditional on covariate-induced splits, adapting classical rank-based hypothesis testing to recursive partitioning (Wright et al., 2016).
• Rank-Conditional Coverage: In high-dimensional multiple testing, confidence coverage is evaluated as a function of estimator rank within a vector of parameters, adjusting classical inference by conditioning on empirical ranking (Morrison et al., 2017).

These divergent settings share the core principle: inference or modeling is not on untransformed data, but on objects whose properties are determined or modulated by ranks, conditioning events, or permutation-invariant statistics.

2. Methodological Frameworks and Estimators

Ranked Set Sampling (RSS) and Distribution Function Estimation

In RSS or JPS, inference on an unknown continuous $X$6 is based on samples indexed by observed ranks. Let $X$7 and $X$8 the stratum-$X$9 empirical CDF. Three key estimators are studied (Duembgen et al., 2013):

Estimator	Definition
Stratified ($U = F_{Y|X}(Y|X)$0)	$U = F_{Y|X}(Y|X)$1
Nonparametric MLE ($U = F_{Y|X}(Y|X)$2)	Maximizes: $U = F_{Y|X}(Y|X)$3; $U = F_{Y|X}(Y|X)$4 is Beta CDF
Moment-based ($U = F_{Y|X}(Y|X)$5)	Solves: $U = F_{Y|X}(Y|X)$6 for $U = F_{Y|X}(Y|X)$7

Where $U = F_{Y|X}(Y|X)$8 is the cumulative distribution of the $U = F_{Y|X}(Y|X)$9th order statistic from a sample of size $F_{Y|X}(y|x)$0 (i.e., Beta$F_{Y|X}(y|x)$1).

Conditional Rank–Rank Regression (CRRR)

Let $F_{Y|X}(y|x)$2 be continuous random variables, $F_{Y|X}(y|x)$3 covariates. Conditional ranks are estimated for each $F_{Y|X}(y|x)$4 as $F_{Y|X}(y|x)$5 (and analogously $F_{Y|X}(y|x)$6 for $F_{Y|X}(y|x)$7), via distribution regression using a link function $F_{Y|X}(y|x)$8 (logit or probit) and basis functions $F_{Y|X}(y|x)$9. CRRR proceeds as an OLS regression or correlation of $Y$0 on $Y$1:

$Y$2

This $Y$3 estimates the average within-$Y$4 Spearman correlation between $Y$5 and $Y$6 (Chernozhukov et al., 2024).

Conditional and Maximally Selected Rank Statistics in Random Forests

In split variable selection for random survival forests, linear rank statistics $Y$7 and their standardized analogues are computed for each covariate $Y$8 and possible split $Y$9, using weights conditional on sample configuration. The maximally selected rank statistic $X$0 across split points $X$1 guards against bias toward high-cardinality covariates (Wright et al., 2016).

3. Asymptotic Theory and Large-Sample Properties

Rank-Based Distribution Estimation

Under mild regularity, as $X$2 (with $X$3 fixed), the process $X$4 for each estimator $X$5 exhibits a uniform linear expansion:

$X$6

with $X$7 independent Brownian bridges (one per rank stratum), and explicit $X$8 weights. In balanced sampling, $X$9 and $Y \mapsto U$0 are asymptotically equivalent, while $Y \mapsto U$1 attains the smallest asymptotic variance (Duembgen et al., 2013).

Conditional Rank–Rank Regression

The CRRR estimator $Y \mapsto U$2 is root-$Y \mapsto U$3 consistent and asymptotically normal:

$Y \mapsto U$4

where $Y \mapsto U$5 is characterized via influence function representation involving the distribution regression estimates (Chernozhukov et al., 2024). Standard errors can be estimated by exchangeable (weighted) bootstrap over the two-stage procedure.

Variable Selection in Random Survival Forests

The maximally selected rank statistic under null is asymptotically Gaussian, permitting analytic or permutation-based $Y \mapsto U$6-value approximations. Careful control via Brownian-bridge, Bonferroni, or multivariate Gaussian approximations reduces selection bias toward variables with many levels without sacrificing consistency (Wright et al., 2016).

4. Confidence Intervals, Bands, and Coverage Properties

Exact and Simultaneous Confidence Intervals

In ranked set data, conditional on the observed rank allocation, the sum $Y \mapsto U$7 ($Y \mapsto U$8 independently) permits construction of non-asymptotic, exact Clopper–Pearson–type confidence intervals for $Y \mapsto U$9. Simultaneous Kolmogorov–Smirnov-type bands are built by simulating the null process under $U \sim \mathrm{Uniform}(0,1)$0 uniform, conditional on ranks (Duembgen et al., 2013).

Rank-Conditional and Coverage-Adjusted Intervals

Rank conditional coverage (RCC) at rank $U \sim \mathrm{Uniform}(0,1)$1 is defined as $U \sim \mathrm{Uniform}(0,1)$2, where $U \sim \mathrm{Uniform}(0,1)$3 indexes estimator ordering. Bootstrapped (parametric or nonparametric) intervals constructed to target RCC achieve nominal coverage uniformly across all ranks, outperforming both marginal and false coverage-statement rate (FCR) controlling methods in addressing the "winner's curse" common in high-dimensional settings (Morrison et al., 2017).

5. Robustness and Sensitivity to Imperfect Ranking

Simulation studies in the context of RSS indicate that estimators react differently to violation of perfect ranking (signal-plus-noise with less than perfect correlation):

• $U \sim \mathrm{Uniform}(0,1)$4 (nonparametric MLE) exhibits increasing bias and mean-squared error under even mild misranking, particularly in the distribution tails.
• $U \sim \mathrm{Uniform}(0,1)$5 (moment) remains nearly unbiased and highly efficient under both perfect and modestly imperfect ranking; it outperforms $U \sim \mathrm{Uniform}(0,1)$6 in robustness.
• $U \sim \mathrm{Uniform}(0,1)$7 (stratified) remains unbiased but is less efficient (higher variance).

A plausible implication is that the moment estimator $U \sim \mathrm{Uniform}(0,1)$8 can be recommended as a compromise between efficiency under perfect conditions and resilience to rank errors (Duembgen et al., 2013).

6. Computational Algorithms and Practical Implementation

Procedures for CRS-based inference are computationally tractable. For RSS distribution estimation, $U \sim \mathrm{Uniform}(0,1)$9 requires only stratum CDFs, $X$0 involves root-finding in a monotone sum over Beta CDFs per $X$1, and $X$2 requires maximizing a strictly concave likelihood, both easily solved by bisection or Newton–Raphson. All estimators are piecewise constant on observed order statistics, allowing precomputation and vectorized algorithms in $X$3 or similar environments (Duembgen et al., 2013).

Distribution regression for CRRR uses binary-response GLMs on fine grids, tail-restricted extrapolations, and standard correlation computation. Bootstrap inference accommodates arbitrary exchangeable weights and is scalable (Chernozhukov et al., 2024).

Maximally selected rank statistics, when deployed in random forests, leverage analytic or fast permutation-based $X$4-value approximations. The "minLau" procedure (minimum of Brownian-bridge and Bonferroni approximations) provides near-unbiased, efficient variable selection even with many candidate splits and large-scale data (Wright et al., 2016).

Bootstrap computation for RCC-controlling intervals is implemented in the R package rcc, providing both parametric (Gaussian) and nonparametric variants. Usage involves resampling, re-ranking, quantile calculation for the error, and interval formation per rank (Morrison et al., 2017).

7. Empirical and Applied Perspectives

CRS and CRRR have demonstrated interpretative value, especially in economic and biomedical science. In intergenerational mobility, CRRR decomposes overall persistence into within-group and between-group components: within-group (conditional rank correlation) and between-group (difference with unconditional rank correlation). An application to Swiss administrative income data demonstrated that within-group persistence explained 62% of total persistence for sons and 52% for daughters (Chernozhukov et al., 2024).

Survival analysis with maximally selected rank statistics yields unbiased variable selection across covariate types and improved predictive performance, as evidenced in simulated and diverse real datasets, including gene expression and GWAS (Wright et al., 2016).

In high-dimensional inference, RCC-based intervals rectify selective inference coverage failures and are computationally practical, with code support in R for a variety of applied contexts (Morrison et al., 2017).

---

References

• Dümbgen & Zamanzade (2018): "Inference on a Distribution Function from Ranked Set Samples" (Duembgen et al., 2013)
• Chernozhukov et al.: "Conditional Rank-Rank Regression" (Chernozhukov et al., 2024)
• Hothorn & Lausen, Genz, et al.: "Unbiased split variable selection for random survival forests using maximally selected rank statistics" (Wright et al., 2016)
• Benjamini & Yekutieli, Weinstein & Reid, et al.: "Rank conditional coverage and confidence intervals in high dimensional problems" (Morrison et al., 2017)