Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multivariate regression”

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 55 records · Page 3

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Unmitigated uncertainties are known to have previously led to failed development programs; in order to combat these uncertainties, risks and their impacts must be understood and handled to ensure program success. In this paper, a probabilistic methodology to handle uncertainties is demonstrated on a three-element Human Landing System (HLS) concept, which allows tracking of current best estimates of the vehicle’s performance and assessment of its robustness against uncertainties. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as to model the HLS architecture. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. Range estimating — a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities — is then adapted with operational parameters as well as vehicle parameters in the DYREQT model to capture mission uncertainty alongside vehicle uncertainty. To perform the range estimation portion of this methodology, the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the set of uncertainty parameters. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo on the surrogate models. This probabilistic methodology was proved to provide insight into the underlying uncertainties of the three-element HLS architecture.

Stephanie Y Zhu↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Unmitigated uncertainties are known to have previously led to failed development programs; in order to combat these uncertainties, risks and their impacts must be understood and handled to ensure program success. In this paper, a probabilistic methodology to handle uncertainties is demonstrated on a three-element Human Landing System (HLS) concept, which allows tracking of current best estimates of the vehicle’s performance and assessment of its robustness against uncertainties. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as to model the HLS architecture. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. Range estimating — a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities — is then adapted with operational parameters as well as vehicle parameters in the DYREQT model to capture mission uncertainty alongside vehicle uncertainty. To perform the range estimation portion of this methodology, the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the set of uncertainty parameters. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo on the surrogate models. This probabilistic methodology was proved to provide insight into the underlying uncertainties of the three-element HLS architecture.

Stephanie Y. Zhu↗

Post-Landing Major Element Quantification Using SuperCam Laser Induced Breakdown Spectroscopy

The SuperCam instrument on the PerseveranceMars 2020 rover uses a pulsed 1064 nm laser to ablate targets at a distance and conduct laser induced breakdown spectroscopy (LIBS) by analyzing the light from the resulting plasma. SuperCam LIBS spectra are preprocessed to remove ambient light, noise, and the continuum signal present in LIBS observations. Prior to quantification, spectra are masked to remove noisier spectrometer regions andspectra are normalized to minimize signal fluctuations and effectsof target distance.In some cases, the spectra are also standardized or binned prior to quantification. To determine quantitative elemental compositionsof diverse geologic materials at Jezero crater, Mars, we use a suite of 1198 laboratory spectra of 334 well-characterized reference samples. The samples were selected to span a wide range of compositions and include typical silicate rocks, pure minerals (e.g.,silicates, sulfates, carbonates, oxides),more unusual compositions (e.g.,Mn oreand sodalite), andreplicates of the sintered SuperCam calibration targets (SCCTs) onboardthe rover. For each major element (SiO2, TiO2, Al2O3, FeOT, MgO, CaO, Na2O, K2O), the database was subdivided into five“folds” with similar distributions of the element of interest. One fold was held out as an independent test set, and the remaining fourfolds were used to optimize multivariate regression models relating the spectrum to the composition. We considered a variety of models, and selected several for further investigation for each element, based primarily on the root mean squared error of prediction (RMSEP) on the test set, when analyzed at 3m. In cases with several models of comparable performance at 3 m, we incorporated the SCCT performance at different distances to choose the preferred model. Shortly after landing on Mars and collecting initial spectra of geologic targets, we selected one model per element. Subsequently, with additional data from geologic targets, some models were revised to ensure results that are more consistent with geochemical constraints. The calibration discussed here is a snapshot of an ongoing effort to deliver the most accurate chemical compositions with SuperCam LIBS.

Mars 2020↗

Distributed Monitoring of the R(sup 2) Statistic for Linear Regression

