Search NASA⌕ Search

SEARCH · Search NASA

Results for “statistical visualization”

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 379 records · Page 21

Near Real-time Data Analysis of Core-collapse Supernova Simulations with Bellerophon

We present an overview of a software system, Bellerophon, built to support a production-level HPC application called CHIMERA, which simulates core-collapse supernova events at the petascale. Developed over the last four years, Bellerophon enables CHIMERA's geographically dispersed team of collaborators to perform data analysis in near real-time. Its n-tier architecture provides an encapsulated, end-to-end software solution that enables the CHIMERA team to quickly and easily access highly customizable animated and static views of results from anywhere in the world via a web- deliverable, cross-platform desktop application. In addition, Bellerophon addresses software engineering tasks for the CHIMERA team by providing an automated mechanism for performing regression testing on a variety of supercomputing platforms. Elements of the team's workflow management needs are met with software tools that dynamically generate code repository statistics, access important online resources, and monitor the current status of several supercomputing resources.

97 MATHEMATICS AND COMPUTING↗

Radiation Dose–Dependent Hippocampal Atrophy Detected With Longitudinal Volumetric Magnetic Resonance Imaging

Purpose: After radiation therapy (RT) to the brain, patients often experience memory impairment, which may be partially mediated by damage to the hippocampus. Hippocampal sparing in RT planning is the subject of recent and ongoing clinical trials. Calculating appropriate hippocampal dose constraints would be improved by efficient in vivo measurements of hippocampal damage. In this study we sought to determine whether brain RT was associated with dose-dependent hippocampal atrophy. Methods and Materials: Hippocampal volume was measured with magnetic resonance imaging (MRI) in 52 patients who underwent fractionated, partial brain RT for primary brain tumors. Study patients had high-resolution, 3-dimensional volumetric MRI before and 1 year after RT. Images were processed using software with clearance from the US Food and Drug Administration and Conformité Européene marking for automated measurement of hippocampal volume. Automated results were inspected visually for accuracy. Tumor and surgical changes were censored. Mean hippocampal dose was tested for correlation with hippocampal atrophy 1 year after RT. Average hippocampal volume change was also calculated for hippocampi receiving high (>40 Gy) or low (<10 Gy) mean RT dose. A multivariate analysis was conducted with linear mixed-effects modeling to evaluate other potential predictors of hippocampal volume change, including patient (random effect), age, hemisphere, sex, seizure history, and baseline volume. Statistical significance was evaluated at α = 0.05. Results: Mean hippocampal dose was significantly correlated with hippocampal volume loss (r=−0.24, P=.03). Mean hippocampal volume was significantly reduced 1 year after high-dose RT (mean −6%, P=.009) but not after low-dose RT. In multivariate analysis, both RT dose and patient age were significant predictors of hippocampal atrophy (P<.01). Conclusions: The hippocampus demonstrates radiation dose–dependent atrophy after treatment for brain tumors. Quantitative MRI is a noninvasive imaging technique capable of measuring radiation effects on intracranial structures. This technique could be investigated as a potential biomarker for development of reliable dose constraints for improved cognitive outcomes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An Extensive Catalog of Early-type Dwarf Galaxies in the Local Universe: Morphology and Environment

We present an extensive catalog of 5405 early-type dwarf (dE) galaxies located in the various environments, i.e., clusters, groups, and fields, of the local universe (z < 0.01). The dEs are selected through visual inspection of the Legacy survey's g–r–z combined tricolor images. The inspected area, covering a total sky area of 7643 deg 2 , encompasses two local clusters (Virgo and Fornax), 265 groups, and the regions around 586 field galaxies of M K < -21 mag. The catalog aims to be one of the most extensive and publicly accessible collections of data on dEs, despite its complex completeness limits, which may not accurately represent its statistical completeness. The strength of the catalog lies in the morphological characteristics, including nucleated, tidal, and ultradiffuse dEs. The two clusters contribute nearly half (2437 out of 5405) of the dEs, and the 265 groups contribute 2103 dEs. There are 864 dEs in 586 fields, i.e., ~1.47 dEs per field. Using a standard definition commonly used in literature, we identify 100 ultradiffuse galaxies (UDGs), which compose ~2% of the dE population. We find that 40% of our sample dEs harbor a central nucleus, and among the UDG population, a majority (79%) are nonnucleated. About 1.3% of dEs suffer from ongoing tidal disturbance by nearby massive galaxies, and only 0.03% show a sign of recent dwarf–dwarf mergers. The association between dEs and their nearest bright neighbor galaxies suggests that dEs are more likely created where their neighbors are non-star-forming ones.

79 ASTRONOMY AND ASTROPHYSICS↗

MALDI Orbitrap Mass Spectrometry Profiling of Dysregulated Sulfoglycosphingolipids in Renal Cell Carcinoma Tissues

Matrix-assisted laser desorption/ionization coupled with Orbitrap mass spectrometry (MALDI-Orbitrap-MS) is used for the clinical study of patients with renal cell carcinoma (RCC), as the most common type of kidney cancer. Significant changes in sulfoglycosphingolipid abundances between tumor and autologous normal kidney tissues are observed. First, sulfoglycosphingolipid species in studied RCC samples are identified using high mass accuracy full scan and tandem mass spectra. Subsequently, optimization, method validation, and statistical evaluation of MALDI-MS data for 158 tissues of 80 patients are discussed. More than 120 sulfoglycosphingolipids containing one to five hexosyl units are identified in human RCC samples based on the systematic study of their fragmentation behavior. Many of them are recorded here for the first time. Multivariate data analysis (MDA) methods, i.e., unsupervised principal component analysis (PCA) and supervised orthogonal partial least square discriminant analysis (OPLS-DA), are used for the visualization of differences between normal and tumor samples to reveal the most up- and downregulated lipids in tumor tissues. Obtained results are closely correlated with MALDI mass spectrometry imaging (MSI) and histologic staining. Important steps of the present MALDI-Orbitrap-MS approach are also discussed, such as the selection of best matrix, correct normalization, validation for semiquantitative study, and problems with possible isobaric interferences on closed masses in full scan mass spectra. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Labile Fe(III) supersaturation controls nucleation and properties of product phases from Fe(II)-catalyzed ferrihydrite transformation

We report Fe(II)-catalyzed ferrihydrite (Fh) transformation to more crystalline iron (oxyhydr)oxide phases is a widely occurring geochemical process which has been extensively studied as a function of Fe(II)/Fh ratios at fixed Fh loadings. However, recent isolation of an intermediate Fe(III) species resulting from Fe(II)-Fh contact that facilitates transformation by dissolution/reprecipitation suggests that the kinetics and properties of product phases will instead depend mostly on its rate of accumulation to a critical concentration, consistent with principles in the classical nucleation theory (CNT). This suggests a dependence both on the loading of Fe(II) on the surface, which controls the rate of labile Fe(III) formation, as well as the available volume of solution, which also impacts how fast it can achieve its critical concentration to nucleate product phases. To specifically examine the latter effect, here we studied transformation of 15 mg Fh in 1 mM FeSO 4 solutions at pH 7.2 in batch suspensions of 30 mL, 150 mL, and 450 mL volumes. Time-dependent concentrations of aqueous Fe(II), surface-associated Fe(II), and resulting labile Fe(III) were monitored along with bulk solids characterization as a function of time. Transmission electron microscopy (TEM) was used to visualize the evolution of phases at identical locations on TEM grids. The collective results show that the rates of Fh loss and emergence of product lepidocrocite (Lp) and goethite (Gt) as well as their phase proportions, nucleation mode and morphological properties depend directly on the rate of accumulation of the labile Fe(III) precursor to its critical concentration, which in our experiments was controlled simply by varying the available volume of solution into which it enters. Statistical analyses of TEM image data suggest that while both heterogeneous and homogeneous nucleation occurred in all experiments, the former was increasingly favored at lower Fh/solution ratio due to its lower nucleation barrier being more favorable at attendant lower supersaturations of Fe(III).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ARCH: Large-scale knowledge graph via aggregated narrative codified health records analysis

