Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonparametric”

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 307 records · Page 17

DESI Survey Validation Spectra Reveal an Increasing Fraction of Recently Quenched Galaxies at z ∼ 1

Abstract We utilize ∼17,000 bright luminous red galaxies (LRGs) from the novel Dark Energy Spectroscopic Instrument Survey Validation spectroscopic sample, leveraging its deep (∼2.5 hr galaxy −1 exposure time) spectra to characterize the contribution of recently quenched galaxies to the massive galaxy population at 0.4 < z < 1.3. We use Prospector to infer nonparametric star formation histories and identify a significant population of recently quenched galaxies that have joined the quiescent population within the past ∼1 Gyr. The highest-redshift subset (277 at z > 1) of our sample of recently quenched galaxies represents the largest spectroscopic sample of post-starburst galaxies at that epoch. At 0.4 < z < 0.8, we measure the number density of quiescent LRGs, finding that recently quenched galaxies constitute a growing fraction of the massive galaxy population with increasing look-back time. Finally, we quantify the importance of this population among massive ( log ( M ⋆ / M ⊙ ) > 11.2) LRGs by measuring the fraction of stellar mass each galaxy formed in the gigayear before observation, f 1 Gyr . Although galaxies with f 1 Gyr > 0.1 are rare at z ∼ 0.4 (≲0.5% of the population), by z ∼ 0.8, they constitute ∼3% of massive galaxies. Relaxing this threshold, we find that galaxies with f 1 Gyr > 5% constitute ∼10% of the massive galaxy population at z ∼ 0.8. We also identify a small but significant sample of galaxies at z = 1.1–1.3 that formed with f 1 Gyr > 50%, implying that they may be analogs to high-redshift quiescent galaxies that formed on similar timescales. Future analysis of this unprecedented sample promises to illuminate the physical mechanisms that drive the quenching of massive galaxies after cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

A consistent dataset for the net income distribution for 190 countries and aggregated to 32 geographical regions from 1958 to 2015

Abstract. Data on income distributions within and across countries are becoming increasingly important for informing analysis of income inequality and understanding the distributional consequences of climate change. While datasets on income distribution collected from household surveys are available for multiple countries, these datasets often do not represent the same concept of inequality (or income concept) and therefore make comparisons across countries, over time and across datasets difficult. Here, we present a consistent dataset of income distributions across 190 countries from 1958 to 2015 measured in terms of net income. We complement the observed values in this dataset with values imputed from a summary measure of the income distribution, specifically the Gini coefficient. For the imputation, we use a recently developed nonparametric principal-component-based approach that shows an excellent fit to data on income distributions compared to other approaches. We also present another version of this dataset aggregated from the country level to 32 geographical regions. Our dataset is developed for the purpose of calibrating models such as integrated human–Earth system models with detailed data on income distributions. This dataset will enable more robust analysis of income distribution at multiple scales. The latest version of our data are available on Zenodo: https://doi.org/10.5281/zenodo.7093997 (Narayan et al., 2022b).

97 MATHEMATICS AND COMPUTING↗

Quantifying wildfire drivers and predictability in boreal peatlands using a two-step error-correcting machine learning framework in TeFire v1.0

Abstract. Wildfires are becoming an increasing challenge to the sustainability of boreal peatland (BP) ecosystems and can alter the stability of boreal carbon storage. However, predicting the occurrence of rare and extreme BP fires proves to be challenging, and gaining a quantitative understanding of the factors, both natural and anthropogenic, inducing BP fires remains elusive. Here, we quantified the predictability of BP fires and their primary controlling factors from 1997 to 2015 using a two-step correcting machine learning (ML) framework that combines multiple ML classifiers, regression models, and an error-correcting technique. We found that (1) the adopted oversampling algorithm effectively addressed the unbalanced data and improved the recall rate by 26.88 %–48.62 % when using multiple datasets, and the error-correcting technique tackled the overestimation of fire sizes during fire seasons; (2) nonparametric models outperformed parametric models in predicting fire occurrences, and the random forest machine learning model performed the best, with the area under the receiver operating characteristic curve ranging from 0.83 to 0.93 across multiple fire datasets; and (3) four sets of factor-control simulations consistently indicated the dominant role of temperature, air dryness, and climate extreme (i.e., frost) for boreal peatland fires, overriding the effects of precipitation, wind speed, and human activities. Our findings demonstrate the efficiency and accuracy of ML techniques in predicting rare and extreme fire events and disentangle the primary factors determining BP fires, which are critical for predicting future fire risks under climate change.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Model Form Uncertainty of Epidemic Forecasting Models from Incidence Data

