A Bayesian approach to quantifying uncertainty from experimental noise in DEER spectroscopy.
Ontology highlight
ABSTRACT:
SUBMITTER: Edwards TH
PROVIDER: S-EPMC4996738 | biostudies-literature | 2016 Sep
REPOSITORIES: biostudies-literature
Publications
A Bayesian approach to quantifying uncertainty from experimental noise in DEER spectroscopy.
Edwards Thomas H TH Stoll Stefan S
Journal of magnetic resonance (San Diego, Calif. : 1997) 20160702
Double Electron-Electron Resonance (DEER) spectroscopy is a solid-state pulse Electron Paramagnetic Resonance (EPR) experiment that measures distances between unpaired electrons, most commonly between protein-bound spin labels separated by 1.5-8nm. From the experimental data, a distance distribution P(r) is extracted using Tikhonov regularization. The disadvantage of this method is that it does not directly provide error bars for the resulting P(r), rendering correct interpretation difficult. He ...[more]
Similar Datasets
Project description:The wiring diagram of the human brain can be described in terms of graph measures that characterize structural regularities. These measures require an estimate of whole-brain structural connectivity for which one may resort to deterministic or thresholded probabilistic streamlining procedures. While these procedures have provided important insights about the characteristics of human brain networks, they ultimately rely on unwarranted assumptions such as those of noise-free data or the use of an arbitrary threshold. Therefore, resulting structural connectivity estimates as well as derived graph measures fail to fully take into account the inherent uncertainty in the structural estimate. In this paper, we illustrate an easy way of obtaining posterior distributions over graph metrics using Bayesian inference. It is shown that this posterior distribution can be used to quantify uncertainty about graph-theoretical measures at the single subject level, thereby providing a more nuanced view of the graph-theoretical properties of human brain connectivity. We refer to this model-based approach to connectivity analysis as Bayesian connectomics.
| S-EPMC4189434 | biostudies-literature
Project description:Double electron-electron resonance (DEER) spectroscopy measures distance distributions between spin labels in proteins, yielding important structural and energetic information about conformational landscapes. Analysis of an experimental DEER signal in terms of a distance distribution is a nontrivial task due to the ill-posed nature of the underlying mathematical inversion problem. This work introduces a Bayesian probabilistic inference approach to analyze DEER data, using a multi-Gauss mixture model for the distance distribution. The method uses Markov chain Monte Carlo (MCMC) sampling to determine a posterior probability distribution over model parameter space. This distribution contains all the information available from the data, including a full quantification of the uncertainty about the parameters. The corresponding uncertainty about the distance distribution is captured via an ensemble of posterior predictive distributions. Several synthetic examples illustrate the method. An experimental example shows the importance of model checking and comparison using residual analysis and Bayes factors. Overall, the Bayesian approach allows for more robust inference about protein conformations from DEER spectroscopy.
| S-EPMC7846514 | biostudies-literature
Project description:Acceptance of novelty depends on the receiver's emotional state. This paper proposes a novel mathematical model for predicting emotions elicited by the novelty of an event under different conditions. It models two emotion dimensions, arousal and valence, and considers different uncertainty levels. A state transition from before experiencing an event to afterwards is assumed, and a Bayesian model estimates a posterior distribution as being proportional to the product of a prior distribution and a likelihood function. Our model uses Kullback-Leibler divergence of the posterior from the prior, which we termed information gain, to represent arousal levels because it corresponds to surprise, a high-arousal emotion, upon experiencing a novel event. Based on Berlyne's hedonic function, we formalized valence as a summation of reward and aversion systems that are modeled as sigmoid functions of information gain. We derived information gain as a function of prediction errors (i.e., differences between the mean of the posterior and the peak likelihood), uncertainty (i.e., variance of the prior that is proportional to prior entropy), and noise (i.e., variance of the likelihood function). This functional model predicted an interaction effect of prediction errors and uncertainty on information gain, which we termed the arousal crossover effect. This effect means that the greater the uncertainty, the greater the information gain for a small prediction error. However, for large prediction errors, greater uncertainty means a smaller information gain. To verify this effect, we conducted an experiment with participants who watched short videos in which different percussion instruments were played. We varied uncertainty levels by using familiar and unfamiliar instruments, and we varied prediction error magnitudes by including congruent or incongruent percussive sounds in the videos. Event-related potential P300 amplitudes and subjective reports of surprise in response to the percussive sounds were used as measures of arousal levels, and the findings supported the hypothesized arousal crossover effect. The concordance between our model's predictions and our experimental results suggests that Bayesian information gain can be decomposed into uncertainty and prediction errors and is a valid measure of emotional arousal. Our model's predictions of arousal may help identify positively accepted novelty.
| S-EPMC6353852 | biostudies-literature
Project description:Model averaging for dichotomous dose-response estimation is preferred to estimate the benchmark dose (BMD) from a single model, but challenges remain regarding implementing these methods for general analyses before model averaging is feasible to use in many risk assessment applications, and there is little work on Bayesian methods that include informative prior information for both the models and the parameters of the constituent models. This article introduces a novel approach that addresses many of the challenges seen while providing a fully Bayesian framework. Furthermore, in contrast to methods that use Monte Carlo Markov Chain, we approximate the posterior density using maximum a posteriori estimation. The approximation allows for an accurate and reproducible estimate while maintaining the speed of maximum likelihood, which is crucial in many applications such as processing massive high throughput data sets. We assess this method by applying it to empirical laboratory dose-response data and measuring the coverage of confidence limits for the BMD. We compare the coverage of this method to that of other approaches using the same set of models. Through the simulation study, the method is shown to be markedly superior to the traditional approach of selecting a single preferred model (e.g., from the U.S. EPA BMD software) for the analysis of dichotomous data and is comparable or superior to the other approaches.
| S-EPMC7722241 | biostudies-literature
Project description:MotivationGenome-wide association studies (GWAS) are used to discover genes underlying complex, heritable disorders for which less powerful study designs have failed in the past. The number of GWAS has skyrocketed recently with findings reported in top journals and the mainstream media. Microarrays are the genotype calling technology of choice in GWAS as they permit exploration of more than a million single nucleotide polymorphisms (SNPs) simultaneously. The starting point for the statistical analyses used by GWAS to determine association between loci and disease is making genotype calls (AA, AB or BB). However, the raw data, microarray probe intensities, are heavily processed before arriving at these calls. Various sophisticated statistical procedures have been proposed for transforming raw data into genotype calls. We find that variability in microarray output quality across different SNPs, different arrays and different sample batches have substantial influence on the accuracy of genotype calls made by existing algorithms. Failure to account for these sources of variability can adversely affect the quality of findings reported by the GWAS.ResultsWe developed a method based on an enhanced version of the multi-level model used by CRLMM version 1. Two key differences are that we now account for variability across batches and improve the call-specific assessment of each call. The new model permits the development of quality metrics for SNPs, samples and batches of samples. Using three independent datasets, we demonstrate that the CRLMM version 2 outperforms CRLMM version 1 and the algorithm provided by Affymetrix, Birdseed. The main advantage of the new approach is that it enables the identification of low-quality SNPs, samples and batches.AvailabilitySoftware implementing of the method described in this article is available as free and open source code in the crlmm R/BioConductor package.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC2804295 | biostudies-literature
Project description:The use of mixed models to determine narrow-sense heritability and related quantities such as SNP heritability has received much recent attention. Less attention has been paid to the inherent variability in these estimates. One approach for quantifying variability in estimates of heritability is a frequentist approach, in which heritability is estimated using maximum likelihood and its variance is quantified through an asymptotic normal approximation. An alternative approach is to quantify the uncertainty in heritability through its Bayesian posterior distribution. In this paper, we develop the latter approach, make it computationally efficient and compare it to the frequentist approach. We show theoretically that, for a sufficiently large sample size and intermediate values of heritability, the two approaches provide similar results. Using the Atherosclerosis Risk in Communities cohort, we show empirically that the two approaches can give different results and that the variance/uncertainty can remain large.
| S-EPMC4521294 | biostudies-literature
Project description:BackgroundAlthough microarrays are analysis tools in biomedical research, they are known to yield noisy output that usually requires experimental confirmation. To tackle this problem, many studies have developed rules for optimizing probe design and devised complex statistical tools to analyze the output. However, less emphasis has been placed on systematically identifying the noise component as part of the experimental procedure. One source of noise is the variance in probe binding, which can be assessed by replicating array probes. The second source is poor probe performance, which can be assessed by calibrating the array based on a dilution series of target molecules. Using model experiments for copy number variation and gene expression measurements, we investigate here a revised design for microarray experiments that addresses both of these sources of variance.ResultsTwo custom arrays were used to evaluate the revised design: one based on 25 mer probes from an Affymetrix design and the other based on 60 mer probes from an Agilent design. To assess experimental variance in probe binding, all probes were replicated ten times. To assess probe performance, the probes were calibrated using a dilution series of target molecules and the signal response was fitted to an adsorption model. We found that significant variance of the signal could be controlled by averaging across probes and removing probes that are nonresponsive or poorly responsive in the calibration experiment. Taking this into account, one can obtain a more reliable signal with the added option of obtaining absolute rather than relative measurements.ConclusionThe assessment of technical variance within the experiments, combined with the calibration of probes allows to remove poorly responding probes and yields more reliable signals for the remaining ones. Once an array is properly calibrated, absolute quantification of signals becomes straight forward, alleviating the need for normalization and reference hybridizations.
| S-EPMC3949741 | biostudies-literature
Project description:BACKGROUND: Locus heterogeneity, wherein a disease can be caused in different individuals by different genes and/or environmental factors, is a ubiquitous feature of complex traits. A Bayesian approach has been proposed to account for variable rates of heterogeneity across families in a parametric linkage analysis setup [Biswas and Lin: J Am Stat Assoc 2006;101:1341-1351]. As with any parametric approach, its application requires specification of the disease model, which limits its practical utility. METHODS: We address this limitation by proposing a Bayesian model averaging (BMA) approach. We consider a finite number of disease models and treat the model as an unknown parameter. In practice, we use simple single-locus disease models as various categories for model. RESULTS: Our simulations as well as analysis of Genetic Analysis Workshop 13 simulated data show that BMA retains at least 80% of the power that is obtained by analyzing under the true disease model. The coverage probability of interval for disease gene is maintained around the nominal level. Finally, we apply BMA to a Late-Onset Alzheimer's Disease dataset and find evidence for linkage on chromosomes 19, 9, and 21. CONCLUSION: We conclude that the BMA approach utilizing simple single-locus models for averaging is effective for mapping heterogeneous traits.
| S-EPMC2889260 | biostudies-literature
Project description:Adaptive radiations and mass extinctions are of critical importance in structuring terrestrial ecosystems. However, the causes and progress of these transitions often remain controversial, in part because of debates surrounding the completeness of the fossil record and biostratigraphy of the relevant fossil-bearing formations. The early-middle Permian, when a substantial faunal turnover in tetrapods coincided with a restructuring of the trophic structure of ecosystems, is such a time. Some have suggested the transition is obscured by a gap in the tetrapod fossil record (Olson's Gap), while others suggest a correlation between North American and Russian tetrapod-bearing formations allows the interval to be documented in detail. The latter biostratigraphic scheme has been used to support a mass extinction at this time (Olson's Extinction). Bayesian tip-dating methods used frequently in phylogenetics are employed to resolve this debate. Bayes factors are used to compare the results of analyses incorporating tip age priors based on different stratigraphic hypotheses, to show which stratigraphic scheme best fits the morphological data and phylogeny. Olson's Gap is rejected, and the veracity of Olson's Extinction is given further support. Tip-dating approaches have great potential to resolve debates surrounding the stratigraphic ages of critical formations where appropriate morphological data is available.
| S-EPMC7341920 | biostudies-literature
Project description:BACKGROUND:Exposure to both occupational and nonoccupational noise is recognized as a risk factor for noise-induced hearing loss (NIHL). Although audiologists routinely inquire regarding history of noise exposure, there are limited tools available for quantifying this history or for identifying those individuals who are at highest risk for NIHL. Identifying those at highest risk would allow hearing conservation activities to be focused on those individuals. PURPOSE:To develop a detailed, task-based questionnaire for quantifying an individual's annual noise exposure (ANE) arising from both occupational and nonoccupational sources (aim 1) and to develop a short screening tool that could be used to identify individuals at high risk of NIHL (aim 2). RESEARCH DESIGN:Review of relevant literature for questionnaire development followed by a cross-sectional descriptive and correlational investigation of the newly developed questionnaire and screening tool. STUDY SAMPLE:One hundred fourteen college freshmen completed the detailed questionnaire for estimating ANE (aim 1) and answered the potential screening questions (aim 2). An additional 59 adults participated in data collection where the accuracy of the screening tool was evaluated (aim 2). DATA COLLECTION AND ANALYSIS:In study aim 1, all participants completed the detailed questionnaire and the potential screening questions. Descriptive statistics were used to quantify participant participation in various noisy activities and their associated ANE estimates. In study aim 2, linear regression techniques were used to identify screening questions that could be used to predict a participant's estimated ANE. Clinical decision theory was then used to assess the accuracy with which the screening tool predicted high and low risk of NIHL in a new group of participants. RESULTS:Responses on the detailed questionnaire indicated that our sample of college freshmen reported high rates of participation in a variety of occupational and nonoccupational activities associated with high sound levels. Although participation rates were high, ANE estimates were below highest-risk levels for many participants because the frequency of participation in these activities was low in many cases. These data illustrate how the Noise Exposure Questionnaire (NEQ) could be used to provide detailed and specific information regarding an individual's exposure to noise. The results of aim 2 suggest that the screening tool, the 1-Minute Noise Screen, can be used to identify those participants with high- and low-risk noise exposure, allowing more in-depth assessment of noise exposure history to be targeted at those most at risk. CONCLUSIONS:The NEQ can be used to estimate an individual's ANE and the 1-Minute Noise Screen can be used to identify those participants at highest risk of NIHL. These tools allow audiologists to focus hearing conservation efforts on those individuals who are most in need of those services.
| S-EPMC5304605 | biostudies-literature