Objective: Electronic health record (EHR) systems contain a wealth of clinical data stored as both codified data and free-text narrative notes (NLP). The complexity of EHR presents challenges in feature representation, information extraction, and uncertainty quantification. Here, to address these challenges, we proposed an efficient Aggregated naRrative Codified Health (ARCH) records analysis to generate a large-scale knowledge graph (KG) for a comprehensive set of EHR codified and narrative features. Methods: Using data from 12.5 million Veterans Affairs patients, ARCH first derives embedding vectors and generates similarities along with associated p-values to measure the strength of relatedness between clinical features with statistical certainty quantification. Next, ARCH performs a sparse embedding regression to remove indirect linkage between features to build a sparse KG. Finally, ARCH was validated on various clinical tasks, including detecting known relationships between entity pairs, predicting drug side effects, disease phenotyping, as well as sub-typing Alzheimer’s disease patients. Results: ARCH produces high-quality clinical embeddings and KG for over 60,000 codified and narrative EHR concepts. The KG and embeddings are visualized in the R-shiny powered web-API.3 ARCH achieved high accuracy in detecting EHR concept relationships, with AUCs of 0.926 (codified) and 0.861 (NLP) for similar EHR concepts, and 0.810 (codified) and 0.843 (NLP) for related pairs. It detected drug side effects with a 0.723 AUC, which improved to 0.826 after fine-tuning. Using both codified and NLP features, the detection power increased significantly. Compared to other methods, ARCH has superior accuracy and enhances weakly supervised phenotyping algorithms’ performance. Notably, it successfully categorized Alzheimer’s patients into two subgroups with varying mortality rates. Conclusion: The proposed ARCH algorithm generates large-scale high-quality semantic representations and knowledge graph for both codified and NLP EHR features, useful for a wide range of predictive modeling tasks.

Electronic health records↗

Dynamics of unstably stratified free shear flows: an experimental investigation of coupled Kelvin–Helmholtz and Rayleigh–Taylor instability

The dynamics of the coupled Kelvin–Helmholtz (KH) and Rayleigh–Taylor (RT) instability (referred to as KHRT instability or KHRTI) is studied using statistically steady experiments performed in a multi-layer gas tunnel. Experiments are performed at four density ratios ranging in Atwood number $A_{t}$ from 0.035 to 0.159, with varying amounts of shear and $\unicode[STIX]{x0394}U/U$ ranging from 0 to 0.48, where $\unicode[STIX]{x0394}U$ is the speed difference between the two flow streams being investigated and $U$ is the mean velocity of these two streams. Three types of diagnostics – back-lit visualization, hot-wire anemometry and particle image velocimetry (PIV) – are employed to obtain the mixing widths, velocity field and density field. The flow is found to be governed by KH dynamics at early times and RT dynamics at late times. This transition from KH-instability-like to RT-instability-like behaviour is quantified using the Richardson number. Transitional Richardson number magnitudes obtained for the present KHRT flows are found to be in the range 0.17–0.56 similar to the critical Richardson numbers for stably stratified free shear flows. Comparing the evolution of density and velocity mixing widths, the density mixing layer is found to be approximately two times as thick as the velocity mixing layer. Scaling of velocity fluctuations is attempted using combinations of KH and RT scales. It is found that the proposed KHRT velocity scale, obtained using the combined mixing-layer growth equation, is appropriate for intermediate stages of the flow when both KH and RT dynamics are comparable. Probability density functions (p.d.f.s) for different fluctuating quantities are presented. Multiple peaks in p.d.f.s are qualitatively explained from the development of coherent KH roll-ups and their subsequent transition into turbulent pockets. The evolution of energy spectra indicates that density fluctuations start to show an inertial subrange from earlier times compared to velocity fluctuations. The spectra exhibit a slightly steeper slope than the Kolmogorov–Obukhov five-thirds law.

Mechanics↗

Multivariate Regression of Pyrotechnic Igniter Output

Multivariate multiple regression models are applied to simplified pyrotechnic igniters for the first time to understand how changes in manufactured parameters can affect the output gas dynamic response and the timing of ignition events. The statistical modeling technique is applied to demonstrate quantification of the effects of a set of independent variables measured in the as-fabricated igniters on a set of responses experimentally measured from the functioned igniters. Two independent process variables were intentionally varied following a full factorial experimental design while several other independent variables varied within their normal manufacturing variability range. The four igniter performance responses consisted of the timing of sequential events during igniter function and visual gas dynamic output in the form of shock wave strength observed with high-speed schlieren imaging. Linear regression models built using the measurements throughout the manufacturing processes and the output variance provide insight into the critical device parameters that dominate performance

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Multiscale Cryo Electron Microscopy Reveals Interfacial Degradation and Stabilization in Battery Electrodes