Objective We present a mathematical framework for non-parametric estimation of the force of infection, together with statistical upper and lower confidence bands. The resulting estimates allow to assess how well simpler models, such as SEIR, fit the observed time series of incidence data.IntroductionUncertainty Quantification (UQ), the ability to quantify the impact of sample-to-sample variations and model misspecification on predictions and forecasts, is a critical aspect of disease surveillance. While quantifying the impact of stochastic uncertainty in the data is well understood, quantifying the impact of model misspecification is significantly harder. For the latter, one needs a "universal model" to which more restrictive parametric models are compared too.MethodsThis talk presents a useful modeling framework for time series of incidence data from contagious diseases that enables one to identify and quantify the impact of model form uncertainty. Specifically, we propose to focus on estimating the timedependent force of infection. The latter is a universal parameters for all contagious disease model. Using a machine learning technique for estimating monotone functions, i.e., isotonic regression and its variants, one can estimate the force ofinfection without addtional assumptions. We note that most contagious disease model do satisfy this monotonicity assumption, due to a combination of factors: depletion of susceptibles, implementation of mitigation strategies, behavior change, etc. Comparing the resulting "non-parametric" estimate with parametric estimates, obtained by fitting an SEIR for example, can reveal model deficiencies and help quantify model form uncertainties.Finally, we discuss how ideas from "strict bound theory" can be used to develop upper and lower uncertainty bands for force of infection that acknowledge the intrinsic stochasticity in the data.ResultsWe demonstrate the application of the methodology to weekly Influenza Like Illness (ILI) incidence data from France andcompare the results to fitted SIR and SEIR models. This comparison can be seen as a nonparametric goodness of fit test, providing one with tools to do simple model selection.ConclusionsWe present a novel and flexible model to statistically describe the force of infection as a function of time. Comparing the fit to incidence data of that model with the fit of simpler parametric models enables the quantification of model form uncertainty and associated model selection.

Hengartner, Nicholas↗

Soda Pop: A Time-Series Clustering, Alarming and Disease Forecasting Application

To introduce Soda Pop, an R/Shiny application designed to be a disease agnostic time-series clustering, alarming, and forecasting tool to assist in disease surveillance “triage, analysis and reporting” workflows within the Biosurveillance Ecosystem (BSVE). In this poster, we highlight the new capabilities that are brought to the BSVE by Soda Pop with an emphasis on the impact of methodological decisions. The Biosurveillance Ecosystem (BSVE) is a biological and chemical threat surveillance system sponsored by the Defense Threat Reduction Agency (DTRA). BSVE is intended to be user-friendly, multi-agency, cooperative, modular and threat agnostic platform for biosurveillance. In BSVE, a web-based workbench presents the analyst with applications (apps) developed by various DTRA funded researchers, which are deployed on-demand in the cloud(e.g., Amazon Web Services). These apps aim to address emerging needs and refine capabilities to enable early warning of chemical and biological threats for multiple users across local, state, and federal agencies. Soda Pop is an app developed by Pacific Northwest National Laboratory (PNNL) to meet the current needs of the BSVE for early warning and detection of disease outbreaks. Aimed for use bya diverse set of analysts, the application is agnostic to data source and spatial scale enabling it to be generalizable across many diseases and locations. To achieve this, we placed a particular emphasis on clustering and alerting of disease signals within Soda Pop without strong prior assumptions on the nature of observed diseased counts. Although designed to be agnostic to the data source, Soda Pop was initially developed and tested on data summarizing Influenza-Like Illness in military hospitals from collaboration with the Armed Forces Health Surveillance Branch. Currently, the data incorporated also includes the CDC’s National Notifiable Diseases Surveillance System(NNDSS) tables and the WHO’s Influenza A/B Influenza Data(Flunet). These data sources are now present in BSVE’s Postgres data storage for direct access. Soda Pop is designed to automate time-series tasks of data summarization, exploration, clustering, alarming and forecasting. Built as an R/Shiny application, Soda Pop is founded on the powerful statistical tool R.. Where applicable, Soda Pop facilitates nonparametric seasonal decomposition of time-series; hierarchical agglomerative clustering across reporting areas and between diseases within reporting areas; and a variety of alarming techniques including Exponential Weighted Moving Average alarms and Early Aberration Detection. Soda Pop embeds these techniques within a user-interface designed to enhance an analyst’s understanding of emerging trends in their data and enables the inclusion of its graphical elements into their dossier for further tracking and reporting. The ultimate goal of this software is to facilitate the discovery of unknown disease signals along with increasing the speed of detection of unusual patterns within these signals. Soda Pop organizes common statistical disease surveillance tasks in a manner integrated with BSVE data source inputs and outputs. The app analyzes time-series disease data and supports a robust set of clustering and alarming routines that avoid strong assumptions on the nature of observed disease counts. This attribute allows for flexibility in the data source, spatial scale, and disease types making it useful to a wide range of analysts.

97 MATHEMATICS AND COMPUTING↗

Code associated with "Assessing the Effect of a Deep-Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights from a Multi-Site US Study"

R code associated with "Assessing the Effect of a Deep-Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights from a Multi-Site US Study" These scripts can be used to analyze the accompanying dataset: https://doi.org/10.5281/zenodo.16620529. (1) resampling_functions.R : Creates functions for resampling data for nonparametric bootstrapping (2) bootstrap_depth_profiles.R : Applies resampling functions to the dataset (3) gapfill_masscalcs.R : Prepares data for equivalent mass calculations (4) get_14C_EM.R : Applies soilR to obtain the 14C end member for mixing model calculations (5) mass_aggregate.R : Performs equivalent mass and mass weighted averaging calculations.

Carbon↗

Quantile-dependent expressivity of plasma adiponectin concentrations may explain its sex-specific heritability, gene-environment interactions, and genotype-specific response to postprandial lipemia

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g. adiponectin) is high or low relative to its distribution. We have previously shown that the heritability ( h 2 ) of adiposity, lipoproteins, postprandial lipemia, pulmonary function, and coffee and alcohol consumption are quantile-specific. Whether adiponectin heritability is quantile specific remains to be determined. Plasma adiponectin concentrations from 4,182 offspring-parent pairs and 1,662 sibships from the Framingham Heart Study were analyzed. Quantile-specific heritability from offspring-parent (β OP , h 2 = 2 β OP /(1 + r spouse )) and full-sib regression slopes ( β FS , h 2 = {(1 + 8r spouse β FS ) 0.05 -1}/(2r spouse )) were robustly estimated by quantile regression with nonparametric significance assigned from 1,000 bootstrap samples. Quantile-specific h 2 (± SE) increased with increasing percentiles of the offspring’s age- and sex-adjusted adiponectin distribution when estimated from β OP ( P trend = 2.2 × 10 -6 ): 0.30 ± 0.03 at the 10th, 0.33 ± 0.04 at the 25th, 0.43 ± 0.04 at the 50th, 0.55 ± 0.05 at the 75th, and 0.57 ± 0.08 at the 90th percentile, and when estimated from β FS (P trend = 7.6 × 10 -7 ): 0.42 ± 0.03 at the 10th, 0.44 ± 0.04 at the 25th, 0.56 ± 0.05 at the 50th, 0.73 ± 0.08 at the 75th, and 0.79 ± 0.11 at the 90th percentile. Consistent with quantile-dependent expressivity, adiponectin’s: (1) heritability was greater in women in accordance with their higher adiponection concentrations; (2) relationships to ADIPOQ polymorphisms were modified by adiposity in accordance with its adiponectin-lowering effect; (3) response to rosiglitazone was predicted by the 45T> G ADIPOQ polymorphism; (4) difference by ADIPOQ haplotypes increased linearly with increasing postprandial adiponectin concentrations. Adiponectin heritability is quantile dependent, which may explain sex-specific heritability, gene-environment and gene-drug interactions, and postprandial response by haplotypes.

