Results for “Bayesian sampling” — Page 4

NASA NTRSDate not supplied

Application of a Bayesian Framework for Plasticity Model Selection

Interpretable Machine Learning (IML) has performed well when tasked with deriving constitutive material models. However, IML has been shown to prefer models that overfit noise in data, which tends to lead to bloat and a decrease in interpretability. Due to these issues, the ability of IML to reliably derive models that fit the data and are both interpretable and generalizable is limited. A method developed recently has shown promise to improve upon traditional IML by using a Bayesian fitness definition for the evolution of free-form models with non-deterministic parameters. This framework was developed for genetic-programming-based symbolic regression(GPSR) and involves model parameter estimation using Sequential Monte Carlo sampling (SMC).The method has demonstrated a reduction in bloat when dealing with noisy data in comparison to conventional GPSR. The results of this framework applied to stress-strain data for copper show models that more effectively predict the experimental data better than was previously shown with GPSR.

plasticity↗

NASA NTRS2004

Perceptual learning through optimization of attentional weighting: human versus optimal Bayesian learner

Human performance in visual detection, discrimination, identification, and search tasks typically improves with practice. Psychophysical studies suggest that perceptual learning is mediated by an enhancement in the coding of the signal, and physiological studies suggest that it might be related to the plasticity in the weighting or selection of sensory units coding task relevant information (learning through attention optimization). We propose an experimental paradigm (optimal perceptual learning paradigm) to systematically study the dynamics of perceptual learning in humans by allowing comparisons to that of an optimal Bayesian algorithm and a number of suboptimal learning models. We measured improvement in human localization (eight-alternative forced-choice with feedback) performance of a target randomly sampled from four elongated Gaussian targets with different orientations and polarities and kept as a target for a block of four trials. The results suggest that the human perceptual learning can occur within a lapse of four trials (<1 min) but that human learning is slower and incomplete with respect to the optimal algorithm (23.3% reduction in human efficiency from the 1st-to-4th learning trials). The greatest improvement in human performance, occurring from the 1st-to-2nd learning trial, was also present in the optimal observer, and, thus reflects a property inherent to the visual task and not a property particular to the human perceptual learning mechanism. One notable source of human inefficiency is that, unlike the ideal observer, human learning relies more heavily on previous decisions than on the provided feedback, resulting in no human learning on trials following a previous incorrect localization decision. Finally, the proposed theory and paradigm provide a flexible framework for future studies to evaluate the optimality of human learning of other visual cues and/or sensory modalities.

Non-NASA Center↗

NASA NTRS2017

OGLE-2013-BLG-1761Lb: A Massive Planet around an MK Dwarf

We report the discovery and the analysis of the planetary microlensing event, OGLE-2013-BLG-1761. There are some degenerate solutions in this event because the planetary anomaly is only sparsely sampled. However, the detailed light curve analysis ruled out all stellar binary models and shows the lens to be a planetary system. There is the so-called close wide degeneracy in the solutions with the planet host mass ratio of q approx.(7.0+/-2.0) x 10(exp -3) and q approx.(8.1+/-2.6) x 10(exp -3) with the projected separation in Einstein radius units of s = 0.95 (close) and s = 1.18(wide), respectively. The microlens parallax effect is not detected, but the finite source effect is detected. Our Bayesian analysis indicates that the lens system is located -D(sub L) = 6.9(+ 1.0 -1.2)kpc away from us and the host star is an M/K dwarf with amass of M(sub L) = 0.33(+ 0.32- 1.9)Stellar Mass orbited by a super-Jupiter mass planet with a mass of m(sub p) = 2.7(+ 2.5 - 1.5) M(sub Jup) at the projected separation of a(sub l) = 1.8(+ 0.5 -0.5)au. The preference of the large lens distance in the Bayesian analysis is due to the relatively large observed source star radius. The distance and other physical parameters may be constrained by the future high-resolution imaging by large ground telescopes or HST. If the estimated lens distance is correct, then this planet provides another sample for testing the claimed deficit of planets in the Galactic bulge.

Hirao, Y.↗

NASA NTRSDate not supplied

Consequences of Asteroid Characterization on the State of Knowledge about Inferred Physical Properties and Impact Risk

