Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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 19 records

LOG-NORMAL DISTRIBUTION OF COSMIC VOIDS IN SIMULATIONS AND MOCKS

Following up on previous studies, we complete here a full analysis of the void size distributions of the Cosmic Void Catalog based on three different simulation and mock catalogs: dark matter (DM), haloes, and galaxies. Based on this analysis, we attempt to answer two questions: Is a three-parameter log-normal distribution a good candidate to satisfy the void size distributions obtained from different types of environments? Is there a direct relation between the shape parameters of the void size distribution and the environmental effects? In an attempt to answer these questions, we find here that all void size distributions of these data samples satisfy the three-parameter log-normal distribution whether the environment is dominated by DM, haloes, or galaxies. In addition, the shape parameters of the three-parameter log-normal void size distribution seem highly affected by environment, particularly existing substructures. Therefore, we show two quantitative relations given by linear equations between the skewness and the maximum tree depth, and between the variance of the void size distribution and the maximum tree depth, directly from the simulated data. In addition to this, we find that the percentage of voids with nonzero central density in the data sets has a critical importance. If the number of voids with nonzero central density reaches ≥3.84% in a simulation/mock sample, then a second population is observed in the void size distributions. This second population emerges as a second peak in the log-normal void size distribution at larger radius.

79 ASTRONOMY AND ASTROPHYSICS↗

On multichannel film dosimetry with channel-independent perturbations