Postprandial↗

Quantile-dependent expressivity of serum C-reactive protein concentrations in family sets

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., C-reactive protein, CRP) is high or low relative to its distribution. We have previously shown that the heritabilities ( h 2 ) of coffee and alcohol consumption, postprandial lipemia, lipoproteins, leptin, adiponectin, adiposity, and pulmonary function are quantile-specific. Whether CRP heritability is quantile-specific is currently unknown. Methods Serum CRP concentrations from 2,036 sibships and 6,144 offspring-parent pairs were analyzed from the Framingham Heart Study. Quantile-specific heritability from full-sib (β FS , h 2 ={(1 + 8r spouse β FS ) 0.5 - 1}/(2r spouse )) and offspring-parent regression slopes (β OP , h 2 = 2β OP /(1 + r spouse )) were estimated robustly by quantile regression with nonparametric significance determined from 1,000 bootstrap samples. Results Quantile-specific h 2 (±SE) increased with increasing percentiles of the offspring’s age- and sex-adjusted CRP distribution when estimated from β OP ( P trend = 0.0004): 0.02 ± 0.01 at the 10th, 0.04 ± 0.01 at the 25th, 0.10 ± 0.02 at the 50th, 0.20 ± 0.05 at the 75th, and 0.33 ± 0.10 at the 90th percentile, and when estimated from β FS ( P trend = 0.0008): 0.03±0.01 at the 10th, 0.06 ± 0.02 at the 25th, 0.14 ± 0.03 at the 50th, 0.24 ± 0.05 at the 75th, and 0.53 ± 0.21 at the 90th percentile. Conclusion Heritability of serum CRP concentration is quantile-specific, which may explain or contribute to the inflated CRP differences between CRP (rs1130864, rs1205, rs1800947, rs2794521, rs3091244), FGB (rs1800787), IL-6 (rs1800795, rs1800796), IL6R (rs8192284), TNF- α (rs1800629) and APOE genotypes following CABG surgery, stroke, TIA, curative esophagectomy, intensive periodontal therapy, or acute exercise; during acute coronary syndrome or Staphylococcus aureus bacteremia; or in patients with chronic rheumatoid arthritis, diabetes, peripheral arterial disease, ankylosing spondylitis, obesity or inflammatory bowel disease or who smoke.

59 BASIC BIOLOGICAL SCIENCES↗

Spirometric traits show quantile-dependent heritability, which may contribute to their gene-environment interactions with smoking and pollution

Background: “Quantile-dependent expressivity” refers to a genetic effect that is dependent upon whether the phenotype (e.g., spirometric data) is high or low relative to its population distribution. Forced vital capacity (FVC), forced expiratory volume in 1 second (FEV 1 ), and the FEV 1 /FVC ratio are moderately heritable spirometric traits. The aim of the analyses is to test whether their heritability (h 2 ) is constant over all quantiles of their distribution. Methods: Quantile regression was applied to the mean age, sex, height and smoking-adjusted spirometric data over multiple visits in 9,993 offspring-parent pairs and 1,930 sibships from the Framingham Heart Study to obtain robust estimates of offspring-parent (β OP ), offspring-midparent (β OM ), and full-sib regression slopes (β FS ). Nonparametric significance levels were obtained from 1,000 bootstrap samples. β OP s were used as simple indicators of quantile-specific heritability (i.e.,h 2 = 2β OP /(1+r spouse ), where r spouse was the correlation between spouses). Results: β OP ± standard error (SE) decreased by 0.0009 ± 0.0003 (P = 0.003) with every one-percent increment in the population distribution of FEV 1 /FVC, i.e., β OP ± SE were: 0.182 ± 0.031, 0.152 ± 0.015; 0.136 ± 0.011; 0.121 ± 0.013; and 0.099 ± 0.013 at the 10th, 25th, 50th, 75th, and 90th percentiles of the FEV 1 /FVC distribution, respectively. These correspond to h 2 ± SEs of 0.350 ± 0.060 at the 10th, 0.292 ± 0.029 at the 25th, 0.262 ± 0.020 at the 50th, 0.234 ± 0.025 at the 75th, and 0.191 ± 0.025 at the 90th percentiles of the FEV 1 /FVC ratio. Maximum mid-expiratory flow (MMEF) h 2 ± SEs increased 0.0025 ± 0.0007 (P = 0.0004) with every one-percent increment in its distribution, i.e.: 0.467 ± 0.046, 0.467 ± 0.033, 0.554 ± 0.038, 0.615 ± 0.042, and 0.675 ± 0.060 at the 10th, 25th, 50th, 75th, and 90th percentiles of its distribution. This was due to forced expiratory flow at 75% of FVC (FEF75%), whose quantile-specific h 2 increased an average of 0.0042 ± 0.0008 for every one-percent increment in its distribution. It is speculated that previously reported gene-environment interactions may be partially attributable to quantile-specific h 2 , i.e., greater heritability in individuals with lower FEV 1 /FVC due to smoking or airborne particles exposure vs. nonsmoking, unexposed individuals. Conclusion: Heritabilities of FEV 1 /FVC, MMEF, and FEF75% from quantile-regression of offspring-parent and sibling spirometric data suggest their quantile-dependent expressivity.

