Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling statistics”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Robust Cross-correlation-based Measurement of Clump Sizes in Galaxies

Stars form in molecular complexes that are visible as giant clouds (∼10{sup 5−6} M{sub ⊙}) in nearby galaxies and as giant clumps (∼10{sup 8−9} M{sub ⊙}) in galaxies at redshifts z≈1–3. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g., Hα) can be recovered statistically using the two-point correlation function (2PCF) if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at ∼20% accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical measurements of Hα clump size agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Representative Day Selection Using Statistical Bootstrapping for Accelerating Annual Distribution Simulations

Capturing technical and economic impacts of solar photovoltaics (PV) and other distributed energy resources (DERs) on electric distribution systems can require high-time resolution (e.g. 1 minute), long-duration (e.g. 1 year) simulations. However, such simulations can be computationally prohibitive, particularly when including complex control schemes in quasi-steady-state time series (QSTS) simulation. Various approaches have been used in the literature to down select representative time segments (e.g. days), but typically these are best suited for lower time resolutions or consider only a single data stream (e.g. PV production) for selection. We present a statistical approach that combines stratified sampling and bootstrapping to select representative days while also providing a simple method to reassemble annual results. We describe the approach in the context of a recent study with a utility partner. This approach enables much faster QSTS analysis by simulating only a subset of days, while maintaining accurate annual estimates.

14 SOLAR ENERGY↗
DOE OSTI2015

Adaptive Bayesian Inference for Prediction

Performing rigorous predictive simulation necessitates that all available data be incorporated into the modeling process in order to accurately quantify uncertainties in predictions of the model. Statistical Bayesian inference is the most commonly used approach for incorporating data into uncertainty estimates. Statistical Bayesian inference uses data and an assumed error model to inform posterior distributions of model input variables and model discrepancies. An explicit characterization of the posterior distribution is not formed, but rather methods such as Markov chain Monte-Carlo (MCMC) are used to draw samples from the posterior. These samples are then used to evaluate statistics, such as integrals, of predicted quantities of interest (QoI). Recently the authors have developed an alternative reformulation of Bayesian inference based upon a measure-theoretic inversion technique. Unlike traditional statistical Bayesian inference, this reformulation allows us to directly interrogate the posterior distribution and, consequently, to generate samples from the posterior without MCMC. In this work we use forward propagation uncertainty quantification techniques to reduce the cost of approximating the posterior distribution. Specifically we investigate goal-oriented adaptive sparse grid approximation methods. These adaptive methods concentrate high-fidelity model evaluations in regions of the parameter space that both significantly contribute to the uncertainty in specified QoI and are informed by the available data. Numerical results are presented to compare our approach with the standard statistical Bayesian approach and to demonstrate the effectiveness of the goal-oriented adaptivity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Simulation Study Using a New Type of Sample Variance

We evaluate with simulated data a new type of sample variance for the characterization of frequency stability. The new statistic (referred to as TOTALVAR and its square root TOTALDEV) is a better predictor of long-term frequency variations than the present sample Allan deviation. The statistical model uses the assumption that a time series of phase or frequency differences is wrapped (periodic) with overall frequency difference removed. We find that the variability at long averaging times is reduced considerably for the five models of power-law noise commonly encountered with frequency standards and oscillators.

Howe, D. A.↗
NASA NTRS1981

Information efficiency of line-scan imaging mechanisms

Information theory is used to formulate a single figure of merit for assessing the performance of line-scan imaging systems as a function of their spatial response (PSF or MTF), sensitivity, and sampling and quantization intervals and of the statistical properties of a random radiance field. Information density and efficiency (i.e., the ratio of information density to data density) tend to be optimum when the MTF and sampling passband of the imaging system are matched to the Wiener spectrum of the radiance field. Computational results for the statistical properties of natural radiance fields and the responses of common line-scan imaging mechanisms indicate that information density and efficiency are not strongly sensitive to variations in typical statistical properties of the radiance field and that the best practically realizable performance is approached when the sampling intervals are approximately 0.5-0.7 times the equivalent diameter of the PSF.

Huck, F. O.↗
NASA NTRS1995

