Search NASA⌕ Search

SEARCH · Search NASA

Results for “robust regression”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantile-specific heritability of sibling leptin concentrations and its implications for gene-environment interactions

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., leptin) is high or low relative to its distribution. Leptin concentrations are strongly related to adiposity, whose heritability is quantile dependent. Whether inheritance of leptin concentrations is quantile dependent, and whether this explains the greater heritability in women than men in accordance with their greater adiposity, and explains other gene-environment interactions, remains to be determined. Therefore, leptin and leptin receptor concentrations from 3068 siblings in 1133 sibships from the Framingham Heart Study Third Generation Cohort were analyzed. Free leptin index (FLI) was calculated as the ratio of leptin to soluble leptin receptor concentrations. Full-sib (β FS ) regression slopes were robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. The analyses showed β FS increased significantly with increasing percentiles of the offspring’s age- and sex-adjusted leptin distribution (P linear = 0.0001), which was accelerated at the higher concentrations (P quadratic = 0.0003). β FS at the 90th percentile (0.418 ± 0.066) was 4.7-fold greater than at the 10th percentile (0.089 ± 0.032, P difference = 3.6 × 10 –6 ). Consistent with quantile-dependent expressivity, the β FS was greater in female sibs, which was attributable to their higher leptin concentrations. Reported gene-environment interactions involving adiposity and LEP , LEPR , MnSOD , PPARγ , PPARγ2 , and IRS-1 polymorphisms were consistent with quantile-dependent expressivity of leptin concentrations. β FS for leptin receptor concentrations and free leptin index also increased significantly with increasing percentiles of their distributions (P linear = 0.04 and P linear = 8.5 × 10 –6 , respectively). In conclusion, inherited genetic and shared environmental effects on leptin concentrations were quantile dependent, which likely explains male–female differences in heritability and some gene-environment interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

An ℓ 0 ℓ 2 -norm regularized regression model for construction of robust cluster expansions in multicomponent systems

In this work we introduce ℓ 0 ℓ 2 -norm regularization and hierarchy constraints into linear regression for the construction of cluster expansions to describe configurational disorder in materials. The approach is implemented through mixed integer quadratic programming (MIQP). The ℓ 2 -norm regularization is used to suppress intrinsic data noise, while the ℓ 0 -norm is used to penalize the number of nonzero elements in the solution. The hierarchy relation between clusters imposes relevant physics and is naturally included by the MIQP paradigm. As such, sparseness and cluster hierarchy can be well optimized to obtain a robust, converged set of effective cluster interactions with improved physical meaning. We demonstrate the effectiveness of ℓ 0 ℓ 2 -norm regularization in two high-component disordered rocksalt cathode material systems, where we compare the cross-validation, convergence speed, and the reproduction of phase diagrams, voltage profiles, and Li-occupancy energies with those of the conventional ℓ 1 -norm regularized cluster expansion models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

How analytic choices can affect the extraction of electromagnetic form factors from elastic electron scattering cross section data

Scientists often try to incorporate prior knowledge into their regression algorithms, such as a particular analytic behavior or a known value at a kinematic endpoint. Unfortunately, there is often no unique way to make use of this prior knowledge, and thus, different analytic choices can lead to very different regression results from the same set of data. In this paper, to illustrate this point in the context of the proton electromagnetic form factors, we use the Mainz elastic data with its 1422 cross section points and 31 normalization parameters. Starting with a complex unbound non-linear regression, we will show how the addition of a single theory-motivated constraint removes an oscillation from the magnetic form factor and shifts the extracted proton charge radius. We then repeat both regressions using the same algorithm, but with a rebinned version of the Mainz dataset. These examples illustrate how analytic choices, such as the function that is being used or even the binning of the data, can dramatically affect the results of a complex regression. These results also demonstrate why it is critical when using regression algorithms to have either a physical model in mind or a firm mathematical basis

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

SU-C-16A-05: OAR Dose Tolerance Recommendations for Prostate and Cervical HDR Brachytherapy: Dose Versus Volume Metrics

