Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bayesian analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Methods for Probabilistic Uncertainty Analysis and Bayesian Analysis with Examples of Statistically Analyzing Data to Revise MMOD Risk Estimates and Compare Models

Probabilistic methods are presented for characterizing and quantifying uncertainties in models and model predictions. Techniques are given for constructing specific uncertainty distributions based on available information. Alternative techniques are given for propagating uncertainties in model inputs to obtain the uncertainty in the model result or prediction. Bayesian techniques are also described for utilizing data and information to update and revise model results and predictions. The focus is on applications with numerous specific examples given. The use of data to revise Micrometeoroid and Orbital Debris (MMOD) risk prediction models are among the examples given.

Risk↗

KMT-2019-BLG-0371 and the Limits of Bayesian Analysis

We show that the perturbation at the peak of the light curve of microlensing event KMT-2019-BLG-0371 is explained by a model with a mass ratio between the host star and planet of q ∼ 0.08. Due to the short event duration (t(sub E) ∼ 6.5 days), the secondary object in this system could potentially be a massive giant planet. A Bayesian analysis shows that the system most likely consists of a host star with a mass M(sub h) = 0.09(+0.14/-0.05) M⨀ and a massive giant planet with a mass = M(sub h) = 7.70(+11.34/-3.90) M(sub Jup). However, the interpretation of the secondary as a planet (i.e., as having M(sub p) < 13M(sub Jup)) rests entirely on the Bayesian analysis. Motivated by this event, we conduct an investigation to determine which constraints meaningfully affect Bayesian analyses for microlensing events. We find that the masses inferred from such a Bayesian analysis are determined almost entirely by the measured value of θ(sub E) and are relatively insensitive to other factors such as the direction of the event (l,b), the lens–source relative proper motion μ(sub rel), or the specific Galactic model prior.

Binary lens microlensing↗

Ockham's razor and Bayesian analysis

'Ockham's razor', the ad hoc principle enjoining the greatest possible simplicity in theoretical explanations, is presently shown to be justifiable as a consequence of Bayesian inference; Bayesian analysis can, moreover, clarify the nature of the 'simplest' hypothesis consistent with the given data. By choosing the prior probabilities of hypotheses, it becomes possible to quantify the scientific judgment that simpler hypotheses are more likely to be correct. Bayesian analysis also shows that a hypothesis with fewer adjustable parameters intrinsically possesses an enhanced posterior probability, due to the clarity of its predictions.

Jefferys, William H.↗

A Bayesian analysis of two probability models describing thunderstorm activity at Cape Kennedy, Florida

A Bayesian analysis of the two discrete probability models, the negative binomial and the modified negative binomial distributions, which have been used to describe thunderstorm activity at Cape Kennedy, Florida, is presented. The Bayesian approach with beta prior distributions is compared to the classical approach which uses a moment method of estimation or a maximum-likelihood method. The accuracy and simplicity of the Bayesian method is demonstrated.

Williford, W. O.↗

Global Bayesian analysis of J/𝜓 photoproduction on proton and lead targets

We perform a global Bayesian analysis of diffractive J/𝜓 production in 𝛾 + 𝑝 and 𝛾 + Pb collisions using a color glass condensate (CGC) based calculation framework. As past calculations have shown that CGC-based models typically overpredict the J/𝜓 production in 𝛾 + Pb collisions at high center of mass energy, we address the question of whether it is possible to describe coherent and incoherent diffractive J/𝜓 data from 𝛾 + 𝑝 collisions at Hadron–Electron Ring Accelerator (HERA) and the LHC, and from 𝛾 + Pb collisions at the LHC simultaneously. Our results indicate that a simultaneous description of 𝛾 + 𝑝 and 𝛾 + Pb data is challenging, with results improving when an overall 𝐾-factor—scaling 𝛾 + 𝑝 and 𝛾 + Pb cross sections to absorb model uncertainties—is introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Bayesian Analysis of SDSS J0914+0853, a Low-mass Dual AGN Candidate

We present the first results from Bayesian AnalYsis of Multiple AGN in X-rays (BAYMAX), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual active galactic nuclei (AGNs) is to use X-ray observations, only sources that are widely separated relative to the instrumentʼs point-spread function are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandraʼs angular resolution (<1″). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ∼ 0"3 from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find a very low probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0"3. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, as well as studying the dual AGN population where past spatial resolution limits have prevented systematic analyses.

Active galaxies↗

Bayesian analysis of (3 +1)⁢D relativistic nuclear dynamics with the RHIC beam energy scan data

This work presents a Bayesian inference study for relativistic heavy-ion collisions in the beam energy scan program at the BNL Relativistic Heavy-Ion Collider. The theoretical model simulates event-by-event (3+1)-dimensional [(3+1)⁢D] collision dynamics using hydrodynamics and hadronic transport theory. We analyze the model's 20-dimensional posterior distributions obtained using three model emulators with different accuracy and demonstrate the essential role of training an accurate model emulator in the Bayesian analysis. Our analysis provides robust constraints on the quark-gluon plasma's transport properties and various aspects of (3+1)⁢D relativistic nuclear dynamics. By running full model simulations with 100 parameter sets sampled from the posterior distribution, we make predictions for p T -differential observables and estimate their systematic theory uncertainty. Here, a sensitivity analysis is performed to elucidate how individual experimental observables respond to different model parameters, providing useful physics insights into the phenomenological model for heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian analysis of the 86 Sr ⁢(𝛼,𝛼) reaction to constrain the 86 Sr ⁢(𝛼,𝑛) cross section at astrophysical energies

