@article{KnoechelKloftHuisinga2018, author = {Kn{\"o}chel, Jane and Kloft, Charlotte and Huisinga, Wilhelm}, title = {Understanding and reducing complex systems pharmacology models based on a novel input-response index}, series = {Journal of pharmacokinetics and pharmacodynamics}, volume = {45}, journal = {Journal of pharmacokinetics and pharmacodynamics}, number = {1}, publisher = {Springer Science + Business Media B.V.}, address = {New York}, issn = {1567-567X}, doi = {10.1007/s10928-017-9561-x}, pages = {139 -- 157}, year = {2018}, abstract = {A growing understanding of complex processes in biology has led to large-scale mechanistic models of pharmacologically relevant processes. These models are increasingly used to study the response of the system to a given input or stimulus, e.g., after drug administration. Understanding the input-response relationship, however, is often a challenging task due to the complexity of the interactions between its constituents as well as the size of the models. An approach that quantifies the importance of the different constituents for a given input-output relationship and allows to reduce the dynamics to its essential features is therefore highly desirable. In this article, we present a novel state- and time-dependent quantity called the input-response index that quantifies the importance of state variables for a given input-response relationship at a particular time. It is based on the concept of time-bounded controllability and observability, and defined with respect to a reference dynamics. In application to the brown snake venom-fibrinogen (Fg) network, the input-response indices give insight into the coordinated action of specific coagulation factors and about those factors that contribute only little to the response. We demonstrate how the indices can be used to reduce large-scale models in a two-step procedure: (i) elimination of states whose dynamics have only minor impact on the input-response relationship, and (ii) proper lumping of the remaining (lower order) model. In application to the brown snake venom-fibrinogen network, this resulted in a reduction from 62 to 8 state variables in the first step, and a further reduction to 5 state variables in the second step. We further illustrate that the sequence, in which a recursive algorithm eliminates and/or lumps state variables, has an impact on the final reduced model. The input-response indices are particularly suited to determine an informed sequence, since they are based on the dynamics of the original system. In summary, the novel measure of importance provides a powerful tool for analysing the complex dynamics of large-scale systems and a means for very efficient model order reduction of nonlinear systems.}, language = {en} } @article{RedmannFreitag2021, author = {Redmann, Martin and Freitag, Melina A.}, title = {Optimization based model order reduction for stochastic systems}, series = {Applied mathematics and computation}, volume = {398}, journal = {Applied mathematics and computation}, publisher = {Elsevier}, address = {New York}, issn = {0096-3003}, doi = {10.1016/j.amc.2020.125783}, pages = {18}, year = {2021}, abstract = {In this paper, we bring together the worlds of model order reduction for stochastic linear systems and H-2-optimal model order reduction for deterministic systems. In particular, we supplement and complete the theory of error bounds for model order reduction of stochastic differential equations. With these error bounds, we establish a link between the output error for stochastic systems (with additive and multiplicative noise) and modified versions of the H-2-norm for both linear and bilinear deterministic systems. When deriving the respective optimality conditions for minimizing the error bounds, we see that model order reduction techniques related to iterative rational Krylov algorithms (IRKA) are very natural and effective methods for reducing the dimension of large-scale stochastic systems with additive and/or multiplicative noise. We apply modified versions of (linear and bilinear) IRKA to stochastic linear systems and show their efficiency in numerical experiments.}, language = {en} }