Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multivariate analysis”

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

Punzi-loss: a non-differentiable metric approximation for sensitivity optimisation in the search for new particles

We present the novel implementation of a non-differentiable metric approximation and a corresponding loss-scheduling aimed at the search for new particles of unknown mass in high energy physics experiments. We call the loss-scheduling, based on the minimisation of a figure-of-merit related function typical of particle physics, a Punzi-loss function, and the neural network that utilises this loss function a Punzi-net. We show that the Punzi-net outperforms standard multivariate analysis techniques and generalises well to mass hypotheses for which it was not trained. This is achieved by training a single classifier that provides a coherent and optimal classification of all signal hypotheses over the whole search space. Our result constitutes a complementary approach to fully differentiable analyses in particle physics. We implemented this work using PyTorch and provide users full access to a public repository containing all the codes and a training example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measuring the electron Yukawa coupling via resonant s-channel Higgs production at FCC-ee

The Future Circular Collider (FCC-ee) offers the unique opportunity of studying the Higgs Yukawa coupling to the electron, \(y_\mathrm {e}\) , via resonant s -channel production, \(\mathrm {e^+e^-}\rightarrow \mathrm {H}\) , in a dedicated run at \(\sqrt{s} = m_\mathrm {H}\) . The signature for direct Higgs production is a small rise in the cross sections for particular final states, consistent with Higgs decays, over the expectations for their occurrence due to Standard Model (SM) background processes involving \(\mathrm {Z}^*\) , \(\gamma ^*\) , or t -channel exchanges alone. Performing such a measurement is remarkably challenging for four main reasons. First, the low value of the e \(^\pm \) mass leads to a tiny \(y_\mathrm {e}\) coupling and correspondingly small cross section: \(\sigma _\mathrm {ee\rightarrow H} \propto m_\mathrm {e}^2 = 0.57\) fb accounting for initial-state \(\gamma \) radiation. Second, the \(\mathrm {e^+e^-}\) beams must be monochromatized such that the spread of their centre-of-mass (c.m.) energy is commensurate with the narrow width of the SM Higgs boson, \(\varGamma _\mathrm {H} = 4.1\) MeV, while keeping large beam luminosities. Third, the Higgs mass must also be known beforehand with a few-MeV accuracy in order to operate the collider at the resonance peak, \(\sqrt{s} = m_\mathrm {H}\) . Last but not least, the cross sections of the background processes are many orders-of-magnitude larger than those of the Higgs decay signals. A preliminary generator-level study of 11 Higgs decay channels using a multivariate analysis, which exploits boosted decision trees to discriminate signal and background events, identifies two final states as the most promising ones in terms of statistical significance: \(\mathrm {H}\rightarrow gg\) and \(\mathrm {H}\rightarrow \mathrm {W}\mathrm {W}^*\!\rightarrow \ell \nu \) + 2 jets. For a benchmark monochromatization with 4.1-MeV c.m. energy spread (leading to \(\sigma _\mathrm {ee\rightarrow H} = 0.28\) fb) and 10 ab \(^{-1}\) of integrated luminosity, a \(1.3\sigma \) signal significance can be reached, corresponding to an upper limit on the e \(^\pm \) Yukawa coupling at 1.6 times the SM value: |y e | < 1.6|y e SM | at 95% confidence level, per FCC-ee interaction point per year. Directions for future improvements of the study are outlined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Association of Serum Triglyceride Levels with Severity in Acute Pancreatitis: Results from an International, Multicenter Cohort Study

Hypertriglyceridemia (HTG) is considered within the top 5 etiologies in acute pancreatitis (AP), but the association of serum triglyceride (TG) levels with the clinical course of AP remains controversial. This study aims to examine the effect of TG levels on severity of AP. Patients were enrolled prospectively through APPRENTICE. High TG levels were defined based on the Endocrine Society Clinical Practice Guidelines. HTG was categorized as mild (serum TG levels 150–199 mg/dL), moderate (200–999 mg/dL), severe (1,000–1,999 mg/dL), and very severe (≥2,000 mg/dL). Severity of AP was based on the revised Atlanta classification criteria. Additionally, early TG levels were measured in 764 subjects and found elevated in 342 (120 with mild; 176, moderate; and 46, severe/very severe HTG). Patients with increased TG levels were younger (age ≥60, 16.7 vs. 30.3%), more likely to be male (66.1 vs. 51.2%), with more frequent alcohol use (62.8 vs. 50.7%), and diabetes mellitus (30.2 vs. 12.3%; all p ≤ 0.005). Severe AP (24.9 vs. 10.0%), ICU admission (32.5 vs. 19.7%), and mortality (5.3 vs. 1.7%; all p ≤ 0.005) were more frequently seen in patients with elevated TG levels. Based on multivariable analysis, elevated TG levels were independently associated with severe AP ( p < 0.05). This large multicenter study confirms that elevated TG levels are associated with severe disease regardless of AP etiology.

