Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “likelihood function”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Joint inference of multiplicative and additive systematics in galaxy density fluctuations and clustering measurements

Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating both additive and multiplicative systematics in galaxy clustering measurements (two-point function) by joint inference of contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with Kilo-Degree Survey-like mock galaxy catalogues and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-percent level and reduces numerous additive and multiplicative systematics from 1.5σ to less than 0.1σ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

New machine learning techniques for simulation-based inference: InferoStatic nets, kernel score estimation, and kernel likelihood ratio estimation

We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential \varphi φ . In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

Kong, Kyoungchul↗
DOE OSTI2023

Codebase release 0.1 for infstat

We propose an intuitive, machine-learning approach to multiparameter inference, dubbed the InferoStatic Networks (ISN) method, to model the score and likelihood ratio estimators in cases when the probability density can be sampled but not computed directly. The ISN uses a backend neural network that models a scalar function called the inferostatic potential \varphi φ . In addition, we introduce new strategies, respectively called Kernel Score Estimation (KSE) and Kernel Likelihood Ratio Estimation (KLRE), to learn the score and the likelihood ratio functions from simulated data. We illustrate the new techniques with some toy examples and compare to existing approaches in the literature. We mention en passant some new loss functions that optimally incorporate latent information from simulations into the training procedure.

Kong, Kyoungchul↗
DOE OSTI2021

Modelling multiple occurrences of activities during a day: An extension of the MDCEV model

The increased interest in time use among transport researchers has led to a search for flexible but tractable models of time use, such as Bhat's Multiple Discrete Continuous Extreme Value (MDCEV) model. MDCEV formulations typically model aggregate time allocation into different activity types during a given period, such as the amount of time spent working and shopping in a day. While these applications provide valuable insights into activity participation, they ignore disaggregate activity-episodes, that is the fact that people might split their total time spent working in multiple separate blocks, with breaks or other activities in between. Insights into this splitting into episodes are necessary for predicting trips and understanding time use satiation. We propose a modified MDCEV model where an activity-episode, rather than an activity type, is the basic choice alternative, using a modified utility function to capture the reduced likelihood of individuals performing a very large number of episodes of the same activity. Results from two large revealed preference datasets exhibit equivalent forecast accuracy between the traditional and proposed approach at an aggregate level, but the latter also provides insights on the number and duration of activity-episodes with significant accuracy.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1968

Concept for simplified serial digital decoder

Modular decoder, which lends itself best to special purpose digital equipment using sequential access memories, decodes the first order Reed-Muller codes. It functions as a maximum-likelihood exhaustive-search decoder and is a modular implementation to accommodate codes of any length.

Green, R. R.↗
NASA NTRS1976

Nonparametric probability density estimation by optimization theoretic techniques

Two nonparametric probability density estimators are considered. The first is the kernel estimator. The problem of choosing the kernel scaling factor based solely on a random sample is addressed. An interactive mode is discussed and an algorithm proposed to choose the scaling factor automatically. The second nonparametric probability estimate uses penalty function techniques with the maximum likelihood criterion. A discrete maximum penalized likelihood estimator is proposed and is shown to be consistent in the mean square error. A numerical implementation technique for the discrete solution is discussed and examples displayed. An extensive simulation study compares the integrated mean square error of the discrete and kernel estimators. The robustness of the discrete estimator is demonstrated graphically.

Scott, D. W.↗
NASA NTRS1981

A complete X-ray sample of the high latitude sky from HEAO-1 A-2: log N lo S and luminosity functions

An experiment was performed in which a complete X-ray survey of the 8.2 steradians of the sky at galactic latitudes where the absolute value of b is 20 deg down to a limiting sensitivity of 3.1 x ten to the minus 11th power ergs/sq cm sec in the 2-10 keV band. Of the 85 detected sources 17 were identified with galactic objects, 61 were identified with extragalactic objects, and 7 remain unidentified. The log N - log S relation for the non-galactic objects is well fit by the Euclidean relationship. The X-ray spectra of these objects were used to construct log N - log S in physical units. The complete sample of identified sources was used to construct X-ray luminosity functions, using the absolute maximum likelihood method, for clusters galaxies and active galactic nuclei.

