Past Networks Seminar

In many applications we are confronted with the following scenario: we observe snapshots of data describing the state of a system at particular times, and based on these observations we want to infer the (dynamical) interactions between the entities we observe. However, often the number of samples we can obtain from such a process are far too few to identify the network exactly. Can we still reliable infer some aspects of the underlying system?
Motivated by this question we consider the following alternative system identification problem: instead of trying to infer the exact network, we aim to recover a (low-dimensional) statistical model of the network based on the observed signals on the nodes.  More concretely, here we focus on observations that consist of snapshots of a diffusive process that evolves over the unknown network. We model the (unobserved) network as generated from an independent draw from a latent stochastic block model (SBM), and our goal is to infer both the partition of the nodes into blocks, as well as the parameters of this SBM. We present simple spectral algorithms that provably solve the partition and parameter inference problems with high-accuracy.

Many soft materials, such as elastomers and hydrogels, are made of long chain molecules crosslinked to form a three-dimensional network. Their mechanical properties depend on network parameters such as chain density, chain length distribution and the functionality of the crosslinks. Understanding the relationships between the topology of polymer networks and their mechanical properties has been a long-standing challenge in polymer physics.

In this work, we focus on so-called “near-ideal” networks, which are produced by the cross-coupling of star-like macromolecules with well-defined chain length. We developed a computational approach based on random discrete networks, according to which the polymer network is represented by an assembly of non-linear springs connected at junction points representing crosslinks. The positions of the crosslink points are determined from the conditions of mechanical equilibrium. Scaling relations for the elastic modulus and maximum extensibility of the network were obtained. Our scaling relations contradict some predictions of classical estimates of rubber elasticity and have implications for the interpretation of experimental data for near-ideal polymer networks.

Reference: G. Alame, L. Brassart. Relative contributions of chain density and topology to the elasticity of two-dimensional polymer networks. Soft Matter 15, 5703 (2019).

In an increasingly urbanized world, where cities are changing continuously, it is essential for policy makers to have access to regularly updated decision-making tools for an effective management of urban areas. An example of these tools is the delineation of cities into functional areas which provides knowledge on high spatial interaction zones and their socioeconomic composition. In this paper, we presented a method for the structural analysis of a city, specifically for the determination of its functional areas, based on communities detection in graphs. The nodes of the graph correspond to geographical units resulting from a cartographic division of the city according to the road network. The edges are weighted using a Gaussian distance-decay function and the amount of spatial interactions between nodes. Our approach optimize the modularity to ensure that the functional areas detected have strong interactions within their borders but lower interactions outside. Moreover, it leverages on POIs' entropy to maintain a good socioeconomic heterogeneity in the detected areas. We conducted experiments using taxi trips and POIs datasets from the city of Porto, as a study case. Trough those experiments, we demonstrate the ability of our method to portray functional areas while including spatial and socioeconomic dynamics.
 

The analysis and characterization of human mobility using population-level mobility models is important for numerous applications, ranging from the estimation of commuter flows to modeling trade flows. However, almost all of these applications have focused on large spatial scales, typically from intra-city level to inter-country level. In this paper, we investigate population-level human mobility models on a much smaller spatial scale by using them to estimate customer mobility flow between supermarket zones. We use anonymized mobility data of customers in supermarkets to calibrate our models and apply variants of the gravity and intervening-opportunities models to fit this mobility flow and estimate the flow on unseen data. We find that a doubly-constrained gravity model can successfully estimate 65-70% of the flow inside supermarkets. We then investigate how to reduce congestion in supermarkets by combining mobility models with queueing networks. We use a simulated-annealing algorithm to find store layouts with lower congestion than the original layout. Our research gives insight both into how customers move in supermarkets and into how retailers can arrange stores to reduce congestion. It also provides a case study of human mobility on small spatial scales.

We investigate a nonlinear version of coevolving voter models, in which both node states and network structure update as a coupled stochastic dynamical process. Most prior work on coevolving voter models has focused on linear update rules with fixed rewiring and adopting probabilities. By contrast, in our nonlinear version, the probability that a node rewires or adopts is a function of how well it "fits in" within its neighborhood. To explore this idea, we incorporate a parameter σ that represents the fraction of neighbors of an updating node that share its opinion state. In an update, with probability σq (for some nonlinearity parameter q), the updating node rewires; with complementary probability 1−σq, the updating node adopts a new opinion state. We study this mechanism using three rewiring schemes: after an updating node deletes a discordant edge, it then either (1) "rewires-to-random" by choosing a new neighbor in a random process; (2) "rewires-to-same" by choosing a new neighbor in a random process from nodes that share its state; or (3) "rewires-to-none" by not rewiring at all (akin to "unfriending" on social media). We compare our nonlinear coevolving model to several existing linear models, and we find in our model that initial network topology can play a larger role in the dynamics, whereas the choice of rewiring mechanism plays a smaller role. A particularly interesting feature of our model is that, under certain conditions, the opinion state that is initially held by a minority of nodes can effectively spread to almost every node in a network if the minority nodes views themselves as the majority. In light of this observation, we relate our results to recent work on the majority illusion in social networks.

Reference: 

Kureh, Yacoub H., and Mason A. Porter. "Fitting In and Breaking Up: A Nonlinear Version of Coevolving Voter Models." arXiv preprint arXiv:1907.11608 (2019).