60 APPLIED LIFE SCIENCES↗

Transcatheter Arterial Chemoembolization for Treatment-Naive Hepatocellular Carcinoma Has Different Treatment Effects Depending on Central or Peripheral Tumor Location

Introduction: The purpose of this study was to evaluate the treatment efficacy of transcatheter arterial chemoembolization (TACE) for treatment-naive hepatocellular carcinoma (HCC) according to tumor location and burden. Methods: Between 2010 and 2019, consecutive patients who underwent TACE as the first treatment were enrolled. Tumors were classified into two categories based on their location, as central or peripheral tumors. Tumors in the central zone, which is within 1 cm of the main trunk or the first branch of the portal vein, were classified as central tumors, while those located in the peripheral zone were classified as peripheral tumors. Patients were grouped according to the HCC location and up-to-7 criteria. Patients with central tumors were classified into the central arm and those with only peripheral tumors were classified into the peripheral arm. Patients within and beyond the up-to-7 criteria were classified into the up-to-7 in and up-to-7 out-groups, respectively. Local recurrence-free survival (LRFS) and progression-free survival (PFS) were compared per nodule (central tumor vs. peripheral tumor) and per patient (central arm vs. peripheral arm), respectively. The prognostic factors of LRFS and PFS were analyzed by univariate and multivariate analyses. Results: A total of 174 treatment-naive patients with 352 HCCs were retrospectively enrolled. Ninety-six patients and 130 lesions were selected by propensity score matching. Median LRFS was longer for peripheral tumors than central tumors (not reached vs. 3.3 months, p < 0.001). Median PFS was 17.1 months (8.3–24.9) in the peripheral arm and up-to-7 in, 7.0 months (3.3–12.7) in the peripheral arm and up-to-7 out, 8.4 months (4.0–12.6) in the central arm and up-to-7 in, and 3.0 months (1.2–4.9) in the central arm and up-to-7 out-groups. The peripheral arm and up-to-7 in-groups had significantly longer PFS than the other three groups (p = 0.013, p = 0.015, p < 0.001, respectively). Multivariate analysis confirmed that the central zone and central arm were associated with high adjusted hazard ratios for tumor recurrence or death (2.87, p < 0.001; 2.89, p < 0.001, respectively). Conclusion: Treatment-naive HCCs in the peripheral zone had a longer LRFS and PFS following TACE compared to those in the central zone.

Gastroenterology & Hepatology↗

MetaboDirect: an analytical pipeline for the processing of FT-ICR MS-based metabolomic data