Physical characteristics of Near-Earth Objects (NEOs) are essential inputs to planetary defense assessments. The size, density, and strength of an NEO are critical inputs to modeling behavior during atmospheric entry as well as assessing the risk of impact. Similarly, knowledge of the physical characteristics of an object are necessary to evaluate the probable result of a mitigation mission. Usually, these attributes cannot be directly measured, but increasingly sophisticated methods have been developed to infer physical properties from related measurements of asteroids, meteors, and/or meteorites. Fortuitously, some of these measurements have been obtained for enough NEOs to elucidate the distribution of values across the sampled population. However, the situation becomes more challenging when considering a specific asteroid, since it is unlikely that all the relevant measurements have been made for any given object. We have developed a Bayesian network that can combine available information about a particular NEO with knowledge of the larger population to infer probabilistic values and uncertainties for physical characteristics of interest. Distributions of asteroid population albedos, taxonomic classes, and macroporosities, along with meteorite density distributions and associations between taxonomic classes and meteorite classes, provide the default distributions for the network’s parameter nodes. The inference network links parameters for each virtual asteroid either deterministically or probabilistically as appropriate, and eliminates any unphysical combinations of parameters. Within the context of planetary defense, our Bayesian network can be used to constrain the ranges of likely impactor properties, which can subsequently reduce the uncertainty in modelling of atmospheric entry, mitigation efficacy, and impact risk assessment. When additional measurements become available for a specific object, the network incorporates those measurements to generate virtual asteroids with property distributions that are consistent with the measurements. We will use the 2023 PDC scenario to demonstrate how the inference network can be combined with plausible characterization measurements to refine the state of knowledge about likely combinations of physical parameters and the resulting impact risk.

risk assessment↗

NASA NTRSDate not supplied

Consequences of Asteroid Characterization on the State of Knowledge about Inferred Physical Properties and Impact Risk

Physical characteristics of Near-Earth Objects (NEOs) are essential inputs to planetary defense assessments. The size, density, and strength of an NEO are critical inputs to modeling behavior during atmospheric entry as well as assessing the risk of impact. Similarly, knowledge of the physical characteristics of an object are necessary to evaluate the probable result of a mitigation mission. Usually, these attributes cannot be directly measured, but increasingly sophisticated methods have been developed to infer physical properties from related measurements of asteroids, meteors, and/or meteorites. Fortuitously, some of these measurements have been obtained for enough NEOs to elucidate the distribution of values across the sampled population. However, the situation becomes more challenging when considering a specific asteroid, since it is unlikely that all the relevant measurements have been made for any given object. We have developed a Bayesian network that can combine available information about a particular NEO with knowledge of the larger population to infer probabilistic values and uncertainties for physical characteristics of interest. Distributions of asteroid population albedos, taxonomic classes, and macroporosities, along with meteorite density distributions and associations between taxonomic classes and meteorite classes, provide the default distributions for the network’s parameter nodes. The inference network links parameters for each virtual asteroid either deterministically or probabilistically as appropriate, and eliminates any unphysical combinations of parameters. Within the context of planetary defense, our Bayesian network can be used to constrain the ranges of likely impactor properties, which can subsequently reduce the uncertainty in modelling of atmospheric entry, mitigation efficacy, and impact risk assessment. When additional measurements become available for a specific object, the network incorporates those measurements to generate virtual asteroids with property distributions that are consistent with the measurements. We will use the 2023 PDC scenario to demonstrate how the inference network can be combined with plausible characterization measurements to refine the state of knowledge about likely combinations of physical parameters and the resulting impact risk.

risk assessment↗

NASA NTRS2012

Metallicities of Emission-Line Galaxies from HST ACS PEARS and HST WFC3 ERS Grism Spectroscopy at 0.6 is less than z is less than 2.4

Galaxies selected on the basis of their emission line strength. show low metallicities, regardless of their redshifts. We conclude this from a sample of faint galaxies at redshifts between 0.6 < z < 2.4, selected by their prominent emission lines in low resolution grism spectra in the optiCa.i with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) and in the near-infrared using Wide-Field Camera 3 (WFC3). Using a sample of 11 emission line galaxies (ELGs) at 0.6 < z < 2.4 with luminosities of -22 approx < MB approx -19 which have [OII], H-Beta, and [OIII] line flux measurements from the combination of two grism spectral surveys, we use the R23 method to derive the gas-phase oxygen abundances: 7.5 <12+log(0/H)<8.5. The galaxy stellar masses are derived using Bayesian based Markov Chain Monte Carlo (pi MC(exp 2)) fitting of their Spectral Energy Distribution (SED), and span the mass range 8.1 < log(M(stellar)/M(solar)) < 10.1. These galaxies show a mass-metal1icity (M-L) and Luminosity-Metallicity (LZ) relation, which is offset by -<J.6 dex in meta1licity at given absolute magnitude and stellar mass relative to the local SDSS galaxies, as well as continuum selected DEEP2 samples at similar redshifts. The emission-line selected galaxies most resemble the local "green peas" galaxies and Lyman-alpha galaxies at z approx = 0.3 and z approx = 2.3 in the M-Z and L-Z relations and their morphologies. The G - M(sub 20) morphology analysis shows that 10 out of 11 show disturbed morphology, even as the star-forming regions are compact. These galaxies may be intrinsically metal poor, being at early stages of formation, or the low metallicities may be due to gas infall and accretion due to mergers.