The vulnerability of a host population to a specific disease measures how likely pathogen introduction will lead to an epidemic outbreak, and how hard it is to contain or eliminate an ongoing one. Predicting vulnerability is thus key to designing risk-reduction strategies that limit disease burden on public health and economic development. To do that, highly-resolved data tracking contacts and mobility of the host population need to integrate into detailed models of disease dynamics. This represents a twofold challenge. Firstly, we need theoretical frameworks that turn data feeds into predictors of epidemic risk, and can identify which of the structural features of the host population drive its vulnerability. Secondly, we need new ways to access, analyze, and share the relevant contact and mobility data: a necessary step to make our predictions realistic and reliable. In my talk, I will address both issues. I will show how to analytically derive the conditions that discriminate between epidemic regime and quick pathogen extinction, by representing empirically measured contacts as time-evolving complex networks. The analytical core of this theory leads to a broad range of applications. At the same time, its data-driven nature prompts context-specific predictions that can inform policymaking, as I will show in two case studies: reorganizing nurse scheduling to reduce the risk of spread of healthcare-associated infections; linking the features of livestock trade movements to the spatial spread of cattle diseases. The latter application is also an example of how limited access and incomplete data collection represent a big hurdle to predictive vulnerability analysis. To overcome this, I will present a collaborative platform for analyzing and comparing trade networks coming from several European countries. Using a bring code to the data approach, our platform surmounts the strict regulations preventing data sharing, and builds an algorithm that predicts vulnerability even in situations when limited data on cattle trade are available. The ultimate goal of all these theoretical and numerical developments is to inform strategies that reduce the vulnerability of the host population by restructuring its contacts. However, such restructuring may entail a feedback effect, acting as selective pressure on the pathogen itself. In the last part of my talk, I will extend the developed formalism to modeling evolutionary pathways that maximize the invasion potential of the pathogen, given the observed host population structure. Specifically, I will link the emergence of exotic replication behaviors in plant-infecting viruses to historical changes in plant distribution patterns.

In recent years, much attention has been given to single-cell RNA sequencing techniques as they allow researchers to examine the functions and relationships of single cells inside a tissue. In this study, we combine single-cell RNA sequencing data with protein–protein interaction networks (PPINs) to detect active modules in cells of different transcriptional states. We achieve this by clustering single-cell RNA sequencing data, constructing node-weighted PPINs, and identifying the maximum-weight connected subgraphs with an exact Steiner-Tree approach. As a case study, we investigate RNA sequencing data from human liver spheroids but the techniques described here are applicable to other organisms and tissues. The benefits of our novel method are two-fold: First, it allows us to identify important proteins (e.g., receptors) which are not detected from a differential gene-expression analysis as they only interact with proteins that are transcribed in higher levels. Second, we find that different transcriptional states have different subnetworks of the PPIN significantly overexpressed. These subnetworks often reflect known biological pathways (e.g., lipid metabolism and stress response) and we obtain a nuanced picture of cellular function as we can associate them with a subset of all analysed cells.

Short- and long-term memories are distinguished by their forgettability. Most of what we perceive and store is lost rather quickly to noise, as new sensations replace older ones, while some memories last for as long as we live. Synaptic dynamics is key to the process of memory storage; in this talk I will discuss a few approaches we have taken to this problem, culminating in a model of synaptic networks containing both cooperative and competitive dynamics. It turns out that the competitionbetween synapses is key to the natural emergence of long-term memory in this model, as in reality.

References
​Mehta, Anita. "Storing and retrieving long-term memories: cooperation and competition in synaptic dynamics." Advances in Physics: X 3.1 (2018): 1480415.

The comparison of graphs is a vitally important, yet difficult task which arises across a number of diverse research areas including biological and social networks. There have been a number of approaches to define graph distance however often these are not metrics (rendering standard data-mining techniques infeasible), or are computationally infeasible for large graphs. In this work, we define a new metric based on the spectrum of the non-backtracking graph operator and show that it can not only be used to compare graphs generated through different mechanisms but can reliably compare graphs of varying size. We observe that the family of Watts-Strogatz graphs lie on a manifold in the non-backtracking spectral embedding and show how this metric can be used in a standard classification problem of empirical graphs.

Complexity Science aims to understand what is that makes some systems to be "more than the sum of their parts". A natural first step to address this issue is to study networks of pairwise interactions, which have been done with great success in many disciplines -- to the extend that many people today identify Complexity Science with network analysis. In contrast, multivariate complexity provides a vast and mostly unexplored territory. As a matter of fact, the "modes of interdependency" that can exist between three or more variables are often nontrivial, poorly understood and, yet, are paramount for our understanding of complex systems in general, and emergence in particular. 
In this talk we present an information-theoretic framework to analyse high-order correlations, i.e. statistical dependencies that exist between groups of variables that cannot be reduced to pairwise interactions. Following the spirit of information theory, our approach is data-driven and model-agnostic, being applicable to discrete, continuous, and categorical data. We review the evolution of related ideas in the context of theoretical neuroscience, and discuss the most prominent extensions of information-theoretic metrics to multivariate settings. Then, we introduce the O-information, a novel metric that quantify various structural (i.e. synchronous) high-order effects. Finally, we provide a critical discussion on the framework of Integrated Information Theory (IIT), which suggests an approach to extend the analysis to dynamical settings. To illustrate the presented methods, we show how the analysis of high-order correlations can reveal critical structures in various scenarios, including cellular automata, Baroque music scores, and various EEG datasets.

References:
[1] F. Rosas, P.A. Mediano, M. Gastpar and H.J. Jensen, ``Quantifying High-order Interdependencies via Multivariate Extensions of the Mutual Information'', submitted to PRE, under review.
https://arxiv.org/abs/1902.11239
[2] F. Rosas, P.A. Mediano, M. Ugarte and H.J. Jensen, ``An information-theoretic approach to self-organisation: Emergence of complex interdependencies in coupled dynamical systems'', in Entropy, vol. 20 no. 10: 793, pp.1-25, Sept. 2018.
https://www.mdpi.com/1099-4300/20/10/793

 

