COVID-19 outbreak in Wuhan demonstrates the limitations of publicly available case numbers for epidemiological modeling

Mathematical models are standard tools for understanding the underlying mechanisms of biological systems. Generally, the parameters of these models are unknown and they need to be inferred from experimental data using statistical methods. Most common measurement techniques only provide relative information about the absolute molecular state and often data is noise-corrupted. Therefore, introducing scaling and noise parameters in the model observables is necessary. Since frequently these parameters are also unknown, the dimensionality of the estimation problem is augmented. Sampling methods are widely used in systems biology to assess parameter and prediction uncertainties. However, the evaluation of sampling methods is usually demanding and often on the border of computational feasibility. Hence, efficient sampling algorithms are required.We propose a marginal sampling scheme for estimating the parameter uncertainties of mechanistic models with relative data. We integrate out the scaling and noise parameters from the original problem, leading to a dimension reduction of the parameter space. Herewith, only reaction rate constants have to be sampled. We find that the marginal sampling scheme retrieves the same parameter probability distributions and outperforms sampling on the full parameter space by substantially increasing the effective sample size and smoothing the transition probability between posterior modes.