Purpose: Different multichannel methods for film dosimetry have been proposed in the literature. Two of them are the weighted mean method and the method put forth byMicke et al. [“Multichannel film dosimetry with nonuniformity correction,” Med. Phys. 38, 2523–2534 (2011)] and Mayer et al. [“Enhanced dosimetry procedures and assessment for EBT2 radiochromic film,” Med. Phys. 39, 2147–2155 (2012)]. The purpose of this work was to compare their results and to develop a generalized channel-independent perturbations framework in which both methods enter as special cases. Methods: Four models of channel-independent perturbations were compared: weighted mean, Micke–Mayer method, uniform distribution, and truncated normal distribution. A closed-form formula to calculate film doses and the associated type B uncertainty for all four models was deduced. To evaluate the models, film dose distributions were compared with planned and measured dose distributions. At the same time, several elements of the dosimetry process were compared: film type EBT2 versus EBT3, different waiting-time windows, reflection mode versus transmission mode scanning, and planned versus measured dose distribution for film calibration and for γ-index analysis. The methods and the models described in this study are publicly accessible through IRISEU. Alpha 1.1 ( http://www.iriseu.com ). IRISEU. is a cloud computing web application for calibration and dosimetry of radiochromic films. Results: The truncated normal distribution model provided the best agreement between film and reference doses, both for calibration and γ-index verification, and proved itself superior to both the weighted mean model, which neglects correlations between the channels, and the Micke–Mayer model, whose accuracy depends on the properties of the sensitometric curves. With respect to the selection of dosimetry protocol, no significant differences were found between transmission and reflection mode scanning, between 75 ± 5 min and 20 ± 1 h waiting-time windows or between employing EBT2 or EBT3 films. Significantly better results were obtained when a measured dose distribution was used instead of a planned one as reference for the calibration, and when a planned dose distribution was used instead of a measured one as evaluation for the γ-analysis. Conclusions: The truncated normal distribution model of channel-independent perturbations was found superior to the other three models under comparison and the authors propose its use for multichannel dosimetry.

61 RADIATION PROTECTION AND DOSIMETRY↗

Quasars with large proper motions: A selection from the LQAC-5 catalogue combined with Gaia EDR3: Focusing on astrometric and photometric properties

Thanks to the cross-identification between the Gaia EDR3 (E-Data Release 3) catalogue and the quasars’ compiled catalogue LQAC-5 (Large Quasar Astrometric Catalogue-5), accurate astrometric parameters as well as photometric measurements were available for a set of 416 113 quasars. After analysing the astrometric and photometric properties of these quasars, we concentrate our study on objects characterized by significant proper motion, which contradicts the postulate that they are fixed in the celestial sphere. We derived statistical properties of the positions, parallaxes, and proper motions of the quasars cross-matched with Gaia EDR3. We carried out a specific investigation of quasars with a proper motion larger than 10 mas yr -1 , and used colour-colour diagrams in mid-infrared to discriminate these quasars from stellar contamination. We discuss, in detail, a few cases where the abnormal proper motion could be caused by misidentifications. We note that 43% of LQAC-5 quasars were not detected by Gaia EDR3, mainly because of the cut-off magnitude at G = 21. The comparison between r(G) (Gaia) and r (SDSS) shows a normal distribution. The normalized distributions of parallaxes and proper motions are identical to that of the Gaia CRF2, as well as the distribution of the G BP -G RP colour index. We have identified 937 LQAC-5 objects (0.27% of the sample) with Gaia EDR3 proper motion which is significant at the 5σ level. We have clearly shown that below z = 0.35, the integrated magnitude G(r) obtained from SDSS spectra is smaller than G (Gaia), leading to the statement that the quasars are generally resolved under this threshold, which should lead to a bad determination of proper motion. Moreover, to discriminate quasars from stellar contamination, we have plotted them in a mid-infrared colour–colour diagram. Finally, we have made a final selection of 41 objects with μ > 10 mas yr -1 , satisfying drastic constraints, which makes them potential candidates for considering their proper motion as real. In this paper we carried out a complete analysis of the quasars of the LQAC-5 which were cross-identified with Gaia EDR3 ; this serves as an alternative and complementary study with respect to Gaia CRF2, involving a different population of quasars. We selected a set of 41 quasars with a proper motion exceeding 10 mas yr -1 , which can be considered as very high for objects which are a priori fixed in the celestial sphere.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-D-BRD-07: Evaluation of the Effectiveness of Statistical Process Control Methods to Detect Systematic Errors For Routine Electron Energy Verification

Purpose: To evaluate the ability of statistical process control methods to detect systematic errors when using a two dimensional (2D) detector array for routine electron beam energy verification. Methods: Electron beam energy constancy was measured using an aluminum wedge and a 2D diode array on four linear accelerators. Process control limits were established. Measurements were recorded in control charts and compared with both calculated process control limits and TG-142 recommended specification limits. The data was tested for normality, process capability and process acceptability. Additional measurements were recorded while systematic errors were intentionally introduced. Systematic errors included shifts in the alignment of the wedge, incorrect orientation of the wedge, and incorrect array calibration. Results: Control limits calculated for each beam were smaller than the recommended specification limits. Process capability and process acceptability ratios were greater than one in all cases. All data was normally distributed. Shifts in the alignment of the wedge were most apparent for low energies. The smallest shift (0.5 mm) was detectable using process control limits in some cases, while the largest shift (2 mm) was detectable using specification limits in only one case. The wedge orientation tested did not affect the measurements as this did not affect the thickness of aluminum over the detectors of interest. Array calibration dependence varied with energy and selected array calibration. 6 MeV was the least sensitive to array calibration selection while 16 MeV was the most sensitive. Conclusion: Statistical process control methods demonstrated that the data distribution was normally distributed, the process was capable of meeting specifications, and that the process was centered within the specification limits. Though not all systematic errors were distinguishable from random errors, process control limits increased the ability to detect systematic errors using routine measurement of electron beam energy constancy.

07 ISOTOPE AND RADIATION SOURCES↗

Sample Generation for Nuclear Data

This report summarizes a NEAMS (Nuclear Energy Advanced Modeling and Simulation) project focused on developing a sampling capability that can handle the challenges of generating samples from nuclear cross-section data. The covariance information between energy groups tends to be very ill-conditioned and thus poses a problem using traditional methods for generated correlated samples. This report outlines a method that addresses the sample generation from cross-section matrices. The treatment allows one to assume the cross sections are distributed with a multivariate normal distribution, lognormal distribution, or truncated normal distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Notes on power of normality tests of error terms in regression models

Normality is one of the basic assumptions in applying statistical procedures. For example in linear regression most of the inferential procedures are based on the assumption of normality, i.e. the disturbance vector is assumed to be normally distributed. Failure to assess non-normality of the error terms may lead to incorrect results of usual statistical inference techniques such as t-test or F-test. Thus, error terms should be normally distributed in order to allow us to make exact inferences. As a consequence, normally distributed stochastic errors are necessary in order to make a not misleading inferences which explains a necessity and importance of robust tests of normality. Therefore, the aim of this contribution is to discuss normality testing of error terms in regression models. In this contribution, we introduce the general RT class of robust tests for normality, and present and discuss the trade-off between power and robustness of selected classical and robust normality tests of error terms in regression models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating the ecological fallacy through sampling distributions constructed from finite populations

Correlation coefficients and linear regression values computed from group averages can differ from correlation coefficients and linear regression values computed using individual scores. This observation known as the ecological fallacy often assumes that all the individual scores are available from a population. In many situations, one must use a sample from the larger population. In such cases, the computed correlation coefficient and linear regression values will depend on the sample that is chosen and the underlying sampling distribution. The sampling distribution of correlation coefficients and linear regression values for group averages will be identical to the sampling distribution for individuals for normally distributed variables for random samples drawn from infinitely large continuous distributions. However, data that is acquired in practice is often acquired when sampling without replacement from a finite population. Our objective is to demonstrate through Monte Carlo simulations that the sampling distributions for correlation and linear regression will also be similar for individuals and group averages when sampling without replacement from normally distributed variables. These simulations suggest that when a random sample from a population is selected, the correlation coefficients and linear regression values computed from individual scores will not be more accurate in estimating the entire population values compared to samples when group averages are used as long as the sample size is the same.

97 MATHEMATICS AND COMPUTING↗

Deep neural network improves the estimation of polygenic risk scores for breast cancer

Polygenic risk scores (PRS) estimate the genetic risk of an individual for a complex disease based on many genetic variants across the whole genome. Here, we compared a series of computational models for estimation of breast cancer PRS. A deep neural network (DNN) was found to outperform alternative machine learning techniques and established statistical algorithms, including BLUP, BayesA, and LDpred. In the test cohort with 50% prevalence, the Area Under the receiver operating characteristic Curve (AUC) were 67.4% for DNN, 64.2% for BLUP, 64.5% for BayesA, and 62.4% for LDpred. BLUP, BayesA, and LPpred all generated PRS that followed a normal distribution in the case population. However, the PRS generated by DNN in the case population followed a bimodal distribution composed of two normal distributions with distinctly different means. This suggests that DNN was able to separate the case population into a high-genetic-risk case subpopulation with an average PRS significantly higher than the control population and a normal-genetic-risk case subpopulation with an average PRS similar to the control population. This allowed DNN to achieve 18.8% recall at 90% precision in the test cohort with 50% prevalence, which can be extrapolated to 65.4% recall at 20% precision in a general population with 12% prevalence. Interpretation of the DNN model identified salient variants that were assigned insignificant p values by association studies, but were important for DNN prediction. These variants may be associated with the phenotype through nonlinear relationships.

59 BASIC BIOLOGICAL SCIENCES↗

Data for "Quantifying the Propagation of Parametric Uncertainty on Flux Balance Analysis"

In the repository are example scripts that perform uncertainty injection and propagation to flux balance analysis with outputs for a small sample size (for demonstration purpose only). For proper analysis, user should download the scripts and run for a large sample size (e.g., 10,000 samples). If you use the scripts, please cite the following Metabolic Engineering article: “Quantifying the propagation of parametric uncertainty on flux balance analysis” (https://doi.org/10.1016/j.ymben.2021.10.012) There are two subdirectories: /uncFBA/uncBiom: injection of normally distributed noise to biomass precursor coeffcients and ATP maintenance (growth-associated ATP maintenance (GAM) and non-growth associated ATP maintenance (NGAM)) /uncFBA/uncRHS: departure from steady-state by adding noise drawn from normal distribution to the RHS terms of mass balance constraints

Metabolomics↗

SU-E-T-586: Optimal Determination of Tolerance Level for Radiation Dose Delivery Verification in An in Vivo Dosimetry System

Purpose: To statistically determine the optimal tolerance level in the verification of delivery dose compared to the planned dose in an in vivo dosimetry system in radiotherapy. Methods: The LANDAUER MicroSTARii dosimetry system with screened nanoDots (optically stimulated luminescence dosimeters) was used for in vivo dose measurements. Ideally, the measured dose should match with the planned dose and falls within a normal distribution. Any deviation from the normal distribution may be redeemed as a mismatch, therefore a potential sign of the dose misadministration. Randomly mis-positioned nanoDots can yield a continuum background distribution. A percentage difference of the measured dose to its corresponding planned dose (ΔD) can be used to analyze combined data sets for different patients. A model of a Gaussian plus a flat function was used to fit the ΔD distribution. Results: Total 434 nanoDot measurements for breast cancer patients were collected across a period of three months. The fit yields a Gaussian mean of 2.9% and a standard deviation (SD) of 5.3%. The observed shift of the mean from zero is attributed to the machine output bias and calibration of the dosimetry system. A pass interval of −2SD to +2SD was applied and a mismatch background was estimated to be 4.8%. With such a tolerance level, one can expect that 99.99% of patients should pass the verification and at most 0.011% might have a potential dose misadministration that may not be detected after 3 times of repeated measurements. After implementation, a number of new start breast cancer patients were monitored and the measured pass rate is consistent with the model prediction. Conclusion: It is feasible to implement an optimal tolerance level in order to maintain a low limit of potential dose misadministration while still to keep a relatively high pass rate in radiotherapy delivery verification.

61 RADIATION PROTECTION AND DOSIMETRY↗

SU-E-T-671: Range-Modulation Effects of Carbon Ion Beams in Lung Tissue

Purpose: When particles traversing inhomogeneous materials like lung they show a characteristic range modulation which cannot be observed in homogeneous materials. It is possible to describe the range modulation by a convolution of an unperturbed Bragg-Curve and a normal distribution. The sigma of the normal distribution is a parameter for the strength of the modulation effect. A new material parameter (modulation power, P-mod) is introduced which is independent of the material thickness. It is defined as the square of sigma divided by the mean water equivalent thickness of the target (µ). Methods: The modulation power of lung tissue was determined by actual Bragg-peak measurements after traversing an ex-vivo porcine lung and by Monte-Carlo simulations with micro-CT data of human lung tissue. The determined modulation powers were used to show the effect of range modulation effects in a simplified treatment situation. A four centimeter spread-out Bragg-peak after traversing eight centimeter of lung tissue was simulated in FLUKA. The SOBP with and without consideration of range modulation effects were compared. Results: As well in the measurements as in the MC simulations range modulation effects of lung tissue were observed. The determined modulation powers showed a great range from 0.05 mm, in the micro-CT data, to 0.7 mm in the lung measurements. The SOBP comparison showed that range modulation effects Result in over- and underdosages at the distal and proximal edge of the SOBP. In the investigated case, the last 0.5 cm of the SOBP showed an underdosage of up to 50% at the distal edge, while 0.5 cm distal to the SOBP an overdosage of up to 50% was observed. Conclusion: Range modulation effects occur in inhomogeneous materials like lung. These modulation effects may Result in clinically relevant over- and underdosages but are currently not considered in commercially available treatment planning systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

New Fitting Formula for Cosmic Nonlinear Density Distribution

We have measured the probability distribution function (PDF) of a cosmic matter density field from a suite of N-body simulations. We propose the generalized normal distribution of version 2 (N{sub v2}) as an alternative fitting formula to the well-known log-normal distribution. We find that N{sub v2} provides a significantly better fit than that of the log-normal distribution for all smoothing radii (2, 5, 10, 25 [Mpc h {sup −1}]) that we studied. The improvement is substantial in the underdense regions. The development of non-Gaussianities in the cosmic matter density field is captured by continuous evolution of the skewness and shift parameters of the N{sub v2} distribution. We present the redshift evolution of these parameters for aforementioned smoothing radii and various background cosmology models. All the PDFs measured from large and high-resolution N-body simulations that we use in this study can be obtained from the web site https://astro.kias.re.kr/jhshin.

79 ASTRONOMY AND ASTROPHYSICS↗

Is the assumption of normality or log-normality for continuous response data critical for benchmark dose estimation?

Continuous responses (e.g. body weight) are widely used in risk assessment for determining the benchmark dose (BMD) which is used to derive a U.S. EPA reference dose. One critical question that is not often addressed in dose–response assessments is whether to model the continuous data as normally or log-normally distributed. Additionally, if lognormality is assumed, and only summarized response data (i.e., mean ± standard deviation) are available as is usual in the peer-reviewed literature, the BMD can only be approximated. In this study, using the “hybrid” method and relative deviation approach, we first evaluate six representative continuous dose–response datasets reporting individual animal responses to investigate the impact on BMD/BMDL estimates of (1) the distribution assumption and (2) the use of summarized versus individual animal data when a log-normal distribution is assumed. We also conduct simulation studies evaluating model fits to various known distributions to investigate whether the distribution assumption has influence on BMD/BMDL estimates. Our results indicate that BMDs estimated using the hybrid method are more sensitive to the distribution assumption than counterpart BMDs estimated using the relative deviation approach. The choice of distribution assumption has limited impact on the BMD/BMDL estimates when the within dose-group variance is small, while the lognormality assumption is a better choice for relative deviation method when data are more skewed because of its appropriateness in describing the relationship between mean and standard deviation. Additionally, the results suggest that the use of summarized data versus individual response data to characterize log-normal distributions has minimal impact on BMD estimates. - Highlights: • We investigate to what extent the distribution assumption can affect BMD estimates. • Both real data analysis and simulation study are conducted. • BMDs estimated using hybrid method are more sensitive to distribution assumption. • Summarized continuous data are adequate for BMD estimation.

60 APPLIED LIFE SCIENCES↗

Energetics and Birth Rates of Supernova Remnants in the Large Magellanic Cloud

Published X-ray emission properties for a sample of 50 supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) are used as input for SNR evolution modeling calculations. The forward shock emission is modeled to obtain the initial explosion energy, age, and circumstellar medium density for each SNR in the sample. The resulting age distribution yields a SNR birthrate of 1/(500 yr) for the LMC. The explosion energy distribution is well fit by a log-normal distribution, with a most-probable explosion energy of 0.5×10{sup 51} erg, with a 1σ dispersion by a factor of 3 in energy. The circumstellar medium density distribution is broader than the explosion energy distribution, with a most-probable density of ∼0.1 cm{sup −3}. The shape of the density distribution can be fit with a log-normal distribution, with incompleteness at high density caused by the shorter evolution times of SNRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Distribution Development for STORM Ingestion Input Parameters

The Sandia-developed Transport of Radioactive Materials (STORM) code suite is used as part of the Radioisotope Power System Launch Safety (RPSLS) program to perform statistical modeling of the consequences due to release of radioactive material given a launch accident. As part of this modeling, STORM samples input parameters from probability distributions with some parameters treated as constants. This report described the work done to convert four of these constant inputs (Consumption Rate, Average Crop Yield, Cropland to Landuse Database Ratio, and Crop Uptake Factor) to sampled values. Consumption rate changed from a constant value of 557.68 kg / yr to a normal distribution with a mean of 102.96 kg / yr and a standard deviation of 2.65 kg / yr. Meanwhile, Average Crop Yield changed from a constant value of 3.783 kg edible / m 2 to a normal distribution with a mean of 3.23 kg edible / m 2 and a standard deviation of 0.442 kg edible / m 2 . The Cropland to Landuse Database ratio changed from a constant value of 0.0996 (9.96%) to a normal distribution with a mean value of 0.0312 (3.12%) and a standard deviation of 0.00292 (0.29%). Finally the crop uptake factor changed from a constant value of 6.37e -4 (Bq crop /kg)/(Bq soil /kg) to a lognormal distribution with a geometric mean value of 3.38e -4 (Bq crop /kg)/(Bq soil /kg) and a standard deviation value of 3.33 (Bq crop /kg)/(Bq soil /kg)

60 APPLIED LIFE SCIENCES↗