Practical Recoverability and Numerical Identifiability of Kinetic Rate Constants from Noisy Time-Activity Curves in a Patient-Derived Renogram Model
Abstract
Close agreement between modeled and observed time-activity curves (TACs) does not necessarily imply accurate recovery of the underlying kinetic parameters. We evaluated this distinction using a partially observed four-compartment renogram model with a latent reservoir and observed left kidney, right kidney, and bladder TACs. Six de-identified renograms provided patient-derived reference parameter vectors for synthetic recovery experiments with added noise, together with analyses of temporal sampling, optimizer behavior, initialization, and numerical identifiability. Median total NRMSE values relative to the clean synthetic TACs were on the order of 10 − 3 , indicating excellent TAC reconstruction. Parameter recovery was substantially less uniform: some rate constants were recovered accurately, whereas others exhibited multiplicative errors of tens to hundreds, repeated boundary contact, or compensation involving the latent reservoir. Repeated optimizations also converged consistently to incorrect parameter vectors while preserving excellent TAC fits. The scaled residual Jacobian was numerically full rank but severely ill-conditioned, and several 95% profile likelihood intervals remained open at the optimization bounds. These findings demonstrate that TAC fit quality alone is insufficient to assess parameter reliability. TAC reconstruction, parameter recovery, numerical rank, boundary behavior, optimizer stability, and profile likelihoods provide complementary evidence and should be interpreted together when evaluating partially observed kinetic models.
Related
- Non-identifiability of autochthonous chikungunya risk from surveillance data, Japan, 2011–2025
- Randomized approaches for identifiability analysis and model reduction in fibrin polymerization kinetics
- Positive Latent Suppression in Logistic Neural Differential Equations: Trajectory Reconstruction, Identifiability and Population Hindcasts
- A Physics-Constrained Marcus-Theory Optimization Framework Restores Parameter Identifiability in Deep-Pocket Oncogenic Drug Kinetics: A KRAS G12C and STAT3 Case Study
- Pediatric Headache: A Practical Approach to Clinical Diagnosis and Management
BibTeX
@article{nguyen2026practical,
title = {{Practical Recoverability and Numerical Identifiability of Kinetic Rate Constants from Noisy Time-Activity Curves in a Patient-Derived Renogram Model}},
author = {Brandon Nguyen MS and Jackson Penning MD, MS and Zuzan Cayci MD and Farhad Jafari PhD},
journal = {Molecular Imaging},
year = {2026},
month = oct,
doi = {10.1177/15353508261492803},
url = {https://doi.org/10.1177/15353508261492803},
note = {CC BY},
}
RIS
TY - JOUR
TI - Practical Recoverability and Numerical Identifiability of Kinetic Rate Constants from Noisy Time-Activity Curves in a Patient-Derived Renogram Model
AU - Brandon Nguyen MS
AU - Jackson Penning MD, MS
AU - Zuzan Cayci MD
AU - Farhad Jafari PhD
PY - 2026
DA - 2026/10/23
JO - Molecular Imaging
DO - 10.1177/15353508261492803
UR - https://doi.org/10.1177/15353508261492803
AB - Close agreement between modeled and observed time-activity curves (TACs) does not necessarily imply accurate recovery of the underlying kinetic parameters. We evaluated this distinction using a partially observed four-compartment renogram model with a latent reservoir and observed left kidney, right kidney, and bladder TACs. Six de-identified renograms provided patient-derived reference parameter vectors for synthetic recovery experiments with added noise, together with analyses of temporal sampling, optimizer behavior, initialization, and numerical identifiability. Median total NRMSE values relative to the clean synthetic TACs were on the order of 10 − 3 , indicating excellent TAC reconstruction. Parameter recovery was substantially less uniform: some rate constants were recovered accurately, whereas others exhibited multiplicative errors of tens to hundreds, repeated boundary contact, or compensation involving the latent reservoir. Repeated optimizations also converged consistently to incorrect parameter vectors while preserving excellent TAC fits. The scaled residual Jacobian was numerically full rank but severely ill-conditioned, and several 95% profile likelihood intervals remained open at the optimization bounds. These findings demonstrate that TAC fit quality alone is insufficient to assess parameter reliability. TAC reconstruction, parameter recovery, numerical rank, boundary behavior, optimizer stability, and profile likelihoods provide complementary evidence and should be interpreted together when evaluating partially observed kinetic models.
SN - 1536-0121
ER -