Hyperphosphatemia in patients with renal failure is associated with increased vascular calcification and mortality. Hemodialysis is a conventional treatment for patients with hyperphosphatemia. Phosphate kinetics during hemodialysis may be described by a diffusion process and modeled by ordinary differential equations. We propose a Bayesian model approach for estimating patient-specific parameters for phosphate kinetics during hemodialysis. The Bayesian approach allows us to both analyze the full parameter space using uncertainty quantification and to compare two types of hemodialysis treatments, the conventional single-pass and the novel multiple-pass treatment. We validate and test our models on synthetic and real data. The results show limited identifiability of the model parameters when only single-pass data are available, and that the Bayesian model greatly reduces the relative standard deviation compared to existing estimates. Moreover, the analysis of the Bayesian models reveal improved estimates with reduced uncertainty when considering consecutive sessions and multiple-pass treatment compared to single-pass treatment.
Keywords: hemodialysis; mathematical modeling; parameter estimation; phosphate kinetics; uncertainty quantification.