Piccinotti, G.↗
NASA NTRS1985

Einstein X-ray survey of the Pleiades - The dependence of X-ray emission on stellar age

The data obtained with two pointed observations of 1 deg by 1 deg fields of the Pleiades region have been analyzed, and the results are presented. The maximum-likelihood X-ray luminosity functions for the Pleiades G and K stars in the cluster are derived, and it is shown that, for the G stars, the Pleiades X-ray luminosity function is significantly brighter than the corresponding function for Hyades G dwarf stars. This finding indicates a dependence of X-ray luminosity on stellar age, which is confirmed by comparison of the same data with median X-ray luminosities of pre-main sequence and local disk population dwarf G stars. It is suggested that the significantly larger number of bright X-ray sources associated with G stars than with K stars, the lack of detection of M stars, and the relatively rapid rotation of the Pleiades K stars can be explained in terms of the onset of internal differential rotation near the convective envelope-radidative core interface after the spin-up phase during evolution to the main sequence.

Micela, G.↗
NASA NTRS1992

X-ray emission at low-mass end of the MS - Results from an extensive Einstein Observatory survey

The 0.16-4.0 keV X-ray emission from K and M stars of luminosity classes IV, V, and VI within 25 parsec of the sun was measured using all available Einstein Observatory IPC data and a critical compilation of catalogued optical data. Fifty four of 88 stars were detected, 70 of 138 M stars with Mv less than 13.4 (corresponding to M6) and 15 or 31 fainter M stars. The surveyed stars were grouped, on the basis of U, V, W space velocity components, into old-disk, young-disk, and halo population stars. Then, a subsample was selected which is statistically representative of the population of K and M stars in the solar neighborhood, on the basis of which unbiased Maximum Likelihood X-ray luminosity functions were constructed for K, early M, and late M stars. The investigation revealed a decrease of X-ray luminosity with increasing stellar age in the range of ages of disk population stars.

Barbera, M.↗
NASA NTRS1993

X-ray emission at the low-mass end - Results from an extensive Einstein Observatory survey

We have used available IPC data and a critical compilation of cataloged optical data to measure the 0.16-3.5 keV X-ray emission from 88 K and 169 M stars of luminosity classes IV, V, and VI within 25 pc from the Sun. The IPC detected 54 out of the 88 K stars, 70 out of the 138 M stars with M(v) less than 13.4, and 15 out of the 31 fainter M stars. We have identified a subsample of surveyed stars that is statistically representative of the population of K and M stars in the solar neighborhood. On the basis of this subsample (1) we have shown the occurrence of a drop in the level of X-ray emission for M stars later than approximately M5; (2) we have built unbiased maximum likelihood X-ray luminosity functions for the K, early M, and late M stars; (3) we have confirmed, both for K and M stars, the decrease of X-ray luminosity with increasing stellar age in the range of ages of disk population stars: and (4) we have shown that no obvious correlation is present between X-ray and bolometric luminosities in the entire representative samples of K and M stars, but only within flare stars which also seem to mark a saturation in X-ray luminosity level.

Barbera, M.↗
DOE OSTI2022

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Measurement of the top quark pole mass using $ \textrm{t}\overline{\textrm{t}} $+jet events in the dilepton final state in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A measurement of the top quark pole mass $\mathcal{m}^\text{pole}_\text{t}$ in events where a top quark-antiquark pair ($\text{t}\bar{\text{t}}$) is produced in association with at least one additional jet ($\text{t}\bar{\text{t}}$ +jet) is presented. This analysis is performed using proton-proton collision data at $\sqrt{s}$ = 13 TeV collected by the CMS experiment at the CERN LHC, corresponding to a total integrated luminosity of 36.3 fb -1 . Events with two opposite-sign leptons in the final state (e + e – , μ + μ – , e ± μ ∓ ) are analyzed. The reconstruction of the main observable and the event classification are optimized using multivariate analysis techniques based on machine learning. The production cross section is measured as a function of the inverse of the invariant mass of the $\text{t}\bar{\text{t}}$ +jet system at the parton level using a maximum likelihood unfolding. Given a reference parton distribution function (PDF), the top quark pole mass is extracted using the theoretical predictions at next-to-leading order. For the ABMP16NLO PDF, this results in $\mathcal{m}^\text{pole}_\text{t}$ = 172.93 ± 1.36 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1967

