Hierarchical Bayesian Joint Modeling of Athlete Performance Trajectories and Latent Peak Ages
Abstract
Determining when athletes reach peak performance is a central question in sports science, with implications for talent identification, training design, and cross-discipline comparisons. Statistical inference is challenging because peak age is not directly observed: competition records are noisy and selective measurements of underlying performance capacity and are often limited to incomplete observation windows. We develop a hierarchical Bayesian joint modeling framework that links a longitudinal model for athlete-specific performance trajectories to a regression submodel for latent peak age. The longitudinal component models observed competition performances as noisy measurements of athlete-specific performance trajectories over age. The peak-age component then treats the maximizer of each trajectory as a latent response in a parametric accelerated failure time (AFT) model, allowing athlete-level covariates to explain earlier or later peak age. The hierarchical structure borrows information across athletes, while joint posterior inference propagates uncertainty from trajectory reconstruction to both individual peak-age estimates and population-level covariate effects. Simulation studies show that, as observation windows become less informative, the joint model preserves substantially smaller regression bias and reflects the information loss through wider uncertainty, whereas pseudo-peak and two-stage alternatives exhibit larger bias. An application to World Athletics ranking data illustrates how the joint model estimates discipline-specific peak-age patterns and covariate effects under incomplete and irregular observation windows. Beyond sports, the proposed framework applies broadly to problems that couple longitudinal trajectories with function-derived latent features.
Related
- Bayesian Hierarchical Modeling Approaches for Combining Information From Multiple Data Sources to Produce Annual Estimates of National Immunization Coverage.
- Bayesian Hierarchical Models for Longitudinal Meta-Analysis in Sports Science: A Practical Tutorial Using brms
- A Bayesian hierarchical model for meta-analysis
- Behavioral Mapping of Hierarchical Dynamical Models Accelerates Bayesian Parameter Inference
- Generative modeling via hierarchical tensor sketching
BibTeX
@misc{peng2026hierarchical,
title = {{Hierarchical Bayesian Joint Modeling of Athlete Performance Trajectories and Latent Peak Ages}},
author = {K. Ken Peng and Tianyu Guan and Xiong Zhao},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.35473},
eprint = {2610.35473},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.35473},
}
RIS
TY - UNPB
TI - Hierarchical Bayesian Joint Modeling of Athlete Performance Trajectories and Latent Peak Ages
AU - K. Ken Peng
AU - Tianyu Guan
AU - Xiong Zhao
PY - 2026
DA - 2026/10/26
PB - arXiv
DO - 10.48550/arxiv.2610.35473
UR - https://arxiv.org/abs/2610.35473
AB - Determining when athletes reach peak performance is a central question in sports science, with implications for talent identification, training design, and cross-discipline comparisons. Statistical inference is challenging because peak age is not directly observed: competition records are noisy and selective measurements of underlying performance capacity and are often limited to incomplete observation windows. We develop a hierarchical Bayesian joint modeling framework that links a longitudinal model for athlete-specific performance trajectories to a regression submodel for latent peak age. The longitudinal component models observed competition performances as noisy measurements of athlete-specific performance trajectories over age. The peak-age component then treats the maximizer of each trajectory as a latent response in a parametric accelerated failure time (AFT) model, allowing athlete-level covariates to explain earlier or later peak age. The hierarchical structure borrows information across athletes, while joint posterior inference propagates uncertainty from trajectory reconstruction to both individual peak-age estimates and population-level covariate effects. Simulation studies show that, as observation windows become less informative, the joint model preserves substantially smaller regression bias and reflects the information loss through wider uncertainty, whereas pseudo-peak and two-stage alternatives exhibit larger bias. An application to World Athletics ranking data illustrates how the joint model estimates discipline-specific peak-age patterns and covariate effects under incomplete and irregular observation windows. Beyond sports, the proposed framework applies broadly to problems that couple longitudinal trajectories with function-derived latent features.
ER -