Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.

Alternative ways to estimate change points in multinomial sequences. An application to an authorship attribution problem

The statistical analysis of literary style is the part of stylometry that compares measurable characteristicsin a text that are rarely controlled by the author, with those in other texts. When thegoal is to settle authorship questions, these characteristics should relate to the author’s style andnot to the genre, epoch or editor, and they should be such that their variation between authors islarger than the variation within comparable texts from the same author.For an overview of the literature on stylometry and some of the techniques involved, see for exampleMosteller and Wallace (1964, 82), Herdan (1964), Morton (1978), Holmes (1985), Oakes (1998) orLebart, Salem and Berry (1998).Tirant lo Blanc, a chivalry book, is the main work in catalan literature and it was hailed to be“the best book of its kind in the world” by Cervantes in Don Quixote. Considered by writterslike Vargas Llosa or Damaso Alonso to be the first modern novel in Europe, it has been translatedseveral times into Spanish, Italian and French, with modern English translations by Rosenthal(1996) and La Fontaine (1993). The main body of this book was written between 1460 and 1465,but it was not printed until 1490.There is an intense and long lasting debate around its authorship sprouting from its first edition,where its introduction states that the whole book is the work of Martorell (1413?-1468), while atthe end it is stated that the last one fourth of the book is by Galba (?-1490), after the death ofMartorell. Some of the authors that support the theory of single authorship are Riquer (1990),Chiner (1993) and Badia (1993), while some of those supporting the double authorship are Riquer(1947), Coromines (1956) and Ferrando (1995). For an overview of this debate, see Riquer (1990).Neither of the two candidate authors left any text comparable to the one under study, and thereforediscriminant analysis can not be used to help classify chapters by author. By using sample textsencompassing about ten percent of the book, and looking at word length and at the use of 44conjunctions, prepositions and articles, Ginebra and Cabos (1998) detect heterogeneities that mightindicate the existence of two authors. By analyzing the diversity of the vocabulary, Riba andGinebra (2000) estimates that stylistic boundary to be near chapter 383.Following the lead of the extensive literature, this paper looks into word length, the use of the mostfrequent words and into the use of vowels in each chapter of the book. Given that the featuresselected are categorical, that leads to three contingency tables of ordered rows and therefore tothree sequences of multinomial observations.Section 2 explores these sequences graphically, observing a clear shift in their distribution. Section 3describes the problem of the estimation of a suden change-point in those sequences, in the followingsections we propose various ways to estimate change-points in multinomial sequences; the methodin section 4 involves fitting models for polytomous data, the one in Section 5 fits gamma modelsonto the sequence of Chi-square distances between each row profiles and the average profile, theone in Section 6 fits models onto the sequence of values taken by the first component of thecorrespondence analysis as well as onto sequences of other summary measures like the averageword length. In Section 7 we fit models onto the marginal binomial sequences to identify thefeatures that distinguish the chapters before and after that boundary. Most methods rely heavilyon the use of generalized linear models

Leontief versus Ghosh: two faces of the same coin

In this paper we explore the sectoral and aggregate implications of some endogeneization rules (i.e. on value-added and final demand) which have been common in the Leontief model and have been recently proposed in the Ghosh model. We detect that these rules may give rise in both models to some allegedly pathological behavior in the sense that sectoral or aggregate output, very often, may not follow the logical and economically expected direct relationship with some underlying endogenous variables—namely, output and value-added in the Ghosh model and output and consumption in the Leontief model. Because of the common mathematical structure, whatever is or seems to be pathological in the Ghosh model also has a symmetric counterpart in the Leontief model. These would not be good news for the inner consistency of these linear models. To avoid such possible inconsistencies, we propose new and simple endogeneization rules that have a sound economic interpretation.

Modeling of hybrid systems by piecewise decomposition

Piecewise linear models systems arise as mathematical models of systems in many practical applications, often from linearization for nonlinear systems. There are two main approaches of dealing with these systems according to their continuous or discrete-time aspects. We propose an approach which is based on the state transformation, more particularly the partition of the phase portrait in different regions where each subregion is modeled as a two-dimensional linear time invariant system. Then the Takagi-Sugeno model, which is a combination of local model is calculated. The simulation results show that the Alpha partition is well-suited for dealing with such a system

Internal test sets studies in a group of antimalarials

Topological indices have been applied to build QSAR models for a set of 20 antimalarial cyclic peroxy cetals. In order to evaluate the reliability of the proposed linear models leave-n-out and Internal Test Sets (ITS) approaches have been considered. The proposed procedure resulted in a robust and consensued prediction equation and here it is shown why it is superior to the employed standard cross-validation algorithms involving multilinear regression models

Emotion and reason in everyday risk perception

