Towards uncertainty quantification of a model for cancer-on-chip experiments

TL;DR

This paper develops a framework for quantifying uncertainty in a cancer-on-chip model using sensitivity analysis, Bayesian parameter estimation, and surrogate models to improve prediction reliability.

Contribution

It introduces a comprehensive uncertainty quantification approach combining sensitivity analysis, Bayesian inference, and surrogate modeling for cancer-on-chip simulations.

Findings

Parameter uncertainty significantly affects model predictions

Bayesian methods effectively estimate parameter distributions

Surrogate models accelerate uncertainty analysis

Abstract

This study is a first step towards using data-informed differential models to predict and control the dynamics of cancer-on-chip experiments. We consider a conceptualized one-dimensional device, containing a cancer and a population of white blood cells. The interaction between the cancer and the population of cells is modeled by a chemotaxis model inspired by Keller-Segel-type equations, which is solved by a Hybridized Discontinuous Galerkin method. Our goal is using (synthetic) data to tune the parameters of the governing equations and to assess the uncertainty on the predictions of the dynamics due to the residual uncertainty on the parameters remaining after the tuning procedure. To this end, we apply techniques from uncertainty quantification for parametric differential models. We first perform a global sensitivity analysis using both Sobol and Morris indices to assess how parameter uncertainty impacts model predictions, and fix the value of parameters with negligible impact. Subsequently, we conduct an inverse uncertainty quantification analysis by Bayesian techniques to compute a data-informed probability distribution of the remaining model parameters. Finally, we carry out a forward uncertainty quantification analysis to compute the impact of the updated (residual) parametric uncertainties on the quantities of interest of the model. The whole procedure is sped up by using surrogate models, based on sparse-grids, to approximate the mapping of the uncertain parameters to the quantities of interest.