Fitting and Selection of Within-Host Models for SARS-CoV-2 Infection
- 1University of Szeged
Abstract
The main question in model selection for longitudinal RT-qPCR data is whether improvements in fit justify the additional parameters of more complex within-host models. We compared a target-cell-limited model (TIV), an eclipse-phase extension (TEIV), and a fractional-order TIV using infection waves reconstructed from SARS-CoV-2 cycle-threshold (Ct) measurements in a National Basketball Association cohort. Each trajectory was fitted on the Ct scale by maximizing a right-censored Gaussian likelihood using hybrid global-local optimization. Model selection was based primarily on the small-sample corrected Akaike information criterion (AICc), computed from the censored likelihood; while robustness to the choice of information criterion was assessed using the Akaike information criterion (AIC) and Bayesian information criterion (BIC). Among 868 evaluable trajectories, the unpenalized likelihood ranking identified TEIV as the clear best-fitting model in 34.0% of cases, fr-TIV in 30.2%, and TIV in 3.3%; the remaining 32.5% were effectively tied. After accounting for model complexity, AICc clearly favored TIV in 97.8% of trajectories and fr-TIV in 0.5%, while 1.7% were effectively tied; TEIV was not a clear winner in any of the trajectories. Thus, for these short, sparsely sampled Ct-only trajectories, TIV provided the best-supported parsimonious effective description. This does not imply that eclipse-phase or memory-like dynamics are biologically absent, but rather that the simpler TIV model is favoured by standard information criteria for the limited information contained in these data.
Related
- Genome-wide within-host evolution during prolonged SARS-CoV-2 infection in a person with advanced HIV infection
- Host metabolic responses to SARS-CoV-2 and influenza viruses: parallels and contrasts
- Hyperinflammation and immunomodulation in influenza, SARS-CoV-2 and syncytial respiratory virus infection in adults: from innate immune responses and biomarkers to clinical and therapeutic implications
- Outcomes of SARS-CoV-2 infection post-solid organ transplantation in the cystic fibrosis population
- Diversification and hybridization within the Spirochona gemmipara complex driven by gammarid host range shifts
BibTeX
@misc{juhasz2026fitting,
title = {{Fitting and Selection of Within-Host Models for SARS-CoV-2 Infection}},
author = {N. Juhasz and F. A. Bartha and G. Rost and B. Das and D. Nagy and D. Selesi and D. Rajter-Ciric},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.15.759195},
url = {https://www.biorxiv.org/content/10.64898/2026.10.15.759195v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Fitting and Selection of Within-Host Models for SARS-CoV-2 Infection
AU - N. Juhasz
AU - F. A. Bartha
AU - G. Rost
AU - B. Das
AU - D. Nagy
AU - D. Selesi
AU - D. Rajter-Ciric
PY - 2026
DA - 2026/10/15
PB - bioRxiv
DO - 10.64898/2026.10.15.759195
UR - https://www.biorxiv.org/content/10.64898/2026.10.15.759195v1
AB - The main question in model selection for longitudinal RT-qPCR data is whether improvements in fit justify the additional parameters of more complex within-host models. We compared a target-cell-limited model (TIV), an eclipse-phase extension (TEIV), and a fractional-order TIV using infection waves reconstructed from SARS-CoV-2 cycle-threshold (Ct) measurements in a National Basketball Association cohort. Each trajectory was fitted on the Ct scale by maximizing a right-censored Gaussian likelihood using hybrid global-local optimization. Model selection was based primarily on the small-sample corrected Akaike information criterion (AICc), computed from the censored likelihood; while robustness to the choice of information criterion was assessed using the Akaike information criterion (AIC) and Bayesian information criterion (BIC). Among 868 evaluable trajectories, the unpenalized likelihood ranking identified TEIV as the clear best-fitting model in 34.0% of cases, fr-TIV in 30.2%, and TIV in 3.3%; the remaining 32.5% were effectively tied. After accounting for model complexity, AICc clearly favored TIV in 97.8% of trajectories and fr-TIV in 0.5%, while 1.7% were effectively tied; TEIV was not a clear winner in any of the trajectories. Thus, for these short, sparsely sampled Ct-only trajectories, TIV provided the best-supported parsimonious effective description. This does not imply that eclipse-phase or memory-like dynamics are biologically absent, but rather that the simpler TIV model is favoured by standard information criteria for the limited information contained in these data.
ER -