The alpha optical model potential (𝛼-OMP ) is a phenomenological approach used to describe elastic scattering where multiple reaction channels are open. It is one of the most critical inputs for the calculation of thermonuclear reaction rates in explosive stellar environments, but uncertainties within the 𝛼-OMP lead to imprecise predictions hindering comparisons between calculations and observations. In order to improve the precision of the 𝛼-OMP, additional nuclear physics data are required. In this paper, a measurement of the 86 Sr (𝛼, 𝛼) elastic scattering cross section at multiple energies is reported. Here, a local optical potential is constructed via a fully Bayesian analysis of the elastic scattering data. The resulting uncertainties on the low-energy cross sections relevant to nuclear astrophysics are then calculated and shown to be on the order of 50%.

59 ≤ A ≤ 89↗

Bayesian Analysis of the Detection Performance of the Lightning Imaging Sensors

Identical Lightning Imaging Sensors aboard the Tropical Rainfall Measuring Mission satellite (TRMM LIS, 1998-2015) and International Space Station (ISS LIS, 2017-present) have collectively provided over two decades of lightning observations over the global tropics, with ISS LIS extending coverage into the mid-latitudes. Quantifying the detection performance of both LIS sensors is a necessary step toward generating a LIS climatological record and accurately combining LIS data with lightning detections from other sensors and networks. We compare lightning observations from both LIS sensors with reference sources including the Geostationary Lightning Mapper (GLM) and ground-based networks operated by Earth Networks (the Earth Networks Total Lightning Network [ENTLN] and Earth Networks Global Lightning Network [ENGLN]) and Vaisala (the National Lightning Detection Network [NLDN] and Global Lightning Dataset [GLD360]). Instead of a relative detection efficiency (RDE) approach that involves assuming perfect performance of the reference sensor, we employ a Bayesian approach to estimate the upper limit of the absolute detection efficiency (ADE) of each sensor being analyzed. The results of this Bayesian analysis illustrate the geographical pattern of ADE as well as its diurnal cycle and yearly evolution, reflecting the growth of the reference networks over time.

Katrina Virts↗

Bayesian Analysis for Risk Assessment of Selected Medical Events in Support of the Integrated Medical Model Effort

The Exploration Medical Capability project is creating a catalog of risk assessments using the Integrated Medical Model (IMM). The IMM is a software-based system intended to assist mission planners in preparing for spaceflight missions by helping them to make informed decisions about medical preparations and supplies needed for combating and treating various medical events using Probabilistic Risk Assessment. The objective is to use statistical analyses to inform the IMM decision tool with estimated probabilities of medical events occurring during an exploration mission. Because data regarding astronaut health are limited, Bayesian statistical analysis is used. Bayesian inference combines prior knowledge, such as data from the general U.S. population, the U.S. Submarine Force, or the analog astronaut population located at the NASA Johnson Space Center, with observed data for the medical condition of interest. The posterior results reflect the best evidence for specific medical events occurring in flight. Bayes theorem provides a formal mechanism for combining available observed data with data from similar studies to support the quantification process. The IMM team performed Bayesian updates on the following medical events: angina, appendicitis, atrial fibrillation, atrial flutter, dental abscess, dental caries, dental periodontal disease, gallstone disease, herpes zoster, renal stones, seizure, and stroke.

Gilkey, Kelly M.↗

Bayesian analysis of nucleon-nucleon scattering data in pionless effective field theory

We perform Bayesian model calibration of two-nucleon (NN) low-energy constants (LECs) appearing in an NN interaction based on pionless effective field theory (πEFT). The calibration is carried out for potentials constructed using naive dimensional analysis in NN relative momenta (p) up to next-to-leading order [NLO, O(p 2 )] and next-to-next-to-next-to-leading order [N3LO, O(p 4 )]. We consider two classes of pionless πEFT potential: one that acts in all partial waves and another that is dominated by s-wave physics. The two classes produce broadly similar results for calibrations to NN data up to E lab = 5 MeV. Our analysis accounts for the correlated uncertainties that arise from the truncation of the pionless πEFT. We simultaneously estimate both the πEFT LECs and the parameters that quantify the truncation error. This permits the first quantitative estimates of the pionless πEFT breakdown scale, Λ b : the 95% intervals are Λ b ∈[50.11,63.03] MeV at NLO and Λ b ∈[72.27,88.54] MeV at N3LO. Furthermore, invoking naive dimensional analysis for the NN potential, therefore, does not lead to consistent results across orders in pionless πEFT. This exemplifies the possible use of Bayesian tools to identify inconsistencies in a proposed EFT power counting.

Bayesian methods↗

Bayesian Analysis of the Cosmic Microwave Background

There is a wealth of cosmological information encoded in the spatial power spectrum of temperature anisotropies of the cosmic microwave background! Experiments designed to map the microwave sky are returning a flood of data (time streams of instrument response as a beam is swept over the sky) at several different frequencies (from 30 to 900 GHz), all with different resolutions and noise properties. The resulting analysis challenge is to estimate, and quantify our uncertainty in, the spatial power spectrum of the cosmic microwave background given the complexities of "missing data", foreground emission, and complicated instrumental noise. Bayesian formulation of this problem allows consistent treatment of many complexities including complicated instrumental noise and foregrounds, and can be numerically implemented with Gibbs sampling. Gibbs sampling has now been validated as an efficient, statistically exact, and practically useful method for low-resolution (as demonstrated on WMAP 1 and 3 year temperature and polarization data). Continuing development for Planck - the goal is to exploit the unique capabilities of Gibbs sampling to directly propagate uncertainties in both foreground and instrument models to total uncertainty in cosmological parameters.

methods - statistical↗

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Composite Bond Line Measurements Based on A Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗

Composite Bond Line Measurements Based on a Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