A Simulated Geochemical Rover Mission to the Taurus-Littrow Valley of the Moon

We test the effectiveness of using an alpha backscatter, alpha-proton, X ray spectrometer on a remotely operated rover to analyze soils and provide geologically useful information about the Moon during a simulated mission to a hypothetical site resembling the Apollo 17 landing site. On the mission, 100 soil samples are "analyzed" for major elements at moderate analytical precision (e.g., typical relative sample standard deviation from counting statistics: Si[11%], Al[18%], Fe[6%], Mg[20%], Ca[5%]). Simulated compositions of soils are generated by combining compositions of components representing the major lithologies occurring at the site in known proportions. Simulated analyses are generated by degrading the simulated compositions according to the expected analytical precision of the analyzer. Compositions obtained from the simulated analyses are modeled by least squares mass balance as mixtures of the components, and the relative proportions of those components as predicted by the model are compared with the actual proportions used to generate the simulated composition. Boundary conditions of the modeling exercise are that all important lithologic components of the regolith are known and are represented by model components, and that the compositions of these components are well known. The effect of having the capability of determining one incompatible element at moderate precision (25%) is compared with the effect of the lack of this capability. We discuss likely limitations and ambiguities that would be encountered, but conclude that much of our knowledge about the Apollo 17 site (based on the return samples) regarding the distribution and relative abundances of lithologies in the regolith could be obtained. This success requires, however, that at least one incompatible element be determined.

Korotev, Randy L.↗
NASA NTRS2000

Statistical Study of Auroral Kilometric Radiation Fine Structure Striations Observed by Polar

We have conducted a statistical survey of a semirandom sample of the auroral kilometric radiation (AKR) data observed by the plasma wave instrument wideband receiver on board the Polar spacecraft. We have determined that AKR fine structure patterns with very narrowband, negative drifting striations occur in approximately 6% of the high-resolution wideband spectrograms when AKR is present. Positive sloping striations are also observed, but at a much lower rate. More than 8200 AKR stripes have been scaled. The stripes are predominantly found in the 40 to 215-kHz frequency range and have a frequency extent of about 4 kHz and a duration of usually less than 2 s. The majority of the stripes have drift rates between -8 and -2 kHz/s, with a peak in the distribution between -6 and -4 kHz/s. There is also a much smaller group of striations with positive drift rates of up to about 5 or 6 kHz/s. We have further investigated the change of drift rate with frequency. Almost all striations are observed in the lowest two frequency bands of the wideband receiver (f < 215 kHz). There is an increase in the statistical drift rate with increasing frequency. The statistical slope of the striations increases with frequency from about -4.4 kHz/s at 75 kHz to about -5.7 kHz/s at 170 kHz. This frequency dependence of the drift rate is consistent, under certain conditions, with a production mechanism stimulated by an upward propagating electromagnetic ion cyclotron wave, as had been suggested earlier. However, such a changing drift rate is also compatible with a stimulated source region that propagates upward along the magnetic field line at the velocity of an ion beam accelerated by a local, upward directed electric field, as is typically observed in the auroral region. An explanation for this association is not apparent at this time.

Menietti, J. D.↗
DOE OSTI2021

A method for predicting failure statistics for steady state elevated temperature structural components

This paper presents the initial development of a high temperature life prediction method that accounts for the variability in the material properties of Grade 91 steel. The method accounts for material variability by fitting a variable 3-parameter Weibull distribution to experimental rupture data and accounts for the variability of creep deformation on the steady-state stresses via a Monte Carlo approach. To ensure reasonable computational times, the model represents the material as an extremely viscous Stokes fluid with a non-Newtonian viscosity, therefore solving the stress relaxation problem with a steady, static, instead of transient, analysis. Furthermore, the complete statistical analysis combines this model for creep deformation with a probabilistic model for creep rupture to evaluate the probability of premature failure for a set of sample problems, comparing the predicted failure statistics to the design life predicted by the ASME Boiler and Pressure Vessel Code rules.

42 ENGINEERING↗
DOE OSTI2019

GMRT 610 MHz observations of galaxy clusters in the ACT equatorial sample