Electrochemical interfaces are dynamic systems, evolving based on their local environment and reactant surface structures. The electrode-electrolyte interface in Li-ion batteries can be protective, limiting parasitic reactions with the electrolyte to passivate the surface [1]. Additionally, this interphase has an impact on the Li-ion transport through that layer based on its composition, bonding environment, and thickness. These parameters are challenging to collect and may vary depending on the electrode surface site investigated relative to its spatial position in a coin cell. This study will detail a multiscale cryogenic electron microscopy approach where millimeter-scale cross-sections through the coin cell batteries were made using a cryogenic stage within a fs-laser plasma focused ion beam (laser PFIB) with complementary energy dispersive X-ray spectroscopy able to detect variations in the composition at electrode interfaces [2]. Microscale cross-sectioning and lamella sample preparation of battery electrodes was conducted at the Center for Integrated Nanotechnologies using a Ga-ion focused ion beam (FIB) with air-free and cryo-transfer [3], followed by nanoscale mapping of composition and bonding within the CEI through cryo-scanning transmission electron microscopy (cryo-STEM) electron energy loss spectroscopy [4]. This multiscale approach enabled identification of millimeter-scale features of a battery stack with visualization of degradation in electrodes such as cracks in cathode particles, gas evolution, and SEI evolution; microscale interfacial characteristics, such as heterogeneity in the SEI or barrier layer and identification of electrolyte networks to the electrode surfaces; and nanoscale measurement of the CEI thickness, mapping of transition metal bonding within the cathode particles to identify loss of active materials, and identification of beneficial electrolyte additives incorporated into the CEI structure. This multiscale approach allows for a statistical understanding of the primary mechanisms and parasitic degradation pathways that impact performance by limiting the ion transport pathways within Li+ batteries.

36 MATERIALS SCIENCE↗

Strain fields in twisted bilayer graphene

Van der Waals heteroepitaxy allows deterministic control over lattice mismatch or azimuthal orientation between atomic layers to produce long-wavelength superlattices. The resulting electronic phases depend critically on the superlattice periodicity and localized structural deformations that introduce disorder and strain. In this study we used Bragg interferometry to capture atomic displacement fields in twisted bilayer graphene with twist angles < 2°. Nanoscale spatial fluctuations in twist angle and uniaxial heterostrain were statistically evaluated, revealing the prevalence of short-range disorder in moiré heterostructures. By quantitatively mapping strain tensor fields, we uncovered two regimes of structural relaxation and disentangled the electronic contributions of constituent rotation modes. Further, we found that applied heterostrain accumulates anisotropically in saddle-point regions, generating distinctive striped strain phases. Our results establish the reconstruction mechanics underpinning the twist-angle-dependent electronic behaviour of twisted bilayer graphene and provide a framework for directly visualizing structural relaxation, disorder and strain in moiré materials.

36 MATERIALS SCIENCE↗

The SRT reconstruction algorithm for semiquantification in PET imaging

Purpose: The spline reconstruction technique (SRT) is a new, fast algorithm based on a novel numerical implementation of an analytic representation of the inverse Radon transform. The mathematical details of this algorithm and comparisons with filtered backprojection were presented earlier in the literature. In this study, the authors present a comparison between SRT and the ordered-subsets expectation–maximization (OSEM) algorithm for determining contrast and semiquantitative indices of {sup 18}F-FDG uptake. Methods: The authors implemented SRT in the software for tomographic image reconstruction (STIR) open-source platform and evaluated this technique using simulated and real sinograms obtained from the GE Discovery ST positron emission tomography/computer tomography scanner. All simulations and reconstructions were performed in STIR. For OSEM, the authors used the clinical protocol of their scanner, namely, 21 subsets and two iterations. The authors also examined images at one, four, six, and ten iterations. For the simulation studies, the authors analyzed an image-quality phantom with cold and hot lesions. Two different versions of the phantom were employed at two different hot-sphere lesion-to-background ratios (LBRs), namely, 2:1 and 4:1. For each noiseless sinogram, 20 Poisson realizations were created at five different noise levels. In addition to making visual comparisons of the reconstructed images, the authors determined contrast and bias as a function of the background image roughness (IR). For the real-data studies, sinograms of an image-quality phantom simulating the human torso were employed. The authors determined contrast and LBR as a function of the background IR. Finally, the authors present plots of contrast as a function of IR after smoothing each reconstructed image with Gaussian filters of six different sizes. Statistical significance was determined by employing the Wilcoxon rank-sum test. Results: In both simulated and real studies, SRT exhibits higher contrast and lower bias than OSEM at the cold lesions. This improvement is achieved at the expense of increasing the noise in the reconstructed images. For the hot lesions, SRT exhibits a small improvement in contrast and LBR over OSEM with 21 subsets and two iterations; however, this improvement is not statistically significant. As the number of iterations increases, the performance of OSEM improves over SRT but again without statistical significance. The curves of contrast and LBR as a function of IR after Gaussian blurring indicate that the advantage of SRT in the cold regions is maintained even after decreasing the noise level by Gaussian blurring. Conclusions: SRT, at the expense of slightly increased noise in the reconstructed images, reconstructs images of higher contrast and lower bias than the clinical protocol of OSEM. This improvement is particularly evident for images involving cold regions. Thus, it appears that SRT should be particularly useful for the quantification of low-count and cold regions.