Coupling dynamics of the states of the nodes of a network to the dynamics of the network topology leads to generic absorbing and fragmentation transitions. The coevolving voter model is a typical system that exhibits such transitions at some critical rewiring. We study the robustness of these transitions under two distinct ways of introducing noise. Noise affecting all the nodes destroys the absorbing-fragmentation transition, giving rise in finite-size systems to two regimes: bimodal magnetization and dynamic fragmentation. Noise targeting a fraction of nodes preserves the transitions but introduces shattered fragmentation with its characteristic fraction of isolated nodes and one or two giant components. Both the lack of absorbing state for homogeneous noise and the shift in the absorbing transition to higher rewiring for targeted noise are supported by analytical approximations.

Paper Link:

https://journals.aps.org/pre/abstract/10.1103/PhysRevE.92.032803

Social surveys are widely used in today's society as a method for obtaining opinions and other information from large groups of people. The questions in social surveys are usually presented in either multiple-choice or free-response formats. Despite their advantages, free-response questions are employed less commonly in large-scale surveys, because in such situations, considerable effort is needed to categorise and summarise the resulting large dataset. This is the so-called coding problem. Here we propose a survey framework in which, respondents not only write down their own opinions, but also input information characterising the similarity between their individual responses and those of other respondents. This is done in much the same way as ``likes" are input in social network services. The information input in this simple procedure constitutes relational data among opinions, which we call the opinion graph. The diversity of typical opinions can be identified as a modular structure of such a graph, and the coding problem is solved through graph clustering in a statistically principled manner. We demonstrate our approach using a poll on the 2016 US presidential election and a survey given to graduates of a particular university.

George Soules [1] introduced a set of vectors $r_1,...,r_N$ with the remarkable property that for any set of ordered numbers $\lambda_1\geq\dots\geq\lambda_N$, the matrix $\sum_n \lambda_nr_nr_n^T$ has nonnegative off-diagonal entries. Later, it was found [2] that there exists a whole class of such vectors - Soules vectors - which are intimately connected to binary rooted trees. In this talk I will describe the construction of Soules vectors starting from a binary rooted tree, and introduce some basic properties. I will also cover a number of applications: the inverse eigenvalue problem, equitable partitions in Laplacian matrices and the eigendecomposition of the Clauset-Moore-Newman hierarchical random graph model.

[1] Soules (1983), Constructing Symmetric Nonnegative Matrices
[2] Elsner, Nabben and Neumann (1998), Orthogonal bases that lead to symmetric nonnegative matrices

Commuting concerns people’s spatial behaviour resulting from the geographic separation of home and workplace and is connected with their willingness to seek economic opportunities outside their place of residence (Rouwendal J., Nijkamp P., 2004). Such opportunities are usually found in the urban areas, so this phenomenon is often a subject of urban studies or research focusing on city centres (Drejerska N., Chrzanowska M., 2014). In literature, commuting patterns are used to determine the boundaries of local and regional labour markets. Furthermore, labour market is one of the most important features for the delimitation of functional regions, as commuting involves not only working outside one’s place of residence but also, among other things, using various services offered there, from shopping to health or cultural services. Taking this into account, it can be stated that commuting is an important characteristic of relations between territories, and these relations form complex networks.

People decide to commute to work for various reasons. Most commuters travel from a small town, village or rural area to a city or town where they have a wider range of employment opportunities. However, people differ in their attitudes toward commuting. While some people find it troublesome, others enjoy their daily travel. There are also people who regard commuting as the necessary condition for supporting themselves and their families. Therefore, commuting is an important factor that should be taken into account in the research on the quality of life and quality of work.

The main goals of this presentation is to identify and analyse relations between communities (municipalities) from the perspective of labour market, especially commuting in the vicinity of Warsaw, Data on the number of commuters come from the Central Statistical Office of Poland and cover the year 2011.

 Bibliography

Drejerska N., Chrzanowska M., 2014: Commuting in the Warsaw suburban area from a spatial perspective – an example of empirical research, Acta Universitatis Lodziensis. Folia Oeconomica 2014, Vol. 6, no 309, pp. 87-96.

Rouwendal J., Nijkamp P., 2004: Living in Two Worlds: A Review of Home-to-Work Decisions, Growth and Change, Volume 35, Issue 3, p. 287.

The advent of social media and the resulting ability to instantaneously communicate ideas and messages to connections worldwide is one of the great consequences arising from the telecommunications revolution over the last century. Individuals do not, however, communicate only upon a single platform; instead there exists a plethora of options available to users, many of whom are active on a number of such media. While each platform offers some unique selling point to attract users, e.g., keeping up to date with friends through messaging and statuses (Facebook), photo sharing (Instagram), seeing information from friends, celebrities and numerous other outlets (Twitter) or keeping track of the career paths of friends and past colleagues (Linkedin), the platforms are all based upon the fundamental mechanisms of connecting with other users and transmitting information to them as a result of this link.

In this talk a model for the spreading of online information or “memes" on multiplex networks is introduced and analyzed using branching-process methods. The model generalizes that of [Gleeson et al., Phys. Rev. X., 2016] in two ways. First, even for a monoplex (single-layer) network, the model is defined for any specific network defined by its adjacency matrix, instead of being restricted to an ensemble of random networks. Second, a multiplex version of the model is introduced to capture the behavior of users who post information from one social media platform to another. In both cases the branching process analysis demonstrates that the dynamical system is, in the limit of low innovation, poised near a critical point, which is known to lead to heavy-tailed distributions of meme popularity similar to those observed in empirical data.

