Markov chain Monte Carlo parameter estimation of the ODE compartmental cell growth model

We use a Markov chain Monte Carlo (MCMC) method to quantify uncertainty in parameters of the heterogeneous linear compartmental model of cell population growth, described by a system of ordinary differential equations. This model allows division number-dependent rates of cell proliferation and death and describes the rate of changes in the numbers of cells having undergone j divisions. The experimental data set specifies the following characteristics of the kinetics of human T lymphocyte proliferation assay in vitro: the total number of live cells and dead but not disintegrated cells and the number of cells divided j times. Our goal is to compare results of the MCMC analysis of the uncertainty in the best-fit parameter estimates with the ones obtained earlier, using the variance-covariance approach, the profile-likelihood based approach and the bootstrap technique. We show that the computed posterior probability density functions are Gaussian for most of the model parameters and they are close to Gaussian ones for other parameters except one. We present posterior uncertainty limits for the model solution and new observations.