Stochastic models of biochemical reaction networks are used for understanding the properties of molecular regulatory circuits in living cells. The state of the cell is defined by the number of copies of each molecular species in the model. The chemical master equation (CME) governs the time evolution of the the probability density function of the often high-dimensional state space. The CME is approximated by a partial differential equation (PDE), the Fokker-Planck equation and solved numerically. Direct solution of the CME rapidly becomes computationally expensive for increasingly complex biological models, since the state space grows exponentially with the number of dimensions. Adaptive numerical methods can be applied in time and space in the PDE framework, and error estimates of the approximate solutions are derived. A method for splitting the CME operator in order to apply the PDE approximation in a subspace of the state space is also developed. The performance is compared to the most widely spread alternative computational method.
Identifer | oai:union.ndltd.org:UPSALLA1/oai:DiVA.org:uu-8293 |
Date | January 2007 |
Creators | Sjöberg, Paul |
Publisher | Uppsala universitet, Avdelningen för teknisk databehandling, Uppsala universitet, Numerisk analys, Uppsala : Acta Universitatis Upsaliensis |
Source Sets | DiVA Archive at Upsalla University |
Language | English |
Detected Language | English |
Type | Doctoral thesis, comprehensive summary, info:eu-repo/semantics/doctoralThesis, text |
Format | application/pdf |
Rights | info:eu-repo/semantics/openAccess |
Relation | Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology, 1651-6214 ; 358 |
Page generated in 0.0079 seconds