[1] J. P. Gleeson et al. “Effects of network structure, competition and memory time on social spreading phenomena”. Physical Review X 6.2 (2016), p. 021019.

[2] J. D. O’Brien et al. "Spreading of memes on multiplex networks." New Journal of Physics 21.2 (2019): 025001.

With the introduction of supermarket loyalty cards in recent decades, there has been an ever-growing body of customer-level shopping data. A natural way to represent this data is with a bipartite network, in which customers are connected to products that they purchased. By predicting likely edges in these networks, one can provide personalised product recommendations to customers.
In this talk, I will first discuss a basic approach for recommendations, based on network community detection, that we have validated on a promotional campaign run by our industrial collaborators. I will then describe a multilayer network model that accounts for the fact that customers tend to buy the same grocery items repeatedly over time. By modelling such correlations explicitly, link-prediction accuracy improves considerably. This approach is also useful in other networks that exhibit significant edge correlations, such as social networks (in which people often have repeated interactions with other people), airline networks (in which popular routes are often served by more than one airline), and biological networks (in which, for example, proteins can interact in multiple ways). 
 

We propose a novel class of network models for temporal dyadic interaction data. Our objective is to capture important features often observed in social interactions: sparsity, degree heterogeneity, community structure and reciprocity. We use mutually-exciting Hawkes processes to model the interactions between each (directed) pair of individuals. The intensity of each process allows interactions to arise as responses to opposite interactions (reciprocity), or due to shared interests between individuals (community structure). For sparsity and degree heterogeneity, we build the non time dependent part of the intensity function on compound random measures following (Todeschini et al., 2016). We conduct experiments on real- world temporal interaction data and show that the proposed model outperforms competing approaches for link prediction, and leads to interpretable parameters.

Link to paper: https://papers.nips.cc/paper/7502-modelling-sparsity-heterogeneity-reci…

Many real-world networks contain small recurring connectivity patterns also known as network motifs. Although network motifs are widely considered to be important structural features of networks that are closely connected to their function methods for characterizing and modelling the local connectivity structure of complex networks remain underdeveloped. In this talk, we will present a non-parametric approach that is based on generative models in which networks are generated by adding not only single edges but also but also copies of larger subgraphs such as triangles to the graph. We show that such models can be formulated in terms of latent states that correspond to subgraph decompositions of the network and derive analytic expressions for the likelihood of such models. Following a Bayesian approach, we present a nonparametric prior for model parameters. Solving the resulting inference problem results in a principled approach for identifying atomic connectivity patterns of networks that do not only identify statistically significant connectivity patterns but also produces a decomposition of the network into such atomic substructures. We tested the presented approach on simulated data for which the algorithm recovers the latent state to a high degree of accuracy. In the case of empirical networks, the method identifies concise sets atomic subgraphs from within thousands of candidates that are plausible and include known atomic substructures.

There have been two main attempts so far to forecast the level of development of artificial intelligence (or ‘computerisation’) over time, Frey and Osborne (2013, 2017) and Manyika et al (2017). Unfortunately, their methodology seems to be flawed. Their results depend upon expert predictions of which occupations will be automatable in 2050, but these predictions are notoriously unreliable. Therefore, we develop an alternative which does not depend upon these expert predictions. We build a dataset of all the start-ups, firms, and university research laboratories working on automating different types of tasks, and use this to build a dynamic network model of them and how they interact. How automatable each type of task is ‘emerges’ from the model. We validate it, predicting the level of development of supervised learning in 2017 using data from the year 2000, and use it to forecast of the automatability of each of these task types from 2018 to 2050. Finally, we discuss extensions for our model; how it could be used to test the impact of public policy decisions or forecast developments in other high-technology industries.

We examine the relationship between social structure and sentiment through the analysis of a large collection of tweets about the Irish Marriage Referendum of 2015. We obtain the sentiment of every tweet with the hashtags #marref and #marriageref that was posted in the days leading to the referendum, and construct networks to aggregate sentiment and use it to study the interactions among users. Our analysis shows that the sentiment of outgoing mention tweets is correlated with the sentiment of incoming mentions, and there are significantly more connections between users with similar sentiment scores than among users with opposite scores in the mention and follower networks. We combine the community structure of the follower and mention networks with the activity level of the users and sentiment scores to find groups that support voting ‘yes’ or ‘no’ in the referendum. There were numerous conversations between users on opposing sides of the debate in the absence of follower connections, which suggests that there were efforts by some users to establish dialogue and debate across ideological divisions. Our analysis shows that social structure can be integrated successfully with sentiment to analyse and understand the disposition of social media users around controversial or polarizing issues. These results have potential applications in the integration of data and metadata to study opinion dynamics, public opinion modelling and polling.

The problem of authorship attribution (AA) involves matching a text of unknown authorship with its creator, found among a pool of candidate authors. In this work, we examine in detail authorship attribution methods that rely on networks of function words to detect an “authorial fingerprint” of literary works. Previous studies interpreted these word adjacency networks (WANs) as Markov chains, giving transition rates between function words, and they compared them using information-theoretic measures. Here, we apply a variety of network flow-based tools, such as role-based similarity and community detection, to perform a direct comparison of the WANs. These tools reveal an interesting relation between communities of function words and grammatical categories. Moreover, we propose two new criteria for attribution based on the comparison of connectivity patterns and the similarity of network partitions. The results are positive, but importantly, we observe that the attribution context is an important limiting factor that is often overlooked in the field's literature. Furthermore, we give important new directions that deserve further consideration.

