Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Cutting-Balloon Angioplasty Versus Balloon Angioplasty as Treatment for Short Atherosclerotic Lesions in the Superficial Femoral Artery: Randomized Controlled Trial

Purpose: To evaluate the treatments of a short-segment atherosclerotic stenosis in the superficial femoral arteries with the cutting balloon angioplasty (CBA) versus conventional balloon angioplasty [percutaneous transluminal angioplasty (PTA)] in a randomized controlled trial. Material and Methods: The study group comprised 60 patients (33 men, 27 women; average age 64 years) with a short ({<=}5 cm) focal SFA de novo atherosclerotic stenosis associated with a history of intermittent claudication or rest pain. The primary end point of this study was the rate of binary restenosis in the treated segment 12 months after the intervention. All patients were evenly randomized to either the PTA or CBA treatment arms. Follow-up angiograms and ankle-brachial index (ABI) measurements were performed after 12 months. The evaluation of the restenosis rates and factors influencing its occurrence were calculated by logistic regression analysis. Results: In the intention-to-treat analysis, restenosis rates after 2-month follow-up were 9 of 30 (30 %) in the PTA group and 4 of 30 (13 %) in the CBA group (p = 0.117). In the actual treatment analysis, after exclusion of patients who required nitinol stent placement for a suboptimal result after angioplasty alone (5 patients in the PTA group and none in the CBA group), restenosis rates were 9 of 25 (36 %) and 4 of 30 (13 %), respectively (p = 0.049). In the intention-to-treat analysis there were also significant differences in ABI values between the PTA and CBA groups at 0.77 {+-} 0.11 versus 0.82 {+-} 0.12, respectively (p = 0.039), at 12 months. Conclusion: Based on the presented results of the trial, CBA seems to be a safer and more effective than PTA for treatment of short atherosclerotic lesions in the superior femoral artery.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Role of Internal Mammary Node Radiation as a Part of Modern Breast Cancer Radiation Therapy: A Systematic Review

Purpose: Despite data from multiple randomized trials, the role of internal mammary lymph node irradiation as a part of regional nodal irradiation (IMLN RT–RNI) remains unanswered. Recent noteworthy data and modern RT techniques might identify a subset of patients who will benefit from IMLN RT–RNI, lending insight into the balance between improved outcomes and acceptable toxicity. We evaluated the current role of IMLN RT–RNI by analyzing randomized, prospective, and retrospective data. Methods and Materials: In accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines, a review of the published data was performed using PubMed to evaluate published studies from 1994 to 2015. The information evaluated included the number of patients, follow-up period, technical aspects of RT, and outcomes (clinical outcomes, complications/toxicity). Results: We included 16 studies (4 randomized, 4 nonrandomized, 7 retrospective, and 1 meta-analysis). Although older randomized trials failed to show differences in clinical outcomes or toxicity with IMLN RT–RNI, recent randomized data suggest the potential for improved outcomes, including overall survival, with IMLN RT–RNI. Furthermore, nonrandomized data have suggested a potential benefit for central tumors with IMLN RT–RNI. Although recent data have suggested a potential increase in pulmonary complications with IMLN RT–RNI with the use of advanced radiation techniques, toxicity rates remain low with limited cardiac toxicity data available. Conclusions: Increasing data from recent randomized trials support the use of IMLN RT–RNI. IMLN RT can be considered based on the inclusion of IMLN RT as a part of RNI in recent trials and the inclusion criteria from IMLN RT–RNI trials and for patients with central or medial tumors and axillary disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Assessment of cardiac single-photon emission computed tomography performance using a scanning linear observer