Xia, Lifang↗

NASA NTRS2012

A Bayesian Framework for Reliability Analysis of Spacecraft Deployments

Deployable subsystems are essential to mission success of most spacecraft. These subsystems enable critical functions including power, communications and thermal control. The loss of any of these functions will generally result in loss of the mission. These subsystems and their components often consist of unique designs and applications for which various standardized data sources are not applicable for estimating reliability and for assessing risks. In this study, a two stage sequential Bayesian framework for reliability estimation of spacecraft deployment was developed for this purpose. This process was then applied to the James Webb Space Telescope (JWST) Sunshield subsystem, a unique design intended for thermal control of the Optical Telescope Element. Initially, detailed studies of NASA deployment history, "heritage information", were conducted, extending over 45 years of spacecraft launches. This information was then coupled to a non-informative prior and a binomial likelihood function to create a posterior distribution for deployments of various subsystems uSing Monte Carlo Markov Chain sampling. Select distributions were then coupled to a subsequent analysis, using test data and anomaly occurrences on successive ground test deployments of scale model test articles of JWST hardware, to update the NASA heritage data. This allowed for a realistic prediction for the reliability of the complex Sunshield deployment, with credibility limits, within this two stage Bayesian framework.

Evans, John W.↗

NASA NTRS2014

The New NASA Orbital Debris Engineering Model ORDEM 3.0

The NASA Orbital Debris Program Office (ODPO) has released its latest Orbital Debris Engineering Model, ORDEM 3.0. It supersedes ORDEM 2000, now referred to as ORDEM 2.0. This newer model encompasses the Earth satellite and debris flux environment from altitudes of low Earth orbit (LEO) through geosynchronous orbit (GEO). Debris sizes of 10 micron through larger than 1 m in non-GEO and 10 cm through larger than 1 m in GEO are available. The inclusive years are 2010 through 2035. The ORDEM model series has always been data driven. ORDEM 3.0 has the benefit of many more hours of data from existing sources and from new sources than past ORDEM versions. The object data range in size from 10 μm to larger than 1 m, and include in situ and remote measurements. The in situ data reveals material characteristics of small particles. Mass densities are grouped in ORDEM 3.0 in terms of 'high-density', represented by 7.9 g/cc, 'medium-density' represented by 2.8 g/cc and 'low-density' represented by 1.4 g/cc. Supporting models have also advanced significantly. The LEO-to-GEO ENvironment Debris model (LEGEND) includes an historical and a future projection component with yearly populations that include launched and maneuvered intact spacecraft and rocket bodies, mission related debris, and explosion and collision event fragments. LEGEND propagates objects with ephemerides and physical characteristics down to 1 mm in size. The full LEGEND yearly population acts as an a priori condition for a Bayesian statistical model. Specific populations are added from sodium potassium droplet releases, recent major accidental and deliberate collisions, and known anomalous debris events. This paper elaborates on the upgrades of this model over previous versions. Sample validation results with remote and in situ measurements are shown, and the consequences of including material density are discussed as it relates to heightened risks to crewed and robotic spacecraft

Krisko, P. H.↗

NASA NTRS1975

A generalized concept for cost-effective structural design

A generalized concept for cost-effective structural design is introduced. It is assumed that decisions affecting the cost effectiveness of aerospace structures fall into three basic categories: design, verification, and operation. Within these basic categories, certain decisions concerning items such as design configuration, safety factors, testing methods, and operational constraints are to be made. All or some of the variables affecting these decisions may be treated probabilistically. Bayesian statistical decision theory is used as the tool for determining the cost optimum decisions. A special case of the general problem is derived herein, and some very useful parametric curves are developed and applied to several sample structures.

Thomas, J. M.↗

NASA NTRSDate not supplied

Effects of Dose Error and Sample Size on Sonic Boom Dose-Response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

NASA NTRS1980

Evaluation of Bayesian Sequential Proportion Estimation Using Analyst Labels