59 BASIC BIOLOGICAL SCIENCES↗

Relationship between vegetation reflectance spectra and soil geochemistry: New data from Catheart Mountain, Maine

Spectral reflectance curves for red spruce and balsam fir in the 350- to 1100-nanometer region showed differences between specimens growing over and away from a concealed copper deposit at Catheart Mountain, Maine. A sample group of 15 anomalous spruce, 9 background spruce, 10 anomalous fir, and 10 background fir was studied. Measurements of copper and molybdenum contents in the supporting soil were used to categorize the trees as anomalous or background. Analyses using parametric and nonparametric statistical tests were performed to establish whether significant differences existed between the anomalous and background groups. At the 95% confidence level, significant differences were found to exist for both species, generally in the chlorophyll band centered at about 550 nanometers and in the region from 700 to 900 nanometers.

Canney, F. C.↗

Permutation codes for sources.

Source encoding techniques based on permutation codes are investigated. For a broad class of distortion measures it is shown that optimum encoding of a source permutation code is easy to instrument even for very long block lengths. Also, the nonparametric nature of permutation encoding is well suited to situations involving unknown source statistics. For the squared-error distortion measure a procedure for generating good permutation codes of a given rate and block length is described. The performance of such codes for a memoryless Gaussian source is compared both with the rate-distortion function bound and with the performance of various quantization schemes. The comparison reveals that permutation codes are asymptotically ideal for small rates and perform as well as the best entropy-coded quantizers presently known for intermediate rates. They can be made to compare favorably at high rates, too, provided the coding delay associated with extremely long block lengths is tolerable.

Berger, T.↗

Nearest Neighbor Algorithms for Pattern Classification

A solution of the discrimination problem is considered by means of the minimum distance classifier, commonly referred to as the nearest neighbor (NN) rule. The NN rule is nonparametric, or distribution free, in the sense that it does not depend on any assumptions about the underlying statistics for its application. The k-NN rule is a procedure that assigns an observation vector z to a category F if most of the k nearby observations x sub i are elements of F. The condensed nearest neighbor (CNN) rule may be used to reduce the size of the training set required categorize The Bayes risk serves merely as a reference-the limit of excellence beyond which it is not possible to go. The NN rule is bounded below by the Bayes risk and above by twice the Bayes risk.

Barrios, J. O.↗

Phase 1 of the earth resources data analysis program

Research completed in the Earth Resources Data Analysis Program is discussed along with recommendations for future study. Projects discussed include use of the Cholesky decomposition in feature selection and classification algorithms; optimal feature selection and extraction, probability density estimation and nonparametric classifiers; use of spatial information in classification; and model for crop row reflectance. The installation of LARSYS on the ICSA's IBM 370/155 is discussed, and a list of technical reports is included.

Source record↗