Although research has documented the importance of emotion in risk perception, little is knownabout its prevalence in everyday life. Using the Experience Sampling Method, 94 part-timestudents were prompted at random via cellular telephones to report on mood state and threeemotions and to assess risk on thirty occasions during their working hours. The emotions valence, arousal, and dominance were measured using self-assessment manikins (Bradley &Lang, 1994). Hierarchical linear models (HLM) revealed that mood state and emotions explainedsignificant variance in risk perception. In addition, valence and arousal accounted for varianceover and above reason (measured by severity and possibility of risks). Six risks were reassessedin a post-experimental session and found to be lower than their real-time counterparts.The study demonstrates the feasibility and value of collecting representative samples of data withsimple technology. Evidence for the statistical consistency of the HLM estimates is provided inan Appendix.

Forecasting with missing data: Application to a real case

This paper presents a comparative analysis of linear and mixed modelsfor short term forecasting of a real data series with a high percentage of missing data. Data are the series of significant wave heights registered at regular periods of three hours by a buoy placed in the Bay of Biscay.The series is interpolated with a linear predictor which minimizes theforecast mean square error. The linear models are seasonal ARIMA models and themixed models have a linear component and a non linear seasonal component.The non linear component is estimated by a non parametric regression of dataversus time. Short term forecasts, no more than two days ahead, are of interestbecause they can be used by the port authorities to notice the fleet.Several models are fitted and compared by their forecasting behavior.

On ignoring scientific evidence: The bumpy road to enlightenment

It is well accepted that people resist evidence that contradicts their beliefs.Moreover, despite their training, many scientists reject results that are inconsistent withtheir theories. This phenomenon is discussed in relation to the field of judgment anddecision making by describing four case studies. These concern findings that clinical judgment is less predictive than actuarial models; simple methods have proven superiorto more theoretically correct methods in times series forecasting; equal weighting ofvariables is often more accurate than using differential weights; and decisions cansometimes be improved by discarding relevant information. All findings relate to theapparently difficult-to-accept idea that simple models can predict complex phenomenabetter than complex ones. It is true that there is a scientific market place for ideas.However, like its economic counterpart, it is subject to inefficiencies (e.g., thinness,asymmetric information, and speculative bubbles). Unfortunately, the market is only correct in the long-run. The road to enlightenment is bumpy.

Individual prediction of automobile bodily injury claims liabilities

[spa] La mayoría de siniestros con daños corporales se liquidan mediante negociación, llegando a juicio menos del 5% de los casos. Una estrategia de negociación bien definida es, por tanto, fundamental para las compañías aseguradoras. En este artículo asumimos que la compensación monetaria concedida en juicio es la máxima cuantía que debería ser ofrecida por el asegurador en el proceso de negociación. Usando una base de datos real, implementamos un modelo log-lineal para estimar la máxima oferta de negociación. Perturbaciones no-esféricas son detectadas. Correlación ocurre cuando más de una siniestro se liquida en la misma sentencia judicial. Heterocedasticidad por grupos se debe a la influencia de la valoración del forense en la indemnización final.

Individual prediction of automobile bodily injury claims liabilities

[spa] La mayoría de siniestros con daños corporales se liquidan mediante negociación, llegando a juicio menos del 5% de los casos. Una estrategia de negociación bien definida es, por tanto, fundamental para las compañías aseguradoras. En este artículo asumimos que la compensación monetaria concedida en juicio es la máxima cuantía que debería ser ofrecida por el asegurador en el proceso de negociación. Usando una base de datos real, implementamos un modelo log-lineal para estimar la máxima oferta de negociación. Perturbaciones no-esféricas son detectadas. Correlación ocurre cuando más de una siniestro se liquida en la misma sentencia judicial. Heterocedasticidad por grupos se debe a la influencia de la valoración del forense en la indemnización final.

Los modelos causales asociados a estructuras de covarianza: una via para la unificación de los dos métodos de la Psicologia cientifica

A través de este articulo he pretendido exponer una introducción a los Modelos Causales desde las perspectivas histórica, metodológica y estadístico-matemtitica. Asimismo, he planteado las posibilidades que tienen estos modelos, asociados a estructuras de covariancia, como via que puede contribuir a superar la dicotomia metodo experimental/mktodo correlacional. Para la exposición de 10s procesos de especificación, identificación, estimación, prueba de ajuste e interpretación de parámetros de un modelo causal he tomado como referencia el modelo LISREL. Sin embargo, en la sistematización de estos procesos no he intentado presentar solamente la simple descripción del modelo, insistiendo, para ello, en la necesidad de integrar este modelo en el contexto del método cientifico.

Brain functional abnormality in schizo-affective disorder: an fMRI study.

Background.Schizo-affective disorder has not been studied to any significant extent using functional imaging. The aim of this study was to examine patterns of brain activation and deactivation in patients meeting strict diagnostic criteria for the disorder. METHOD: Thirty-two patients meeting research diagnostic criteria (RDC) for schizo-affective disorder (16 schizomanic and 16 schizodepressive) and 32 matched healthy controls underwent functional magnetic resonance imaging (fMRI) during performance of the n-back task. Linear models were used to obtain maps of activations and deactivations in the groups. RESULTS: Controls showed activation in a network of frontal and other areas and also deactivation in the medial frontal cortex, the precuneus and the parietal cortex. Schizo-affective patients activated significantly less in prefrontal, parietal and temporal regions than the controls, and also showed failure of deactivation in the medial frontal cortex. When task performance was controlled for, the reduced activation in the dorsolateral prefrontal cortex (DLPFC) and the failure of deactivation of the medial frontal cortex remained significant. CONCLUSIONS: Schizo-affective disorder shows a similar pattern of reduced frontal activation to schizophrenia. The disorder is also characterized by failure of deactivation suggestive of default mode network dysfunction.

Brain functional abnormality in schizo-affective disorder: an fMRI study.

Background.Schizo-affective disorder has not been studied to any significant extent using functional imaging. The aim of this study was to examine patterns of brain activation and deactivation in patients meeting strict diagnostic criteria for the disorder. METHOD: Thirty-two patients meeting research diagnostic criteria (RDC) for schizo-affective disorder (16 schizomanic and 16 schizodepressive) and 32 matched healthy controls underwent functional magnetic resonance imaging (fMRI) during performance of the n-back task. Linear models were used to obtain maps of activations and deactivations in the groups. RESULTS: Controls showed activation in a network of frontal and other areas and also deactivation in the medial frontal cortex, the precuneus and the parietal cortex. Schizo-affective patients activated significantly less in prefrontal, parietal and temporal regions than the controls, and also showed failure of deactivation in the medial frontal cortex. When task performance was controlled for, the reduced activation in the dorsolateral prefrontal cortex (DLPFC) and the failure of deactivation of the medial frontal cortex remained significant. CONCLUSIONS: Schizo-affective disorder shows a similar pattern of reduced frontal activation to schizophrenia. The disorder is also characterized by failure of deactivation suggestive of default mode network dysfunction.

Brain functional abnormality in schizo-affective disorder: an fMRI study.

Background.Schizo-affective disorder has not been studied to any significant extent using functional imaging. The aim of this study was to examine patterns of brain activation and deactivation in patients meeting strict diagnostic criteria for the disorder. METHOD: Thirty-two patients meeting research diagnostic criteria (RDC) for schizo-affective disorder (16 schizomanic and 16 schizodepressive) and 32 matched healthy controls underwent functional magnetic resonance imaging (fMRI) during performance of the n-back task. Linear models were used to obtain maps of activations and deactivations in the groups. RESULTS: Controls showed activation in a network of frontal and other areas and also deactivation in the medial frontal cortex, the precuneus and the parietal cortex. Schizo-affective patients activated significantly less in prefrontal, parietal and temporal regions than the controls, and also showed failure of deactivation in the medial frontal cortex. When task performance was controlled for, the reduced activation in the dorsolateral prefrontal cortex (DLPFC) and the failure of deactivation of the medial frontal cortex remained significant. CONCLUSIONS: Schizo-affective disorder shows a similar pattern of reduced frontal activation to schizophrenia. The disorder is also characterized by failure of deactivation suggestive of default mode network dysfunction.

Quantitative developmental response to the length of exposure to long photoperiod in wheat and barley

The present study was designed to analyse the effect of the length of exposure to a long photoperiod imposed c. 3 weeks after sowing in spring wheat (cv. UQ189) and barley (cv. Arapiles) to (i) establish whether the response to the number of cycles of exposure is quantitative or qualitative, (ii) determine the existence of a commitment to particular stages well before the stage has been observable, and (iii) study the interrelationships between the effects on final leaf number and phyllochron when the stimulus is provided several days after seedling emergence. Both wheat and barley seemed to respond quantitatively to the number of long-day cycles they were exposed to. However, wheat showed a requirement of approximately 4 long-day cycles to be able to produce a significant response in time to heading. The barley cultivar used in the study was responsive to the minimum length of exposure. The response to extended photoperiod cycles during the stem elongation phase was due to the ‘ memory’ photoperiod effects being related, in the case of wheat, to the fact that the pre-terminal spikelet appearance phase saturated its photoperiod response well before that stage was reached. Therefore, the commitment to the terminal spikelet appearance in wheat may be reached well before this stage could be recognized. As the response in duration to heading exceeded that of the final leaf number, and the stem elongation phase responded to memory effects of photoperiod, the phyllochron of both cereals was responsive to the treatments accelerating the average phyllochron when exposed to longer periods of long days. The response in average phyllochron was due to a switch from bi-linear to linear models of leaf number v. time when the conditions were increasingly inductive, with the phyllochron of the initial (6–8) leaves being similar for all treatments (within each species), and from then on increased.