60 APPLIED LIFE SCIENCES↗

Gamma-ray imaging and holdup assays of 235-F PuFF cells 1 & 2

Savannah River National Laboratory (SRNL) Nuclear Measurements (L4120) was tasked with performing enhanced characterization of the holdup in the PuFF shielded cells. Assays were performed in accordance with L16.1-ADS-2460 using two high-resolution gamma-ray detectors. The first detector, an In Situ Object Counting System (ISOCS)-characterized detector, was used in conjunction with the ISOCS Geometry Composer software to quantify grams of holdup. The second detector, a Germanium Gamma-ray Imager (GeGI), was used to visualize the location and relative intensity of the holdup in the cells. Carts and collimators were specially designed to perform optimum assays of the cells. Thick, pencil-beam tungsten collimators were fabricated to allow for extremely precise targeting of items of interest inside the cells. Carts were designed with a wide range of motion to position and align the detectors. A total of 24 measurements were made, each typically 24 hours or longer to provide sufficient statistical precision. This report presents the results of the enhanced characterization for cells 1 and 2. The measured gram values agree very well with results from the 2014 study. In addition, images were created using both the 2014 data and the new GeGI data. The GeGI images of the cells walls reveal significant Pu-238 holdup on the surface of the walls in cells 1 and 2. Additionally, holdup is visible in the two pass-throughs from cell 1 to the wing cabinets. This report documents the final element (exterior measurements coupled with gamma-ray imaging and modeling) of the enhanced characterization of cells 1-5 (East Cell Line).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Overview of the Tolerance Limit Calculations with Application to TSURFER

To establish confidence in the results of computerized physics models, a key regulatory requirement is to develop a scientifically defendable process. The methods employed for confidence, characterization, and consolidation, or C3, are statistically involved and are often accessible only to avid statisticians. This manuscript serves as a pedagogical presentation of the C3 process to all stakeholders—including researchers, industrial practitioners, and regulators—to impart an intuitive understanding of the key concepts and mathematical methods entailed by C3. The primary focus is on calculation of tolerance limits, which is the overall goal of the C3 process. Tolerance limits encode the confidence in the calculation results as communicated to the regulator. Understanding the C3 process is especially critical today, as the nuclear industry is considering more innovative ways to assess new technologies, including new reactor and fuel concepts, via an integrated approach that optimally combines modeling and simulation and minimal targeted validation experiments. This manuscript employs intuitive, analytical, numerical, and visual representations to explain how tolerance limits may be calculated for a wide range of configurations, and it also describes how their values may be interpreted. Various verification tests have been developed to test the calculated tolerance limits and to help delineate their values. The manuscript demonstrates the calculation of tolerance limits for TSURFER, a computer code developed by the Oak Ridge National Laboratory for criticality safety applications. The goal is to evaluate the tolerance limit for TSURFER-determined criticality biases to support the determination of upper, subcritical limits for regulatory purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

UrbanScaping: Community Spatial Data Visualization & Analytics