Purpose: HDR brachytherapy consensus dose tolerance recommendations for organs at risk (OARs) remain widely debated. Prospective trials reporting metrics must be sufficiently data-dense to assess adverse affects and identify optimally predictive tolerances. We explore the tradeoffs between reporting dose-metrics versus volume-metrics and the potential impact on trial outcome analysis and tolerance recommendations. Methods: We analyzed 26 prostate patients receiving 15 Gy HDR single-fraction brachytherapy boost to 45 Gy external beam radiation therapy and 28 cervical patients receiving 28 Gy HDR brachytherapy monotherapy in 4 fractions using 2 implants. For each OAR structure, a robust linear regression fit was performed for the dose-metrics as a function of the volume-metrics. The plan quality information provided by recommended dose-metric and volume-metric values were compared. Results: For prostate rectal dose, D2cc and V75 lie close to the regression line, indicating they are similarly informative. Two outliers for prostate urethral dose are substantially different from the remaining cohort in terms of D0.1cc and V75, but not D1cc, suggesting the choice of reporting dose metric is essential. For prostate bladder and cervical bladder, rectum, and bowel, dose outliers are more apparent via V75 than recommended dose-metrics. This suggests that for prostate bladder dose and all cervical OAR doses, the recommended volume-metrics may be better predictors of clinical outcome than dose-metrics. Conclusion: For plan acceptance criteria, dose and volume-metrics are reciprocally equivalent. However, reporting dosemetrics or volume-metrics alone provides substantially different information. Our results suggest that volume-metrics may be more sensitive to differences in planned dose, and if one metric must be chosen, volumemetrics are preferable. However, reporting discrete DVH points severely limits the ability to identify planning tolerances most predictive of adverse effects. Thus, we recommend that full OAR DVH reporting be required for future prospective trials.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1993

A Look at a Day of Data from the Topex/Poseidon GPS Receiver

Launched on August 10, 1992, the Topex/Poseidon satellite is a joint project of NASA and the French space agency CNES. The pupose of this Earth orbiting altimetric satellite (for brevity, referred to as Topex) is to enable monitoring of the ocean topography to higher levels of accuracy than previously achieved.

robust regression spline graphics↗
NASA NTRS2019

Radar Altimetry as a Proxy for Determining Terrestrial Water Storage Variability in Tropical Basins

The Gravity Recovery and Climate Experiment (GRACE) mission has provided us with unforeseen information on terrestrial water-storage (TWS) variability, contributing to our understanding of global hydrological processes, including hydrological extreme events and anthropogenic impacts on water storage. Attempts to decompose GRACE-based TWS signals into its different water storage layers, i.e., surface water storage (SWS), soil moisture, groundwater and snow, have shown that SWS is a principal component, particularly in the tropics, where major rivers flow over arid regions at high latitudes. Here, we demonstrate that water levels, measured with radar altimeters at a limited number of locations, can be used to reconstruct gridded GRACE-based TWS signals in the Amazon basin, at spatial resolutions ranging from 0.5 to 3°, with mean absolute errors (MAE) as low as 2.5 cm and correlations as high as 0.98. We show that, at 3° spatial resolution, spatially-distributed TWS time series can be precisely reconstructed with as few as 41 water-level time series located within the basin. The proposed approach is competitive when compared to existing TWS estimates derived from physically based and computationally expensive methods. Also, a validation experiment indicates that TWS estimates can be extrapolated to periods beyond that of the model regression with low errors. The approach is robust, based on regression models and interpolation techniques, and offers a new possibility to reproduce spatially and temporally distributed TWS that could be used to fill inter-mission gaps and to extend GRACE-based TWS time series beyond its timespan.

Terrestrial water storage↗
NASA NTRS1984

Use of Thematic Mapper for water quality assessment

The evaluation of simulated TM data obtained on an ER-2 aircraft at twenty-five predesignated sample sites for mapping water quality factors such as conductivity, pH, suspended solids, turbidity, temperature, and depth, is discussed. Using a multiple regression for the seven TM bands, an equation is developed for the suspended solids. TM bands 1, 2, 3, 4, and 6 are used with logarithm conductivity in a multiple regression. The assessment of regression equations for a high coefficient of determination (R-squared) and statistical significance is considered. Confidence intervals about the mean regression point are calculated in order to assess the robustness of the regressions used for mapping conductivity, turbidity, and suspended solids, and by regressing random subsamples of sites and comparing the resultant range of R-squared, cross validation is conducted.

Horn, E. M.↗
DOE OSTI2021

Quantile-dependent expressivity of serum C-reactive protein concentrations in family sets

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., C-reactive protein, CRP) is high or low relative to its distribution. We have previously shown that the heritabilities ( h 2 ) of coffee and alcohol consumption, postprandial lipemia, lipoproteins, leptin, adiponectin, adiposity, and pulmonary function are quantile-specific. Whether CRP heritability is quantile-specific is currently unknown. Methods Serum CRP concentrations from 2,036 sibships and 6,144 offspring-parent pairs were analyzed from the Framingham Heart Study. Quantile-specific heritability from full-sib (β FS , h 2 ={(1 + 8r spouse β FS ) 0.5 - 1}/(2r spouse )) and offspring-parent regression slopes (β OP , h 2 = 2β OP /(1 + r spouse )) were estimated robustly by quantile regression with nonparametric significance determined from 1,000 bootstrap samples. Results Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s age- and sex-adjusted CRP distribution when estimated from β OP ( P trend = 0.0004): 0.02 ± 0.01 at the 10th, 0.04 ± 0.01 at the 25th, 0.10 ± 0.02 at the 50th, 0.20 ± 0.05 at the 75th, and 0.33 ± 0.10 at the 90th percentile, and when estimated from β FS ( P trend = 0.0008): 0.03±0.01 at the 10th, 0.06 ± 0.02 at the 25th, 0.14 ± 0.03 at the 50th, 0.24 ± 0.05 at the 75th, and 0.53 ± 0.21 at the 90th percentile. Conclusion Heritability of serum CRP concentration is quantile-specific, which may explain or contribute to the inflated CRP differences between CRP (rs1130864, rs1205, rs1800947, rs2794521, rs3091244), FGB (rs1800787), IL-6 (rs1800795, rs1800796), IL6R (rs8192284), TNF- α (rs1800629) and APOE genotypes following CABG surgery, stroke, TIA, curative esophagectomy, intensive periodontal therapy, or acute exercise; during acute coronary syndrome or Staphylococcus aureus bacteremia; or in patients with chronic rheumatoid arthritis, diabetes, peripheral arterial disease, ankylosing spondylitis, obesity or inflammatory bowel disease or who smoke.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2010

Model Robust Calibration: Method and Application to Electronically-Scanned Pressure Transducers

This article presents the application of a recently developed statistical regression method to the controlled instrument calibration problem. The statistical method of Model Robust Regression (MRR), developed by Mays, Birch, and Starnes, is shown to improve instrument calibration by reducing the reliance of the calibration on a predetermined parametric (e.g. polynomial, exponential, logarithmic) model. This is accomplished by allowing fits from the predetermined parametric model to be augmented by a certain portion of a fit to the residuals from the initial regression using a nonparametric (locally parametric) regression technique. The method is demonstrated for the absolute scale calibration of silicon-based pressure transducers.

Walker, Eric L.↗
DOE OSTI2026

Poisson Log-Normal Process for Count Data Prediction

Modeling count data is important in physics and other scientific disciplines, where measurements often involve discrete, non-negative quantities such as photon or neutrino detection events. Traditional parametric approaches can be trained to generate integer-count predictions but may struggle with capturing complex, non-linear dependencies often observed in the data. Gaussian process (GP) regression provides a robust non-parametric alternative to modeling continuous data; however, it cannot generate integer outputs. We propose the Poisson Log-Normal (PoLoN) process, a framework that employs GP to model Poisson log-rates. As in GP regression, our approach relies on the correlations between data points captured via GP kernel structure rather than explicit functional parameterizations. We demonstrate that the PoLoN predictive distribution is Poisson-LogNormal and provide an algorithm for optimizing kernel hyperparameters. Furthermore, we adapt the PoLoN approach to the problem of detecting weak localized signals superimposed on a smoothly varying background - a task of considerable interest in many areas of science and engineering. Our framework allows us to predict the strength, location and width of the detected signals. We evaluate PoLoN's performance using both synthetic and real-world datasets, including the open dataset from CERN which was used to detect the Higgs boson at the Large Hadron Collider. Our results indicate that the PoLoN process can be used as a non-parametric alternative for analyzing, predicting, and extracting signals from integer-valued data.

