Showing 1–27 of 27 results for author: Louzoun, Y

IRMA: Iterative Repair for graph MAtching

Authors: Barak Babayov, Yoram Louzoun

Abstract: The alignment of two similar graphs from different domains is a well-studied problem. In many practical usages, there is no reliable information or labels over the vertices or edges, leaving structural similarity as the only information available to match such a graph. In such cases, one often assumes a small amount of already aligned vertices -- called a seed. Current state-of-the-art scalable se… ▽ More The alignment of two similar graphs from different domains is a well-studied problem. In many practical usages, there is no reliable information or labels over the vertices or edges, leaving structural similarity as the only information available to match such a graph. In such cases, one often assumes a small amount of already aligned vertices -- called a seed. Current state-of-the-art scalable seeded alignment algorithms are based on percolation. Namely, aligned vertices are used to align their neighbors and gradually percolate in parallel in both graphs. However, percolation-based graph alignment algorithms are still limited in scale-free degree distributions. We here propose `IRMA' -- Iterative Repair for graph MAtching to show that the accuracy of percolation-based algorithms can be improved in real-world graphs with a limited additional computational cost, and with lower run time when used in a parallel version. IRMA starts by creating a primary alignment using an existing percolation algorithm, then it iteratively repairs the mistakes in the previous alignment steps. We prove that IRMA improves on single-iteration algorithms. We then numerically show that it is significantly better than all state-of-the-art seeded graph alignment algorithms on the graphs that they tested. In scale-free networks, many vertices have a very low degree. Such vertices have a high probability of erroneous alignments. We show that combining iterations with high recall but low precision in the alignment leads in the long run to higher recall and precision for the entire alignment. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Family based HLA imputation and optimization of haplo-identical transplants

Authors: Zuriya Ansbacher-Feldman, Sapir Israeli, Martin Maiers, Loren Gragert, Dianne De Santis, Moshe Israeli, Yoram Louzoun

Abstract: Recently, haplo-identical transplantation with multiple HLA mismatches has become a viable option for system cell transplants. Haplotype sharing detection requires imputation of donor and recipient. We show that even in high-resolution typing when all alleles are known, there is a 15% error rate in haplotype phasing, and even more in low resolution typings. Similarly, in related donors, parents ha… ▽ More Recently, haplo-identical transplantation with multiple HLA mismatches has become a viable option for system cell transplants. Haplotype sharing detection requires imputation of donor and recipient. We show that even in high-resolution typing when all alleles are known, there is a 15% error rate in haplotype phasing, and even more in low resolution typings. Similarly, in related donors, parents haplotypes should be imputed to determine what haplotype each child inherited. We propose GRAMM (GRaph bAsed FaMilly iMputation) to phase alleles in family pedigree HLA typing data, and in mother-cord blood unit pairs. We show that GRAMM has practically no phasing errors when pedigree data are available. We apply GRAMM to simulations with different typing resolutions as well as paired cord-mother typings, and show very high phasing accuracy, and improved alleles imputation accuracy. We use GRAMM to detect recombination events and show that the rate of falsely detected recombination events (False Positive Rate) in simulations is very low. We then apply recombination detection to typed families to estimate the recombination rate in Israeli and Australian population datasets. The estimated recombination rate has an upper bound of 10-20% per family (1-4% per individual). GRAMM is available at: https://gramm.math.biu.ac.il/. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 23 pages (20 + 3 of Supp. Mat.), 9 figures (6 + 3 of Supp. Mat.) Submitted to NAR (Nucleic Acids Research)

BFS based distributed algorithm for parallel local directed sub-graph enumeration

Authors: Itay Levinas, Roy Scherz, Yoram Louzoun

Abstract: Estimating the frequency of sub-graphs is of importance for many tasks, including sub-graph isomorphism, kernel-based anomaly detection, and network structure analysis. While multiple algorithms were proposed for full enumeration or sampling-based estimates, these methods fail in very large graphs. Recent advances in parallelization allow for estimates of total sub-graphs counts in very large grap… ▽ More Estimating the frequency of sub-graphs is of importance for many tasks, including sub-graph isomorphism, kernel-based anomaly detection, and network structure analysis. While multiple algorithms were proposed for full enumeration or sampling-based estimates, these methods fail in very large graphs. Recent advances in parallelization allow for estimates of total sub-graphs counts in very large graphs. The task of counting the frequency of each sub-graph associated with each vertex also received excellent solutions for undirected graphs. However, there is currently no good solution for very large directed graphs. We here propose VDMC (Vertex specific Distributed Motif Counting) -- a fully distributed algorithm to optimally count all the 3 and 4 vertices connected directed graphs (sub-graph motifs) associated with each vertex of a graph. VDMC counts each motif only once and its efficacy is linear in the number of counted motifs. It is fully parallelized to be efficient in GPU-based computation. VDMC is based on three main elements: 1) Ordering the vertices and only counting motifs containing increasing order vertices, 2) sub-ordering motifs based on the average length of the BFS composing the motif, and 3) removing isomorphisms only once for the entire graph. We here compare VDMC to analytical estimates of the expected number of motifs and show its accuracy. VDMC is available as a highly efficient CPU and GPU code with a novel data structure for efficient graph manipulation. We show the efficacy of VDMC and real-world graphs. VDMC allows for the precise analysis of sub-graph frequency around each vertex in large graphs and opens the way for the extension of methods until now limited to graphs of thousands of edges to graphs with millions of edges and above. GIT: https://github.com/louzounlab/graph-measures △ Less

Submitted 6 January, 2022; originally announced January 2022.

MSC Class: 05C30; 05C75; 05C85

Planted Dense Subgraphs in Dense Random Graphs Can Be Recovered using Graph-based Machine Learning

Authors: Itay Levinas, Yoram Louzoun

Abstract: Multiple methods of finding the vertices belonging to a planted dense subgraph in a random dense $G(n, p)$ graph have been proposed, with an emphasis on planted cliques. Such methods can identify the planted subgraph in polynomial time, but are all limited to several subgraph structures. Here, we present PYGON, a graph neural network-based algorithm, which is insensitive to the structure of the pl… ▽ More Multiple methods of finding the vertices belonging to a planted dense subgraph in a random dense $G(n, p)$ graph have been proposed, with an emphasis on planted cliques. Such methods can identify the planted subgraph in polynomial time, but are all limited to several subgraph structures. Here, we present PYGON, a graph neural network-based algorithm, which is insensitive to the structure of the planted subgraph. This is the first algorithm that uses advanced learning tools for recovering dense subgraphs. We show that PYGON can recover cliques of sizes $Θ\left(\sqrt{n}\right)$, where $n$ is the size of the background graph, comparable with the state of the art. We also show that the same algorithm can recover multiple other planted subgraphs of size $Θ\left(\sqrt{n}\right)$, in both directed and undirected graphs. We suggest a conjecture that no polynomial time PAC-learning algorithm can detect planted dense subgraphs with size smaller than $O\left(\sqrt{n}\right)$, even if in principle one could find dense subgraphs of logarithmic size. △ Less

Submitted 26 November, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

ACM Class: G.2.2; I.2.1

Journal ref: Journal of Artificial Intelligence Research 75 (2022) 541-568

Exposing individual differences through network topology

Authors: Yuval Samoilov-Katz, Yoram Louzoun, Lev Muchnik, Adam Zaidel

Abstract: Social animals, including humans, have a broad range of personality traits, which can be used to predict individual behavioral responses and decisions. Current methods to quantify individual personality traits in humans rely on self-report questionnaires, which require time and effort to collect and rely on active cooperation. However, personality differences naturally manifest in social interacti… ▽ More Social animals, including humans, have a broad range of personality traits, which can be used to predict individual behavioral responses and decisions. Current methods to quantify individual personality traits in humans rely on self-report questionnaires, which require time and effort to collect and rely on active cooperation. However, personality differences naturally manifest in social interactions such as online social networks. Here, we demonstrate that the topology of an online social network can be used to characterize the personality traits of its members. We analyzed the directed social graph formed by the users of the LiveJournal (LJ) blogging platform. Individual users personality traits, inferred from their self-reported domains of interest (DOIs), were associated with their network measures. Empirical clustering of DOIs by topological similarity exposed two main self-emergent DOI groups that were in alignment with the personality meta-traits plasticity and stability. Closeness, a global topological measure of network centrality, was significantly higher for bloggers associated with plasticity (vs. stability). A local network motif (a triad of 3 connected bloggers) that correlated with closeness also separated the personality meta-traits. Finally, topology-based classification of DOIs (without analyzing the content of the blogs) attained > 70% accuracy (average AUC of the test-set). These results indicate that personality traits are evident and detectable in network topology. This has serious implications for user privacy. But, if used responsibly, network identification of personality traits could aid in early identification of health-related risks, at the population level. △ Less

Submitted 27 March, 2022; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: Peer-review of this preprint exposed criticisms of the methods and analyses that we believe to be valid. These question the results and conclusions. Therefore, we would like to withdraw the preprint

High fraction of silent recombination in a finite population two-locus neutral birth-death-mutation model

Authors: Akiva Bruno Melka, Yoram Louzoun

Abstract: A precise estimate of allele and haplotype polymorphism is of great interest in theoretical population genetics, but also has practical applications, such as bone marrow registries management. Allele polymorphism is driven mainly by point mutations, while haplotype polymorphism is also affected by recombination. Current estimates treat recombination as mutations in an infinite site model. We here… ▽ More A precise estimate of allele and haplotype polymorphism is of great interest in theoretical population genetics, but also has practical applications, such as bone marrow registries management. Allele polymorphism is driven mainly by point mutations, while haplotype polymorphism is also affected by recombination. Current estimates treat recombination as mutations in an infinite site model. We here show that even in the simple case of two loci in a haploid individual, for a finite population, most recombination events produce existing haplotypes, and as such are silent. Silent recombination considerably reduces the total number of haplotypes expected from the infinite site model for populations that are not much larger than one over the mutation rate. Moreover, in contrast with mutations, the number of haplotypes does not grow linearly with the population size. We hence propose a more accurate estimate of the total number of haplotypes that takes into account silent recombination. We study large-scale Human Leukocyte Antigen (HLA) haplotype frequencies from human populations to show that the current estimated recombination rate in the HLA region is underestimated. △ Less

Submitted 10 July, 2022; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 12 pages, 7 figures

Journal ref: Phys. Rev. E 106, 024409 (2022)

Quadratic GCN for Graph Classification

Authors: Omer Nagar, Shoval Frydman, Ori Hochman, Yoram Louzoun

Abstract: Graph Convolutional Networks (GCNs) have been extensively used to classify vertices in graphs and have been shown to outperform other vertex classification methods. GCNs have been extended to graph classification tasks (GCT). In GCT, graphs with different numbers of edges and vertices belong to different classes, and one attempts to predict the graph class. GCN based GCT have mostly used pooling a… ▽ More Graph Convolutional Networks (GCNs) have been extensively used to classify vertices in graphs and have been shown to outperform other vertex classification methods. GCNs have been extended to graph classification tasks (GCT). In GCT, graphs with different numbers of edges and vertices belong to different classes, and one attempts to predict the graph class. GCN based GCT have mostly used pooling and attention-based models. The accuracy of existing GCT methods is still limited. We here propose a novel solution combining GCN, methods from knowledge graphs, and a new self-regularized activation function to significantly improve the accuracy of the GCN based GCT. We present quadratic GCN (QGCN) - A GCN formalism with a quadratic layer. Such a layer produces an output with fixed dimensions, independent of the graph vertex number. We applied this method to a wide range of graph classification problems, and show that when using a self regularized activation function, QGCN outperforms the state of the art methods for all graph classification tasks tested with or without external input on each graph. The code for QGCN is available at: https://github.com/Unknown-Data/QGCN . △ Less

Submitted 14 April, 2021; originally announced April 2021.

Self-driven criticality in a stochastic epidemic model

Authors: Gil Ariel, Yoram Louzoun

Abstract: We present a generic epidemic model with stochastic parameters, in which the dynamics self-organize to a critical state with suppressed exponential growth. More precisely, the dynamics evolve into a quasi-steady-state, where the effective reproduction rate fluctuates close to the critical value one, as observed for different epidemics. The main assumptions underlying the model are that the rate at… ▽ More We present a generic epidemic model with stochastic parameters, in which the dynamics self-organize to a critical state with suppressed exponential growth. More precisely, the dynamics evolve into a quasi-steady-state, where the effective reproduction rate fluctuates close to the critical value one, as observed for different epidemics. The main assumptions underlying the model are that the rate at which each individual becomes infected changes stochastically in time with a heavy-tailed steady state. The critical regime is characterized by an extremely long duration of the epidemic. Its stability is analyzed both numerically and analytically. △ Less

Submitted 7 June, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Journal ref: Phys. Rev. E 103, 062303 (2021)

Initial growth rates of malware epidemics fail to predict their reach

Authors: Lev Muchnik, Elad Yom-Tov, Nir Levy, Amir Rubin, Yoram Louzoun

Abstract: Empirical studies show that epidemiological models based on an epidemic's initial spread rate often fail to predict the true scale of that epidemic. Most epidemics with a rapid early rise die out before affecting a significant fraction of the population, whereas the early pace of some pandemics is rather modest. Recent models suggest that this could be due to the heterogeneity of the target popula… ▽ More Empirical studies show that epidemiological models based on an epidemic's initial spread rate often fail to predict the true scale of that epidemic. Most epidemics with a rapid early rise die out before affecting a significant fraction of the population, whereas the early pace of some pandemics is rather modest. Recent models suggest that this could be due to the heterogeneity of the target population's susceptibility. We study a computer malware ecosystem exhibiting spread mechanisms resembling those of biological systems while offering details unavailable for human epidemics. Rather than comparing models, we directly estimate reach from a new and vastly more complete data from a parallel domain, that offers superior details and insight as concerns biological outbreaks. We find a highly heterogeneous distribution of computer susceptibilities, with nearly all outbreaks initially over-affecting the tail of the distribution, then collapsing quickly once this tail is depleted. This mechanism restricts the correlation between an epidemic's initial growth rate and its total reach, thus preventing the majority of epidemics, including initially fast-growing outbreaks, from reaching a macroscopic fraction of the population. The few pervasive malwares distinguish themselves early on via the following key trait: they avoid infecting the tail, while preferentially targeting computers unaffected by typical malware. △ Less

Submitted 4 May, 2021; v1 submitted 2 August, 2020; originally announced August 2020.

Explicit Gradient Learning

Authors: Mor Sinay, Elad Sarafian, Yoram Louzoun, Noa Agmon, Sarit Kraus

Abstract: Black-Box Optimization (BBO) methods can find optimal policies for systems that interact with complex environments with no analytical representation. As such, they are of interest in many Artificial Intelligence (AI) domains. Yet classical BBO methods fall short in high-dimensional non-convex problems. They are thus often overlooked in real-world AI tasks. Here we present a BBO method, termed Expl… ▽ More Black-Box Optimization (BBO) methods can find optimal policies for systems that interact with complex environments with no analytical representation. As such, they are of interest in many Artificial Intelligence (AI) domains. Yet classical BBO methods fall short in high-dimensional non-convex problems. They are thus often overlooked in real-world AI tasks. Here we present a BBO method, termed Explicit Gradient Learning (EGL), that is designed to optimize high-dimensional ill-behaved functions. We derive EGL by finding weak-spots in methods that fit the objective function with a parametric Neural Network (NN) model and obtain the gradient signal by calculating the parametric gradient. Instead of fitting the function, EGL trains a NN to estimate the objective gradient directly. We prove the convergence of EGL in convex optimization and its robustness in the optimization of integrable functions. We evaluate EGL and achieve state-of-the-art results in two challenging problems: (1) the COCO test suite against an assortment of standard BBO methods; and (2) in a high-dimensional non-convex image generation task. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Evaluation of the number of undiagnosed infected in an outbreak using source of infection measurements

Authors: Akiva B. Melka, Yoram Louzoun

Abstract: In times of outbreaks, an essential requirement for better monitoring is the evaluation of the number of undiagnosed infected individuals. An accurate estimate of this fraction is crucial for the assessment of the situation and the establishment of protective measures. In most current studies using epidemics models, the total number of infected is either approximated by the number of diagnosed ind… ▽ More In times of outbreaks, an essential requirement for better monitoring is the evaluation of the number of undiagnosed infected individuals. An accurate estimate of this fraction is crucial for the assessment of the situation and the establishment of protective measures. In most current studies using epidemics models, the total number of infected is either approximated by the number of diagnosed individuals or is dependent on the model parameters and assumptions, which are often debated. We here study the relationship between the fraction of diagnosed infected out of all infected, and the fraction of infected with known contaminator out of all diagnosed infected. We show that those two are approximately the same in exponential models and across most models currently used in the study of epidemics, independently of the model parameters. As an application, we compute an estimate of the effective number of infected by the SARS-CoV-2 virus in various countries. △ Less

Submitted 16 March, 2021; v1 submitted 9 June, 2020; originally announced June 2020.

Journal ref: Sci Rep 11, 3601 (2021)

Invasion rate versus diversity in population dynamics with catastrophes

Authors: Akiva Bruno Melka, Navit Dori, Yoram Louzoun

Abstract: A key question in the current diversity crisis is how diversity has been maintained throughout evolution and how to preserve it. Modern coexistence theories suggest that a high invasion rate of rare new types is directly related to diversity. We show that adding almost any mechanism of catastrophes to a stochastic birth, death, and mutation process with limited carrying capacity induces a novel ph… ▽ More A key question in the current diversity crisis is how diversity has been maintained throughout evolution and how to preserve it. Modern coexistence theories suggest that a high invasion rate of rare new types is directly related to diversity. We show that adding almost any mechanism of catastrophes to a stochastic birth, death, and mutation process with limited carrying capacity induces a novel phase transition characterized by a positive invasion rate but a low diversity. In this phase, new types emerge and grow rapidly, but the resulting growth of very large types decreases diversity. This model also resolves two major drawbacks of neutral evolution models: their failure to explain balancing selection without resorting to fitness differences and the unrealistic time required for the creation of the observed large types. We test this model on a classical case of genetic polymorphism: the HLA locus. △ Less

Submitted 18 April, 2020; v1 submitted 12 January, 2020; originally announced January 2020.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 124, 158301 (2020)

Topological based classification using graph convolutional networks

Authors: Roy Abel, Idan Benami, Yoram Louzoun

Abstract: In colored graphs, node classes are often associated with either their neighbors class or with information not incorporated in the graph associated with each node. We here propose that node classes are also associated with topological features of the nodes. We use this association to improve Graph machine learning in general and specifically, Graph Convolutional Networks (GCN). First, we show th… ▽ More In colored graphs, node classes are often associated with either their neighbors class or with information not incorporated in the graph associated with each node. We here propose that node classes are also associated with topological features of the nodes. We use this association to improve Graph machine learning in general and specifically, Graph Convolutional Networks (GCN). First, we show that even in the absence of any external information on nodes, a good accuracy can be obtained on the prediction of the node class using either topological features, or using the neighbors class as an input to a GCN. This accuracy is slightly less than the one that can be obtained using content based GCN. Secondly, we show that explicitly adding the topology as an input to the GCN does not improve the accuracy when combined with external information on nodes. However, adding an additional adjacency matrix with edges between distant nodes with similar topology to the GCN does significantly improve its accuracy, leading to results better than all state of the art methods in multiple datasets. △ Less

Submitted 26 October, 2019; originally announced November 2019.

Comments: arXiv admin note: text overlap with arXiv:1904.07787

Estimating differential entropy using recursive copula splitting

Authors: Gil Ariel, Yoram Louzoun

Abstract: A method for estimating the Shannon differential entropy of multidimensional random variables using independent samples is described. The method is based on decomposing the distribution into a product of the marginal distributions and the joint dependency, also known as the copula. The entropy of marginals is estimated using one-dimensional methods. The entropy of the copula, which always has a co… ▽ More A method for estimating the Shannon differential entropy of multidimensional random variables using independent samples is described. The method is based on decomposing the distribution into a product of the marginal distributions and the joint dependency, also known as the copula. The entropy of marginals is estimated using one-dimensional methods. The entropy of the copula, which always has a compact support, is estimated recursively by splitting the data along statistically dependent dimensions. Numerical examples demonstrate that the method is accurate for distributions with compact and non-compact supports, which is imperative when the support is not known or of mixed type (in different dimensions). At high dimensions (larger than 20), our method is not only more accurate, but also significantly more efficient than existing approaches. △ Less

Submitted 9 January, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Regional based query in graph active learning

Authors: Roy Abel, Yoram Louzoun

Abstract: Graph convolution networks (GCN) have emerged as the leading method to classify node classes in networks, and have reached the highest accuracy in multiple node classification tasks. In the absence of available tagged samples, active learning methods have been developed to obtain the highest accuracy using the minimal number of queries to an oracle. The current best active learning methods use the… ▽ More Graph convolution networks (GCN) have emerged as the leading method to classify node classes in networks, and have reached the highest accuracy in multiple node classification tasks. In the absence of available tagged samples, active learning methods have been developed to obtain the highest accuracy using the minimal number of queries to an oracle. The current best active learning methods use the sample class uncertainty as selection criteria. However, in graph based classification, the class of each node is often related to the class of its neighbors. As such, the uncertainty in the class of a node's neighbor may be a more appropriate selection criterion. We here propose two such criteria, one extending the classical uncertainty measure, and the other extending the page-rank algorithm. We show that the latter is optimal when the fraction of tagged nodes is low, and when this fraction grows to one over the average degree, the regional uncertainty performs better than all existing methods. While we have tested this methods on graphs, such methods can be extended to any classification problem, where a distance metrics can be defined between the input samples. All the code used can be accessed at : https://github.com/louzounlab/graph-al All the datasets used can be accessed at : https://github.com/louzounlab/DataSets △ Less

Submitted 20 June, 2019; originally announced June 2019.

Comments: 9 pages, 7 figures

ACM Class: E.1; I.2.8; I.5.m; I.2.m

Topological based classification of paper domains using graph convolutional networks

Authors: Idan Benami, Keren Cohen, Oved Nagar, Yoram Louzoun

Abstract: The main approaches for node classification in graphs are information propagation and the association of the class of the node with external information. State of the art methods merge these approaches through Graph Convolutional Networks. We here use the association of topological features of the nodes with their class to predict this class. Moreover, combining topological information with inform… ▽ More The main approaches for node classification in graphs are information propagation and the association of the class of the node with external information. State of the art methods merge these approaches through Graph Convolutional Networks. We here use the association of topological features of the nodes with their class to predict this class. Moreover, combining topological information with information propagation improves classification accuracy on the standard CiteSeer and Cora paper classification task. Topological features and information propagation produce results almost as good as text-based classification, without no textual or content information. We propose to represent the topology and information propagation through a GCN with the neighboring training node classification as an input and the current node classification as output. Such a formalism outperforms state of the art methods. △ Less

Submitted 10 April, 2019; originally announced April 2019.

A combined network and machine learning approaches for product market forecasting

Authors: Jingfang Fan, Keren Cohen, Louis M. Shekhtman, Sibo Liu, Jun Meng, Yoram Louzoun, Shlomo Havlin

Abstract: Sustainable financial markets play an important role in the functioning of human society. Still, the detection and prediction of risk in financial markets remain challenging and draw much attention from the scientific community. Here we develop a new approach based on combined network theory and machine learning to study the structure and operations of financial product markets. Our network links… ▽ More Sustainable financial markets play an important role in the functioning of human society. Still, the detection and prediction of risk in financial markets remain challenging and draw much attention from the scientific community. Here we develop a new approach based on combined network theory and machine learning to study the structure and operations of financial product markets. Our network links are based on the similarity of firms' products and are constructed using the Securities Exchange Commission (SEC) filings of US listed firms. We find that several features in our network can serve as good precursors of financial market risks. We then combine the network topology and machine learning methods to predict both successful and failed firms. We find that the forecasts made using our method are much better than other well-known regression techniques. The framework presented here not only facilitates the prediction of financial markets but also provides insight and demonstrate the power of combining network theory and machine learning. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Family-size variability grows with collapse rate in Birth-Death-Catastrophe model

Authors: N. Dori, H. Behar, H. Brot, Y. Louzoun

Abstract: Forest-fire and avalanche models support the notion that frequent catastrophes prevent the growth of very large populations and as such prevent rare large-scale catastrophes. We show that this notion is not universal. A new model class leads to a paradigm shift in the influence of catastrophes on the family-size distribution of sub-populations. We study a simple population dynamics model where ind… ▽ More Forest-fire and avalanche models support the notion that frequent catastrophes prevent the growth of very large populations and as such prevent rare large-scale catastrophes. We show that this notion is not universal. A new model class leads to a paradigm shift in the influence of catastrophes on the family-size distribution of sub-populations. We study a simple population dynamics model where individuals, as well as a whole family, may die with a constant probability, accompanied by a logistic population growth model. We compute the characteristics of the family-size distribution in steady-state and the phase diagram of the steady-state distribution, and show that the family and catastrophe size variances increase with the catastrophe frequency, which is the opposite of common intuition. Frequent catastrophes are balanced by a larger net-growth rate in surviving families, leading to the exponential growth of these families. When the catastrophe rate is further increased, a second phase transition to extinction occurs, when the rate of new families creations is lower than their destruction rate by catastrophes. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 6 pages 4 figures, plus Supplemental Material 8 pages 6 figues

Journal ref: Phys. Rev. E 98, 012416 (2018)

Non-invasive biomarkers of fetal brain development reflecting prenatal stress: an integrative multi-scale multi-species perspective on data collection and analysis

Authors: Martin G. Frasch, Silvia Lobmaier, Tamara Stampalija, Paula Desplats, María Eugenia Pallarés, Verónica Pastor, Marcela Brocco, Hau-tieng Wu, Jay Schulkin, Christophe Herry, Andrew Seely, Gerlinde A. S. Metz, Yoram Louzoun, Marta Antonelli

Abstract: Prenatal stress (PS) impacts early postnatal behavioural and cognitive development. This process of 'fetal programming' is mediated by the effects of the prenatal experience on the developing hypothalamic-pituitary-adrenal (HPA) axis and autonomic nervous system (ANS). The HPA axis is a dynamic system regulating homeostasis, especially the stress response, and is highly sensitive to adverse early… ▽ More Prenatal stress (PS) impacts early postnatal behavioural and cognitive development. This process of 'fetal programming' is mediated by the effects of the prenatal experience on the developing hypothalamic-pituitary-adrenal (HPA) axis and autonomic nervous system (ANS). The HPA axis is a dynamic system regulating homeostasis, especially the stress response, and is highly sensitive to adverse early life experiences. We review the evidence for the effects of PS on fetal programming of the HPA axis and the ANS. We derive a multi-scale multi-species approach to devising preclinical and clinical studies to identify early non-invasively available pre- and postnatal biomarkers of these programming effects. Such approach would identify adverse postnatal brain developmental trajectories, a prerequisite for designing therapeutic interventions. The multiple scales include the biomarkers reflecting changes in the brain epigenome, metabolome, microbiome and the ANS activity gauged via an array of advanced non-invasively obtainable properties of fetal heart rate fluctuations. The proposed framework has the potential to reveal mechanistic links between maternal stress during pregnancy and changes across these physiological scales. Such biomarkers may hence be useful as early and non-invasive predictors of neurodevelopmental trajectories influenced by the PS. We conclude that studies into PS effects must be conducted on multiple scales derived from concerted observations in multiple animal models and human cohorts performed in an interactive and iterative manner and deploying machine learning for data synthesis, identification and validation of the best non-invasive biomarkers. △ Less

Submitted 31 December, 2017; originally announced January 2018.

Comments: Focused review, 13 pages, 5 figures

Volatility driven market in a generalised Lotka Voltera formalism

Authors: Yoram Louzoun, Sorin Solomon

Abstract: The Generalized Lotka Voltera (GLV) formalism has been introduced in order to explain the power law distributions in the individual wealth (w_i (t)) (Pareto law) and financial markets returns (fluctuations) (r) as a result of the auto-catalytic (multiplicative random) character of the individual capital dynamics. As long as the multiplicative random factor (lambda) is extracted from the same p… ▽ More The Generalized Lotka Voltera (GLV) formalism has been introduced in order to explain the power law distributions in the individual wealth (w_i (t)) (Pareto law) and financial markets returns (fluctuations) (r) as a result of the auto-catalytic (multiplicative random) character of the individual capital dynamics. As long as the multiplicative random factor (lambda) is extracted from the same probability distribution for all the individuals, the exponent of the power laws turns out to be independent on the time variations of the average (<lambda >). This explains also the stability over the past century of experimentally measured Pareto exponent. In contrast to the scaling properties of the single time ("unconditional") probability distributions, the (auto-)correlations between observables measured at different times are not correctly reproduced by the original GLV, if the variance (sigma^2) of $λ$ is time independent. In the GLV formalism the volatility (r^2) auto-correlations decay exponentially while the measurements in real markets indicate a power law with a very small exponent. We show in the present paper that by making the variance of the individual wealth changes (sigma^2) a function of the market volatility <r^2>, one correctly reproduces the market volatility long range correlations. Moreover, we show that this non-trivial feedback loop between the market price volatility and the variance of the investors wealth leads to non-trivial patterns in the overall market trends. If the feedback is too strong, it may even endanger the market stability. △ Less

Submitted 26 September, 2001; originally announced September 2001.

Comments: 8 figures 17 pages(including figures) to be published in physica A

Destruction of CD4 T Lymphocytes Alone Cannot Account for their Long-term Decrease in AIDS

Authors: Yoram Louzoun, Irun. R. Cohen, Henri Atlan

Abstract: Following previous models describing a quasi steady state (QSS) for the evolution of HIV infection and AIDS, we have developed a larger formalism simulating the long-term evolution of the QSS We show that the long-term evolution of AIDS cannot be explained by the destruction alone of CD4 T cells, either directly or indirectly. The destruction of CD4 T cells can lead only to a QSS with a lower co… ▽ More Following previous models describing a quasi steady state (QSS) for the evolution of HIV infection and AIDS, we have developed a larger formalism simulating the long-term evolution of the QSS We show that the long-term evolution of AIDS cannot be explained by the destruction alone of CD4 T cells, either directly or indirectly. The destruction of CD4 T cells can lead only to a QSS with a lower concentration of CD4 T cells, but CD4 destruction cannot generate the sustained long-term decrease in T cells leading to AIDS. We here suggest some workable explanations. △ Less

Submitted 7 August, 2000; originally announced August 2000.

Comments: 5 pages 4 figures

The Emergence of Spatial Complexity in the immune System

Authors: Yoram Louzoun, Sorin Solomon, Henri Atlan, Irun R. Cohen

Abstract: Biological systems, unlike physical or chemical systems, are characterized by the very inhomogeneous distribution of their components. The immune system, in particular, is notable for self-organizing its structure. Classically, the dynamics of natural systems have been described using differential equations. But, differential equation models fail to account for the emergence of large-scale inhom… ▽ More Biological systems, unlike physical or chemical systems, are characterized by the very inhomogeneous distribution of their components. The immune system, in particular, is notable for self-organizing its structure. Classically, the dynamics of natural systems have been described using differential equations. But, differential equation models fail to account for the emergence of large-scale inhomogeneities and for the influence of inhomogeneity on the overall dynamics of biological systems. Here, we show that a microscopic simulation methodology enables us to model the emergence of large-scale objects and to extend the scope of mathematical modeling in biology. We take a simple example from immunology and illustrate that the methods of classical differential equations and microscopic simulation generate contradictory results. Microscopic simulations generate a more faithful approximation of the reality of the immune system. △ Less

Submitted 8 August, 2000; originally announced August 2000.

Adaptation of Autocatalytic Fluctuations to Diffusive Noise

Authors: N. M. Shnerb, E. Bettelheim, Y. Louzoun, O. Agam, S. Solomon

Abstract: Evolution of a system of diffusing and proliferating mortal reactants is analyzed in the presence of randomly moving catalysts. While the continuum description of the problem predicts reactant extinction as the average growth rate becomes negative, growth rate fluctuations induced by the discrete nature of the agents are shown to allow for an active phase, where reactants proliferate as their sp… ▽ More Evolution of a system of diffusing and proliferating mortal reactants is analyzed in the presence of randomly moving catalysts. While the continuum description of the problem predicts reactant extinction as the average growth rate becomes negative, growth rate fluctuations induced by the discrete nature of the agents are shown to allow for an active phase, where reactants proliferate as their spatial configuration adapts to the fluctuations of the catalysts density. The model is explored by employing field theoretical techniques, numerical simulations and strong coupling analysis. For d<=2, the system is shown to exhibits an active phase at any growth rate, while for d>2 a kinetic phase transition is predicted. The applicability of this model as a prototype for a host of phenomena which exhibit self organization is discussed. △ Less

Submitted 6 July, 2000; originally announced July 2000.

Comments: 6 pages 6 figure

Journal ref: Phys. Rev. E 63, 21103 (2001)

Microscopic Discrete Proliferating Components Cause the Self-Organized Emergence of Macroscopic Adaptive Features in Biological Systems

Authors: Yoram Louzoun, Sorin Solomon, Henri Atlan, Irun. R. Cohen

Abstract: We study the emergence of the collective spatio-temporal macroscopic properties of the immune system, by representing individually the elementary interactions between its microscopic components (antibodies, antigens, cytokines). The results of this detailed explicit analysis are compared with the traditional procedure of averaging over individuals and representing the collective dynamics in term… ▽ More We study the emergence of the collective spatio-temporal macroscopic properties of the immune system, by representing individually the elementary interactions between its microscopic components (antibodies, antigens, cytokines). The results of this detailed explicit analysis are compared with the traditional procedure of averaging over individuals and representing the collective dynamics in terms of densities that obey partial differential equations (PDE). The simulations show even for very simple elementary reactions the spontaneous emergence of localized complex structures. In turn the effective dynamics of these structures affects the average behavior of the system in a very decisive way: systems which would according to the differential equations approximation die, display in reality a very lively behavior. Our conclusions are supported both by explicit microscopic simulations and by analytic calculations. As the optimal method we propose a mixture of microscopic simulation (MS) systems describing each reaction separately, and average methods describing the average behavior of the agents. △ Less

Submitted 25 June, 2000; originally announced June 2000.

Comments: 35 pages, 14 pages

HIV time hierarchy: Winning the war while, loosing all the battles

Authors: Uri Hershberg, Yoram Louzoun, Henri Atlan, Sorin Solomon

Abstract: AIDS is the pandemic of our era. A disease that scares us not only because it is fatal but also because its insidious time course makes us all potential carriers long before it hands us our heads in a basket. The strange three stage dynamics of aids is also one of the major puzzles in describing the disease theoretically (1). Aids starts, like most diseases, in a peak of virus expression (2,3),… ▽ More AIDS is the pandemic of our era. A disease that scares us not only because it is fatal but also because its insidious time course makes us all potential carriers long before it hands us our heads in a basket. The strange three stage dynamics of aids is also one of the major puzzles in describing the disease theoretically (1). Aids starts, like most diseases, in a peak of virus expression (2,3), which is practically wiped out by the immune system. However it then remains in the body at a low level of expression until later (some time years later) when there is an outbreak of the disease which terminally cripples the immune system causing death from various common pathogens. In this paper we show, using a microscopic simulation, that the time course of AIDS is determined by the interactions of the virus and the immune cells in the shape space of antigens and that it is the virus's ability to move more rapidly in this space (it's high mutability) that causes the time course and eventual 'victory' of the disease. These results open the way for further experimental and therapeutic conclusions in the ongoing battle with the HIV epidemic. △ Less

Submitted 23 August, 2000; v1 submitted 15 June, 2000; originally announced June 2000.

Comments: 16 pages, TeX 5 eps figures, 1 epsi figure

Modeling the influence of TH1 and TH2 type cells in autoimmune diseases

Authors: Y. Louzoun, H. Atlan, I. R. Cohen

Abstract: A sharp TH1/TH2 dichotomy has often been used to define the effects of cytokines on autoimmune diseases. However contradictory results in recent research indicate that the situation may be more complex. We build here a simple mathematical model aimed at settling the contradictions. The model is based on a neural network paradigm, and is applied using Partial Differential Equations (PDE). We show… ▽ More A sharp TH1/TH2 dichotomy has often been used to define the effects of cytokines on autoimmune diseases. However contradictory results in recent research indicate that the situation may be more complex. We build here a simple mathematical model aimed at settling the contradictions. The model is based on a neural network paradigm, and is applied using Partial Differential Equations (PDE). We show here that a TH1/TH2 paradigm is only an external view of a complex multivariate system. △ Less

Submitted 19 June, 2000; originally announced June 2000.

Comments: 24 pages, 8 figures, 1 table

MSC Class: 92B08

The Importance of Being Discrete - Life Always Wins on the Surface

Authors: Nadav M. Shnerb, Yoram Louzoun, Eldad Bettelheim, Sorin Solomon

Abstract: Many systems in chemistry, biology, finance and social sciences present emerging features which are not easy to guess from the elementary interactions of their microscopic individual components. In the past, the macroscopic behavior of such systems was modeled by assuming that the collective dynamics of microscopic components can be effectively described collectively by equations acting on spati… ▽ More Many systems in chemistry, biology, finance and social sciences present emerging features which are not easy to guess from the elementary interactions of their microscopic individual components. In the past, the macroscopic behavior of such systems was modeled by assuming that the collective dynamics of microscopic components can be effectively described collectively by equations acting on spatially continuous density distributions. It turns out that quite contrary, taking into account the actual individual/discrete character of the microscopic components of these systems is crucial for explaining their macroscopic behavior. In fact, we find that in conditions in which the continuum approach would predict the extinction of all the population (respectively the vanishing of the invested capital or of the concentration of a chemical substance, etc), the microscopic granularity insures the emergence of macroscopic localized sub-populations with collective adaptive properties which allow their survival and development. In particular it is found that in 2 dimensions "life" (the localized proliferating phase) always prevails. △ Less

Submitted 29 December, 1999; originally announced December 1999.

Comments: 4 pages, 4 figures

Journal ref: P NATL ACAD SCI USA 97: (19) 10322-10324 SEP 12 2000