Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models, Statistical”

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 73 records · Page 4

SU-C-BRD-01: A Statistical Modeling Method for Quality Control of Intensity- Modulated Radiation Therapy Planning

Purpose: To apply a statistical modeling approach, threshold modeling (TM), for quality control of intensity-modulated radiation therapy (IMRT) treatment plans. Methods: A quantitative measure, which was the weighted sum of violations of dose/dose-volume constraints, was first developed to represent the quality of each IMRT plan. Threshold modeling approach, which is is an extension of extreme value theory in statistics and is an effect way to model extreme values, was then applied to analyze the quality of the plans summarized by our quantitative measures. Our approach modeled the plans generated by planners as a series of independent and identically distributed random variables and described the behaviors of them if the plan quality was controlled below certain threshold. We tested our approach with five locally advanced head and neck cancer patients retrospectively. Two statistics were incorporated for numerical analysis: probability of quality improvement (PQI) of the plans and expected amount of improvement on the quantitative measure (EQI). Results: After clinical planners generated 15 plans for each patient, we applied our approach to obtain the PQI and EQI as if planners would generate additional 15 plans. For two of the patients, the PQI was significantly higher than the other three (0.17 and 0.18 comparing to 0.08, 0.01 and 0.01). The actual percentage of the additional 15 plans that outperformed the best of initial 15 plans was 20% and 27% comparing to 11%, 0% and 0%. EQI for the two potential patients were 34.5 and 32.9 and the rest of three patients were 9.9, 1.4 and 6.6. The actual improvements obtained were 28.3 and 20.5 comparing to 6.2, 0 and 0. Conclusion: TM is capable of reliably identifying the potential quality improvement of IMRT plans. It provides clinicians an effective tool to assess the trade-off between extra planning effort and achievable plan quality. This work was supported in part by NIH/NCI grant CA130814.

60 APPLIED LIFE SCIENCES↗

Statistical Modeling of CTV Motion and Deformation for IMRT of Early-Stage Rectal Cancer

Purpose: To derive and validate a statistical model of motion and deformation for the clinical target volume (CTV) of early-stage rectal cancer patients. Methods and Materials: For 16 patients, 4 to 5 magnetic resonance images (MRI) were acquired before each fraction was administered. The CTV was delineated on each MRI. Using a leave-one-out methodology, we constructed a population-based principal component analysis (PCA) model of the CTV motion and deformation of 15 patients, and we tested the model on the left-out patient. The modeling error was calculated as the amount of the CTV motion-deformation of the left-out-patient that could not be explained by the PCA model. Next, the PCA model was used to construct a PCA target volume (PCA-TV) by accumulating motion-deformations simulated by the model. A PCA planning target volume (PTV) was generated by expanding the PCA-TV by uniform margins. The PCA-PTV was compared with uniform and nonuniform CTV-to-PTV margins. To allow comparison, geometric margins were determined to ensure adequate coverage, and the volume difference between the PTV and the daily CTV (CTV-to-PTV volume) was calculated. Results: The modeling error ranged from 0.9 ± 0.5 to 2.9 ± 2.1 mm, corresponding to a reduction of the CTV motion-deformation between 6% and 60% (average, 23% ± 11%). The reduction correlated with the magnitude of the CTV motion-deformation (P<.001, R=0.66). The PCA-TV and the CTV required 2-mm and 7-mm uniform margins, respectively. The nonuniform CTV-to-PTV margins were 4 mm in the left, right, inferior, superior, and posterior directions and 8 mm in the anterior direction. Compared to uniform and nonuniform CTV-to-PTV margins, the PCA-based PTV significantly decreased (P<.001) the average CTV-to-PTV volume by 128 ± 20 mL (49% ± 4%) and by 35 ± 6 mL (20% ± 3.5%), respectively. Conclusions: The CTV motion-deformation of a new patient can be explained by a population-based PCA model. A PCA model-generated PTV significantly improved sparing of organs at risk compared to uniform and nonuniform CTV-to-PTV margins.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A statistical model for radar images of agricultural scenes

The presently derived and validated statistical model for radar images containing many different homogeneous fields predicts the probability density functions of radar images of entire agricultural scenes, thereby allowing histograms of large scenes composed of a variety of crops to be described. Seasat-A SAR images of agricultural scenes are accurately predicted by the model on the basis of three assumptions: each field has the same SNR, all target classes cover approximately the same area, and the true reflectivity characterizing each individual target class is a uniformly distributed random variable. The model is expected to be useful in the design of data processing algorithms and for scene analysis using radar images.

Frost, V. S.↗

Statistical model based iterative reconstruction (MBIR) in clinical CT systems: Experimental assessment of noise performance

Purpose: To reduce radiation dose in CT imaging, the statistical model based iterative reconstruction (MBIR) method has been introduced for clinical use. Based on the principle of MBIR and its nonlinear nature, the noise performance of MBIR is expected to be different from that of the well-understood filtered backprojection (FBP) reconstruction method. The purpose of this work is to experimentally assess the unique noise characteristics of MBIR using a state-of-the-art clinical CT system. Methods: Three physical phantoms, including a water cylinder and two pediatric head phantoms, were scanned in axial scanning mode using a 64-slice CT scanner (Discovery CT750 HD, GE Healthcare, Waukesha, WI) at seven different mAs levels (5, 12.5, 25, 50, 100, 200, 300). At each mAs level, each phantom was repeatedly scanned 50 times to generate an image ensemble for noise analysis. Both the FBP method with a standard kernel and the MBIR method (Veo{sup ®}, GE Healthcare, Waukesha, WI) were used for CT image reconstruction. Three-dimensional (3D) noise power spectrum (NPS), two-dimensional (2D) NPS, and zero-dimensional NPS (noise variance) were assessed both globally and locally. Noise magnitude, noise spatial correlation, noise spatial uniformity and their dose dependence were examined for the two reconstruction methods. Results: (1) At each dose level and at each frequency, the magnitude of the NPS of MBIR was smaller than that of FBP. (2) While the shape of the NPS of FBP was dose-independent, the shape of the NPS of MBIR was strongly dose-dependent; lower dose lead to a “redder” NPS with a lower mean frequency value. (3) The noise standard deviation (σ) of MBIR and dose were found to be related through a power law of σ ∝ (dose){sup −β} with the component β ≈ 0.25, which violated the classical σ ∝ (dose){sup −0.5} power law in FBP. (4) With MBIR, noise reduction was most prominent for thin image slices. (5) MBIR lead to better noise spatial uniformity when compared with FBP. (6) A composite image generated from two MBIR images acquired at two different dose levels (D1 and D2) demonstrated lower noise than that of an image acquired at a dose level of D1+D2. Conclusions: The noise characteristics of the MBIR method are significantly different from those of the FBP method. The well known tradeoff relationship between CT image noise and radiation dose has been modified by MBIR to establish a more gradual dependence of noise on dose. Additionally, some other CT noise properties that had been well understood based on the linear system theory have also been altered by MBIR. Clinical CT scan protocols that had been optimized based on the classical CT noise properties need to be carefully re-evaluated for systems equipped with MBIR in order to maximize the method's potential clinical benefits in dose reduction and/or in CT image quality improvement.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Statistical Modeling of Landscape Pattern and the Effects of Coarse Spatial Resolution

Analysis of classified satellite imagery was conducted to characterize errors in estimates of area based on coarse resolution satellite imagery which are due to distortions in sizes of small fragments, and to explore the feasibility of correcting for these errors using a statistical modeling approach. Sizes of bodies of open water on European Remote-Sensing (ERS-1) synthetic aperture radar (SAR) and fire scars on Landsat Multispectral Scanner (MSS) imagery were measured using geographic information system (GIS) software. Statistical analysis of the smaller scars and ponds as observed with this imagery of relatively fine resolution demonstrated that the distribution of the sizes could be modeled by either of two types of statistical distributions-a power distribution related to fractal processes or a simple exponential distribution. A comparison of the size distribution of small burn scars as observed with Landsat with the distribution observed with the much coarser Advanced Very High Resolution Radiometer (AVHRR) imagery revealed effects of coarse resolution, each of which were associated with a particular size range, and whose net effect accounted for the bias in total burn area due to coarse spatial resolution. Size dependent area correction procedures are therefore feasible. A statistical modeling approach to area estimation, based models of the underlying and observed size distributions of fragments in the scene, may lead to a procedure which requires little or no ancillary data.

Hlavka, C. A.↗

Most theoretical approaches used in nuclear astrophysics to model the nucleosynthesis of heavy elements incorporate the so-called statistical model in order to describe the excitation and decay properties of atomic nuclei. One of the basic assumptions of this model is the validity of the Brink–Axel hypothesis and the related concept of so-called photon strength functions to describe γ-ray transition probabilities. We present a novel experimental approach that allows for the first time to experimentally determine the photon strength function simultaneously in two independent ways by a unique combination of quasi-monochromatic photon beams and a newly implemented γ–γ coincidence setup. This technique does not assume a priori the validity of the Brink–Axel hypothesis and sets a benchmark in terms of the detection sensitivity for measuring decay properties of photo-excited states below the neutron separation energy. The data for the spherical off-shell nucleus 128 Te were obtained for γ-ray beam-energy settings between 3 MeV and 9 MeV in steps of 130 keV for the lower beam energies and in steps of up to 280 keV for the highest beam settings. We present a quantitative analysis on the consistency of the derived photon strength function with the Brink–Axel hypothesis. The data clearly demonstrate a discrepancy of up to a factor of two between the photon strength functions extracted from the photoabsorption and photon emission process, respectively. In addition, we observe that the photon strength functions are not independent of the excitation energy, as usually assumed. Thus, we conclude, that the Brink–Axel hypothesis is not strictly fulfilled in the excitation-energy region below the neutron separation threshold (S n = 8.78 MeV) for the studied case of 128 Te.

Not Available

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enantioselective Intramolecular Allylic Substitution via Synergistic Palladium/Chiral Phosphoric Acid Catalysis: Insight into Stereoinduction through Statistical Modeling

The mode of asymmetric induction in an enantioselective intramolecular allylic substitution reaction catalyzed by a combination of palladium and a chiral phosphoric acid was investigated by a combined experimental and statistical modeling approach. Experiments to probe nonlinear effects, the reactivity of deuterium-labeled substrates, and control experiments revealed that nucleophilic attack to the π-allylpalladium intermediate is the enantio-determining step, in which the chiral phosphate anion is involved in stereoinduction. Here, using multivariable linear regression analysis, we determined that multiple noncovalent interactions with the chiral environment of the phosphate anion are integral to enantiocontrol in the transition state. The synthetic protocol to form chiral pyrrolidines was further applied to the asymmetric construction of C-O bonds at fully substituted carbon centers in the synthesis of chiral 2,2-disubstituted benzomorpholines.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Statistical Modeling Efforts for Headspace Gas

The purpose of this document is to describe the statistical modeling effort for gas concentrations in WIPP storage containers. The concentration (in ppm) of CO 2 in the headspace volume of standard waste box (SWB) 68685 is shown. A Bayesian approach and an adaptive Metropolis-Hastings algorithm were used.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Comparing Two Statistical Models for Low Boom Dose-response Relationships with Correlated Responses

This study compares two statistical modeling approaches to construct summary dose-response curves for low boom community noise surveys. NASA field survey data were used that consist of multiple responses from a survey participant. These data require an approach that accounts for the correlation among repeated annoyance observations from the same participant. A multilevel model accounts for the correlation by allowing estimated parameters to vary with each survey participant. On the other hand, a population average model utilizes generalized estimating equations and accounts for the correlation via a user-specified within-subject correlation structure. A visual comparison of the dose-response curves for these two methods reveals similar results. When comparing specific points along the summary curves, the multilevel model yields more precise confidence bounds than the population average model. The similarity between the summary curves derived from each model lends validity to both approaches for approximating a population representative summary curve.

X-59↗

Comparison of Two Statistical Models for Low Boom Dose-response Relationships with Correlated Responses