Abstract We present Giant Metrewave Radio Telescope (GMRT) 610 MHz observations of 14 Atacama Cosmology Telescope (ACT) clusters, including new data for nine. The sample includes 73 per cent of ACT equatorial clusters with M500 > 5 × 1014 M⊙. We detect diffuse emission in three of these (27$^{+20}_{-14}$ per cent): we detect a radio minihalo in ACT-CL J0022.2–0036 at z = 0.8, making it the highest redshift minihalo known; we detect potential radio relic emission in ACT-CL J0014.9–0057 (z = 0.533); and we confirm the presence of a radio halo in low-mass cluster ACT-CL J0256.5+0006, with flux density S610 = 6.3 ± 0.4 mJy. We also detect residual diffuse emission in ACT-CL J0045.9–0152 (z = 0.545), which we cannot conclusively classify. For systems lacking diffuse radio emission, we determine radio halo upper limits in two ways and find via survival analysis that these limits do not significantly affect radio power scaling relations. Several clusters with no diffuse emission detection are known or suspected mergers, based on archival X-ray and/or optical measures; given the limited sensitivity of our observations, deeper observations of these disturbed systems are required in order to rule out the presence of diffuse emission consistent with known scaling relations. In parallel with our diffuse emission results, we present catalogues of individual radio sources, including a few interesting extended sources. Our study represents the first step towards probing the occurrence of diffuse emission in high-redshift (z ≳ 0.5) clusters, and serves as a pilot for statistical studies of larger cluster samples with the new radio telescopes available in the pre-SKA era.

Knowles, Kenda↗
NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the respective standard approach. Of the methods tested, standard ComBat and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

GeneLab↗
NASA NTRSDate not supplied

Combining RNA-SEQ Datasets from NASA GENELAB: An Evaluation of Correction Methods

Background: Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets. Methods: In the present study, NASA GeneLab-hosted RNAseq datasets from mouse liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, the median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC to determine statistical differences between datasets before and after correction, Principal Component Analysis to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. Results: The results showed that the reference-based approach introduced several additional (and likely artificial) differentially expressed genes when compared with the respective standard approach. Conclusions: Of the methods tested, standard ComBat_seq and DESeq2 were identified as the most robust correction methods for combining spaceflight mouse liver RNAseq datasets hosted on GeneLab.

Finsam Samson↗
NASA NTRS1993

Statistics of arcs in clusters of galaxies

Samples of gravitational lens events in clusters show many large arcs compared to arclets, relative to what can be obtained by idealized singular lens models. We describe the probability of image magnification for point sources and for simple but more realistic gravitational lensing models that include a finite core size and an ellipticity. In addition, we explore the changes in the probability distribution of image magnifications, distortions, and angular extents for sources of different sizes as the parameters of the lenses are varied. A finite core in spherically symmetric lens models introduces a discontinuity in the probability distribution at which the relative number of highly magnified images is increased. In elliptical lenses, this discontinuity and its effect are replaced by a continuous increase in the probability of obtaining high-magnification images relative to singular spherically symmetric models. We also find that the finite size of the source causes a further increase in the expected number of images just below the maximum possible magnification.

Bergmann, Anton G.↗
DOE OSTI2022

Application of Linear Additive Conditions for Near-Infrared Diffuse Reflectance Absorption Spectroscopy

Determining the homogeneity of material mixing in real time during product processing is critical for quality control. According to the Kubelka–Munk (K-M) function of diffuse reflectance absorption spectrum, absorbance (A) is approximately linear with the content of the components when the sample scattering coefficient (S) is in a certain range. The S is determined by the particle size of powder samples. Therefore, this study determined particle size ranges that satisfy linear additivity in near-infrared diffuse reflectance spectroscopy (NIRDRS). Thus, the proposed NIRDRS analysis technique can be used to determine the homogeneity of material mixes or analyze the percentages of the components in the mixture. In this study, vitamin B3 and vitamin C were used for preparing mixed samples with varying percentages. The experimental results revealed that linear additivity is satisfied when the powder particle size is in the range of less than 280, 280–450, and 450–900 μm. When the confidence level is 0.01, the actual mixed spectra are not significantly different from the “simulated mixed spectra” constructed by linear addition, with their relative deviations less than 1.08%. The absolute errors of the actual and analytic percentages were within 2.98% for each component in the mixtures. The above conclusions also hold for sorghum, which has a complex material composition. Statistical models cannot analyze the percentages of components in the mixture. In contrast, linear addition and direct calibration approach avoids the use of a large number of samples for statistical modeling and analyze the percentages of mixed samples. Meanwhile, it can be used to discriminate and analyze the material mixing uniformity by building a mechanistic model.