Background: Microbiomes are now recognized as the main drivers of ecosystem function ranging from the oceans and soils to humans and bioreactors. However, a grand challenge in microbiome science is to characterize and quantify the chemical currencies of organic matter (i.e., metabolites) that microbes respond to and alter. Critical to this has been the development of Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS), which has drastically increased molecular characterization of complex organic matter samples, but challenges users with hundreds of millions of data points where readily available, user-friendly, and customizable software tools are lacking. Results: Here, we build on years of analytical experience with diverse sample types to develop MetaboDirect, an open-source, command-line-based pipeline for the analysis (e.g., chemodiversity analysis, multivariate statistics), visualization (e.g., Van Krevelen diagrams, elemental and molecular class composition plots), and presentation of direct injection high-resolution FT-ICR MS data sets after molecular formula assignment has been performed. When compared to other available FT-ICR MS software, MetaboDirect is superior in that it requires a single line of code to launch a fully automated framework for the generation and visualization of a wide range of plots, with minimal coding experience required. Among the tools evaluated, MetaboDirect is also uniquely able to automatically generate biochemical transformation networks (ab initio) based on mass differences (mass difference network-based approach) that provide an experimental assessment of metabolite connections within a given sample or a complex metabolic system, thereby providing important information about the nature of the samples and the set of microbial reactions or pathways that gave rise to them. Finally, for more experienced users, MetaboDirect allows users to customize plots, outputs, and analyses. Conclusion: Application of MetaboDirect to FT-ICR MS-based metabolomic data sets from a marine phage-bacterial infection experiment and a Sphagnum leachate microbiome incubation experiment showcase the exploration capabilities of the pipeline that will enable the research community to evaluate and interpret their data in greater depth and in less time. It will further advance our knowledge of how microbial communities influence and are influenced by the chemical makeup of the surrounding system. The source code and User’s guide of MetaboDirect are freely available through (https://github.com/Coayala/MetaboDirect) and (https://metabodirect.readthedocs.io/en/latest/), respectively.

54 ENVIRONMENTAL SCIENCES↗

National variability in Americans’ COVID-19 protective behaviors: Implications for vaccine roll-out

Protective behaviors such as mask wearing and physical distancing are critical to slow the spread of COVID-19, even in the context of vaccine scale-up. Understanding the variation in self-reported COVID-19 protective behaviors is critical to developing public health messaging. The purpose of the study is to provide nationally representative estimates of five self-reported COVID-19 protective behaviors and correlates of such behaviors. In this cross-sectional survey study of US adults, surveys were administered via internet and telephone. Adults were surveyed from April 30-May 4, 2020, a time of peaking COVID-19 incidence within the US. Participants were recruited from the probability-based AmeriSpeak® national panel. Brief surveys were completed by 994 adults, with 73.0% of respondents reported mask wearing, 82.7% reported physical distancing, 75.1% reported crowd avoidance, 89.8% reported increased hand-washing, and 7.7% reported having prior COVID-19 testing. Multivariate analysis (p critical value .05) indicates that women were more likely to report protective behaviors than men, as were those over age 60. Respondents who self-identified as having low incomes, histories of criminal justice involvement, and Republican Party affiliation, were less likely to report four protective behaviors, though Republicans and individuals with criminal justice histories were more likely to report having received COVID-19 testing. The majority of Americans engaged in COVID-19 protective behaviors, with low-income Americans, those with histories of criminal justice involvement, and self-identified Republicans less likely to engage in these preventive behaviors. Culturally competent public health messaging and interventions might focus on these latter groups to prevent future infections. These findings will remain highly relevant even with vaccines widely available, given the complementarities between vaccines and protective behaviors, as well as the many challenges in delivering vaccines.

60 APPLIED LIFE SCIENCES↗

Association of Serum Bile Acid and Unsaturated Fatty Acid Profiles with the Risk of Diabetic Retinopathy in Type 2 Diabetic Patients

Aim: We aimed to identify the ability of serum bile acids (BAs) and unsaturated fatty acids (UFAs) profiles to predict the development of diabetic retinopathy (DR) in type 2 diabetes mellitus (T2DM) patients. Methods: We first used univariate and multivariate analysis to compare 15 serum BA and 11 UFA levels in healthy control (HC) group (n = 82), T2DM patients with DR (n = 58) and T2DM patients without DR (n = 60). Forty T2DM patients were considered for validation. Then, the receiver operating characteristic curve (ROC) and decision curve analysis were used to assess the diagnostic value and clinical benefit of serum biomarkers alone, clinical variables alone or in combination, and the area under the curve (AUC), integrated discrimination improvement (IDI), and net reclassification improvement (NRI) were used to further assess whether the addition of biomarkers significantly improved the predictive ability of the model. Results: Orthogonal partial least squares-discriminant analysis (OPLS-DA) of serum BAs and UFAs separated the three cohorts including HC, T2DM patients with or without DR. The difference in serum BA and UFA profiles of T2DM patients with or without DR was mainly manifested in the three metabolites of taurolithocholic acid (TLCA), tauroursodeoxycholic acid (TUDCA) and arachidonic acid (AA). Together, they had an AUC of 0.785 (0.918 for validation cohort) for predicting DR in T2DM patients. After adjusting for numerous confounding factors, TLCA, TUDCA, and AA were independent predictors that differentiated T2DM with or without DR. The results of AUC, IDI, and NRI demonstrated that adding these three biomarkers to a model with clinical variables statistically increased their predictive value and were replicated in our independent validation cohort. Conclusion: These findings highlight the association of three metabolites, TLCA, TUDCA and AA, with DR and may indicate their potential value in the pathogenesis of DR.

60 APPLIED LIFE SCIENCES↗

Within and Among Fish Species Differences in Simulated Turbine Blade Strike Mortality: Limits on the Use of Surrogacy for Untested Species

Use of surrogacy remains a useful method for prioritizing research on representatives of at-risk groups of fishes, yet quantifiable evidence in support of its use is generally not available. Blade strike impact represents one of the most traumatic stressors experienced by fish during non-volitional movements through hydropower turbines. Here, we use data generated from laboratory trials on blade strike impact experiments to directly test use of surrogacy for salmonid and clupeid fishes. Results of logistic regression indicated that a -taxonomic (genus) variable was not a significant predictor of mortality among large rainbow trout and brook trout. Similar results were found for young-of-the-year shad species, but genus-level taxonomy was a significant predictor of mortality while species was not. Multivariate analysis of morphometric data showed that shad clustered together based on similarities in fish shape which was also closely associated with genus. Logistic regression including size as a major covariate suggested total fish length was not a significant predictor of mortality, yet dose–response data suggest differential susceptibility to lower strike velocities. We suggest that use of surrogacy among species is justifiable but should be avoided within a species since the effects of size remain unclear.

59 BASIC BIOLOGICAL SCIENCES↗

An evaluation of prognostic factors, oncologic outcomes, and management for primary and recurrent squamous cell carcinoma of the vulva

Objective: To evaluate prognostic factors, outcomes, and management patterns of patients treated for squamous cell carcinoma of the vulva. Methods: One hundred sixty-four women were retrospectively identified with primary squamous cell carcinoma of the vulva treated at our institution between 1/1996–12/2018. Descriptive statistics were performed on patient, tumor, and treatment characteristics. The χ 2 tests and t-tests were used to compare categorical variables and continuous variables, respectively. Recurrence free survival (RFS), overall survival (OS), and disease-specific survival (DSS) were analyzed with Kaplan-Meier estimates, the log-rank test, and Cox proportional hazards. Results: Median follow-up was 52.5 months. Five-year RFS was 67.9%, 60.0%, 42.1%, and 20.0% for stage I–IV, respectively. Five-year DSS was 86.2%, 81.6%, 65.0%, and 42.9% for stage I–IV, respectively. On multivariate analysis, positive margins predicted overall RFS (hazard ratio [HR]=3.55; 95% confidence interval [CI]=1.18–10.73; p=0.025), while presence of lichen sclerosus on pathology (HR=2.78; 95% CI=1.30–5.91; p=0.008) predicted local RFS. OS was predicted by nodal involvement (HR=2.51; 95% CI=1.02–6.13; p=0.043) and positive margins (HR=5.19; 95% CI=2.03–13.26; p=0.001). Adjuvant radiotherapy significantly improved RFS (p=0.016) and DSS (p=0.012) in node-positive patients. Median survival after treatment of local, groin, and pelvic/distant recurrence was 52, 8, and 5 months, respectively. Conclusion: For primary treatment, more conservative surgical approaches can be considered with escalation of treatment in patients with concurrent precursor lesions, positive margins, and/or nodal involvement. Further studies are warranted to improve risk stratification in order to optimize treatment paradigms for vulvar cancer patients.

60 APPLIED LIFE SCIENCES↗

MSW Variability Mapping and Conversion to Biofuel

MSW (Municipal Solid Waste) is a form of biomass which consists of categorized components of waste/trash. The general categories are paper, yard trash, construction & debris, appliances, tires, glass, metals, aluminum & steel cans, plastics, organics, inorganics, and HHW (Household Hazardous Waste). This project focuses on the factors within a region or population that contribute to variability in the composition of MSW and in turn MSW’s convertibility to biofuel. A list of contributors was determined (Social Vulnerability Index, Access to Public Transportation, Racial Distribution, GDP, Personal Income) and then JMP was used to perform a Multivariate analysis to determine correlations and a Partial Least-Squares regression to determine Variable Importance Plots for each MSW category. In addition to data analysis, the convertibility of MSW to biofuel was studied via microwave pyrolysis system in order to separate and characterize the various gaseous and bio-oil products.

09 BIOMASS FUELS↗

Carbon fiber classification using raman spectroscopy

Carbon fiber characterization processes are described that include multi-condition Raman spectroscopy-based examination combined with multivariate data analyses. Methods are a nondestructive material characterization approach that can provide predictions as to carbon fiber bulk physical properties, as well as identification of unknown carbon fiber materials for quality control purposes. The framework of the multivariate analysis methods includes a principal component-based identification protocol including comparison of Raman spectral data from an unknown carbon fiber with a data library of multiple principal component spaces.

Houk, Amanda L.↗

Using Boosted Decision Trees to Select High Quality Measurements in the Mu2e Experiment at Fermilab

This thesis presents the implementation and evaluation of a Boosted Decision Tree (BDT) model to improve the selection of high-quality track measurements in the Mu2e experiment at Fermilab. The Mu2e experiment is a high-energy physics experiments seeking to observe a rare theoretical physics process known as Charged Lepton Flavor Violation. A significant challenge faced by the Mu2e experiment are so-called background events, which are events whose data mimics that of the rare physics process the experiment seeks to observe. Without a mechanism to reduce background, it would be impossible to know whether Charged Lepton Flavor Violation occurred or not. To this end, high-quality track measurements must be distinguished from low-quality track measurements. A track can be conceived of as the reconstructed path of a particle that traveled through the Mu2e detector. In addition to other data, data about such tracks is stored using a C++-based framework, specific to the domain of high-energy physics, known as ROOT. A boosted decision tree model was trained using ROOT’s Toolkit For Multivariate Analysis by leveraging variables ancillary to track quality. In evaluation, the BDT achieves a ROC-AUC of 0.927 in discriminating good-quality tracks from poor-quality tracks. Such a score is indicative of both strong discrimination and strong generalization. Subsequently, it is shown that applying a BDT-based quality cut to the distribution of particle momenta significantly enhances the signal-to-background distinction for signal electrons, paving the way for improved sensitivity to Charged Lepton Flavor Violation.

Mullany, Brendan T. [Drew U.] (ORCID:0009000818888↗

Spatio-temporal multivariate cluster evolution analysis for detecting and tracking climate impacts

Recent years have seen a growing concern about climate change and its impacts. While Earth System Models (ESMs) can be invaluable tools for studying the impacts of climate change, the complex coupling processes encoded in ESMs and the large amounts of data produced by these models, together with the high internal variability of the Earth system, can obscure important source-to-impact relationships. Here, this paper presents a novel and efficient unsupervised data-driven approach for detecting statistically-significant impacts and tracing spatio-temporal source-impact pathways in the climate through a unique combination of ideas from anomaly detection, clustering and Natural Language Processing (NLP). Using as an exemplar the 1991 eruption of Mount Pinatubo in the Philippines, we demonstrate that the proposed approach is capable of detecting known post-eruption impacts/events. We additionally describe a methodology for extracting meaningful sequences of post-eruption impacts/events by using NLP to efficiently mine frequent multivariate cluster evolutions, which can be used to confirm or discover the chain of physical processes between a climate source and its impact(s).

Anomaly detection↗

A State-Space Model of an Inverter-Based Microgrid for Multivariable Feedback Control Analysis and Design

In this work, a synchronous model for grid-connected and islanded microgrids is presented. The grid-connected model is based on the premise that the reference frame is synchronized with the AC bus. The quadrature component of the AC bus voltage can be cancelled, which allows to express output power as a linear equation for nominal values in the AC bus amplitude voltage. The model for the islanded microgrid is developed by integrating all the inverter dynamics using a state-space model for the load currents. This model is presented in a comprehensive way such that it could be scalable to any number of inverter-based generators using inductor–capacitor–inductor (LCL) output filters. The use of these models allows designers to assess microgrid stability and robustness using modern control methods such as eigenvalue analysis and singular value diagrams. Both models were tested and validated in an experimental setup to demonstrate their accuracy in describing microgrid dynamics. In addition, three scenarios are presented: non-controlled model, Linear-Quadratic Integrator (LQI) power control, and Power-Voltage (PQ/Vdq) droop–boost controller. Experimental results demonstrate the effectiveness of the control strategies and the accuracy of the models to describe microgrid dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Inorganic characterization of switchgrass biomass using laser-induced breakdown spectroscopy

The inorganic characterization of 74 samples of switchgrass using laser-induced breakdown spectroscopy (LIBS) was undertaken. Determination of ash and inorganic elements content in biomass materials is vital for feedstock screening for bioconversion processes. Hierarchical models using principal component analysis (PCA) and partial least square analysis (PLS) were used to determine the presence of specific elemental micronutrients that are important in determining plant health for robust biomass production. LIBS uses a 532 nm laser with 45 mJ of laser power to excite the samples of switchgrass plant material and the emission of all the elements present in the plant samples were recorded in single spectra with a wide wavelength range of 200–800 nm. The results were compared to the laboratory standard technique, e.g., ICP-OES technique, to determine the true values for major micronutrients such as, silicon (Si), potassium (K), calcium (Ca), magnesium (Mg), phosphorus (P), and sulfur (S). Overall, our objectives were: 1) To determine the spectral features of switchgrass containing different amounts of these elements and 2) To examine the viability of this technique for determining the quality of the feedstock in terms of its inorganic composition. Cross-validation results showed that the broad-based model developed is promising for inorganics prediction in switchgrass. The LIBS validation prediction for the micronutrient elements mentioned here have been obtained. The regression coefficients for Si, were obtained to be 0.995, 0.994 for calibration and validation respectively, in case of Ca the regression coefficients were, 0.994 and 0.992 for calibration and validation. Similarly, in the case of Mg and K these were calculated to be 0.992 and 0.985, and 0.994 and 0.993 respectively. The regression coefficients are not as good as those for the elements mentioned, in case of the two elements S and P. They are 0.957, and 0.878, and 0.952 and 0.894 respectively for calibration, validation for the two elements. This demonstrates that LIBS-based techniques are inherently well suited for diverse environmental applications. Furthermore, LIBS along with PLS model can show capability in determining the viability of switchgrass as a biomass in the production of biofuels and survivability of switchgrass in processes associated with climate change. LIBS can help determining which switchgrass would be appropriate for a specific conversion process that favors low ash content overall or low value of specific inorganics.

