Dataset Information

Identifying the sources of structural sensitivity in partially specified biological models.

ABSTRACT: Biological systems are characterised by a high degree of uncertainty and complexity, which implies that exact mathematical equations to describe biological processes cannot generally be justified. Moreover, models can exhibit sensitivity to the precise formulations of their component functions-a property known as structural sensitivity. Structural sensitivity can be revealed and quantified by considering partially specified models with uncertain functions, but this goes beyond well-established, parameter-based sensitivity analysis, and currently presents a mathematical challenge. Here we build upon previous work in this direction by addressing the crucial question of identifying the processes which act as the major sources of model uncertainty and those which are less influential. To achieve this goal, we introduce two related concepts: (1) the gradient of structural sensitivity, accounting for errors made in specifying unknown functions, and (2) the partial degree of sensitivity with respect to each function, a global measure of the uncertainty due to possible variation of the given function while the others are kept fixed. We propose an iterative framework of experiments and analysis to inform a heuristic reduction of structural sensitivity in a model. To demonstrate the framework introduced, we investigate the sources of structural sensitivity in a tritrophic food chain model.

SUBMITTER: Adamson MW

PROVIDER: S-EPMC7547730 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identifying the sources of structural sensitivity in partially specified biological models.

Adamson Matthew W MW Morozov Andrew Yu AY

Scientific reports 20201009 1

Biological systems are characterised by a high degree of uncertainty and complexity, which implies that exact mathematical equations to describe biological processes cannot generally be justified. Moreover, models can exhibit sensitivity to the precise formulations of their component functions-a property known as structural sensitivity. Structural sensitivity can be revealed and quantified by considering partially specified models with uncertain functions, but this goes beyond well-established, ...[more]

PMID: 33037267

Similar Datasets

Project description:Statistical models built using different data sources and methods can exhibit conflicting patterns. We used the northern stock of black sea bass (Centropristis striata) as a case study to assess the impacts of using different fisheries data sources and laboratory-derived physiological metrics in the development of thermal habitat models for marine fishes. We constructed thermal habitat models using generalized additive models (GAMs) based on various fisheries datasets as input, including the NOAA Northeast Fisheries Science Center (NEFSC) bottom trawl surveys, various inshore fisheries-independent trawl surveys (state waters), NEFSC fisheries-dependent observer data, and laboratory-based physiological metrics. We compared each model's GAM response curve and coupled them to historical ocean conditions in the U.S. Northeast Shelf using bias-corrected ocean temperature output from a regional ocean model. Thermal habitat models based on shelf-wide data (NEFSC fisheries-dependent observer data and fisheries-independent spring and fall surveys) explained the most variation in black sea bass presence/absence data at ~15% deviance explained. Models based on a narrower range of sampled thermal habitat from inshore survey data in the Northeast Area Monitoring and Assessment Program (NEAMAP) and the geographically isolated Long Island Sound data performed poorly. All models had similar lower thermal limits around 8.5℃, but thermal optima, when present, ranged from 16.7 to 24.8℃. The GAMs could reliably predict habitat from years excluded from model training, but due to strong seasonal temperature fluctuations in the region, could not be used to predict habitat in seasons excluded from training. We conclude that survey data source can greatly impact development and interpretation of thermal habitat models for marine fishes. We suggest that model development be based on data sources that sample the widest range of ocean temperature and physical habitat throughout multiple seasons when possible, and encourage thorough consideration of how data gaps may influence model uncertainty.

Project description:The inference of reaction rate parameters in biochemical network models from time series concentration data is a central task in computational systems biology. Under the assumption of well mixed conditions the network dynamics are typically described by the chemical master equation, the Fokker Planck equation, the linear noise approximation or the macroscopic rate equation. The inverse problem of estimating the parameters of the underlying network model can be approached in deterministic and stochastic ways, and available methods often compare individual or mean concentration traces obtained from experiments with theoretical model predictions when maximizing likelihoods, minimizing regularized least squares functionals, approximating posterior distributions or sequentially processing the data. In this article we assume that the biological reaction network can be observed at least partially and repeatedly over time such that sample moments of species molecule numbers for various time points can be calculated from the data. Based on the chemical master equation we furthermore derive closed systems of parameter dependent nonlinear ordinary differential equations that predict the time evolution of the statistical moments. For inferring the reaction rate parameters we suggest to not only compare the sample mean with the theoretical mean prediction but also to take the residual of higher order moments explicitly into account. Cost functions that involve residuals of higher order moments may form landscapes in the parameter space that have more pronounced curvatures at the minimizer and hence may weaken or even overcome parameter sloppiness and uncertainty. As a consequence both deterministic and stochastic parameter inference algorithms may be improved with respect to accuracy and efficiency. We demonstrate the potential of moment fitting for parameter inference by means of illustrative stochastic biological models from the literature and address topics for future research.

Dataset Information

Identifying the sources of structural sensitivity in partially specified biological models.

Publications

Identifying the sources of structural sensitivity in partially specified biological models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets