Dataset Information

Hybrid algorithms for generating optimal designs for discriminating multiple nonlinear models under various error distributional assumptions.

ABSTRACT: Finding a model-based optimal design that can optimally discriminate among a class of plausible models is a difficult task because the design criterion is non-differentiable and requires 2 or more layers of nested optimization. We propose hybrid algorithms based on particle swarm optimization (PSO) to solve such optimization problems, including cases when the optimal design is singular, the mean response of some models are not fully specified and problems that involve 4 layers of nested optimization. Using several classical examples, we show that the proposed PSO-based algorithms are not models or criteria specific, and with a few repeated runs, can produce either an optimal design or a highly efficient design. They are also generally faster than the current algorithms, which are generally slow and work for only specific models or discriminating criteria. As an application, we apply our techniques to find optimal discriminating designs for a dose-response study in toxicology with 5 possible models and compare their performances with traditional and a recently proposed algorithm. In the supplementary material, we provide a R package to generate different types of discriminating designs and evaluate efficiencies of competing designs so that the user can implement an informed design.

SUBMITTER: Chen RB

PROVIDER: S-EPMC7535070 | biostudies-literature | 2020

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Hybrid algorithms for generating optimal designs for discriminating multiple nonlinear models under various error distributional assumptions.

Chen Ray-Bing RB Chen Ping-Yang PY Hsu Cheng-Lin CL Wong Weng Kee WK

PloS one 20201005 10

Finding a model-based optimal design that can optimally discriminate among a class of plausible models is a difficult task because the design criterion is non-differentiable and requires 2 or more layers of nested optimization. We propose hybrid algorithms based on particle swarm optimization (PSO) to solve such optimization problems, including cases when the optimal design is singular, the mean response of some models are not fully specified and problems that involve 4 layers of nested optimiza ...[more]

PMID: 33017415

Similar Datasets

Project description:ObjectivesIntra-uterine growth retardation is often of unknown origin, and is of great interest as a "Fetal Origin of Adult Disease" has been now well recognized. We built a benchmark based upon a previously analysed data set related to Intrauterine Growth Retardation with 46 subjects described by 14 variables, related with the insulin-like growth factor system and pro-inflammatory cytokines, namely interleukin-6 and tumor necrosis factor-?.Design and methodsWe used new algorithms for optimal information sorting based on the combination of two neural network algorithms: Auto-contractive Map and Activation and Competition System. Auto-Contractive Map spatializes the relationships among variables or records by constructing a suitable embedding space where 'closeness' among variables or records reflects accurately their associations. The Activation and Competition System algorithm instead works as a dynamic non linear associative memory on the weight matrices of other algorithms, and is able to produce a prototypical variable profile of a given target.ResultsClassical statistical analysis, proved to be unable to distinguish intrauterine growth retardation from appropriate-for-gestational age (AGA) subjects due to the high non-linearity of underlying functions. Auto-contractive map succeeded in clustering and differentiating completely the conditions under study, while Activation and Competition System allowed to develop the profile of variables which discriminated the two conditions under study better than any other previous form of attempt. In particular, Activation and Competition System showed that ppropriateness for gestational age was explained by IGF-2 relative gene expression, and by IGFBP-2 and TNF-? placental contents. IUGR instead was explained by IGF-I, IGFBP-1, IGFBP-2 and IL-6 gene expression in placenta.ConclusionThis further analysis provided further insight into the placental key-players of fetal growth within the insulin-like growth factor and cytokine systems. Our previous published analysis could identify only which variables were predictive of fetal growth in general, and identified only some relationships.

Project description:Gene-set analysis (GSA) is a standard procedure for exploring potential biological functions of a group of genes. The development of its methodology has been an active research topic in recent decades. Many GSA methods, when newly proposed, rely on simulation studies to evaluate their performance with an implicit assumption that the multivariate expression values are normally distributed. This assumption is commonly adopted in GSAs, particularly those in the group of functional class scoring (FCS) methods. The validity of the normality assumption, however, has been disputed in several studies, yet no systematic analysis has been carried out to assess the effect of this distributional assumption. Our goal in this study is not to propose a new GSA method but to first examine if the multi-dimensional gene expression data in gene sets follow a multivariate normal (MVN) distribution. Six statistical methods in three categories of MVN tests were considered and applied to a total of 24 RNA data sets. These RNA values were collected from cancer patients as well as normal subjects, and the values were derived from microarray experiments, RNA sequencing, and single-cell RNA sequencing. Our first finding suggests that the MVN assumption is not always satisfied. This assumption does not hold true in many applications tested here. In the second part of this research, we evaluated the influence of non-normality on the statistical power of current FCS methods, both parametric and nonparametric ones. Specifically, the scenario of mixture distributions representing more than one population for the RNA values was considered. This second investigation demonstrates that the non-normality distribution of the RNA values causes a loss in the statistical power of these GSA tests, especially when subtypes exist. Among the FCS GSA tools examined here and among the scenarios studied in this research, the N-statistics outperform the others. Based on the results from these two investigations, we conclude that the assumption of MVN should be used with caution when evaluating new GSA tools, since this assumption cannot be guaranteed and violation may lead to spurious results, loss of power, and incorrect comparison between methods. If a newly proposed GSA tool is to be evaluated, we recommend the incorporation of a wide range of multivariate non-normal distributions or sampling from large databases if available.

Project description:IntroductionThere is a growing need in clinical research domains for direct comparability between amyloid-beta (Aβ) Positron Emission Tomography (PET) measures obtained via different radiotracers and processing methodologies. Previous efforts to provide a common measurement scale fail to account for non-linearities between measurement scales that can arise from these differences. We introduce a new application of distribution mapping, based on well established statistical orthodoxy, that we call Nonlinear Distribution Mapping (NoDiM). NoDiM uses cumulative distribution functions to derive mappings between Aβ-PET measurements from different tracers and processing streams that align data based on their location in their respective distributions.MethodsUtilizing large datasets of Florbetapir (FBP) from the Alzheimer's Disease Neuroimaging Initiative (n = 349 female (%) = 53) and Pittsburgh Compound B (PiB) from the Harvard Aging Brain Study (n = 305 female (%) = 59.3) and the Australian Imaging, Biomarker & Lifestyle Flagship Study of Ageing (n = 184 female (%) = 53.3), we fit explicit mathematical models of a mixture of two normal distributions, with parameter estimates from Gaussian Mixture Models, to each tracer's empirical data. We demonstrate the accuracy of these fits, and then show the ability of NoDiM to transform FBP measurements into PiB-like units.ResultsA mixture of two normal distributions fit both the FBP and PiB empirical data and provides a strong basis for derivation of a transfer function. Transforming Aβ-PET data with NoDiM results in FBP and PiB distributions that are closely aligned throughout their entire range, while a linear transformation does not. Additionally the NoDiM transform better matches true positive and false positive profiles across tracers.DiscussionThe NoDiM transformation provides a useful alternative to the linear mapping advocated in the Centiloid project, and provides improved correspondence between measurements from different tracers across the range of observed values. This improved alignment enables disparate measures to be merged on to continuous scale, and better enables the use of uniform thresholds across tracers.

Dataset Information

Hybrid algorithms for generating optimal designs for discriminating multiple nonlinear models under various error distributional assumptions.

Publications

Hybrid algorithms for generating optimal designs for discriminating multiple nonlinear models under various error distributional assumptions.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets