Fast variance reduction for steady-state simulation and sensitivity analysis of stochastic chemical systems using
Andreas Milias-Argeitis1, John Lygeros2, Mustafa Khammash1
1Department of Biosystems Science and Engineering, ETH Zurich, 4058 Basel, Switzerland.
This study introduces a new variance reduction algorithm for estimating steady-state quantities in stochastic chemical kinetics. The novel method, inspired by queueing theory, efficiently calculates parametric sensitivities.
Area of Science:
- Computational Chemistry
- Chemical Kinetics
- Stochastic Systems
Background:
- Estimating steady-state quantities in stochastic chemical kinetics is crucial but often analytically intractable.
- Computational methods are typically required for these estimations.
- Existing methods may face challenges with efficiency and accuracy.
Purpose of the Study:
- To introduce a novel variance reduction algorithm for stochastic chemical kinetics.
- To improve the computational efficiency of estimating steady-state quantities.
- To assess the algorithm's performance in calculating steady-state parametric sensitivities.
Main Methods:
- Development of a new variance reduction algorithm inspired by queueing theory and shadow functions.
- Application of the algorithm to systems of stochastic chemical kinetics.
- Numerical evaluation using two distinct examples.
Main Results:
- The proposed algorithm demonstrates efficiency in estimating steady-state parametric sensitivities.
- The method's performance was compared favorably against other existing estimation techniques.
- Successful application to numerical examples validates the algorithm's utility.
Conclusions:
- The novel variance reduction algorithm offers an efficient computational approach for stochastic chemical kinetics.
- This method provides a valuable tool for analyzing complex chemical systems.
- Further research can explore its application to a broader range of kinetic models.
Related Concept Videos
Determination of Multiple Dosing Parameters: Steady-State, Minimum and Maximum Concentrations
Parameters Affecting Nonlinear Elimination: Zero-Order Input, First-Order Absorption and Two-Compartment Model
When a drug is administered through a constant intravenous infusion and eliminated via nonlinear pharmacokinetics, it follows zero-order input. For example, oral drugs undergo first-order absorption upon administration and are eliminated through nonlinear pharmacokinetics.
In the case of subcutaneously administered drugs,...
Mechanistic Models: Compartment Models in Algorithms for Numerical Problem Solving
In individual population analyses, different algorithms are employed, such as Cauchy's method, which uses a...
Testing a Claim about Standard Deviation
The hypothesis testing for the claim of population standard deviation (or variance) requires the data and samples to be random and unbiased. The population distribution also must be normal. There is no specific requirement on the sample size as the estimation is based on the chi-square distribution.
As a first step, the hypothesis (null and alternative) concerning the claim about...
Mechanistic Models: Compartment Models in Individual and Population Analysis
One-Compartment Open Model: Wagner-Nelson and Loo Riegelman Method for ka Estimation
On...