Purpose: Single-photon emission computed tomography (SPECT) is widely used to detect myocardial ischemia and myocardial infarction. It is important to assess and compare different SPECT system designs in order to achieve the highest detectability of cardiac defects. Methods: Whitaker et al.'s study ['Estimating random signal parameters from noisy images with nuisance parameters: linear and scanning-linear methods,' Opt. Express 16(11), 8150-8173 (2008)] on the scanning linear observer (SLO) shows that the SLO can be used to estimate the location and size of signals. One major advantage of the SLO is that it can be used with projection data rather than with reconstruction data. Thus, this observer model assesses the overall hardware performance independent of any reconstruction algorithm. In addition, the computation time of image quality studies is significantly reduced. In this study, three systems based on the design of the GE cadmium zinc telluride-based dedicated cardiac SPECT camera Discovery 530c were assessed. This design, which is officially named the Alcyone Technology: Discovery NM 530c, was commercialized in August, 2009. The three systems, GE27, GE19, and GE13, contain 27, 19, and 13 detectors, respectively. Clinically, a human heart can be virtually segmented into three coronary artery territories: the left-anterior descending artery, left-circumflex artery, and right coronary artery. One of the most important functions of a cardiac SPECT system is to produce images from which a radiologist can accurately predict in which territory the defect exists [http://www.asnc.org/media/PDFs/PPReporting081511.pdf, Guideline from American Society of Nuclear Cardiology]. A good estimation of the extent of the defect from the projection images is also very helpful for determining the seriousness of the myocardial ischemia. In this study, both the location and extent of defects were estimated by the SLO, and the system performance was assessed by localization receiver operating characteristic (LROC) [P. Khurd and G. Gindi, 'Decision strategies maximizing the area under the LROC curve,' Proc. SPIE 5749, 150-161 (2005)] or estimation receiver operating characteristic (EROC) [E. Clarkson, 'Estimation receiver operating characteristic curve and ideal observers for combined detection/estimation tasks,' J. Opt. Soc. Am. A 24, B91-B98 (2007)] curves. Results: The area under the LROC/EROC curve (AULC/AUEC) and the true positive fraction (TPF) at a specific false positive fraction (FPF) can be treated as the figures of merit. For radii estimation with a 1 mm tolerance, the AUEC values of the GE27, GE19, and GE13 systems are 0.8545, 0.8488, and 0.8329, and the TPF at FPF = 5% are 77.1%, 76.46%, and 73.55%, respectively. The assessment of all three systems revealed that the GE19 system yields estimated information and cardiac defect detectability very close to those of the GE27 system while using eight fewer detectors. Thus, 30% of the expensive detector units can be removed with confidence. Conclusions: As the results show, a combination of the SLO and LROC/EROC curves can determine the configuration that yields the most relevant estimation/detection information. Thus, this is a useful method for assessing cardiac SPECT systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Fast Extraction and Characterization of Fundamental Frequency Events from a Large PMU Dataset using Big Data Analytics

A novel method for fast extraction of fundamental frequency events (FFE) based on measurements of frequency and rate of change of frequency by Phasor Measurement Units (PMU) is introduced. The method is designed to work with exceptionally large historical PMU datasets. Statistical analysis was used to extract the features and train Random Forest and Catboost classifiers. The method is capable of fast extraction of FFE from a historical dataset containing measurements from hundreds of PMUs captured over multiple years. The reported accuracy of the best algorithm for classification expressed as Area Under the receiver operating Characteristic curve reaches 0.98, which was obtained in out-of-sample evaluations on 109 system-wide events over 2 years observed at 43 PMUs. Then Minimum Volume Enclosing Ellipsoid Algorithm was used to further analyze the events. 93.72% events were correctly characterized, where average duration of the event as seen by the PMU was 9.93 sec.

Baembitov, Rashid↗
DOE OSTI2014

Extremely small test cell structure for resistive random access memory element with removable bottom electrode

We established a method of preparing an extremely small memory cell by fabricating a resistive random access memory (ReRAM) structure on the tip of a cantilever of an atomic force microscope. This structure has the high robustness against the drift of the cantilever, and the effective cell size was estimated to be less than 10 nm in diameter due to the electric field concentration at the tip of the cantilever, which was confirmed using electric field simulation. The proposed structure, which has a removable bottom electrode, enables not only the preparation of a tiny ReRAM structure but also the performance of unique experiments, by making the most of its high robustness against the drift of the cantilever.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Deep Learning Classification of Cheatgrass Invasion in the Western United States Using Biophysical and Remote Sensing Data

Cheatgrass (Bromus tectorum) invasion is driving an emerging cycle of increased fire frequency and irreversible loss of wildlife habitat in the western US. Yet, detailed spatial information about its occurrence is still lacking for much of its presumably invaded range. Deep learning (DL) has demonstrated success for remote sensing applications but is less tested on more challenging tasks like identifying biological invasions using sub-pixel phenomena. We compare two DL architectures and the more conventional Random Forest and Logistic Regression methods to improve upon a previous effort to map cheatgrass occurrence at >2% canopy cover. High-dimensional sets of biophysical, MODIS, and Landsat-7 ETM+ predictor variables are also compared to evaluate different multi-modal data strategies. All model configurations improved results relative to the case study and accuracy generally improved by combining data from both sensors with biophysical data. Cheatgrass occurrence is mapped at 30 m ground sample distance (GSD) with an estimated 78.1% accuracy, compared to 250-m GSD and 71% map accuracy in the case study. Furthermore, DL is shown to be competitive with well-established machine learning methods in a limited data regime, suggesting it can be an effective tool for mapping biological invasions and more broadly for multi-modal remote sensing applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Spinbox: tools for many-body quantum systems in a Monte Carlo context

Spinbox is a piece of software that facilitates quantum mechanical calculations relevant to Monte Carlo simulation of atomic nuclei. At the front lines of research on the nuclear many-body problem are a large number of supercomputer-scale simulation codes. These codes produce valuable results but can be hard to understand, especially for those without intimate knowledge of the relevant theoretical methods. Thus, tools that fill pedagogical roles are extremely valuable. Spinbox makes it easy for one to replicate and analyze the computational processes relevant to a Quantum Monte Carlo (QMC) simulation that may be difficult to understand/debug/analyze due to the scale of the corresponding simulation software. Spinbox is written in Python using other state-of-the-art Python modules for numerical calculations. While a number of Python libraries exist that are suited to general quantum many-body calculations, the motivation of Spinbox is quite particular. In Diffusion Monte Carlo methods (DMC, GFMC, AFDMC), the central calculation is the imaginary-time propagation of individual samples of the many-body wavefunction. Although quantum wavefunctions generally must be described by a probability distribution over a basis, DMC imbues particles (within one sample) with classical spatial coordinates. This method is unusual, so other Python packages are typically not set up to do this easily. Furthermore, the software has built-in options for nuclear systems assuming isospin symmetry, which can be set up with other libraries but is a nontrivial process to do so. Features: - numerical representation of samples of the many-body wavefunctions, including tensor-product states (used in AFDMC) - numerical representation of many-body operators, including tensor-product operators: general, spin, imaginary-time propagation, etc. - the correct associated arithmetic and algebra, implemented as class methods - classes for representing realistic nuclear two- and three-body Hamiltonians (e.g. Argonne V18, Illinois NNN) - large-scale parallel integration over random variables, crucial for the AFDMC method My goal is to make this package open source so that anyone may use it and contribute to it, particularly other researchers doing AFDMC calculations

Fox, Jordan↗
DOE OSTI2020

Interpreting High-resolution Spectroscopy of Exoplanets using Cross-correlations and Supervised Machine Learning

We present a new method for performing atmospheric retrieval on ground-based, high-resolution data of exoplanets. Our method combines cross-correlation functions with a random forest, a supervised machine-learning technique, to overcome challenges associated with high-resolution data. A series of cross-correlation functions are concatenated to give a “CCF-sequence” for each model atmosphere, which reduces the dimensionality by a factor of ∼100. The random forest, trained on our grid of ∼65,000 models, provides a likelihood-free method of retrieval. The precomputed grid spans 31 values of both temperature and metallicity, and incorporates a realistic noise model. We apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b and obtain a metallicity consistent with solar (logM = − 0.2 ± 0.2). Our retrieved transit chord temperature (T=6000{sub −200}{sup +0}K) is unreliable as strong ion lines lie outside of the extent of the training set, which we interpret as being indicative of missing physics in our atmospheric model. We compare our method to traditional nested sampling, as well as other machine-learning techniques, such as Bayesian neural networks. We demonstrate that the likelihood-free aspect of the random forest makes it more robust than nested sampling to different error distributions, and that the Bayesian neural network we tested is unable to reproduce complex posteriors. We also address the claim in Cobb et al. 2019 that our random forest retrieval technique can be overconfident but incorrect. We show that this is an artifact of the training set, rather than of the machine-learning method, and that the posteriors agree with those obtained using nested sampling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Online updating of context-aware landmark detectors for prostate localization in daily treatment CT images

Purpose: In image guided radiation therapy, it is crucial to fast and accurately localize the prostate in the daily treatment images. To this end, the authors propose an online update scheme for landmark-guided prostate segmentation, which can fully exploit valuable patient-specific information contained in the previous treatment images and can achieve improved performance in landmark detection and prostate segmentation. Methods: To localize the prostate in the daily treatment images, the authors first automatically detect six anatomical landmarks on the prostate boundary by adopting a context-aware landmark detection method. Specifically, in this method, a two-layer regression forest is trained as a detector for each target landmark. Once all the newly detected landmarks from new treatment images are reviewed or adjusted (if necessary) by clinicians, they are further included into the training pool as new patient-specific information to update all the two-layer regression forests for the next treatment day. As more and more treatment images of the current patient are acquired, the two-layer regression forests can be continually updated by incorporating the patient-specific information into the training procedure. After all target landmarks are detected, a multiatlas random sample consensus (multiatlas RANSAC) method is used to segment the entire prostate by fusing multiple previously segmented prostates of the current patient after they are aligned to the current treatment image. Subsequently, the segmented prostate of the current treatment image is again reviewed (or even adjusted if needed) by clinicians before including it as a new shape example into the prostate shape dataset for helping localize the entire prostate in the next treatment image. Results: The experimental results on 330 images of 24 patients show the effectiveness of the authors’ proposed online update scheme in improving the accuracies of both landmark detection and prostate segmentation. Besides, compared to the other state-of-the-art prostate segmentation methods, the authors’ method achieves the best performance. Conclusions: By appropriate use of valuable patient-specific information contained in the previous treatment images, the authors’ proposed online update scheme can obtain satisfactory results for both landmark detection and prostate segmentation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Evaluating disease surveillance strategies for early outbreak detection in contact networks with varying community structure

Disease surveillance systems allow public health agencies to respond to emerging diseases before they become widespread. Developing such systems requires identifying optimal ways to monitor in the context of an epidemic outbreak; this problem is known as sensor selection. Contact networks represent the dynamics of interaction in a population and are used to model how a disease spreads in a population and to explore strategies of sensor selection. We evaluated five sensor selection strategies on their ability to provide an early warning of a COVID-like outbreak in synthetic contact networks encapsulated in four network scenarios. Three of these scenarios assessed different aspects of community structure. The fourth scenario employed a contact network representing the population and interactions of 6.8 million people in New York City, constructed from an agent-based simulation using census and transportation data. This scenario exemplifies how sensor selection strategies may perform in a real-world, urban context. Our findings suggest that the choice of the optimal strategy depends heavily on the community structure of the network. Strategies that select highly connected nodes or maximize network coverage are the optimal surveillance strategy for outbreak detection in many network community structures. However, a naive implementation of these strategies may fail to provide an early warning at all—including in the New York City scenario. Moreover, these methods are impractical for real-world use as they require knowledge of the underlying contact network. Instead, a selection strategy that starts with a set of random nodes and then performs a random walk through a chain of neighbors reliably provides early warnings without requiring prior knowledge of the network. We find this method, called “random chain”, to be the most pragmatic for implementation in a real-world disease surveillance context.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Novel 2D velocity estimation method for large transient events in plasmas