Computer program reduces calculation time of normal response functions

FORTRAN 2 computer program rapidly calculates parameters of maximum likelihood estimates from sensitivity experiment data populations. The program uses the Newton-Raphson iterative procedure to calculate the mean and standard deviation of portions of the cumulative normal response function.

Alexander, M. J.↗
NASA NTRSDate not supplied

Problem Complexity and LLM: H-M Team Reliability in Challenging Environments

In traditional human-machine operations, the functional decomposition of actions and responsibilities among various agents is assigned a priori. For instance, in current air traffic operations, although assisted by software, human pilots have the ultimate control of aircraft. Multiagent human-machine and machine-machine systems will face problems of varying and potentially unpredicted complexity in future challenging scenarios of planetary, en route, and orbital activities. Hence, it is important to enable dynamic transfer of decision-making to appropriate team members, human or machine, depending on which agent is best equipped to solve that specific problem on a time budget. In this paper, we consider aspects of problem-solving and its modeling that affect the outcomes of decision-making as a function of solution quality and the likelihood of solving the problem on a required time budget. We focus on Large Language Models (LLM) as potential machine teammates and conclude that practical, predictive modeling of their performance, at the current stage of their development, is infeasible. Simple examples help us illustrate that current LLM will require fundamental advancements to provide reliable support in team decision-making, especially in safety-critical and time-critical domains. The study is not meant to diminish the value of the remarkable capabilities of LLM, but rather to gain a better understanding of the technology’s appropriate use and the needed additions.

function allocation↗
DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗
DOE OSTI2021

A selection function toolbox for subsets of astronomical catalogues

ABSTRACT Large catalogues are ubiquitous throughout astronomy, but most scientific analyses are carried out on smaller samples selected from these catalogues by chosen cuts on catalogued quantities. The selection function of that scientific sample – the probability that a star in the catalogue will satisfy these cuts and so make it into the sample – is thus unique to each scientific analysis. We have created a general framework that can flexibly estimate the selection function of a sample drawn from a catalogue in terms of position, magnitude and colour. Our method is unique in using the binomial likelihood and accounting for correlations in the selection function across position, magnitude, and colour using Gaussian processes and spherical harmonics. We have created a new open-source python package selectionfunctiontoolbox that implements this framework and used it to make three different estimates of the APOGEE DR16 red giant sample selection function, as a subset of 2MASS, with each estimate using an increasing amount of technical knowledge of the APOGEE targeting. In a companion paper, we applied our methodology to derive estimates of the astrometric and spectroscopic selection functions of Gaia EDR3. Our framework will make it trivial for astrophysicists to estimate the selection function that they should be using with the custom sample of stars that they have chosen to answer their scientific question.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Magnetic character of the low-energy enhancement in 70 Zn

The transition of an atomic nucleus from an excited state to a state of lower energy is typically accompanied by the emission of photons. The likelihood of photon emission is a function of the photon energy and of properties of the initial and final nuclear states and is known as the photon (or γ-ray) strength function (γSF). Of the features present in the γSF, an enhancement in the low-energy region has been observed in some nuclei, but the electromagnetic nature of this enhancement remains an open question. This low-energy enhancement in the γSF significantly impacts our understanding of how elements are created in stars and the structure of the nucleus. Here we present the results from an experiment on the γSF of 70 Zn that show that the enhancement in the γSF of this nucleus is of magnetic nature. This observation helps address a long-open question about γSFs in atomic nuclei and directly impacts our predictive capabilities in nuclear science.

Physics - Nuclear physics and radiation physics↗