Bayesian inference using qualitative observations of underlying continuous variables.
Ontology highlight
ABSTRACT: Motivation
Recent work has demonstrated the feasibility of using non-numerical, qualitative data to parameterize mathematical models. However, uncertainty quantification (UQ) of such parameterized models has remained challenging because of a lack of a statistical interpretation of the objective functions used in optimization.Results
We formulated likelihood functions suitable for performing Bayesian UQ using qualitative observations of underlying continuous variables or a combination of qualitative and quantitative data. To demonstrate the resulting UQ capabilities, we analyzed a published model for immunoglobulin E (IgE) receptor signaling using synthetic qualitative and quantitative datasets. Remarkably, estimates of parameter values derived from the qualitative data were nearly as consistent with the assumed ground-truth parameter values as estimates derived from the lower throughput quantitative data. These results provide further motivation for leveraging qualitative data in biological modeling.Availability and implementation
The likelihood functions presented here are implemented in a new release of PyBioNetFit, an open-source application for analyzing Systems Biology Markup Language- and BioNetGen Language-formatted models, available online at www.github.com/lanl/PyBNF.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Mitra ED
PROVIDER: S-EPMC7214020 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Selenium (Se) is an essential dietary element for humans and animals, and the atmosphere is an important source of Se to soils. However, estimates of global atmospheric Se fluxes are highly uncertain. To constrain these uncertainties, we use a global model of atmospheric Se cycling and a database of more than 600 sites where Se in aerosol has been measured. Applying Bayesian inference techniques, we determine the probability distributions of global Se emissions from the four major sources: anthropogenic activities, volcanoes, marine biosphere, and terrestrial biosphere. Between 29 and 36 Gg of Se are emitted to the atmosphere every year, doubling previous estimates of emissions. Using emission parameters optimized by aerosol network measurements, our model shows good agreement with the aerosol Se observations (R2 = 0.66), as well as with independent aerosol (0.59) and wet deposition measurements (0.57). Both model and measurements show a decline in Se over North America in the last two decades because of changes in technology and energy policy. Our results highlight the role of the ocean as a net atmospheric Se sink, with around 7 Gg yr-1 of Se transferred from land through the atmosphere. The constrained Se emissions represent a substantial step forward in understanding the global Se cycle.
| S-EPMC7301612 | biostudies-literature
Project description:When presented with complex rhythmic auditory stimuli, humans are able to track underlying temporal structure (e.g., a "beat"), both covertly and with their movements. This capacity goes far beyond that of a simple entrained oscillator, drawing on contextual and enculturated timing expectations and adjusting rapidly to perturbations in event timing, phase, and tempo. Previous modeling work has described how entrainment to rhythms may be shaped by event timing expectations, but sheds little light on any underlying computational principles that could unify the phenomenon of expectation-based entrainment with other brain processes. Inspired by the predictive processing framework, we propose that the problem of rhythm tracking is naturally characterized as a problem of continuously estimating an underlying phase and tempo based on precise event times and their correspondence to timing expectations. We present two inference problems formalizing this insight: PIPPET (Phase Inference from Point Process Event Timing) and PATIPPET (Phase and Tempo Inference). Variational solutions to these inference problems resemble previous "Dynamic Attending" models of perceptual entrainment, but introduce new terms representing the dynamics of uncertainty and the influence of expectations in the absence of sensory events. These terms allow us to model multiple characteristics of covert and motor human rhythm tracking not addressed by other models, including sensitivity of error corrections to inter-event interval and perceived tempo changes induced by event omissions. We show that positing these novel influences in human entrainment yields a range of testable behavioral predictions. Guided by recent neurophysiological observations, we attempt to align the phase inference framework with a specific brain implementation. We also explore the potential of this normative framework to guide the interpretation of experimental data and serve as building blocks for even richer predictive processing and active inference models of timing.
| S-EPMC8216548 | biostudies-literature
Project description:[This corrects the article DOI: 10.1371/journal.pcbi.1009025.].
| S-EPMC8668110 | biostudies-literature
Project description:Detailed demographic data on wild Asian elephants have been difficult to collect due to habitat characteristics of much of the species' remaining range. Such data, however, are critical for understanding and modeling population processes in this endangered species. We present data from six years of an ongoing study of Asian elephants (Elephas maximus) in Uda Walawe National Park, Sri Lanka. This relatively undisturbed population numbering over one thousand elephants is individually monitored, providing cohort-based information on mortality and reproduction. Reproduction was seasonal, such that most births occurred during the long inter-monsoon dry season and peaked in May. During the study, the average age at first reproduction was 13.4 years and the 50(th) percentile inter-birth interval was approximately 6 years. Birth sex ratios did not deviate significantly from parity. Fecundity was relatively stable throughout the observed reproductive life of an individual (ages 11-60), averaging between 0.13-0.17 female offspring per individual per year. Mortalities and injuries based on carcasses and disappearances showed that males were significantly more likely than females to be killed or injured through anthropogenic activity. Overall, however, most observed injuries did not appear to be fatal. This population exhibits higher fecundity and density relative to published estimates on other Asian elephant populations, possibly enhanced by present range constriction. Understanding the factors responsible for these demographic dynamics can shed insight on the future needs of this elephant population, with probable parallels to other populations in similar settings.
| S-EPMC3869725 | biostudies-literature
Project description:In genomic studies with both genotypes and gene or protein expression profile available, causal effects of gene or protein on clinical outcomes can be inferred through using genetic variants as instrumental variables (IVs). The goal of introducing IV is to remove the effects of unobserved factors that may confound the relationship between the biomarkers and the outcome. A valid inference under the IV framework requires pairwise associations and pathway exclusivity. Among these assumptions, the IV expression association needs to be strong for the casual effect estimates to be unbiased. However, a small number of single nucleotide polymorphisms (SNPs) often provide limited explanation of the variability in the gene or protein expression and can only serve as weak IVs. In this study, we propose to replace SNPs with haplotypes as IVs to increase the variant-expression association and thus improve the casual effect inference of the expression. In the classical two-stage procedure, we developed a haplotype regression model combined with a model selection procedure to identify optimal instruments. The performance of the new method was evaluated through simulations and compared with the IV approaches using observed multiple SNPs. Our results showed the gain of power to detect a causal effect of gene or protein on the outcome using haplotypes compared with using only observed SNPs, under either complete or missing genotype scenarios. We applied our proposed method to a study of the effect of interleukin-1 beta (IL-1?) protein expression on the 90-day survival following sepsis and found that overly expressed IL-1? is likely to increase mortality.
| S-EPMC5625343 | biostudies-literature
Project description:We study the application of a widely used ordinal regression model, the cumulative probability model (CPM), for continuous outcomes. Such models are attractive for the analysis of continuous response variables because they are invariant to any monotonic transformation of the outcome and because they directly model the cumulative distribution function from which summaries such as expectations and quantiles can easily be derived. Such models can also readily handle mixed type distributions. We describe the motivation, estimation, inference, model assumptions, and diagnostics. We demonstrate that CPMs applied to continuous outcomes are semiparametric transformation models. Extensive simulations are performed to investigate the finite sample performance of these models. We find that properly specified CPMs generally have good finite sample performance with moderate sample sizes, but that bias may occur when the sample size is small. Cumulative probability models are fairly robust to minor or moderate link function misspecification in our simulations. For certain purposes, the CPMs are more efficient than other models. We illustrate their application, with model diagnostics, in a study of the treatment of HIV. CD4 cell count and viral load 6 months after the initiation of antiretroviral therapy are modeled using CPMs; both variables typically require transformations, and viral load has a large proportion of measurements below a detection limit.
| S-EPMC5675816 | biostudies-literature
Project description:The study of animal movement is challenging because movement is a process modulated by many factors acting at different spatial and temporal scales. In order to describe and analyse animal movement, several models have been proposed which differ primarily in the temporal conceptualization, namely continuous and discrete time formulations. Naturally, animal movement occurs in continuous time but we tend to observe it at fixed time intervals. To account for the temporal mismatch between observations and movement decisions, we used a state-space model where movement decisions (steps and turns) are made in continuous time. That is, at any time there is a non-zero probability of making a change in movement direction. The movement process is then observed at regular time intervals. As the likelihood function of this state-space model turned out to be intractable yet simulating data is straightforward, we conduct inference using different variations of Approximate Bayesian Computation (ABC). We explore the applicability of this approach as a function of the discrepancy between the temporal scale of the observations and that of the movement process in a simulation study. Simulation results suggest that the model parameters can be recovered if the observation time scale is moderately close to the average time between changes in movement direction. Good estimates were obtained when the scale of observation was up to five times that of the scale of changes in direction. We demonstrate the application of this model to a trajectory of a sheep that was reconstructed in high resolution using information from magnetometer and GPS devices. The state-space model used here allowed us to connect the scales of the observations and movement decisions in an intuitive and easy to interpret way. Our findings underscore the idea that the time scale at which animal movement decisions are made needs to be considered when designing data collection protocols. In principle, ABC methods allow to make inferences about movement processes defined in continuous time but in terms of easily interpreted steps and turns.
| S-EPMC7020826 | biostudies-literature
Project description:This paper addresses the problem of creating probabilistic brain atlases from manually labeled training data. Probabilistic atlases are typically constructed by counting the relative frequency of occurrence of labels in corresponding locations across the training images. However, such an "averaging" approach generalizes poorly to unseen cases when the number of training images is limited, and provides no principled way of aligning the training datasets using deformable registration. In this paper, we generalize the generative image model implicitly underlying standard "average" atlases, using mesh-based representations endowed with an explicit deformation model. Bayesian inference is used to infer the optimal model parameters from the training data, leading to a simultaneous group-wise registration and atlas estimation scheme that encompasses standard averaging as a special case. We also use Bayesian inference to compare alternative atlas models in light of the training data, and show how this leads to a data compression problem that is intuitive to interpret and computationally feasible. Using this technique, we automatically determine the optimal amount of spatial blurring, the best deformation field flexibility, and the most compact mesh representation. We demonstrate, using 2-D training datasets, that the resulting models are better at capturing the structure in the training data than conventional probabilistic atlases. We also present experiments of the proposed atlas construction technique in 3-D, and show the resulting atlases' potential in fully-automated, pulse sequence-adaptive segmentation of 36 neuroanatomical structures in brain MRI scans.
| S-EPMC3274721 | biostudies-literature
Project description:In this paper, we present the implementation of two types of Bayesian inference problems to demonstrate the potential of building probabilistic algorithms in hardware using single set of building blocks with the ability to perform these computations in real time. The first implementation, referred to as the BEAST (Bayesian Estimation and Stochastic Tracker), demonstrates a simple problem where an observer uses an underlying Hidden Markov Model (HMM) to track a target in one dimension. In this implementation, sensors make noisy observations of the target position at discrete time steps. The tracker learns the transition model for target movement, and the observation model for the noisy sensors, and uses these to estimate the target position by solving the Bayesian recursive equation online. We show the tracking performance of the system and demonstrate how it can learn the observation model, the transition model, and the external distractor (noise) probability interfering with the observations. In the second implementation, referred to as the Bayesian INference in DAG (BIND), we show how inference can be performed in a Directed Acyclic Graph (DAG) using stochastic circuits. We show how these building blocks can be easily implemented using simple digital logic gates. An advantage of the stochastic electronic implementation is that it is robust to certain types of noise, which may become an issue in integrated circuit (IC) technology with feature sizes in the order of tens of nanometers due to their low noise margin, the effect of high-energy cosmic rays and the low supply voltage. In our framework, the flipping of random individual bits would not affect the system performance because information is encoded in a bit stream.
| S-EPMC4796016 | biostudies-literature
Project description:Haemagglutination inhibition (HI) assays are typically used for comparing and characterizing influenza viruses. Data obtained from the assays (titres) are used quantitatively to determine antigenic differences between influenza strains. However, the use of these titres has been criticized as they sometimes fail to capture accurate antigenic differences between strains. Our previous analytical work revealed how antigenic and non-antigenic variables contribute to the titres. Building on this previous work, we have developed a Bayesian method for decoupling antigenic and non-antigenic contributions to the titres in this paper. We apply this method to a compendium of HI titres of influenza A (H3N2) viruses curated from 1968 to 2016. Remarkably, the results of this fit indicate that the non-antigenic variable, which is inversely correlated with viral avidity for the red blood cells used in HI assays, oscillates during the course of influenza virus evolution, with a period that corresponds roughly to the timescale on which antigenic variants replace each other. Together, the results suggest that the new Bayesian method is applicable to the analysis of long-term dynamics of both antigenic and non-antigenic properties of influenza virus.
| S-EPMC6083687 | biostudies-literature