The problem of monitoring a multivariate linear regression model is relevant in studying the evolving relationship between a set of input variables (features) and one or more dependent target variables. This problem becomes challenging for large scale data in a distributed computing environment when only a subset of instances is available at individual nodes and the local data changes frequently. Data centralization and periodic model recomputation can add high overhead to tasks like anomaly detection in such dynamic settings. Therefore, the goal is to develop techniques for monitoring and updating the model over the union of all nodes data in a communication-efficient fashion. Correctness guarantees on such techniques are also often highly desirable, especially in safety-critical application scenarios. In this paper we develop DReMo a distributed algorithm with very low resource overhead, for monitoring the quality of a regression model in terms of its coefficient of determination (R2 statistic). When the nodes collectively determine that R2 has dropped below a fixed threshold, the linear regression model is recomputed via a network-wide convergecast and the updated model is broadcast back to all nodes. We show empirically, using both synthetic and real data, that our proposed method is highly communication-efficient and scalable, and also provide theoretical guarantees on correctness.

Bhaduri, Kanishka↗

Expanding the application of soil moisture monitoring systems through regression-based transformation

Relative to other geophysical variables, soil moisture (SM) estimates derived from land surface models (LSMs) and land data assimilation systems (LDAS) are difficult to transfer between platforms and applications. This difficulty stems from the highly model-dependent nature of LSM SM estimates and differences in the vertical support of discretized SM values. As a result, operational SM estimates generated by one LSM (or LDAS) cannot generally be directly applied to a hydrologic monitoring or forecast system designed around a second LSM. This lack of transferability is particularly problematic for LDAS applications, where the time, expertise, and computational resources required to generate an operational LDAS analysis cannot be practically duplicated for every LSM-specific application. Here, we develop a set of simple regression tools for translating SM estimates between LSMs and multiple LDAS analyses. Results demonstrate that simple multivariate linear regression - utilizing independent variables based on multi-layer and temporally lagged SM estimates - can significantly improve upon baseline transformation approaches using direct percentile matching. The proposed regression approaches are effective for both the LSM-to-LSM and LDAS-to-LDAS transformation of multi-layer SM percentiles. Application of this approach will expand the utility of existing, high-quality (but LSM-specific) operational sources of SM information like the NASA Soil Moisture Active Passive Level-4 Soil Moisture product.

Soil Moisture↗

Determining the response of sea level to atmospheric pressure forcing using TOPEX/POSEIDON data

The static response of sea level to the forcing of atmospheric pressure, the so-called inverted barometer (IB) effect, is investigated using TOPEX/POSEIDON data. This response, characterized by the rise and fall of sea level to compensate for the change of atmospheric pressure at a rate of -1 cm/mbar, is not associated with any ocean currents and hence is normally treated as an error to be removed from sea level observation. Linear regression and spectral transfer function analyses are applied to sea level and pressure to examine the validity of the IB effect. In regions outside the tropics, the regression coefficient is found to be consistently close to the theoretical value except for the regions of western boundary currents, where the mesoscale variability interferes with the IB effect. The spectral transfer function shows near IB response at periods of 30 degrees is -0.84 +/- 0.29 cm/mbar (1 standard deviation). The deviation from = 1 cm /mbar is shown to be caused primarily by the effect of wind forcing on sea level, based on multivariate linear regression model involving both pressure and wind forcing. The regression coefficient for pressure resulting from the multivariate analysis is -0.96 +/- 0.32 cm/mbar. In the tropics the multivariate analysis fails because sea level in the tropics is primarily responding to remote wind forcing. However, after removing from the data the wind-forced sea level estimated by a dynamic model of the tropical Pacific, the pressure regression coefficient improves from -1.22 +/- 0.69 cm/mbar to -0.99 +/- 0.46 cm/mbar, clearly revealing an IB response. The result of the study suggests that with a proper removal of the effect of wind forcing the IB effect is valid in most of the open ocean at periods longer than 20 days and spatial scales larger than 500 km.

Fu, Lee-Lueng↗

Analysis of Multivariate Experimental Data Using A Simplified Regression Model Search Algorithm

