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: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: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:BackgroundThe apnea-hypopnea index (AHI) (or one of its derivatives) is the primary clinical metric for characterizing sleep disordered breathing-the value of which with respect to a threshold determines severity of diagnosis and eligibility for treatment reimbursement. The index value, however, is taken as a perfect point estimate, with no measure of statistical uncertainty. Thus, current practice does not robustly account for variability in diagnosis/eligibility due to chance. In this paper, we quantify the statistical uncertainty associated with respiratory event indices for sleep disordered breathing and the effect of uncertainty on treatment eligibility.MethodsWe develop an empirical estimate of uncertainty using a non-parametric bootstrap on the interevent times, as well as a theoretical Poisson estimate reflecting the current formulation of the AHI. We then apply these methods to estimate AHI uncertainty for 2049 subjects (954/1095 M/F, age: mean 69 ± 9.1) from the Multi-Ethnic Study of Atherosclerosis (MESA).Results and conclusionsThe mean 95% empirical confidence interval width was 11.500 ± 6.208 events per hour and the mean 95% theoretical Poisson confidence interval width was 5.998 ± 2.897 events per hour, suggesting that uncertainty is likely a major confounding factor within the current diagnostic framework. Of the 278 subjects in the symptomatic population (ESS>10), 27% (76/278) had uncertain diagnoses given the 95% empirical confidence interval. Of the 2049 subjects in the full population, 43% (880/2049) had uncertain diagnoses given the 95% empirical confidence interval. The inclusion of subjects with uncertain diagnoses increases the number of eligible patients by 21.3% for the symptomatic population and by 84.8% for the full population. The exclusion of subjects with uncertain diagnoses given the 95% empirical confidence interval decreases the number of eligible patients by 12.4% for the symptomatic population and by 34.8% for full population. Additional analyses suggest that it is practically infeasible to gain diagnostic statistical significance through additional testing for a broad range of borderline cases. Overall, these results suggest that AHI uncertainty is a vital additional piece of information that would greatly benefit clinical practice, and that the inclusion of uncertainty in epidemiological analysis might help improve the ability for researchers to robustly link AHI with co-morbidities and long-term outcomes.
| S-EPMC6938549 | 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: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
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