59 BASIC BIOLOGICAL SCIENCES↗

Transfer Learning-Based Independent Component Analysis

Understanding the underlying component structure is crucial for multivariate signal analysis. Among all the techniques that try to learn the latent structure, independent component analysis (ICA) is one of the most important and popular methods, which aims to extract independent components from multivariate signals and enables further analysis. For example, in electroencephalogram (EEG) analysis, artifacts filtering and disease detection are conducted based on the independent components of the signals. One critical challenge in existing ICA approaches is that the component extraction accuracy may degrade when the available data of a unit are limited. To address this issue, this paper proposes a transfer learning-based ICA method by innovatively transferring component distribution from a source domain, so that accurate component extraction results can be achieved even when only limited data are available in the target domain. To the best of our knowledge, this is the first work that leverages transfer learning to improve ICA accuracy with limited available data. In particular, we first extract all the independent components from the source domain by maximizing the log-likelihood function with a Newton-like method on a smooth manifold. Then for the target domain, the component with the largest negentropy is extracted in each round. To effectively leverage the knowledge from the source domain and to prevent the negative transfer, we try to find a component in the source domain that matches the component we are extracting. The probability density function of the matched component will then be used to improve the component extraction accuracy if such matched component can be found; otherwise, no knowledge will be transferred. Finally, numerical simulations and a case study with electrocardiogram (ECG) data are conducted, showing the effectiveness of the proposed method in transferring knowledge and reducing negative transfer.

42 ENGINEERING↗

A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study

High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Closure to “Coincidence Risk Analysis of Floods Using Multivariate Copulas: Case Study of Jinsha River and Min River, China” by Yang Peng, Yulong Shi, Hongxiang Yan, Kai Chen, and Jipeng Zhang

We would like to thank the discusser for his comments on the coincidence flood risk modeling in Jinsha River and Min River, China. The discusser’s comment brings focus to the discussion of extreme-value selection methodology and historical flood information. In flood risk analysis, we acknowledge that there is always an alternative modeling approach that could be used to generate desired prediction results. In this comment, the discusser made two valuable points regarding the original paper, and the following two paragraphs provide our responses.

42 ENGINEERING↗