Evaluating the electrification potential of buildings through retrofitting is crucial for reducing carbon emissions and the carbon footprint of built environments. This study leverages the Automatic Building Energy Modeling (AutoBEM) software, integrating the Model America database to create an urban context-based spatial analysis platform for community engagement and development. We selected Camp Hill Borough, PA, as a case study to analyze building-specific energy performance and evaluate the electrification potential of each building by switching to different Heating, Ventilation, and Air Conditioning (HVAC) systems and measurement components. The simulation results generated by the workflow provide retrofitting suggestions to help mitigate the carbon footprint as well as energy saving statistics of buildings. Additionally, the developed web-based interface serves as a community engagement platform, allowing residents to provide feedback and further develop interactive communication protocols. The outcomes of this project offer a baseline for community electrification planning and contribute to the design of low-carbon communities.

Chowdhury, Shovan [ORNL]↗

Electronic Production Control System News Note

The CMD subprogram develops and deploys multi-application weapon component manufacturing capabilities needed to replace sunset technologies, upgrade existing technologies, and introduce new technologies that support the nuclear weapons stockpile. Sandia National Laboratories (SNL) is currently developing an Electronic Production Control System (EPCS) for ASIC production in the back end of line (BEOL) flow. A pre-production version of EPCS has been released and is operating as intended. We are currently processing all ASIC products, including CVV validation ASICs. The BEOL flow includes wafer probe, packaging, quality inspection, initial electrical test, dynamic burn-in, final electrical test, quality inspection, final visual inspection, radiation testing, SNL acceptance, DOE SSO acceptance, and qualification testing steps. EPCS will execute this production workflow (initiate, control, and electronically track each product); it will implement defect prevention strategies; it will create an electronic inventory control system that is accountable for all ASIC material from inception to delivery; it will electronically collect and store any and all required quality evidence and most of the electronic BEOL data; and it will create a management reporting module with production statistics and timelines. The defect prevention strategies include the use of 2D marks on our ASICs with 2D scanners so that there will be no manual entry of part numbers (scanners), will not allow testing the wrong parts in Test and Burn-in. Other defect prevention strategies will negate the possibility of neither using the wrong test or burn-in program for any part nor the use of BEOL machines or boards that are out of calibration. In summary, EPCS is being implemented and it will greatly improve the product acceptance through full electronic data collection and quality control and defect prevention through bar codes, scanners, tracking training of personnel, and tracking tool availability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The snapshot Hubble U-band cluster survey (SHUCS). II. The star cluster population of NGC 2997

We study the star cluster population of NGC 2997, a giant spiral galaxy located at 9.5 Mpc and targeted by the Snapshot Hubble U-band Cluster Survey (SHUCS). Combining our U-band imaging from SHUCS with archival BVI imaging from the Hubble Space Telescope, we select a high confidence sample of clusters in the circumnuclear ring and disk through a combination of automatic detection procedures and visual inspection. The cluster luminosity functions in all four filters can be approximated by power laws with indices of –1.7 to –2.3. Some deviations from pure power-law shape are observed, hinting at the presence of a high-mass truncation in the cluster mass function. However, upon inspection of the cluster mass function, we find it is consistent with a pure power law of index –2.2 ± 0.2 despite a slight bend at ∼2.5 × 10{sup 4} M {sub ☉}. No statistically significant truncation is observed. From the cluster age distributions, we find a low rate of disruption (ζ ∼ –0.1) in both the disk and circumnuclear ring. Finally, we estimate the cluster formation efficiency (Γ) over the last 100 Myr in each region, finding 7% ± 2% for the disk, 12% ± 4% for the circumnuclear ring, and 10% ± 3% for the entire UBVI footprint. This study highlights the need for wide-field UBVI coverage of galaxies to study cluster populations in detail, though a small sample of clusters can provide significant insight into the characteristics of the population.

79 ASTRONOMY AND ASTROPHYSICS↗

Predictive Big Data Analytics: A Study of Parkinson’s Disease Using Large, Complex, Heterogeneous, Incongruent, Multi-Source and Incomplete Observations