This study compares two statistical models to construct summary dose-response curves for low boom community noise surveys. Data from two NASA field surveys are used that consist of multiple responses per survey participant. These data require an approach that accounts for the correlation among repeated annoyance observations from the same participant. A multilevel model accounts for the correlation by allowing estimated parameters to vary with each survey participant. On the other hand, a population average model utilizes generalized estimating equations and accounts for the correlation via a userspecified within-subject correlation structure. A visual comparison of the dose-response curves for these two methods reveals similar results. When comparing specific points along the summary curves, the multilevel model yields more precise confidence bounds than the population average model. The similarity between the summary curves derived from each model lends validity to both approaches for approximating a population representative summary curve, though modeling assumptions may lend favor to the multilevel logistic modeling approach over the population average model.

dose-response↗

A Statistical Model of Tropical Cyclone Tracks in the Western North Pacific with ENSO-Dependent Cyclogenesis

A new statistical model for western North Pacific Ocean tropical cyclone genesis and tracks is developed and applied to estimate regionally resolved tropical cyclone landfall rates along the coasts of the Asian mainland, Japan, and the Philippines. The model is constructed on International Best Track Archive for Climate Stewardship (IBTrACS) 1945-2007 historical data for the western North Pacific. The model is evaluated in several ways, including comparing the stochastic spread in simulated landfall rates with historic landfall rates. Although certain biases have been detected, overall the model performs well on the diagnostic tests, for example, reproducing well the geographic distribution of landfall rates. Western North Pacific cyclogenesis is influenced by El Nino-Southern Oscillation (ENSO). This dependence is incorporated in the model s genesis component to project the ENSO-genesis dependence onto landfall rates. There is a pronounced shift southeastward in cyclogenesis and a small but significant reduction in basinwide annual counts with increasing ENSO index value. On almost all regions of coast, landfall rates are significantly higher in a negative ENSO state (La Nina).

Yonekura, Emmi↗

Cellular Statistical Models of Broken Cloud Fields. Part IV: Effects of Pixel Size on Idealized Satellite Observations

In the fourth part of our “Cellular Statistical Models of Broken Cloud Fields” series we use the binary Markov processes framework for quantitative investigation of the effects of low resolution of idealized satellite observations on the statistics of the retrieved cloud masks. We assume that the cloud fields are Markovian and are characterized by the “actual” cloud fraction (CF) and scale length. We use two different models of observations: a simple discrete-point sampling and a more realistic “pixel” protocol. The latter is characterized by a state attribution function (SAF) which has the meaning of the probability that the pixel with a certain CF is declared cloudy in the observed cloud mask. The stochasticity of the SAF means that the cloud/clear attribution is not ideal and can be affected by external or unknown factors. We show that the observed cloud masks can be accurately described as Markov chains of pixels and use the master-matrix formalism (introduced in Part III of the series) for analytical computation of their parameters: the “observed” CF and scale length. This procedure allows us to establish a quantitative relationship (which is pixel-size dependent) between the actual and the observed cloud-field statistics. The feasibility of restoring the former from the latter is considered. The adequacy of our analytical approach to idealized observations is evaluated using numerical simulations. Comparison of the observed parameters of the simulated datasets with their theoretical expectations showed an agreement within 0.005 for the CF, while for the scale length it is within 1% in the sampling case and within 4% in the pixel case.

Mikhail D. Alexandrov↗

Markovian Statistical Model of Cloud Optical Thickness. Part I: Theory and Examples

We present a generalization of the binary-value Markovian model previously used for statistical characterization of cloud masks to a continuous-value model describing 1D fields of cloud optical thickness (COT). This model has simple functional expressions and is specified by four parameters: the cloud fraction, the autocorrelation (scale) length, and the two parameters of the normalized probability density function of (non-zero) COT values (this PDF is assumed to have gamma-distribution form). Cloud masks derived from this model by separation between the values above and below some threshold in COT appear to have the same statistical properties as in binary-value model described in our previous publications. We demonstrate the ability of our model to generate examples of various cloud-field types by using it to statistically imitate actual cloud observations made by the Research Scanning Polarimeter (RSP) during two field experiments.

binary-value Markovian model↗