Abstract
Dynamical models are widely used in systems biology to describe biological processes ranging from single cell transcription of genes to the tissue scale formation of gradients for cell guidance. One of the key issues for this class of models is the estimation of kinetic parameters from given measurement data, the so called parameter estimation. Measurement noise and the limited amount of data, give rise to uncertainty in estimates which can be captured in a probability density over the parameter space. Unfortunately, studying this probability density, using e.g. Markov chain Monte-Carlo, is often computationally demanding as it requires the repeated simulation of the underlying model. In the case of highly complex models, such as PDE models, this can render the study intractable. In this paper, we will present novel methods for analysis of such probability densities using networks of radial basis functions. We employed lattice generation algorithms, adaptive interacting particle sampling schemes as well as classical sampling schemes for the generation of approximation nodes coupled to the respective weighting scheme and compared their efficiency on different application examples. Our analysis showed that the novel method can yield an expected L 2 approximation error in marginals that is several orders of magnitude lower compared to classical approximations. This allows for a drastic reduction of the number of model evaluations. This facilitates the analysis of uncertainty for problems with high computational complexity. Finally, we successfully applied our method to a complex partial differential equation model for guided cell migration of dendritic cells.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Robert, C., Casella, G.: Monte Carlo Statistical Methods. Springer, New York (2004)
Haario, H., Laine, M., Mira, A., Saksman, E.: DRAM: Efficient Adaptive MCMC. Statistics and Computing 16(4), 339–354 (2006)
Schmidl, D., Czado, C., Hug, S., Theis, F.J.: A Vine-copula Based Adaptive MCMC Sampler for Efficient Inference of Dynamical Systems. Bayesian Analysis 8(1), 1–22 (2013)
MacKay, D.: Information Theory, Inference and Learning Algorithms (2003)
Jarner, S., Roberts, G.: Convergence of Heavy-tailed Monte Carlo Markov Chain Algorithms. Scandinavian Journal of Statistics 34(1994), 781–815 (2007)
Fasshauer, G.: Meshfree Approximation Methods with MATLAB. World Scientific (2007)
Sarra, S.A., Kansa, E.J.: Multiquadric Radial Basis Function Approximation Methods for the Numerical Solution of Partial Differential Equations. Advances in Computational Mechanics 2 (2009)
Müller, J.: Surrogate Model Algorithms for Computationally Expensive Black-Box Global Optimization Problems. PhD thesis (2012)
Potts, D., Steidl, G.: Fast Summation at Nonequispaced Knots by NFFT. SIAM Journal on Scientific Computing 24(6), 2013–2037 (2003)
Silverman, B.: Density Estimation for Statistics and Data Analysis. Chapman & Hall/CRC (1986)
Ramsay, P.H., Scott, D.W.: Multivariate Density Estimation, Theory, Practice, and Visualization. Wiley (1993)
Wu, Z., Schaback, R.: Local Error Estimates for Radial Basis Function Interpolation of Scattered Data. IMA Journal of Numerical Analysis 13(1), 1–15 (1993)
Fröhlich, F.: Approximation and Analysis of Probability Densities using Radial Basis Functions. Master’s thesis, Technische Universität München, Germany
Bungartz, H.J., Griebel, M.: Sparse grids. Acta Numerica 13, 147 (2004)
Torquato, S.: Reformulation of the Covering and Quantizer Problems as Ground States of Interacting Particles. Physical Review E 82(5), 1–52 (2010)
Reboux, S., Schrader, B., Sbalzarini, I.F.: A Self-Organizing Lagrangian Particle Method for Adaptive Resolution Advection-Diffusion Simulations. Journal of Computational Physics 231(9), 3623–3646 (2012)
Hock, S., Hasenauer, J., Theis, F.J.: Modeling of 2D Diffusion Processes based on Microscopy Data: Parameter Estimation and Practical Identifiability Analysis. BMC Bioinformatics 14(suppl.7) (2013)
Blizniouk, N., Ruppert, D., Shoemaker, C., Regis, R.: Bayesian Calibration of Computationally Expensive Models using Optimization and Radial Basis Function Approximation. Journal of Computational and Graphical Statistics 17(2) (2008)
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2014 Springer International Publishing Switzerland
About this paper
Cite this paper
Fröhlich, F., Hross, S., Theis, F.J., Hasenauer, J. (2014). Radial Basis Function Approximations of Bayesian Parameter Posterior Densities for Uncertainty Analysis. In: Mendes, P., Dada, J.O., Smallbone, K. (eds) Computational Methods in Systems Biology. CMSB 2014. Lecture Notes in Computer Science(), vol 8859. Springer, Cham. https://doi.org/10.1007/978-3-319-12982-2_6
Download citation
DOI: https://doi.org/10.1007/978-3-319-12982-2_6
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-12981-5
Online ISBN: 978-3-319-12982-2
eBook Packages: Computer ScienceComputer Science (R0)