Background A unique archive of Big Data on Parkinson’s Disease is collected, managed and disseminated by the Parkinson’s Progression Markers Initiative (PPMI). The integration of such complex and heterogeneous Big Data from multiple sources offers unparalleled opportunities to study the early stages of prevalent neurodegenerative processes, track their progression and quickly identify the efficacies of alternative treatments. Many previous human and animal studies have examined the relationship of Parkinson’s disease (PD) risk to trauma, genetics, environment, co-morbidities, or life style. The defining characteristics of Big Data–large size, incongruency, incompleteness, complexity, multiplicity of scales, and heterogeneity of information-generating sources–all pose challenges to the classical techniques for data management, processing, visualization and interpretation. We propose, implement, test and validate complementary model-based and model-free approaches for PD classification and prediction. To explore PD risk using Big Data methodology, we jointly processed complex PPMI imaging, genetics, clinical and demographic data. Methods and Findings Collective representation of the multi-source data facilitates the aggregation and harmonization of complex data elements. This enables joint modeling of the complete data, leading to the development of Big Data analytics, predictive synthesis, and statistical validation. Using heterogeneous PPMI data, we developed a comprehensive protocol for end-to-end data characterization, manipulation, processing, cleaning, analysis and validation. Specifically, we (i) introduce methods for rebalancing imbalanced cohorts, (ii) utilize a wide spectrum of classification methods to generate consistent and powerful phenotypic predictions, and (iii) generate reproducible machine-learning based classification that enables the reporting of model parameters and diagnostic forecasting based on new data. We evaluated several complementary model-based predictive approaches, which failed to generate accurate and reliable diagnostic predictions. However, the results of several machine-learning based classification methods indicated significant power to predict Parkinson’s disease in the PPMI subjects (consistent accuracy, sensitivity, and specificity exceeding 96%, confirmed using statistical n-fold cross-validation). Clinical (e.g., Unified Parkinson's Disease Rating Scale (UPDRS) scores), demographic (e.g., age), genetics (e.g., rs34637584, chr12), and derived neuroimaging biomarker (e.g., cerebellum shape index) data all contributed to the predictive analytics and diagnostic forecasting. Conclusions Model-free Big Data machine learning-based classification methods (e.g., adaptive boosting, support vector machines) can outperform model-based techniques in terms of predictive precision and reliability (e.g., forecasting patient diagnosis). We observed that statistical rebalancing of cohort sizes yields better discrimination of group differences, specifically for predictive analytics based on heterogeneous and incomplete PPMI data. UPDRS scores play a critical role in predicting diagnosis, which is expected based on the clinical definition of Parkinson’s disease. Even without longitudinal UPDRS data, however, the accuracy of model-free machine learning based classification is over 80%. The methods, software and protocols developed here are openly shared and can be employed to study other neurodegenerative disorders (e.g., Alzheimer’s, Huntington’s, amyotrophic lateral sclerosis), as well as for other predictive Big Data analytics applications.

Science & Technology - Other Topics↗

A Catalog of Merging Dwarf Galaxies in the Local Universe

We introduce the largest publicly available catalog of interacting dwarf galaxies. It includes 177 nearby merging dwarf galaxies of stellar mass M * < 10 10 M ⊙ and redshifts z < 0.02. These galaxies are selected by visual inspection of publicly available archival imaging from two wide-field optical surveys (SDSS-III and the Legacy Survey), and they possess low-surface-brightness features that are likely the result of an interaction between dwarf galaxies. We list UV and optical photometric data that we use to estimate stellar masses and star formation rates. Currently, the study of interacting dwarf galaxies has largely been done on an individual basis, and lacks a sufficiently large catalog to give statistics on the properties of interacting dwarf galaxies, and their role in the evolution of low-mass galaxies. We anticipate that this public catalog can be used as a reference sample to investigate the effects of the tidal interaction on the evolution of star formation, and the morphology/structure of dwarf galaxies. Our sample is overwhelmingly dominated by star-forming galaxies, and they are generally found significantly below the red sequence in the color–magnitude relation. The number of early-type galaxies is only 3 out of 177. We classify them, according to observed low-surface-brightness features, into various categories including shells, stellar streams, loops, antennae, or simply interacting. We determine that dwarf–dwarf interactions tend to prefer the low-density environment. Only 41 out of the 177 candidate dwarf–dwarf interaction systems have giant neighbors within a sky-projected distance of 700 kpc and a line-of-sight radial velocity range ±700 km s –1 , and compared to the LMC–SMC, they are generally located at much larger sky-projected distances from their nearest giant neighbors.

79 ASTRONOMY AND ASTROPHYSICS↗