Joint work with Zsolt Vizi (Bolyai Institute, University of Szeged, Hungary), Istvan Kiss (Department
of Mathematics, University of Sussex, United Kingdom)

Pairwise models have been proven to be a flexible framework for analytical approximations
of stochastic epidemic processes on networks that are in many situations much more accurate
than mean field compartmental models. The non-Markovian aspects of disease transmission
are undoubtedly important, but very challenging to incorporate them into both numerical
stochastic simulations and analytical investigations. Here we present a generalization of
pairwise models to non-Markovian epidemics on networks. For the case of infectious periods
of fixed length, the resulting pairwise model is a system of delay differential equations, which
shows excellent agreement with results based on the explicit stochastic simulations. For more
general distribution classes (uniform, gamma, lognormal etc.) the resulting models are PDEs
that can be transformed into systems of integro-differential equations. We derive pairwise
reproduction numbers and relations for the final epidemic size, and initiate a systematic
study of the impact of the shape of the particular distributions of recovery times on how
the time evolution of the disease dynamics play out.

Crime is a major risk to society’s well-being, particularly in cities, and yet the scientific literature lacks a comprehensive statistical characterization of crime that could uncover some of the mechanisms behind such pervasive social phenomenon. Evidence of nonlinear scaling of urban indicators in cities, such as wages and serious crime, has motivated the understanding of cities as complex systems—a perspective that offers insights into resources limits and sustainability, but usually without examining the details of indicators. Notably, since the nineteenth century, criminal activities have been known not to occur uniformly within a city. Crime concentrates in such way that most of the offenses take place in few regions of the city. However, though this concentration is confirmed by different studies, the absence of broad examinations of the characteristics of crime concentration hinders not only the comprehension of crime dynamics but also the proposal of sounding counter-measures. Here, we developed a framework to characterize crime concentration which splits cities into regions with the same population size. We used disaggregated criminal data from 25 locations in the U.S. and the U.K. which include offenses in places spanning from 2 to 15 years of data. Our results confirmed that crime concentrates regardless of city and revealed that the level of concentration does not scale with city size. We found that distribution of crime in a city can be approximated by a power-law distribution with exponent α that depends on the type of crime. In particular, our results showed that thefts tend to concentrate more than robberies, and robberies more than burglaries. Though criminal activities present regularities of concentration, we found that criminal ranks have the tendency to change continuously over time. Such features support the perspective of crime as a complex system which demands analyses and evolving urban policies covering the city as a whole. 

Disease transmission and rumour spreading are ubiquitous in social and technological networks. In this talk, we will present our last results on the modelling of rumour and disease spreading in multilayer networks.  We will derive analytical expressions for the epidemic threshold of the susceptible-infected-susceptible (SIS) and susceptible-infected-recovered dynamics, as well as upper and lower bounds for the disease prevalence in the steady state for the SIS scenario. Using the quasistationary state method, we numerically show the existence of disease localization and the emergence of two or more susceptibility peaks in a multiplex network. Moreover, we will introduce a model of epidemic spreading with awareness, where the disease and information are propagated in different layers with different time scales. We will show that the time scale determines whether the information awareness is beneficial or not to the disease spreading. 

Mechanical percolation is a phenomenon in materials processing wherein ‘filler’ rod-like particles are incorporated into polymeric materials to enhance the composite’s mechanical properties. Experiments have well-characterized a nonlinear phase transition from floppy to rigid behavior at a threshold filler concentration, but the underlying mechanism is not well understood. We develop and utilize an iterative graph compression algorithm to demonstrate that this experimental phenomenon coincides with the formation of a spatially extending set of mutually rigid rods (‘rigidity percolation’). First, we verify the efficacy of this method in two-dimensional fiber systems (intersecting line segments), then moving to the more interesting and mechanically representative problem of three-dimensional fiber systems (cylinders). We show that, when the fibers are uniformly distributed both spatially and orientationally, the onset of rigidity percolation appears to co-occur with a mean field prediction that is applicable across a wide range of aspect ratios.

There is a mismatch between levels of crime and its fear and often, cities might see an increase or a decrease in crime over time while the fear of crime remains unchanged. A model that considers fear of crime as an opinion shared by simulated individuals on a network will be presented, and the impact that different distributions of crime have on the fear experienced by the population will be explored. Results show that the dynamics of the fear is sensitive to the distribution of crime and that there is a phase transition for high levels of concentration of crime.

Matrix completion is an area of great mathematical interest and has numerous applications, including recommender systems for e-commerce. The recommender problem can be viewed as follows: given a database where rows are users and and columns are products, with entries indicating user preferences, fill in the entries so as to be able to recommend new products based on the preferences of other users. Viewing the interactions between user and product as links in a bipartite graph, the problem is equivalent to approximating a partially observed graph using clusters. We propose a divide and conquer algorithm inspired by the work of [1], who use recursive rank-1 approximation. We make the case for using an LP rank-1 approximation, similar to that of [2] by a showing that it guarantees a 2-approximation to the optimal, even in the case of missing data. We explore our algorithm's performance for different test cases.

[1]  Shen, B.H., Ji, S. and Ye, J., 2009, June. Mining discrete patterns via binary matrix factorization. In Proceedings of the 15th ACM SIGKDD international conference on Knowledge discovery and data mining (pp. 757-766). ACM.

[2] Koyutürk, M. and Grama, A., 2003, August. PROXIMUS: a framework for analyzing very high dimensional discrete-attributed datasets. In Proceedings of the ninth ACM SIGKDD international conference on Knowledge discovery and data mining (pp. 147-156). ACM.
 