Dynamics of fast transient events are challenging to be analyzed with high time resolution. Such events can occur in fusion plasmas such as the filaments during edge-localized modes (ELMs). Here, we present a robust method—the spatial displacement estimation—for estimating the displacements of structures with fast dynamics from high spatial and time resolution imaging diagnostics [e.g., gas-puff imaging (GPI)] with sampling time temporal resolution. First, a background suppression method is shown, which suppresses the slowly time-evolving and spatially non-uniform background in the signal. In the second step, a two-dimensional polynomial trend subtraction method is presented to tackle the remaining polynomial order trend in the signal. After performing these pre-processing steps, the spatial displacement of the propagating structure is estimated from the two-dimensional spatial cross-correlation coefficient function calculated between consecutive frames. The method is tested for its robustness and accuracy by simulated Gaussian events and spatially displaced random noise. An example application of the method is presented on propagating ELM filaments measured by the GPI system on the National Spherical Torus Experiment spherical tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

UNDERSTANDING THE SEMI-PROBABILISTIC APPROACHES IN STRUCTURAL RELIABILITY USED TO SET DESIGN RELIABILITY TARGETS FOR GRAPHITE COMPONENTS USING ASME BPVC METHODS

Graphite is a quasi-brittle material, resulting in random variability in tensile strength distributions. To account for the random variability in strength, HHA-3000 of the ASME BPVC provides two semi-probabilistic methods for qualifying nuclear graphite components in the design stage, the simplified and full assessments. The full and simplified assessments apply statistical methods to engineering-based design problems. This is often referred to as reliability-based design. Reliability-based design (RBD) is a method to develop reliable designs by accounting for uncertainties and result in small chances of failure when also considering safety factors. RBDs provide reliability targets using semi-probabilistic approaches. RBD is implemented in ASME BPVC HHA-3000 for nuclear graphite components, but is not specific to that application. There has been much confusion around the methods implemented in ASME BPVC HHA-3000 for qualifying nuclear graphite components. To address the confusion, this paper takes a hierarchical approach. First, the general RBD framework is presented. Then, the semi-probabilistic methods and the underlying assumptions implemented in the assessments are presented. The semi-probabilistic methods are separated from the engineering modifications that have been made to the assessments. After building the framework and underlying assumptions, the specific methods in the full and simplified assessments are explained in three steps: inputs, methods, outputs. The methods are applied to an H-451 reflector block. Tensile strength properties for other graphite grades are provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes

This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large-scale sampling-based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right-hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. Here, we demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallel forward model problem solver, for large-scale 3D MLMC simulations with up to 1.9·109 unknowns.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Phase III Study of Radiation Therapy With or Without Cis-Platinum in Patients With Unresectable Squamous or Undifferentiated Carcinoma of the Head and Neck: An Intergroup Trial of the Eastern Cooperative Oncology Group (E2382)