The author has identified the following significant results. A total of ten Large Area Crop Inventory Experiment Phase 3 blind sites and analyst-interpreter labels were used in a study to compare proportional estimates obtained by the Bayes sequential procedure with estimates obtained from simple random sampling and from Procedure 1. The analyst error rate using the Bayes technique was shown to be no greater than that for the simple random sampling. Also, the segment proportion estimates produced using this technique had smaller bias and mean squared errors than the estimates produced using either simple random sampling or Procedure 1.

Lennington, R. K.↗

NASA NTRS2020

Bayesian Statistical Models for Community Annoyance Survey Data

This paper demonstrates the use of two Bayesian statistical models to analyze single-event sonic boom exposure and human annoyance data from community response surveys. Each model is fit to data from a NASA pilot study.Unlike many community noise surveys, this study used a panel sample to collect multiple observations per participant instead of a single observation. Thus, a multilevel (also known as hierarchical or mixed-effects) model is used to account for the within-subject correlation in the panel sample data. This paper describes a multilevel logistic regression model and a multilevel ordinal regression model. The paper also proposes a method for calculating a summary dose-response curve from the multilevel models that represents the population. The two models’ summary dose-response curves are visually similar. However, their estimates differ when calculating the noise dose at a fixed percent highly annoyed.

Musical instruments↗

NASA NTRS2021

New giant planet beyond the snow line for an extended MOA exoplanet microlens sample

Characterizing a planet detected by microlensing is hard if the planetary signal is weak or the lens-source relative trajectory is far from caustics. However, statistical analyses of planet demography must include those planets to accurately determine occurrence rates. As part of a systematic modelling effort in the context of a >10-yr retrospective analysis of MOA’s survey observations to build an extended MOA statistical sample, we analyse the light curve of the planetary microlensing event MOA-2014-BLG-472. This event provides weak constraints on the physical parameters of the lens, as a result of a planetary anomaly occurring at low magnification in the light curve. We use a Bayesian analysis to estimate the properties of the planet, based on a refined Galactic model and the assumption that all Milky Way’s stars have an equal planet-hosting probability. We find that a lens consisting of a 1.9(+2.2,−1.2)M(J) giant planet orbiting a 0.31(+0.36,−0.19)Mꙩ host at a projected separation of 0.75±0.24au is consistent with the observations and is most likely, based on the Galactic priors. The lens most probably lies in the Galactic bulge, at 7.2(+0.6,−1.7)kpc from Earth. The accurate measurement of the measured planet-to-host star mass ratio will be included in the next statistical analysis of cold planet demography detected by microlensing.

Clément Ranc↗

NASA NTRS2014

A Practical Philosophy of Complex Climate Modelling

We give an overview of the practice of developing and using complex climate models, as seen from experiences in a major climate modelling center and through participation in the Coupled Model Intercomparison Project (CMIP).We discuss the construction and calibration of models; their evaluation, especially through use of out-of-sample tests; and their exploitation in multi-model ensembles to identify biases and make predictions. We stress that adequacy or utility of climate models is best assessed via their skill against more naive predictions. The framework we use for making inferences about reality using simulations is naturally Bayesian (in an informal sense), and has many points of contact with more familiar examples of scientific epistemology. While the use of complex simulations in science is a development that changes much in how science is done in practice, we argue that the concepts being applied fit very much into traditional practices of the scientific method, albeit those more often associated with laboratory work.

complex simulation↗

NASA NTRS2012

The Rise and Fall of Star Formation Histories of Blue Galaxies at Redshifts 0.2 < z < 1.4

Popular cosmological scenarios predict that galaxies form hierarchically from the merger of many progenitor, each with their own unique star formation history (SFH). We use the approach recently developed by Pacifici et al. to constrain the SFHs of 4517 blue (presumably star-forming) galaxies with spectroscopic redshifts in the range O.2 < z < 1:4 from the All-Wavelength Extended Groth Strip International Survey (AEGIS). This consists in the Bayesian analysis of the observed galaxy spectral ' energy distributions with a comprehensive library of synthetic spectra assembled using state-of-the-art models of star formation and chemical enrichment histories, stellar population synthesis, nebular emission and attenuation by dust. We constrain the SFH of each galaxy in our sample by comparing the observed fluxes in the B, R,l and K(sub s) bands and rest-frame optical emission-line luminosities with those of one million model spectral energy distributions. We explore the dependence of the resulting SFH on galaxy stellar mass and redshift. We find that the average SFHs of high-mass galaxies rise and fall in a roughly symmetric bell-shaped manner, while those of low-mass galaxies rise progressively in time, consistent with the typically stronger activity of star formation in low-mass compared to high-mass galaxies. For galaxies of all masses, the star formation activity rises more rapidly at high than at low redshift. These findings imply that the standard approximation of exponentially declining SFHs wIdely used to interpret observed galaxy spectral energy distributions is not appropriate to constrain the physical parameters of star-forming galaxies at intermediate redshifts.

Pacifici, Camilla↗

NASA NTRSDate not supplied

Effects of Dose Error and Sample Size on Sonic Boom Dose-response Curves

NASA will soon be collecting noise-annoyance community survey data as the X-59 aircraft flies supersonically over several communities in the USA. Sparse measurements of the X-59 sonic thumps will be used together with physics-based simulations to estimate noise doses at survey participant locations. These dose estimates have associated error that affects the accuracy of modeled dose-response curves, which can result in misestimation of annoyance. The precision in dose-response curves is also a consideration in selecting the number of survey participants. To enable pretest studies of dose error and precision, simulated dose-response data were generated based on NASA’s Quiet Supersonic Flights 2018 test. The data included various degrees of dose error and sample size. Frequentist multilevel logistic regression models were fit to the true and perturbed dose-response data. Simple proportional relationships were identified between the model parameters and the perturbation standard deviation. The summary dose-response curves illustrate the impact on accuracy if dose error is not accounted for in the model. The precision in the dose-response curves is also shown as the number of participants and degree of participation is varied. Finally, sampling variability is illustrated by showing the dose-response curves for several replicates with random draws of participants and errors.

X-59↗

NASA NTRS2012

Accelerated Aging in Electrolytic Capacitors for Prognostics

The focus of this work is the analysis of different degradation phenomena based on thermal overstress and electrical overstress accelerated aging systems and the use of accelerated aging techniques for prognostics algorithm development. Results on thermal overstress and electrical overstress experiments are presented. In addition, preliminary results toward the development of physics-based degradation models are presented focusing on the electrolyte evaporation failure mechanism. An empirical degradation model based on percentage capacitance loss under electrical overstress is presented and used in: (i) a Bayesian-based implementation of model-based prognostics using a discrete Kalman filter for health state estimation, and (ii) a dynamic system representation of the degradation model for forecasting and remaining useful life (RUL) estimation. A leave-one-out validation methodology is used to assess the validity of the methodology under the small sample size constrain. The results observed on the RUL estimation are consistent through the validation tests comparing relative accuracy and prediction error. It has been observed that the inaccuracy of the model to represent the change in degradation behavior observed at the end of the test data is consistent throughout the validation tests, indicating the need of a more detailed degradation model or the use of an algorithm that could estimate model parameters on-line. Based on the observed degradation process under different stress intensity with rest periods, the need for more sophisticated degradation models is further supported. The current degradation model does not represent the capacitance recovery over rest periods following an accelerated aging stress period.

Celaya, Jose R.↗

NASA NTRS2016

Predictable Patterns in Planetary Transit Timing Variations and Transit Duration Variations Due to Exomoons

We present new ways to identify single and multiple moons around extrasolar planets using planetary transit timing variations (TTVs) and transit duration variations (TDVs). For planets with one moon, measurements from successive transits exhibit a hitherto undescribed pattern in the TTV-TDV diagram, originating from the stroboscopic sampling of the planet's orbit around the planet-moon barycenter. This pattern is fully determined and analytically predictable after three consecutive transits. The more measurements become available, the more the TTV-TDV diagram approaches an ellipse. For planets with multiple moons in orbital mean motion resonance (MMR), like the Galilean moon system, the pattern is much more complex and addressed numerically in this report. Exomoons in MMR can also form closed, predictable TTV-TDV figures, as long as the drift of the moons' pericenters is suciently slow.We find that MMR exomoons produce loops in the TTV-TDV diagram and that the number of these loops is equal to the order of the MMR, or the largest integer in the MMR ratio.We use a Bayesian model and Monte Carlo simulations to test the discoverability of exomoons using TTV-TDV diagrams with current and near-future technology. In a blind test, two of us (BP, DA) successfully retrieved a large moon from simulated TTV-TDV by co-authors MH and RH, which resembled data from a known Kepler planet candidate. Single exomoons with a 10 percent moon-to-planet mass ratio, like to Pluto-Charon binary, can be detectable in the archival data of the Kepler primary mission. Multi-exomoon systems, however, require either larger telescopes or brighter target stars. Complementary detection methods invoking a moon's own photometric transit or its orbital sampling effect can be used for validation or falsification. A combination of TESS, CHEOPS, and PLATO data would offer a compelling opportunity for an exomoon discovery around a bright star.

stroboscopic sampling↗