Feng, Zhiyue↗
NASA NTRSDate not supplied

Evaluation of Correction Methods for NASA GeneLab Transcriptomic Datasets

Conducting space biology experiments aboard the International Space Station, particularly those utilizing complex model organisms like mice, is expensive and difficult due to limited crew availability, hardware, and space. As a result, sample numbers from these studies are low, reducing the statistical power of any one experiment. Aggregating spaceflight datasets serves as a method to increase sample numbers, allowing for novel insights through bioinformatic analysis of ‘omics data from merged datasets. However, aggregating datasets can introduce unwanted variation including 1) differences in sample handling, processing, and sequencing platforms between datasets (technical variation) as well as 2) differences in experimental design between datasets such as sex or age of the model organism used. In the present study, NASA GeneLab-hosted RNAseq datasets from rodent liver tissues were used to evaluate several statistical methods to correct for this unwanted variation through two approaches, reference-based and standard. The following correction algorithms were applied with (reference-based) and/or without (standard) considering Universal Mouse RNA Reference samples: ComBat and ComBat_seq from the SVA package, median polish, empirical Bayes, and ANOVA-based algorithms from the MBatch package, and negative binomial regression normalization in the DESeq2 package. For each approach, after the correction algorithm was applied, differential gene expression (DGE) analysis of flight and ground control samples was performed with the combined data. The robustness of each tool was evaluated using BatchQC, to determine statistical differences between datasets before and after correction, Principal Component Analysis, to evaluate global gene expression in samples before and after correction, and by comparing DGE analysis of individual datasets and combined datasets before and after correction. The results showed that the reference-based approach introduced several additional (and likely artificial) DEGs when compared with the standard approach. Thus, the most robust standard correction will be implemented in the GeneLab Visualization 2.0 platform when datasets are combined.

GeneLab, RNA-seq, Batch Correction↗
DOE OSTI2024

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Mechanical properties of silicate glasses exposed to a low-Earth orbit

The effects of a 5.8 year exposure to low earth orbit environment upon the mechanical properties of commercial optical fused silica, low iron soda-lime-silica, Pyrex 7740, Vycor 7913, BK-7, and the glass ceramic Zerodur were examined. Mechanical testing employed the ASTM-F-394 piston on 3-ball method in a liquid nitrogen environment. Samples were exposed on the Long Duration Exposure Facility (LDEF) in two locations. Impacts were observed on all specimens except Vycor. Weibull analysis as well as a standard statistical evaluation were conducted. The Weibull analysis revealed no differences between control samples and the two exposed samples. We thus concluded that radiation components of the Earth orbital environment did not degrade the mechanical strength of the samples examined within the limits of experimental error. The upper bound of strength degradation for meteorite impacted samples based upon statistical analysis and observation was 50 percent.

Wiedlocher, David E.↗
NASA NTRS2013

Ensemble Data Assimilation Without Ensembles: Methodology and Application to Ocean Data Assimilation

Two methods to estimate background error covariances for data assimilation are introduced. While both share properties with the ensemble Kalman filter (EnKF), they differ from it in that they do not require the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The first method is referred-to as SAFE (Space Adaptive Forecast error Estimation) because it estimates error covariances from the spatial distribution of model variables within a single state vector. It can thus be thought of as sampling an ensemble in space. The second method, named FAST (Flow Adaptive error Statistics from a Time series), constructs an ensemble sampled from a moving window along a model trajectory. The underlying assumption in these methods is that forecast errors in data assimilation are primarily phase errors in space and/or time.

Data Assimilation↗