Search NASA⌕ Search

SEARCH · Search NASA

Results for “Partial correlation graph”

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.

9 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Information-incorporated gene network construction with FDR control

Abstract Motivation Large-scale gene expression studies allow gene network construction to uncover associations among genes. To study direct associations among genes, partial correlation-based networks are preferred over marginal correlations. However, FDR control for partial correlation-based network construction is not well-studied. In addition, currently available partial correlation-based methods cannot take existing biological knowledge to help network construction while controlling FDR. Results In this paper, we propose a method called Partial Correlation Graph with Information Incorporation (PCGII). PCGII estimates partial correlations between each pair of genes by regularized node-wise regression that can incorporate prior knowledge while controlling the effects of all other genes. It handles high-dimensional data where the number of genes can be much larger than the sample size and controls FDR at the same time. We compare PCGII with several existing approaches through extensive simulation studies and demonstrate that PCGII has better FDR control and higher power. We apply PCGII to a plant gene expression dataset where it recovers confirmed regulatory relationships and a hub node, as well as several direct associations that shed light on potential functional relationships in the system. We also introduce a method to supplement observed data with a pseudogene to apply PCGII when no prior information is available, which also allows checking FDR control and power for real data analysis. Availability and implementation R package is freely available for download at https://cran.r-project.org/package=PCGII.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

A convex pseudolikelihood framework for high dimensional partial correlation estimation with convergence guarantees

Sparse high dimensional graphical model selection is a topic of much interest in modern day statistics. A popular approach is to apply l1-penalties to either parametric likelihoods, or regularized regression/pseudolikelihoods, with the latter having the distinct advantage that they do not explicitly assume Gaussianity. As none of the popular methods proposed for solving pseudolikelihood-based objective functions have provable convergence guarantees, it is not clear whether corresponding estimators exist or are even computable, or if they actually yield correct partial correlation graphs. We propose a new pseudolikelihood-based graphical model selection method that aims to overcome some of the shortcomings of current methods, but at the same time retain all their respective strengths. In particular, we introduce a novel framework that leads to a convex formulation of the partial covariance regression graph problem, resulting in an objective function comprised of quadratic forms. The objective is then optimized via a coordinate wise approach. The specific functional form of the objective function facilitates rigorous convergence analysis leading to convergence guarantees; an important property that cannot be established by using standard results, when the dimension is larger than the sample size, as is often the case in high dimensional applications. These convergence guarantees ensure that estimators are well defined under very general conditions and are always computable. In addition, the approach yields estimators that have good large sample properties and also respect symmetry. Furthermore, application to simulated and real data, timing comparisons and numerical convergence is demonstrated. We also present a novel unifying framework that places all graphical pseudolikelihood methods as special cases of a more general formulation, leading to important insights.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

First-principles treatment of Mott insulators: linearized QSGW+DMFT approach

The theoretical understanding of emergent phenomena in quantum materials is one of the greatest challenges in condensed matter physics. In contrast to simple materials such as noble metals and semiconductors, macroscopic properties of quantum materials cannot be predicted by the properties of individual electrons. One of the examples of scientific importance is strongly correlated electron system. Neither localized nor itinerant behaviors of electrons in partially filled 3d, 4f, and 5f orbitals give rise to rich physics such as Mott insulators, high-temperature superconductors, and superior thermoelectricity, but hinder quantitative understanding of low-lying excitation spectrum. Here we present a new first-principles approach to strongly correlated solids. It is based on a combination of the quasiparticle self-consistent GW approximation and the dynamical mean-field theory. The sole input in this method is the projector to the set of correlated orbitals for which all local Feynman graphs are being evaluated. For that purpose, we choose very localized quasiatomic orbitals spanning large energy window, which contains most strongly hybridized bands, as well as upper and lower Hubbard bands. The self-consistency is carried out on the Matsubara axis. This method enables the first-principles study of Mott insulators in both their paramagnetic and antiferromagnetic phases. We illustrate the method on the archetypical charge transfer correlated insulators La 2 CuO 4 and NiO, and obtain spectral properties and magnetic moments in good agreement with experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2016

SU-D-201-05: On the Automatic Recognition of Patient Safety Hazards in a Radiotherapy Setup Using a Novel 3D Camera System and a Deep Learning Framework

Purpose: Patient safety hazards such as a wrong patient/site getting treated can lead to catastrophic results. The purpose of this project is to automatically detect potential patient safety hazards during the radiotherapy setup and alert the therapist before the treatment is initiated. Methods: We employed a set of co-located and co-registered 3D cameras placed inside the treatment room. Each camera provided a point-cloud of fraxels (fragment pixels with 3D depth information). Each of the cameras were calibrated using a custom-built calibration target to provide 3D information with less than 2 mm error in the 500 mm neighborhood around the isocenter. To identify potential patient safety hazards, the treatment room components and the patient’s body needed to be identified and tracked in real-time. For feature recognition purposes, we used a graph-cut based feature recognition with principal component analysis (PCA) based feature-to-object correlation to segment the objects in real-time. Changes in the object’s position were tracked using the CamShift algorithm. The 3D object information was then stored for each classified object (e.g. gantry, couch). A deep learning framework was then used to analyze all the classified objects in both 2D and 3D and was then used to fine-tune a convolutional network for object recognition. The number of network layers were optimized to identify the tracked objects with >95% accuracy. Results: Our systematic analyses showed that, the system was effectively able to recognize wrong patient setups and wrong patient accessories. The combined usage of 2D camera information (color + depth) enabled a topology-preserving approach to verify patient safety hazards in an automatic manner and even in scenarios where the depth information is partially available. Conclusion: By utilizing the 3D cameras inside the treatment room and a deep learning based image classification, potential patient safety hazards can be effectively avoided.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2001

Ceiling Fires Studied to Simulate Low-Gravity Fires

A unique new way to study low-gravity flames in normal gravity has been developed. To study flame structure and extinction characteristics in low-stretch environments, a normal gravity low-stretch diffusion flame was generated using a cylindrical PMMA sample of varying large radii, as shown in the photograph. These experiments have demonstrated that low-gravity flame characteristics can be generated in normal gravity through the proper use of scaling. On the basis of this work, it is feasible to apply this concept toward the development of an Earth-bound method of evaluating material flammability in various gravitational environments from normal gravity to microgravity, including the effects of partial gravity low-stretch rates such as those found on the Moon (1/6g) or Mars (1/3g). During these experiments, the surface regression rates for PMMA were measured for the first time over the full range of flammability in air, from blowoff at high stretch, to quenching at low stretch, as plotted in the graph. The solid line drawn through the central portion of the data (3<a<100 sec-1) has a slope of unity, which indicates regression is proportional to stretch. The figure coordinates assume that the values of stretch are equivalent, whether derived from forced stretch or from buoyant stretch. The excellent correlation of the regression-rate data over the two-order-of-magnitude variation of stretch shows the reasonableness of this assumption.

Olson, Sandra L.↗
DOE OSTI2018

Socioeconomic and particulate air pollution correlates of heart disease risk

Highlights: • Heart disease risk varies with age, sex, smoking, income, education, and particulate (PM2.5) air pollution. • How these risk factors interact is unclear. • We examine correlates of self-reported heart disease risk in a large data set. • Income is negatively associated with both PM2.5 exposure and heart disease risk. • Income is an important confounder of air pollution-heart disease associations. • Nonparametric methods can quantify associations and suggest causal interpretations. How does risk of heart disease depend on age, sex, smoking, income, education, marital status, and outdoor concentrations of fine particulate matter (PM2.5)? We join data available from the Centers for Disease Control and Prevention (CDC) Behavioral Risk Factor Surveillance (BRFSS) System for years 2008–2012 to US Environmental Protection Agency (EPA) data on county-specific concentrations of fine particulate matter (PM2.5) to quantify associations among these variables and to explore possible causal interpretations. Low income is identified as a direct cause of increased heart disease risk in this data set. The effect depends on age and sex: it is most pronounced for men under age 70 and for women under age 80. Income is significantly associated with all of the other variables examined and confounds the association between PM2.5 and heart disease risk. This association is significant in regression models that exclude income, but not in regression models that include it, both in the data set as a whole and in the subset of observations with PM2.5 < 15 μg/m{sup 3}. Causal directed acyclic graph (DAG) models and non-parametric model ensemble partial dependence plots confirm that higher incomes reduce heart disease risk, consistent with previous observations of socioeconomic gradients in health risks. They support interpretation of this as a robust causal relation apparent in non-parametric analyses, and hence independent of any specific parametric modeling assumptions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Eigenvector centrality for characterization of protein allosteric pathways

Determining the principal energy-transfer pathways responsible for allosteric communication in biomolecules remains challenging, partially due to the intrinsic complexity of the systems and the lack of effective characterization methods. In this work, we introduce the eigenvector centrality metric based on mutual information to elucidate allosteric mechanisms that regulate enzymatic activity. Moreover, we propose a strategy to characterize the range of correlations that underlie the allosteric processes. We use the V-type allosteric enzyme imidazole glycerol phosphate synthase (IGPS) to test the proposed methodology. The eigenvector centrality method identifies key amino acid residues of IGPS with high susceptibility to effector binding. The findings are validated by solution NMR measurements yielding important biological insights, including direct experimental evidence for interdomain motion, the central role played by helix h α 1 , and the short-range nature of correlations responsible for the allosteric mechanism. Beyond insights on IGPS allosteric pathways and the nature of residues that could be targeted by therapeutic drugs or site-directed mutagenesis, the reported findings demonstrate the eigenvector centrality analysis as a general cost-effective methodology to gain fundamental understanding of allosteric mechanisms at the molecular level.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