In this work, study the mean first saturation time (MFST), a generalization to the mean first passage time, on networks and show an application to the 2015 Burundi refugee crisis. The MFST between a sink node j, with capacity s, and source node i, with n random walkers, is the average number of time steps that it takes for at least s of the random walkers to reach a sink node j. The same concept, under the name of extreme events, has been studied in previous work for degree biased-random walks [2]. We expand the literature by exploring the behaviour of the MFST for node-biased random walks [1] in Erdős–Rényi random graph and geographical networks. Furthermore, we apply MFST framework to study the distribution of refugees in camps for the 2015 Burundi refugee crisis. For this last application, we use the geographical network of the Burundi conflict zone in 2015 [3]. In this network, nodes are cities or refugee camps, and edges denote the distance between them. We model refugees as random walkers who are biased towards the refugee camps which can hold s_j people. To determine the source nodes (i) and the initial number of random walkers (n), we use data on where the conflicts happened and the number of refugees that arrive at any camp under a two-month period after the start of the conflict [3]. With such information, we divide the early stage of the Burundi 2015 conflict into two waves of refugees. Using the first wave of refugees we calibrate the biased parameter β of the random walk to best match the distribution of refugees on the camps. Then, we test the prediction of the distribution of refugees in camps for the second wave using the same biased parameters. Our results show that the biased random walk can capture, to some extent, the distribution of refugees in different camps. Finally, we test the probability of saturation for various camps. Our model suggests the saturation of one or two camps (Nakivale and Nyarugusu) when in reality only Nyarugusu camp saturated.

[1] Sood, Vishal, and Peter Grassberger. ”Localization transition of biased random walks on random
networks.” Physical review letters 99.9 (2007): 098701.
[2] Kishore, Vimal, M. S. Santhanam, and R. E. Amritkar. ”Extreme event-size fluctuations in biased
random walks on networks.” arXiv preprint arXiv:1112.2112 (2011).
[3] Suleimenova, Diana, David Bell, and Derek Groen. ”A generalized simulation development approach
for predicting refugee destinations.” Scientific reports 7.1 (2017): 13377.

Graphs are a central object of study in various scientific fields, such as discrete mathematics, theoretical computer science and network science. These graphs are typically studied using combinatorial, algebraic or probabilistic methods, each of which highlights the properties of graphs in a unique way. I will discuss a novel approach to study graphs: the simplex geometry (a simplex is a generalized triangle). This perspective, proposed by Miroslav Fiedler, introduces techniques from (simplex) geometry into the field of graph theory and conversely, via an exact correspondence. We introduce the graph-simplex correspondence, identify a number of basic connections between graph characteristics and simplex properties, and suggest some applications as example.

Reference: https://arxiv.org/abs/1807.06475
 

Many protein interaction databases provide confidence scores based on the experimental evidence underpinning each in- teraction. The databases recommend that protein interac- tion networks (PINs) are built by thresholding on these scores. We demonstrate that varying the score threshold can re- sult in PINs with significantly different topologies. We ar- gue that if a node metric is to be useful for extracting bio- logical signal, it should induce similar node rankings across PINs obtained at different thresholds. We propose three measures—rank continuity, identifiability, and instability— to test for threshold robustness. We apply these to a set of twenty-five metrics of which we identify four: number of edges in the step-1 ego network, the leave-one-out dif- ference in average redundancy, average number of edges in the step-1 ego network, and natural connectivity, as robust across medium-high confidence thresholds. Our measures show good agreement across PINs from different species and data sources. However, analysis of synthetically gen- erated scored networks shows that robustness results are context-specific, and depend both on network topology and on how scores are placed across network edges. 

Fake accounts detection and users’ polarization are two very well known topics concerning the social media sphere, that have been extensively discussed and analyzed, both in the academic literature and in everyday life. Social bots are autonomous accounts that are explicitly created to increase the number of followers of a target user, in order to inflate its visibility and consensus in a social media context. For this reason, a great variety of methods for their detection have been proposed and tested. Polarisation, also known as confirmation bias, is instead the common tendency to look for information that confirms one's preexisting beliefs, while ignoring opposite ones. Within this environment, groups of individuals characterized by the same system of beliefs are very likely to form. In the present talk we will first review part of the literature discussing both these topics. Then we will focus on a new dataset collecting tweets from the last Italian parliament elections in 2018 and some preliminary results will be discussed.

Neural networks underpin both biological intelligence and modern AI systems, yet there is relatively little theory for how the observed behavior of these networks arises. Even the connectivity of neurons within the brain remains largely unknown, and popular deep learning algorithms lack theoretical justification or reliability guarantees.  In this talk, we consider paths towards a more rigorous understanding of neural networks. We characterize and, where possible, prove essential properties of neural algorithms: expressivity, learning, and robustness. We show how observed emergent behavior can arise from network dynamics, and we develop algorithms for learning more about the network structure of the brain.

Cascading phenomena are prevalent in natural and social-technical complex networks. We study the persistent cascade-recovery dynamics on random networks which are robust against small trigger but may collapse for larger one. It is observed that depending on the relative intensity of triggering and recovery, the network belongs one of the two dynamical phases: collapsing or active phase. We devise an analytical framework which characterizes not only the critical behaviour but also the temporal evolution of network activity in both phases. Results from agent-based simulations show good agreement with theoretical calculations. This work is an important attempt in understanding networked systems gradually evolving into a state of critical transition, with many potential applications.
 