A new regression model search algorithm was developed in 2011 that may be used to analyze both general multivariate experimental data sets and wind tunnel strain-gage balance calibration data. The new algorithm is a simplified version of a more complex search algorithm that was originally developed at the NASA Ames Balance Calibration Laboratory. The new algorithm has the advantage that it needs only about one tenth of the original algorithm's CPU time for the completion of a search. In addition, extensive testing showed that the prediction accuracy of math models obtained from the simplified algorithm is similar to the prediction accuracy of math models obtained from the original algorithm. The simplified algorithm, however, cannot guarantee that search constraints related to a set of statistical quality requirements are always satisfied in the optimized regression models. Therefore, the simplified search algorithm is not intended to replace the original search algorithm. Instead, it may be used to generate an alternate optimized regression model of experimental data whenever the application of the original search algorithm either fails or requires too much CPU time. Data from a machine calibration of NASA's MK40 force balance is used to illustrate the application of the new regression model search algorithm.

multivariate experimental data↗

The saturated activities of Na-22, Mn-54, and Co-56 and the depth of sampling of soils

A multivariate linear regression study shows that the saturated activities of Na-22, Co-56, and Mn-54 are functions of the chemical composition and the reciprocal of the cube root of the sample weight. A procedure for taking the temporal variation in these short-lived radionuclides into account is described. The regression results are used to estimate the average depths of sampling of soil material. The data indicate that all 28 samples were in the same conditions of exposure when recovered as they had been for the last approximately 10 years. The saturation activity of Mn-54 is a very weak function of time.

Keith, J. E.↗

Synthesis of rotor test data for real-time simulation

A mathematical model of a hingeless tilting rotor is presented. The model was obtained by a systematic curve fit procedure applied to an extensive set of model scale wind tunnel data. The math model equations were used in a real time flight simulation model of a hingeless tilt rotor XV-15 to assess changes in flying qualities compared to those obtained using a previous rotor model. Extensive plots of the rotor derivatives are given. Discussions of attempts to apply multivariable linear regression technqiues to the data and the use of an analytical rotor representation are included.

Mcveigh, M. A.↗

Soil temperature extrema recovery rates after precipitation cooling

From a one dimensional view of temperature alone variations at the Earth's surface manifest themselves in two cyclic patterns of diurnal and annual periods, due principally to the effects of diurnal and seasonal changes in solar heating as well as gains and losses of available moisture. Beside these two well known cyclic patterns, a third cycle has been identified which occurs in values of diurnal maxima and minima soil temperature extrema at 10 cm depth usually over a mesoscale period of roughly 3 to 14 days. This mesoscale period cycle starts with precipitation cooling of soil and is followed by a power curve temperature recovery. The temperature recovery clearly depends on solar heating of the soil with an increased soil moisture content from precipitation combined with evaporation cooling at soil temperatures lowered by precipitation cooling, but is quite regular and universal for vastly different geographical locations, and soil types and structures. The regularity of the power curve recovery allows a predictive model approach over the recovery period. Multivariable linear regression models alloy predictions of both the power of the temperature recovery curve as well as the total temperature recovery amplitude of the mesoscale temperature recovery, from data available one day after the temperature recovery begins.

Welker, J. E.↗

Regional climate change predictions from the Goddard Institute for Space Studies high resolution GCM

A new diagnostic tool is developed for examining relationships between the synoptic scale circulation and regional temperature distributions in GCMs. The 4 x 5 deg GISS GCM is shown to produce accurate simulations of the variance in the synoptic scale sea level pressure distribution over the U.S. An analysis of the observational data set from the National Meteorological Center (NMC) also shows a strong relationship between the synoptic circulation and grid point temperatures. This relationship is demonstrated by deriving transfer functions between a time-series of circulation parameters and temperatures at individual grid points. The circulation parameters are derived using rotated principal components analysis, and the temperature transfer functions are based on multivariate polynomial regression models. The application of these transfer functions to the GCM circulation indicates that there is considerable spatial bias present in the GCM temperature distributions. The transfer functions are also used to indicate the possible changes in U.S. regional temperatures that could result from differences in synoptic scale circulation between a 1XCO2 and a 2xCO2 climate, using a doubled CO2 version of the same GISS GCM.

Crane, Robert G.↗

Modeling and control for closed environment plant production systems

A computer program was developed to study multiple crop production and control in controlled environment plant production systems. The program simulates crop growth and development under nominal and off-nominal environments. Time-series crop models for wheat (Triticum aestivum), soybean (Glycine max), and white potato (Solanum tuberosum) are integrated with a model-based predictive controller. The controller evaluates and compensates for effects of environmental disturbances on crop production scheduling. The crop models consist of a set of nonlinear polynomial equations, six for each crop, developed using multivariate polynomial regression (MPR). Simulated data from DSSAT crop models, previously modified for crop production in controlled environments with hydroponics under elevated atmospheric carbon dioxide concentration, were used for the MPR fitting. The model-based predictive controller adjusts light intensity, air temperature, and carbon dioxide concentration set points in response to environmental perturbations. Control signals are determined from minimization of a cost function, which is based on the weighted control effort and squared-error between the system response and desired reference signal.

NASA Discipline Life Support Systems↗

Usefulness of microvolt T-wave alternans for prediction of ventricular tachyarrhythmic events in patients with dilated cardiomyopathy: results from a prospective observational study

OBJECTIVES: This study was designed to evaluate the ability of microvolt-level T-wave alternans (MTWA) to identify prospectively patients with idiopathic dilated cardiomyopathy (DCM) at risk of ventricular tachyarrhythmic events and to compare its predictive accuracy with that of conventional risk stratifiers. BACKGROUND: Patients with DCM are at increased risk of sudden death from ventricular tachyarrhythmias. At present, there are no established methods of assessing this risk. METHODS: A total of 137 patients with DCM underwent risk stratification through assessment of MTWA, left ventricular ejection fraction, baroreflex sensitivity (BRS), heart rate variability, presence of nonsustained ventricular tachycardia (VT), signal-averaged electrocardiogram, and presence of intraventricular conduction defect. The study end point was either sudden death, resuscitated ventricular fibrillation, or documented hemodynamically unstable VT. RESULTS: During an average follow-up of 14 +/- 6 months, MTWA and BRS were significant univariate predictors of ventricular tachyarrhythmic events (p < 0.035 and p < 0.015, respectively). Multivariate Cox regression analysis revealed that only MTWA was a significant predictor. CONCLUSIONS: Microvolt-level T-wave alternans is a powerful independent predictor of ventricular tachyarrhythmic events in patients with DCM.

Evaluation Studies↗

T wave alternans as a predictor of recurrent ventricular tachyarrhythmias in ICD recipients: prospective comparison with conventional risk markers

INTRODUCTION: The current standard for arrhythmic risk stratification is electrophysiologic (EP) testing, which, due to its invasive nature, is limited to patients already known to be at high risk. A number of noninvasive tests, such as determination of left ventricular ejection fraction (LVEF) or heart rate variability, have been evaluated as additional risk stratifiers. Microvolt T wave alternans (TWA) is a promising new risk marker. Prospective evaluation of noninvasive risk markers in low- or moderate-risk populations requires studies involving very large numbers of patients, and in such studies, documentation of the occurrence of ventricular tachyarrhythmias is difficult. In the present study, we identified a high-risk population, recipients of an implantable cardioverter defibrillator (ICD), and prospectively compared microvolt TWA with invasive EP testing and other risk markers with respect to their ability to predict recurrence of ventricular tachyarrhythmias as documented by ICD electrograms. METHODS AND RESULTS: Ninety-five patients with a history of ventricular tachyarrhythmias undergoing implantation of an ICD underwent EP testing, assessment of TWA, as well as determination of LVEF, baroreflex sensitivity, signal-averaged ECG, analysis of 24-hour Holter monitoring, and QT dispersion from the 12-lead surface ECG. The endpoint of the study was first appropriate ICD therapy for electrogram-documented ventricular fibrillation or tachycardia during follow-up. Kaplan-Meier survival analysis revealed that TWA (P < 0.006) and LVEF (P < 0.04) were the only significant univariate risk stratifiers. EP testing was not statistically significant (P < 0.2). Multivariate Cox regression analysis revealed that TWA was the only statistically significant independent risk factor. CONCLUSIONS: Measurement of microvolt TWA compared favorably with both invasive EP testing and other currently used noninvasive risk assessment methods in predicting recurrence of ventricular tachyarrhythmias in ICD recipients. This study suggests that TWA might also be a powerful tool for risk stratification in low- or moderate-risk patients, and needs to be prospectively evaluated in such populations.

Non-NASA Center↗

Diurnal Variability of Vertical Structure from a TRMM Passive Microwave "Virtual Radar" Retrieval

Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Boccippio, Dennis J.↗

An Update on Mortality in the U.S. Astronaut Corps: 1959-2009

Although it has now been over 50 years since mankind first ventured into space, the long-term health impacts of human space flight remain largely unknown. Identifying factors that affect survival and prognosis among those who participate in space flight is vitally important, as the era of commercial space flight approaches and NASA prepares for missions to Mars. The Longitudinal Study of Astronaut Health is a prospective study designed to examine trends in astronaut morbidity and mortality. The purpose of this analysis was to describe and explore predictors of overall and cause-specific mortality among individuals selected for the U.S. astronaut corps. All U.S. astronauts (n=321), regardless of flight status, were included in this analysis. Death certificate searches were conducted to ascertain vital status and cause of death through April 2009. Data were collected from medical records and lifestyle questionnaires. Multivariable Cox regression modeling was used to calculate the mortality hazard associated with embarking on space flight, adjusted for sex, race, and age at selection. Between 1959 and 2009, there were 39 (12.1%) deaths. Of these deaths, 18 (42.2%) were due to occupational accidents; 7 (17.9%) were due to other accidents; 6 (15.4%) were attributable to cancer; 6 (15.4%) resulted from cardiovascular/circulatory diseases; and 2 (5.1%) were from other causes. Participation in space flight did not significantly increase mortality hazard over time (adjusted hazard ratio=0.57; 95% confidence interval=0.26-1.26. Because our results are based on a small sample size, future research that includes payload specialists, other space flight participants, and international crew members is warranted to maximize statistical power.

Amirian, E.↗

Physical Interpretation of the Correlation Between Multi-Angle Spectral Data and Canopy Height

Recent empirical studies have shown that multi-angle spectral data can be useful for predicting canopy height, but the physical reason for this correlation was not understood. We follow the concept of canopy spectral invariants, specifically escape probability, to gain insight into the observed correlation. Airborne Multi-Angle Imaging Spectrometer (AirMISR) and airborne Laser Vegetation Imaging Sensor (LVIS) data acquired during a NASA Terrestrial Ecology Program aircraft campaign underlie our analysis. Two multivariate linear regression models were developed to estimate LVIS height measures from 28 AirMISR multi-angle spectral reflectances and from the spectrally invariant escape probability at 7 AirMISR view angles. Both models achieved nearly the same accuracy, suggesting that canopy spectral invariant theory can explain the observed correlation. We hypothesize that the escape probability is sensitive to the aspect ratio (crown diameter to crown height). The multi-angle spectral data alone therefore may not provide enough information to retrieve canopy height globally

Schull, M. A.↗

Analysis of Multivariate Experimental Data Using A Simplified Regression Model Search Algorithm

A new regression model search algorithm was developed that may be applied to both general multivariate experimental data sets and wind tunnel strain-gage balance calibration data. The algorithm is a simplified version of a more complex algorithm that was originally developed for the NASA Ames Balance Calibration Laboratory. The new algorithm performs regression model term reduction to prevent overfitting of data. It has the advantage that it needs only about one tenth of the original algorithm's CPU time for the completion of a regression model search. In addition, extensive testing showed that the prediction accuracy of math models obtained from the simplified algorithm is similar to the prediction accuracy of math models obtained from the original algorithm. The simplified algorithm, however, cannot guarantee that search constraints related to a set of statistical quality requirements are always satisfied in the optimized regression model. Therefore, the simplified algorithm is not intended to replace the original algorithm. Instead, it may be used to generate an alternate optimized regression model of experimental data whenever the application of the original search algorithm fails or requires too much CPU time. Data from a machine calibration of NASA's MK40 force balance is used to illustrate the application of the new search algorithm.

Ulbrich, Norbert M.↗