Purpose: The Head and Neck Intergroup conducted a Phase III randomized trial to determine whether the addition weekly cisplatin to daily radiation therapy (RT) would improve survival in patients with unresectable squamous cell head-and-neck carcinoma. Methods and Materials: Eligible patients were randomized to RT (70 Gy at 1.8-2 Gy/day) or to the identical RT with weekly cisplatin dosed at 20 mg/m{sup 2}. Failure-free survival (FFS) and overall survival (OS) curves were estimated with the Kaplan-Meier method and compared with the log rank test. Results: Between 1982 and 1987, 371 patients were accrued, and 308 patients were found eligible for analysis. Median follow-up was 62 months. The median FFS was 6.5 and 7.2 months for the RT and RT + cisplatin groups, respectively (p = 0.30). The p value for the treatment difference was p = 0.096 in multivariate modeling of FFS (compared to a p = 0.30 in univariate analysis). Expected acute toxicities were significantly increased with the addition of cisplatin except for in-field RT toxicities. Late toxicities were not significantly different except for significantly more esophageal (9% vs. 3%, p = 0.03) and laryngeal (11% vs. 4%, p = 0.05) late toxicities in the RT + cisplatin group. Conclusion: The addition of concurrent weekly cisplatin at 20 mg/m{sup 2} to daily radiation did not improve survival, although there was evidence of activity. Low-dose weekly cisplatin seems to have modest tumor radiosensitization but can increase the risk of late swallowing complications.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

From micro-correlations to macro-correlations

Random vectors with a symmetric correlation structure share a common value of pair-wise correlation between their different components. The symmetric correlation structure appears in a multitude of settings, e.g. mixture models. In a mixture model the components of the random vector are drawn independently from a general probability distribution that is determined by an underlying parameter, and the parameter itself is randomized. In this paper we study the overall correlation of high-dimensional random vectors with a symmetric correlation structure. Considering such a random vector, and terming its pair-wise correlation “micro-correlation”, we use an asymptotic analysis to derive the random vector’s “macro-correlation” : a score that takes values in the unit interval, and that quantifies the random vector’s overall correlation. The method of obtaining macro-correlations from micro-correlations is then applied to a diverse collection of frameworks that demonstrate the method’s wide applicability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nonlinear dynamics and quantum chaos of a family of kicked p -spin models

Herein we introduce kicked p-spin models describing a family of transverse Ising-like models for an ensemble of spin-1/2 particles with all-to-all p-body interaction terms occurring periodically in time as delta-kicks. This is the natural generalization of the well-studied quantum kicked top (p = 2) [Haake, Kus', and Scharf, Z. Phys. B 65, 381 (1987)]. We fully characterize the classical nonlinear dynamics of these models, including the transition to global Hamiltonian chaos. The classical analysis allows us to build a classification for this family of models, distinguishing between p = 2 and p > 2, and between models with odd and even p's. Quantum chaos in these models is characterized in both kinematic and dynamic signatures. For the latter, we show numerically that the growth rate of the out-of-time-order correlator is dictated by the classical Lyapunov exponent. Finally, we argue that the classification of these models constructed in the classical system applies to the quantum system as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Single Cell RNA-Seq and Machine Learning Reveal Novel Subpopulations in Low-Grade Inflammatory Monocytes With Unique Regulatory Circuits

Subclinical doses of LPS (SD-LPS) are known to cause low-grade inflammatory activation of monocytes, which could lead to inflammatory diseases including atherosclerosis and metabolic syndrome. Sodium 4-phenylbutyrate is a potential therapeutic compound which can reduce the inflammation caused by SD-LPS. To understand the gene regulatory networks of these processes, we have generated scRNA-seq data from mouse monocytes treated with these compounds and identified 11 novel cell clusters. We have developed a machine learning method to integrate scRNA-seq, ATAC-seq, and binding motifs to characterize gene regulatory networks underlying these cell clusters. Using guided regularized random forest and feature selection, our method achieved high performance and outperformed a traditional enrichment-based method in selecting candidate regulatory genes. Our method is particularly efficient in selecting a few candidate genes to explain observed expression pattern. In particular, among 531 candidate TFs, our method achieves an auROC of 0.961 with only 10 motifs. Finally, we found two novel subpopulations of monocyte cells in response to SD-LPS and we confirmed our analysis using independent flow cytometry experiments. Our results suggest that our new machine learning method can select candidate regulatory genes as potential targets for developing new therapeutics against low grade inflammation.

60 APPLIED LIFE SCIENCES↗