Several scholars of evolutionary biology have suggested that functional redundancy (also known as "biological degener-
acy") is important for robustness of biological networks. Structural redundancy indicates the existence of structurally
similar subsystems that can perform the same function. Functional redundancy indicates the existence of structurally
di erent subsystems that can perform the same function. For networks with Ornstein--Uhlenbeck dynamics, Tononi et al.
[Proc. Natl. Acad. Sci. U.S.A. 96, 3257{3262 (1999)] proposed measures of structural and functional redundancy that are
based on mutual information between subnetworks. For a network of n vertices, an exact computation of these quantities
requires O(n!) time. We derive expansions for these measures that one can compute in O(n3) time. We use the expan-
sions to compare the contributions of di erent types of motifs to a network's functional redundancy.

Wikipedia has more than 40 million articles in 280 languages. It represents a decent coverage of human knowledge.
Even with its biases it can tell us a lot about what's important for people. London has an article in 238 languages and
Swansea has in 73 languages. Is London more "culturally" important than Swansea? Probably. 
We use this information and look at various factors that could help us model "cultural" importance of a city and hence
try to find the driving force behind sister city relationships.
We also look at creating cultural maps of different cities, finding the artsy/hipster, academic, political neighbourhoods of a city.

With the advent of large-scale data and the concurrent development of robust scientific tools to analyze them, important discoveries are being made in a wider range of scientific disciplines than ever before. A field of research that has gained substantial attention recently is the analytical, large-scale study of human behavior, where many analytical and statistical techniques are applied to various behavioral data from online social media, markets, and mobile communication, enabling meaningful strides in understanding the complex patterns of humans and their social actions.

The importance of such research originates from the social nature of humans, an essential human nature that clearly needs to be understood to ultimately understand ourselves. Another essential human nature is that they are creative beings, continually expressing inspirations or emotions in various physical forms such as a picture, sound, or writing. As we are successfully probing the social behaviors humans through science and novel data, it is natural and potentially enlightening to pursue an understanding of the creative nature of humans in an analogous way. Further, what makes such research even more potentially beneficial is that human creativity has always been in an interplay of mutual influence with the scientific and technological advances, being supplied with new tools and media for creation, and in return providing valuable scientific insights.

In this talk I will present two recent ongoing works on the mathematical analysis of color contrast in painting and measuring novelty in piano music.

The story of human progress is often described as a succession of ‘eras’ or ‘ages’ that are characterised by their most dominant technologies (e.g., the bronze age, the industrial revolution or the information age). In modern times, the fast pace of technological progress has accelerated the succession of eras. In addition, the increasing complexity of inventions has made the task of determining when eras begin and end more challenging, as eras are less about the dominance of a single technology and more about the way in which different technologies are combined. We present a data-driven method to determine and uncover technological eras based on networks and patent classification data. We construct temporal networks of technologies that co-appear in patents. By analyzing the evolution of the core-periphery structure and centrality time-series in these networks, we identify periods of time dominated by technological combinations which we identify as distinct ‘eras’. We test the performance of our method using a database of patents in Great Britain spanning a century, and identify five distinct eras.

Protein interaction networks (PINs) allow the representation and analysis of biological processes in cells. Because cells are dynamic and adaptive, these processes change over time. Thus far, research has focused either on the static PIN analysis or the temporal nature of gene expression. By analysing temporal PINs using multilayer networks, we want to link these efforts. The analysis of temporal PINs gives insights into how proteins, individually and in their entirety, change their biological functions. We present a general procedure that integrates temporal gene expression information with a monolayer PIN to a temporal PIN and allows the detection of modular structure using multilayer modularity maximisation.

The structure of the state of art of scientific research is an important object of study motivated by the understanding of how research evolves and how new fields of study stem from existing research. In the last years complex networks tools contributed to provide insights on the structure of research, through the study of collaboration, citation and co-occurrence networks, in particular keyword co-occurrence networks proved useful to provide maps of knowledge inside a scientific domain. The network approach focuses on pairwise relationships, often compressing multidimensional data structures and inevitably losing information. In this paper we propose to adopt a simplicial complex approach to co-occurrence relations, providing a natural framework for the study of higher-order relations in the space of scientific knowledge. Using topological methods we explore the shape of concepts in mathematical research, focusing on homological cycles, regions with low connectivity in the simplicial structure, and we discuss their role in the understanding of the evolution of scientific research. In addition, we map authors’ contribution to the conceptual space, and explore their role in the formation of homological cycles.

Authors: Daniele Cassese, Vsevolod Salnikov, Renaud Lambiotte
 

We study how the synchronizability of a diffusive network increases (or decreases) when we add some links in its underlying graph. This is of interest in the context of power grids where people want to prevent from having blackouts, or for neural networks where synchronization is responsible of many diseases such as Parkinson. Based on spectral properties for Laplacian matrices, we show some classification results obtained (with Tiago Pereira and Philipp Pade) with respect to the effects of these links.
 

Assortative mixing in networks is the tendency for nodes with the same attributes, or metadata, to link to each other. It is a property often found in social networks manifesting as a higher tendency of links occurring between people with the same age, race, or political belief. Quantifying the level of assortativity or disassortativity (the preference of linking to nodes with different attributes) can shed light on the factors involved in the formation of links and contagion processes in complex networks. It is common practice to measure the level of assortativity according to the assortativity coefficient, or modularity in the case of discrete-valued metadata. This global value is the average level of assortativity across the network and may not be a representative statistic when mixing patterns are heterogeneous. For example, a social network spanning the globe may exhibit local differences in mixing patterns as a consequence of differences in cultural norms. Here, we introduce an approach to localise this global measure so that we can describe the assortativity, across multiple scales, at the node level. Consequently we are able to capture and qualitatively evaluate the distribution of mixing patterns in the network. We find that for many real-world networks the distribution of assortativity is skewed, overdispersed and multimodal. Our method provides a clearer lens through which we can more closely examine mixing patterns in networks.

Link to arxiv paper:  https://arxiv.org/abs/1708.01236

Online discussions are the essence of many social platforms on the Internet. Discussion platforms are receiving increasing interest because of their potential to become deliberative spaces. Although previous studies have proposed approaches to measure online deliberation using the complexity of discussion networks as a proxy, little research has focused on how these networks are affected by changes of platform features.

In this talk, we will focus on how interfaces might influence the network structures of discussions using techniques like interrupted time series analysis and regression discontinuity design. Futhermore, we will review and extend state-of-the-art generative models of discussion threads to explain better the structure and growth of online discussions.
 

Systemic risk arises as a multi-layer network phenomenon. Layers represent direct financial exposures of various types, including interbank liabilities, derivative or foreign exchange exposures. Another network layer of systemic risk emerges through common asset holdings of financial institutions. Strongly overlapping portfolios lead to similar exposures that are caused by price movements of the underlying financial assets. Based on the knowledge of portfolio holdings of financial agents we quantify systemic risk of overlapping portfolios. We present an optimization procedure, where we minimize the systemic risk in a given financial market by optimally rearranging overlapping portfolio networks, under the constraints that the expected returns and risks of the individual portfolios are unchanged. We explicitly demonstrate the power of the method on the overlapping portfolio network of sovereign exposure between major European banks by using data from the European Banking Authority stress test of 2016. We show that systemic-risk-efficient allocations are accessible by the optimization. In the case of sovereign exposure, systemic risk can be reduced by more than a factor of two, without any detrimental effects for the individual banks. These results are confirmed by a simple simulation of fire sales in the government bond market. In particular we show that the contagion probability is reduced dramatically in the optimized network.
 

Many studies of digital communication, in particular of Twitter, use natural language processing (NLP) to find topics, assess sentiment, and describe user behaviour.
In finding topics often the relationships between users who participate in the topic are neglected.
We propose a novel method of describing and classifying online conversations using only the structure of the underlying temporal network and not the content of individual messages.
This method utilises all available information in the temporal network (no aggregation), combining both topological and temporal structure using temporal motifs and inter-event times.
This allows us to describe the behaviour of individuals and collectives over time and examine the structure of conversation over multiple timescales.
 

Current technological progress has raised concerns about automation of tasks performed by workers resulting in job losses. Previous studies have used machine learning techniques to compute the automation probability of occupations and thus, studied the impact of automation on employment. However, such studies do not consider second-order effects, for example, an occupation with low automation probability can have a  surplus of labor supply due to similar occupations being automated. In this work, we study such second-order effects of automation using a network approach.  In our network – the Job Space – occupations are nodes and edges link occupations which share a significant amount of work activities. By mapping employment, automation probabilities into the network, and considering the movement of workers, we show that an occupation’s position in the network may be crucial to determining its employment future.

A great deal of effort has gone into trying to model social influence --- including the spread of behavior, norms, and ideas --- on networks. Most models of social influence tend to assume that individuals react to changes in the states of their neighbors without any time delay, but this is often not true in social contexts, where (for various reasons) different agents can have different response times. To examine such situations, we introduce the idea of a timer into threshold models of social influence. The presence of timers on nodes delays the adoption --- i.e., change of state --- of each agent, which in turn delays the adoptions of its neighbors. With a homogeneous-distributed timer, in which all nodes exhibit the same amount of delay, adoption delays are also homogeneous, so the adoption order of nodes remains the same. However, heterogeneously-distributed timers can change the adoption order of nodes and hence the "adoption paths" through which state changes spread in a network. Using a threshold model of social contagions, we illustrate that heterogeneous timers can either accelerate or decelerate the spread of adoptions compared to an analogous situation with homogeneous timers, and we investigate the relationship of such acceleration or deceleration with respect to timer distribution and network structure. We derive an analytical approximation for the temporal evolution of the fraction of adopters by modifying a pair approximation of the Watts threshold model, and we find good agreement with numerical computations. We also examine our new timer model on networks constructed from empirical data.

Link to arxiv paper: https://arxiv.org/abs/1706.04252

Temporal networks are increasingly being used to model the interactions of complex systems. 
Most studies require the temporal aggregation of edges (or events) into discrete time steps to perform analysis.
In this article we describe a static, behavioural representation of a temporal network, the temporal event graph (TEG).
The TEG describes the temporal network in terms of both inter-event time and two-event temporal motifs.
By considering the distributions of these quantities in unison we provide a new method to characterise the behaviour of individuals and collectives in temporal networks as well as providing a natural decomposition of the network.
We illustrate the utility of the TEG by providing examples on both synthetic and real temporal networks.

Identifying clusters or "communities" of densely connected nodes in networks is an active area of research, with relevance to many applications. Recent advances in the field have focused especially on temporal, multiplex, and other kinds of multilayer networks.

One method for detecting communities in multilayer networks is to maximise a generalised version of an objective function known as modularity. Writing down multilayer modularity requires the specification of two types of resolution parameters, and choosing appropriate values is crucial for uncovering meaningful community structure. In the simplest case, there are just two parameters, one controlling the sizes of detected communities, and the other influencing how much communities change from layer to layer. By establishing an equivalence between modularity optimisation and a multilayer maximum-likelihood approach to community detection, we are able to determine statistically optimal values for these two parameters. 

When applied to existing multilayer benchmarks, our optimized approach performs significantly better than using parameter choices guided by heuristics. We also apply the method to supermarket data, revealing changes in consumer behaviour over time.