Saha, Anushka [Rutgers U., Piscataway]↗
DOE OSTI2026

Data driven investigation to understand the influence of total solids on biological biogas upgrading

In situ biogas upgrading achieves CO 2 conversion to CH 4 via hydrogenotrophic methanogenesis; however, gas-liquid mass transfer constraints limit the upgrading performance. Recognizing that optimization studies often underrepresent the effects of total solids (TS) and organic loading rate (OLR), this study undertook a holistic, statistics driven assessment of operating conditions for in situ H 2 assisted biogas upgrading, centering the analysis on TS and OLR. A dataset of 31 studies was compiled and comprised 99 observations. A rigorous analytical framework was employed, combining data standardization, fixed- and random-effects (REML) weighted regressions with cluster-robust errors, stratified analyses, and machine learning. Mixed-effects meta regression indicated that TS was the main factor explaining differences of methane fraction (CH 4 %) when considering the between studies heterogeneity. Focusing on a near-stoichiometric subset (H 2 /CO 2 ≈ 4:1), TS remained significant. Stratified results showed a stronger negative relationship between TS and CH 4 % in UASB reactors than in CSTRs, with a negative effect under mesophilic conditions and no significant effect under thermophilic conditions. A Random Forest model corroborated the statistical findings, consistently ranking H 2 /CO 2 ratio, OLR, TS, and hydrogen injection rate (HIR) as the most influential predictors. These findings delineate trends across increasing TS levels, particularly between 1% and 10%, and provide preliminary insights for TS above 15% in in situ biogas upgrading. They further provide insights for the influence of TS by reactor type and temperature, thereby advancing the evidence base for implementing biological CO 2 conversion to CH 4 in practice.

In situ biogas upgrading↗
DOE OSTI2016

AGR-1, AGR-2 and AGR-3/4 Dimensional Change Data Analysis

A series of Advanced Gas Reactor (AGR) experiments have been completed in the Advanced Test Reactor at Idaho National Laboratory in support of qualification and development of tristructural isotropic fuel. Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder. These capsules are instrumented with thermocouples embedded in the graphite, enabling temperature control. The fuel compacts are composed of fuel particles surrounded by a graphitic A3 matrix material. Dimensional change in AGR fuel compacts is vital because the swelling or shrinkage affects the size of the gas gaps that are used to control temperatures. Analysis of dimensional change in the AGR fuel compacts is needed to establish the variables directly relating to compact shrinkage. The variables initially identified for consideration were matrix density, compact density, fuel packing fraction, uranium loading, fuel particle diameter, cumulative fast neutron fluence, and volume average time average fuel temperature. In addition to the data available from the AGR experiments, the analysis included specimens formed from the same A3 matrix material used in Advanced Graphite Creep (AGC) experiments, which provide graphite creep data during irradiation for design and licensing purposes. The primary purpose of including the AGC specimens was to encompass dimensional behavior at zero packing fraction, zero uranium loading, and zero particle diameter. All possible combinations of first-order variable regressions were considered in the analysis. The study focused on identifying the best regression models for percent change in diameter, length, and volume. Bootstrap analysis was used to ensure the resulting regression models were robust and well-performing. The variables identified as very significant in predicting change in one or more dimensions (diameter, length, and volume) are volume average time average temperature, fast fluence, compact density, packing fraction and fuel particle diameter. Due to the presence of confounding effects between several variables, interpretation of these results is equivocal; the use of multiple statistical tests provides additional confidence in the conclusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