When Does a Pathway Score Preserve a Mediated Effect? Aggregation Geometry, Minimax Adequacy, and Gap Inference for High-Dimensional Mediation
- 1Michigan State University
Abstract
High-dimensional mediators are routinely collapsed into low-dimensional scores before mediation analysis, but aggregation can change the mediated estimand rather than merely reduce noise. We characterize when a linear score preserves the full mediated effect. Exact preservation for every outcome model occurs precisely when the score contains the covariance-adjusted exposure-predictive mediator direction, the classical sufficient reduction for exposure under a Gaussian inverse-regression model. We introduce a distribution-free minimax adequacy loss equal to the fraction of exposure-relevant mediator information discarded by the score and exactly linked to overlap recovered by aggregation. A split-sample inverse-Wishart correction yields unbiased estimators of retained and full exposure information in moderate dimensions. Under spherically symmetric outcome-effect benchmarks, relative aggregation error has a Cauchy location-scale profile, and paired influence-function inference targets the realized gap for an observed outcome. In simulations, paired interval coverage ranged from 93.3% to 97.5%. In an outcome-free TCGA lung adenocarcinoma audit of 1,314 smoking-responsive genes, the pooled gene-mean loss was 0.811 (descriptive 95% gene-bootstrap interval 0.797 to 0.823), versus a dimension-only random-subspace benchmark of 0.913; retaining promoter and non-promoter means separately reduced pooled loss to 0.625. Independent half-sample screening was unstable, but top-100 and top-250 cross-screen pooled gene-mean losses ranged from 0.827 to 0.932. The framework separates structural adequacy, benchmark fragility, and outcome-specific discrepancy.
Related
- An M-learner approach for heterogeneous mediation analysis with high-dimensional omics mediators
- Flagella are required to activate expression of aggregation factors necessary for T6SS-mediated competition in host-like conditions
- Forest-weighted non-linear causal inference for high-dimensional and multimodal time series
- Hormone oscillations preserve cellular responsiveness to future physiological demands
- Patient-Derived Melanoma Organoids Preserve Tumor-Immune Heterogeneity and Reveal Context-Dependent Responses to Immune Checkpoint Blockade
BibTeX
@misc{hait2026when,
title = {{When Does a Pathway Score Preserve a Mediated Effect? Aggregation Geometry, Minimax Adequacy, and Gap Inference for High-Dimensional Mediation}},
author = {S. Hait},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.02.756279},
url = {https://www.biorxiv.org/content/10.64898/2026.10.02.756279v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - When Does a Pathway Score Preserve a Mediated Effect? Aggregation Geometry, Minimax Adequacy, and Gap Inference for High-Dimensional Mediation
AU - S. Hait
PY - 2026
DA - 2026/10/02
PB - bioRxiv
DO - 10.64898/2026.10.02.756279
UR - https://www.biorxiv.org/content/10.64898/2026.10.02.756279v1
AB - High-dimensional mediators are routinely collapsed into low-dimensional scores before mediation analysis, but aggregation can change the mediated estimand rather than merely reduce noise. We characterize when a linear score preserves the full mediated effect. Exact preservation for every outcome model occurs precisely when the score contains the covariance-adjusted exposure-predictive mediator direction, the classical sufficient reduction for exposure under a Gaussian inverse-regression model. We introduce a distribution-free minimax adequacy loss equal to the fraction of exposure-relevant mediator information discarded by the score and exactly linked to overlap recovered by aggregation. A split-sample inverse-Wishart correction yields unbiased estimators of retained and full exposure information in moderate dimensions. Under spherically symmetric outcome-effect benchmarks, relative aggregation error has a Cauchy location-scale profile, and paired influence-function inference targets the realized gap for an observed outcome. In simulations, paired interval coverage ranged from 93.3% to 97.5%. In an outcome-free TCGA lung adenocarcinoma audit of 1,314 smoking-responsive genes, the pooled gene-mean loss was 0.811 (descriptive 95% gene-bootstrap interval 0.797 to 0.823), versus a dimension-only random-subspace benchmark of 0.913; retaining promoter and non-promoter means separately reduced pooled loss to 0.625. Independent half-sample screening was unstable, but top-100 and top-250 cross-screen pooled gene-mean losses ranged from 0.827 to 0.932. The framework separates structural adequacy, benchmark fragility, and outcome-specific discrepancy.
ER -