A framework for quantifying uncertainty in DFT energy corrections.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang A
PROVIDER: S-EPMC8322326 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
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: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:In this article, we present a nanoelectromechanical system (NEMS) designed to detect changes in the Casimir energy. The Casimir effect is a result of the appearance of quantum fluctuations in an electromagnetic vacuum. Previous experiments have used nano- or microscale parallel plate capacitors to detect the Casimir force by measuring the small attractive force these fluctuations exert between the two surfaces. In this new set of experiments, we aim to directly detect the shifts in the Casimir energy in a vacuum due to the presence of the metallic parallel plates, one of which is a superconductor. A change in the Casimir energy of this configuration is predicted to shift the superconducting transition temperature (T c) because of the interaction between it and the superconducting condensation energy. In our experiment, we take a superconducting film, carefully measure its transition temperature, bring a conducting plate close to the film, create a Casimir cavity, and then measure the transition temperature again. The expected shifts are smaller than the normal shifts one sees in cycling superconducting films to cryogenic temperatures, so using a NEMS resonator in situ is the only practical way to obtain accurate, reproducible data. Using a thin Pb film and opposing Au surface, we observe no shift in T c >12 µK down to a minimum spacing of ~70 nm at zero applied magnetic field.
| S-EPMC7767790 | biostudies-literature
Project description:High-level assessments of climate change impacts aggregate multiple perils into a common framework. This requires incorporating multiple dimensions of uncertainty. Here we propose a methodology to transparently assess these uncertainties within the 'Reasons for Concern' framework, using extreme heat as a case study. We quantitatively discriminate multiple dimensions of uncertainty, including future vulnerability and exposure to changing climate hazards. High risks from extreme heat materialise after 1.5-2 °C and very high risks between 2-3.5 °C of warming. Risks emerge earlier if global assessments were based on national risk thresholds, underscoring the need for stringent mitigation to limit future extreme heat risks.
| S-EPMC8655081 | biostudies-literature
Project description:Researchers have long been interested in the evolution of culture and the ways in which change in cultural systems can be reconstructed and tracked. Within the realm of language, these questions are increasingly investigated with Bayesian phylogenetic methods. However, such work in cultural phylogenetics could be improved by more explicit quantification of reconstruction and transition probabilities. We apply such methods to numerals in the languages of Australia. As a large phylogeny with almost universal 'low-limit' systems, Australian languages are ideal for investigating numeral change over time. We reconstruct the most likely extent of the system at the root and use that information to explore the ways numerals evolve. We show that these systems do not increment serially, but most commonly vary their upper limits between 3 and 5. While there is evidence for rapid system elaboration beyond the lower limits, languages lose numerals as well as gain them. We investigate the ways larger numerals build on smaller bases, and show that there is a general tendency to both gain and replace 4 by combining 2 + 2 (rather than inventing a new unanalysable word 'four'). We develop a series of methods for quantifying and visualizing the results.
| S-EPMC4614748 | biostudies-literature
Project description:Understanding disease transmission in the workplace is essential for protecting workers. To model disease outbreaks, the small populations in many workplaces require that stochastic effects are considered, which results in higher uncertainty. The aim of this study was to quantify and interpret the uncertainty inherent in such circumstances. We assessed how uncertainty of an outbreak in workplaces depends on i) the infection dynamics in the community, ii) the workforce size, iii) spatial structure in the workplace, iv) heterogeneity in susceptibility of workers, and v) heterogeneity in infectiousness of workers. To address these questions, we developed a multiscale model: A deterministic model to predict community transmission, and a stochastic model to predict workplace transmission. We extended this basic workplace model to allow for spatial structure, and heterogeneity in susceptibility and infectiousness in workers. We found a non-monotonic relationship between the workplace transmission rate and the coefficient of variation (CV), which we use as a measure of uncertainty. Increasing community transmission, workforce size and heterogeneity in susceptibility decreased the CV. Conversely, increasing the level of spatial structure and heterogeneity in infectiousness increased the CV. However, when the model predicts bimodal distributions, for example when community transmission is low and workplace transmission is high, the CV fails to capture this uncertainty. Overall, our work informs modellers and policy makers on how model complexity impacts outbreak uncertainty. In particular: workforce size, community and workplace transmission, spatial structure and individual heterogeneity contribute in a specific and individual manner to the predicted workplace outbreak size distribution.
| S-EPMC10834607 | biostudies-literature
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:In recent years, a number of methods have been proposed to estimate the times at which a neuron spikes on the basis of calcium imaging data. However, quantifying the uncertainty associated with these estimated spikes remains an open problem. We consider a simple and well-studied model for calcium imaging data, which states that calcium decays exponentially in the absence of a spike, and instantaneously increases when a spike occurs. We wish to test the null hypothesis that the neuron did not spike-i.e., that there was no increase in calcium-at a particular timepoint at which a spike was estimated. In this setting, classical hypothesis tests lead to inflated Type I error, because the spike was estimated on the same data used for testing. To overcome this problem, we propose a selective inference approach. We describe an efficient algorithm to compute finite-sample $p$-values that control selective Type I error, and confidence intervals with correct selective coverage, for spikes estimated using a recent proposal from the literature. We apply our proposal in simulation and on calcium imaging data from the $\texttt{spikefinder}$ challenge.
| S-EPMC10449000 | biostudies-literature
Project description:Long-term streamflow datasets inevitably include gaps, which must be filled to allow estimates of runoff and ultimately catchment water budgets. Uncertainty introduced by filling gaps in discharge records is rarely, if ever, reported. We characterized the uncertainty due to streamflow gaps in a reference watershed at the Hubbard Brook Experimental Forest (HBEF) from 1996 to 2009 by simulating artificial gaps of varying duration and flow rate, with the objective of quantifying their contribution to uncertainty in annual streamflow. Gaps were filled using an ensemble of regressions relating discharge from nearby streams, and the predicted flow was compared to the actual flow. Differences between the predicted and actual runoff increased with both gap length and flow rate, averaging 2.8% of the runoff during the gap. At the HBEF, the sum of gaps averaged 22 days per year, with the lowest and highest annual uncertainties due to gaps ranging from 1.5 mm (95% confidence interval surrounding mean runoff) to 21.1 mm. As a percentage of annual runoff, uncertainty due to gap filling ranged from 0.2-2.1%, depending on the year. Uncertainty in annual runoff due to gaps was small at the HBEF, where infilling models are based on multiple similar catchments in close proximity to the catchment of interest. The method demonstrated here can be used to quantify uncertainty due to gaps in any long-term streamflow data set, regardless of the gap-filling model applied.
| S-EPMC7380281 | biostudies-literature
Project description:Short-duration precipitation extremes (PE) increase at a rate of around 7%/K explained by the Clausius-Clapeyron relationship. Previous studies show uncertainty in the extreme precipitation-temperature relationship (scaling) due to various thermodynamic/dynamic factors. Here, we show that uncertainty may arise from the choice of data and methods. Using hourly precipitation (PPT) and daily dewpoint temperature (DPT) across 2,905 locations over the United States, we found higher scaling for quality-controlled data, all locations showing positive (median 6.2%/K) scaling, as compared to raw data showing positive (median 5.3%/K) scaling over 97.5% of locations. We found higher scaling for higher measurement precision of PPT (0.25 mm: median 7.8%/K; 2.54 mm: median 6.6%/K). The method that removes seasonality in PPT and DPT gives higher (with seasonality: median 6.2%/K; without seasonality: median 7.2%/K) scaling. Our results demonstrate the importance of quality-controlled, high-precision observations and robust methods in estimating accurate scaling for a better understanding of PE change with warming.
| S-EPMC9285755 | biostudies-literature