Source author record

Michael Small

Michael Small appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2026arXiv

Dynamics, Complexity and Time Series Analysis

The aim of this text is to provide a linguistically accessible, but comprehensive introduction into a variety of topics in dynamical systems and its applications. Whilst preliminary knowledge of dynamical systems is useful, it is not essential and readers are only assumed to have familiarity with foundational undergraduate mathematics topics of calculus, linear algebra and rudimentary statistics. A variety of extended topics on recent publications and research activities in the field have been included in the last four chapters, which the interested reader may use as an introduction into further reading. A collection of exercises and questions both theoretical and computational are also included in this text.

preprint2020arXiv

Modelling remote epidemic transmission in Western Australia and implications for pandemic response

We develop an agent-based model of disease transmission in remote communities in Western Australia. Despite extreme isolation, we show that the movement of people amongst a large number of small but isolated communities has the effect of causing transmission to spread quickly. Significant movement between remote communities, and regional and urban centres allows for infection to quickly spread to and then among these remote communities. Our conclusions are based on two characteristic features of remote communities in Western Australia: (1) high mobility of people amongst these communities, and (2) relatively high proportion of travellers from very small communities to major population centres. In models of infection initiated in the state capital, Perth, these remote communities are collectively and uniquely vulnerable. Our model and analysis does not account for possibly heightened impact due to preexisting conditions, such additional assumptions would only make the projections of this model more dire. We advocate stringent monitoring and control of movement to prevent significant impact on the indigenous population of Western Australia.

preprint2020arXiv

Modelling strong control measures for epidemic propagation with networks -- A COVID-19 case study

We show that precise knowledge of epidemic transmission parameters is not required to build an informative model of the spread of disease. We propose a detailed model of the topology of the contact network under various external control regimes and demonstrate that this is sufficient to capture the salient dynamical characteristics and to inform decisions. Contact between individuals in the community is characterised by a contact graph, the structure of that contact graph is selected to mimic community control measures. Our model of city-level transmission of an infectious agent (SEIR model) characterises spread via a (a) scale-free contact network (no control); (b) a random graph (elimination of mass gatherings); and (c) small world lattice (partial to full lockdown -- "social" distancing). This model exhibits good qualitative agreement between simulation and data from the 2020 pandemic spread of coronavirus. Estimates of the relevant rate parameters of the SEIR model are obtained and we demonstrate the robustness of our model predictions under uncertainty of those estimates. The social context and utility of this work is identified, contributing to a highly effective pandemic response in Western Australia.

preprint2020arXiv

Navigating differential structures in complex networks

Structural changes in a network representation of a system (e.g.,different experimental conditions, time evolution), can provide insight on its organization, function and on how it responds to external perturbations. The deeper understanding of how gene networks cope with diseases and treatments is maybe the most incisive demonstration of the gains obtained through this differential network analysis point-of-view, which lead to an explosion of new numeric techniques in the last decade. However, {\it where} to focus ones attention, or how to navigate through the differential structures can be overwhelming even for few experimental conditions. In this paper, we propose a theory and a methodological implementation for the characterization of shared "structural roles" of nodes simultaneously within and between networks, whose outcome is a highly {\em interpretable} map. The main features and accuracy are investigated with numerical benchmarks generated by a stochastic block model. Results show that it can provide nuanced and interpretable information in scenarios with very different (i) community sizes and (ii) total number of communities, and (iii) even for a large number of 100 networks been compared (e.g., for 100 different experimental conditions). Then, we show evidence that the strength of the method is its "story-telling"-like characterization of the information encoded in a set of networks, which can be used to pinpoint unexpected differential structures, leading to further investigations and providing new insights. We provide an illustrative, exploratory analysis of four gene co-expression networks from two cell types $\times$ two treatments (interferon-$β$ stimulated or control). The method proposed here allowed us to elaborate and test a set of very specific hypotheses related to {\em unique} and {\em subtle} nuances of the structural differences between these networks.

preprint2020arXiv

West Australian Pandemic Response: The Black Swan of Black Swans

The COVID-19 Pandemic has been described as the global challenge of our time, an enormous human tragedy with dramatic economic impacts. This paper describes the response and expected recovery process for Western Australia, where a rapid and effective response was implemented. This has enabled an early transition into an expected recovery both in health and economic terms. The positive lessons learned from this experience are documented as they emerge in order to support other states and nations as they address this issue globally in the near-term and consider enduring improvements for the longer term. While the authors have personal experience in the WA context, wider observations across Australia and selected international benchmarks are also included. Key lessons include the importance of good health advice in Australia's interest; timely, synchronized and aligned action at all levels of government; a program of well communicated, aligned health and economic measures which support all in society allowing a very high level of appropriate community behaviour, ensuring the health system was not overloaded; innovation in telehealth, testing, pandemic modelling, and integrated operations which also allowed essential industries to continue; and strong border and travel controls with highly effective isolation preventing community spread, ultimately enabling rapid elimination of the disease from the hospital system. In combination, these demonstrate that in the case of Western Australia the result of first eliminating the disease from the community, and then reopening the economy progressively at a strong pace, has enabled a world leading outcome in both in health and economic terms. The lessons from this experience are widely applicable, shareable both as supporting service to other regions and through knowledge transfer.

preprint2019arXiv

Fast Automatic Detection of Geological Boundaries from Multivariate Log Data Using Recurrence

Manual interpretation of data collected from drill holes for mineral or oil and gas exploration is time-consuming and subjective. Identification of geological boundaries and distinctive rock physical property domains is the first step of interpretation. We introduce a multivariate technique, that can identify geological boundaries from petrophysical or geochemical data. The method is based on time-series techniques that have been adapted to be applicable for detecting transitions in geological spatial data. This method allows for the use of multiple variables in detecting different lithological layers. Additionally, it reconstructs the phase space of a single drill-hole or well to be applicable for further investigations across other holes or wells. The computationally cheap method shows efficiency and accuracy in detecting boundaries between lithological layers, which we demonstrate using examples from mineral exploration boreholes and an offshore gas exploration well.

preprint2016arXiv

Impact of asymptomatic infection on coupled disease-behavior dynamics in complex networks

Studies on how to model the interplay between diseases and behavioral responses (so-called coupled disease-behavior interaction) have attracted increasing attention. Owing to the lack of obvious clinical evidence of diseases, or the incomplete information related to the disease, the risks of infection cannot be perceived and may lead to inappropriate behavioral responses. Therefore, how to quantitatively analyze the impacts of asymptomatic infection on the interplay between diseases and behavioral responses is of particular importance. In this Letter, under the complex network framework, we study the coupled disease-behavior interaction model by dividing infectious individuals into two states: U-state (without evident clinical symptoms, labelled as U) and I-state (with evident clinical symptoms, labelled as I). A susceptible individual can be infected by U- or I-nodes, however, since the U-nodes cannot be easily observed, susceptible individuals take behavioral responses \emph{only} when they contact I-nodes. The mechanism is considered in the improved Susceptible-Infected-Susceptible (SIS) model and the improved Susceptible-Infected-Recovered (SIR) model, respectively. Then, one of the most concerned problems in spreading dynamics: the epidemic thresholds for the two models are given by two methods. The analytic results \emph{quantitatively} describe the influence of different factors, such as asymptomatic infection, the awareness rate, the network structure, and so forth, on the epidemic thresholds. Moreover, because of the irreversible process of the SIR model, the suppression effect of the improved SIR model is weaker than the improved SIS model.

preprint2016arXiv

Navigation by anomalous random walks on complex networks

Anomalous random walks having long-range jumps are a critical branch of dynamical processes on networks, which can model a number of search and transport processes. However, traditional measurements based on mean first passage time are not useful as they fail to characterize the cost associated with each jump. Here we introduce a new concept of mean first traverse distance (MFTD) to characterize anomalous random walks that represents the expected traverse distance taken by walkers searching from source node to target node, and we provide a procedure for calculating the MFTD between two nodes. We use Levy walks on networks as an example, and demonstrate that the proposed approach can unravel the interplay between diffusion dynamics of Levy walks and the underlying network structure. Interestingly, applying our framework to the famous PageRank search, we can explain why its damping factor empirically chosen to be around 0.85. The framework for analyzing anomalous random walks on complex networks offers a new useful paradigm to understand the dynamics of anomalous diffusion processes, and provides a unified scheme to characterize search and transport processes on networks.

preprint2016arXiv

Prevention of infectious diseases by public vaccination and individual protection

In the face of serious infectious diseases, governments endeavour to implement containment measures such as public vaccination at a macroscopic level. Meanwhile, individuals tend to protect themselves by avoiding contacts with infections at a microscopic level. However, a comprehensive understanding of how such combined strategy influences epidemic dynamics is still lacking. We study a susceptible-infected-susceptible epidemic model with imperfect vaccination on dynamic contact networks, where the macroscopic intervention is represented by random vaccination of the population and the microscopic protection is characterised by susceptible individuals rewiring contacts from infective neighbours. In particular, the model is formulated both in populations without and then with demographic effects. Using the pairwise approximation and the probability generating function approach, we investigate both dynamics of the epidemic and the underlying network. For populations without demography, the emerging degree correlations, bistable states, and oscillations demonstrate the combined effects of the public vaccination program and individual protective behavior. Compared to either strategy in isolation, the combination of public vaccination and individual protection is more effective in preventing and controlling the spread of infectious diseases by increasing both the invasion threshold and the persistence threshold. For populations with additional demographic factors, the integration between vaccination intervention and individual rewiring may promote epidemic spreading due to the birth effect. Moreover, the degree distributions of both networks in the steady state is closely related to the degree distribution of newborns, which leads to uncorrelated connectivity. All the results demonstrate the importance of both local protection and global intervention, as well as the demographic effects.

preprint2015arXiv

Preferential imitation of vaccinating behavior can invalidate the targeted subsidy on complex network

We consider the effect of inducement to vaccinate during the spread of an infectious disease on complex networks. Suppose that public resources are finite and that only a small proportion of individuals can be vaccinated freely (complete subsidy), for the remainder of the population vaccination is a voluntary behavior --- and each vaccinated individual carries a perceived cost. We ask whether the classical targeted subsidy strategy is definitely better than the random strategy: does targeting subsidy at individuals perceived to be with the greatest risk actually help? With these questions, we propose a model to investigate the \emph{interaction effects} of the subsidy policies and individuals responses when facing subsidy policies on the epidemic dynamics on complex networks. In the model, a small proportion of individuals are freely vaccinated according to either the targeted or random subsidy policy, the remainder choose to vaccinate (or not) based on voluntary principle and update their vaccination decision via an imitation rule. Our findings show that the targeted strategy is only advantageous when individuals prefer to imitate the subsidized individuals' strategy. Otherwise, the effect of the targeted policy is worse than the random immunization, since individuals preferentially select non-subsidized individuals as the imitation objects. More importantly, we find that under the targeted subsidy policy, increasing the proportion of subsidized individuals may increase the final epidemic size. We further define social cost as the sum of the costs of vaccination and infection, and study how each of the two policies affect the social cost. Our result shows that there exist some optimal intermediate regions leading to the minimal social cost.

preprint2015arXiv

Time lagged ordinal partition networks for capturing dynamics of continuous dynamical systems

We investigate a generalised version of the recently proposed ordinal partition time series to network transformation algorithm. Firstly we introduce a fixed time lag for the elements of each partition that is selected using techniques from traditional time delay embedding. The resulting partitions define regions in the embedding phase space that are mapped to nodes in the network space. Edges are allocated between nodes based on temporal succession thus creating a Markov chain representation of the time series. We then apply this new transformation algorithm to time series generated by the Rössler system and find that periodic dynamics translate to ring structures whereas chaotic time series translate to band or tube-like structures -- thereby indicating that our algorithm generates networks whose structure is sensitive to system dynamics. Furthermore we demonstrate that simple network measures including the mean out degree and variance of out degrees can track changes in the dynamical behaviour in a manner comparable to the largest Lyapunov exponent. We also apply the same analysis to experimental time series generated by a diode resonator circuit and show that the network size, mean shortest path length and network diameter are highly sensitive to the interior crisis captured in this particular data set.

preprint2014arXiv

Basin of Attraction Determines Hysteresis in Explosive Synchronization

Spontaneous explosive emergent behavior takes place in heterogeneous networks when the frequencies of the nodes are positively correlated to the node degree. A central feature of such explosive transitions is a hysteretic behavior at the transition to synchronization. We unravel the underlying mechanisms and show that the dynamical origin of the hysteresis is a change of basin of attraction of the synchronization state. Our findings hold for heterogeneous networks with star graph motifs such as scale free networks, and hence reveal how microscopic network parameters such as node degree and frequency affect the global network properties and can be used for network design and control.

preprint2014arXiv

Exactly scale-free scale-free networks

Many complex natural and physical systems exhibit patterns of interconnection that conform, approximately, to a network structure referred to as scale-free. Preferential attachment is one of many algorithms that have been introduced to model the growth and structure of scale-free networks. With so many different models of scale-free networks it is unclear what properties of scale-free networks are typical, and what properties are peculiarities of a particular growth or construction process. We propose a simple maximum entropy process which provides the best representation of what are typical properties of scale-free networks, and provides a standard against which real and algorithmically generated networks can be compared. As an example we consider preferential attachment and find that this particular growth model does not yield typical realizations of scale-free networks. In particular, the widely discussed "fragility" of scale-free networks is actually found to be due to the peculiar "hub-centric" structure of preferential attachment networks. We provide a method to generate or remove this latent hub-centric bias --- thereby demonstrating exactly which features of preferential attachment networks are atypical of the broader class of scale-free networks. We are also able to statistically demonstrate whether real networks are typical realizations of scale-free networks, or networks with that particular degree distribution; using a new surrogate generation method for complex networks, exactly analogous the the widely used surrogate tests of nonlinear time series analysis.

preprint2014arXiv

Expected degree of finite preferential attachment networks

We provide an analytic expression for the quantity described in the title. Namely, we perform a preferential attachment growth process to generate a scale-free network. At each stage we add a new node with $m$ new links. Let $k$ denote the degree of a node, and $N$ the number of nodes in the network. The degree distribution is assumed to converge to a power-law (for $k\geq m$) of the form $k^{-γ}$ and we obtain an exact implicit relationship for $γ$, $m$ and $N$. We verify this with numerical calculations over several orders of magnitude. Although this expression is exact, it provides only an implicit expression for $γ(m)$. Nonetheless, we provide a reasonable guess as to the form of this curve and perform curve fitting to estimate the parameters of that curve --- demonstrating excellent agreement between numerical fit, theory, and simulation.

preprint2014arXiv

Exploration of period-doubling cascade route to chaos with complex network based time series construction

In this work, the topologies of networks constructed from time series from an underlying system undergo a period doubling cascade have been explored by means of the prevalence of different motifs using an efficient computational motif detection algorithm. By doing this we adopt a refinement based on the $k$ nearest neighbor recurrence-based network has been proposed. We demonstrate that the refinement of network construction together with the study of prevalence of different motifs allows a full explosion of the evolving period doubling cascade route to chaos in both discrete and continuous dynamical systems. Further, this links the phase space time series topologies to the corresponding network topologies, and thus helps to understand the empirical "superfamily" phenomenon, as shown by Xu.

preprint2014arXiv

Super-star networks: Growing optimal scale-free networks via likelihood

Preferential attachment --- by which new nodes attach to existing nodes with probability proportional to the existing nodes' degree --- has become the standard growth model for scale-free networks, where the asymptotic probability of a node having degree $k$ is proportional to $k^{-γ}$. However, the motivation for this model is entirely ad hoc. We use exact likelihood arguments and show that the optimal way to build a scale-free network is to attach most new links to nodes of low degree. Curiously, this leads to a scale-free networks with a single dominant hub: a star-like structure we call a super-star network. Asymptotically, the optimal strategy is to attach each new node to one of the nodes of degree $k$ with probability proportional to $\frac{1}{N+ζ(γ)(k+1)^γ}$ (in a $N$ node network) --- a stronger bias toward high degree nodes than exhibited by standard preferential attachment. Our algorithm generates optimally scale-free networks (the super-star networks) as well as randomly sampling the space of all scale-free networks with a given degree exponent $γ$. We generate viable realisation with finite $N$ for $1\ll γ<2$ as well as $γ>2$. We observe an apparently discontinuous transition at $γ\approx 2$ between so-called super-star networks and more tree-like realisations. Gradually increasing $γ$ further leads to re-emergence of a super-star hub. To quantify these structural features we derive a new analytic expression for the expected degree exponent of a pure preferential attachment process, and introduce alternative measures of network entropy. Our approach is generic and may also be applied to an arbitrary degree distribution.

preprint2013arXiv

A surrogate for networks -- How scale-free is my scale-free network?

Complex networks are now being studied in a wide range of disciplines across science and technology. In this paper we propose a method by which one can probe the properties of experimentally obtained network data. Rather than just measuring properties of a network inferred from data, we aim to ask how typical is that network? What properties of the observed network are typical of all such scale free networks, and which are peculiar? To do this we propose a series of methods that can be used to generate statistically likely complex networks which are both similar to the observed data and also consistent with an underlying null-hypothesis -- for example a particular degree distribution. There is a direct analogy between the approach we propose here and the surrogate data methods applied to nonlinear time series data.

preprint2013arXiv

Complex Network Approach to the Statistical Features of the Sunspot Series

Complex network approaches have been recently developed as an alternative framework to study the statistical features of time-series data. We perform a visibility-graph analysis on both the daily and monthly sunspot series. Based on the data, we propose two ways to construct the network: one is from the original observable measurements and the other is from a negative-inverse-transformed series. The degree distribution of the derived networks for the strong maxima has clear non-Gaussian properties, while the degree distribution for minima is bimodal. The long-term variation of the cycles is reflected by hubs in the network which span relatively large time intervals. Based on standard network structural measures, we propose to characterize the long-term correlations by waiting times between two subsequent events. The persistence range of the solar cycles has been identified over 15\,--\,1000 days by a power-law regime with scaling exponent $γ= 2.04$ of the occurrence time of the two subsequent and successive strong minima. In contrast, a persistent trend is not present in the maximal numbers, although maxima do have significant deviations from an exponential form. Our results suggest some new insights for evaluating existing models. The power-law regime suggested by the waiting times does indicate that there are some level of predictable patterns in the minima.

preprint2013arXiv

Temporal prediction of epidemic patterns in community networks

Most previous studies of epidemic dynamics on complex networks suppose that the disease will eventually stabilize at either a disease-free state or an endemic one. In reality, however, some epidemics always exhibit sporadic and recurrent behaviour in one region because of the invasion from an endemic population elsewhere. In this paper we address this issue and study a susceptible-infected-susceptible epidemiological model on a network consisting of two communities, where the disease is endemic in one community but alternates between outbreaks and extinctions in the other. We provide a detailed characterization of the temporal dynamics of epidemic patterns in the latter community. In particular, we investigate the time duration of both outbreak and extinction, and the time interval between two consecutive inter-community infections, as well as their frequency distributions. Based on the mean-field theory, we theoretically analyze these three timescales and their dependence on the average node degree of each community, the transmission parameters, and the number of intercommunity links, which are in good agreement with simulations, except when the probability of overlaps between successive outbreaks is too large. These findings aid us in better understanding the bursty nature of disease spreading in a local community, and thereby suggesting effective time-dependent control strategies.

preprint2013arXiv

What exactly are the properties of scale-free and other networks?

The concept of scale-free networks has been widely applied across natural and physical sciences. Many claims are made about the properties of these networks, even though the concept of scale-free is often vaguely defined. We present tools and procedures to analyse the statistical properties of networks defined by arbitrary degree distributions and other constraints. Doing so reveals the highly likely properties, and some unrecognised richness, of scale-free networks, and casts doubt on some previously claimed properties being due to a scale-free characteristic.

preprint2012arXiv

Generating self-organizing collective behavior using separation dynamics from experimental data

Mathematical models for systems of interacting agents using simple local rules have been proposed and shown to exhibit emergent swarming behavior. Most of these models are constructed by intuition or manual observations of real phenomena, and later tuned or verified to simulate desired dynamics. In contrast to this approach, we propose using a model that attempts to follow an averaged rule of the essential distance-dependent collective behavior of real pigeon flocks, which was abstracted from experimental data. By using a simple model to follow the behavioral tendencies of real data, we show that our model can exhibit emergent self-organizing dynamics such as flocking, pattern formation, and counter-rotating vortices. The range of behaviors observed in our simulations are richer than the standard models of collective dynamics, and should thereby give potential for new models of complex behavior.

preprint2012arXiv

Multiscale characterization of recurrence-based phase space networks constructed from time series

Recently, a framework for analyzing time series by constructing an associated complex network has attracted significant research interest. One of the advantages of the complex network method for studying time series is that complex network theory provides a tool to describe either important nodes, or structures that exist in the networks, at different topological scale. This can then provide distinct information for time series of different dynamical systems. In this paper, we systematically investigate the recurrence-based phase space network of order $k$ that has previously been used to specify different types of dynamics in terms of the motif ranking from a different perspective. Globally, we find that the network size scales with different scale exponents and the degree distribution follows a quasi-symmetric bell shape around the value of $2k$ with different values of degree variance from periodic to chaotic Rössler systems. Local network properties such as the vertex degree, the clustering coefficients and betweenness centrality are found to be sensitive to the local stability of the orbits and hence contain complementary information.

preprint2012arXiv

Neuronal avalanches of a self-organized neural network with active-neuron-dominant structure

Neuronal avalanche is a spontaneous neuronal activity which obeys a power-law distribution of population event sizes with an exponent of -3/2. It has been observed in the superficial layers of cortex both \emph{in vivo} and \emph{in vitro}. In this paper we analyze the information transmission of a novel self-organized neural network with active-neuron-dominant structure. Neuronal avalanches can be observed in this network with appropriate input intensity. We find that the process of network learning via spike-timing dependent plasticity dramatically increases the complexity of network structure, which is finally self-organized to be active-neuron-dominant connectivity. Both the entropy of activity patterns and the complexity of their resulting post-synaptic inputs are maximized when the network dynamics are propagated as neuronal avalanches. This emergent topology is beneficial for information transmission with high efficiency and also could be responsible for the large information capacity of this network compared with alternative archetypal networks with different neural connectivity.

preprint2012arXiv

Pairwise interaction pattern in the weighted communication network

Although recent studies show that both topological structures and human dynamics can strongly affect information spreading on social networks, the complicated interplay of the two significant factors has not yet been clearly described. In this work, we find a strong pairwise interaction based on analyzing the weighted network generated by the short message communication dataset within a Chinese tele-communication provider. The pairwise interaction bridges the network topological structure and human interaction dynamics, which can promote local information spreading between pairs of communication partners and in contrast can also suppress global information (e.g., rumor) cascade and spreading. In addition, the pairwise interaction is the basic pattern of group conversations and it can greatly reduce the waiting time of communication events between a pair of intimate friends. Our findings are also helpful for communication operators to design novel tariff strategies and optimize their communication services.

preprint2012arXiv

Predicting the outcome of roulette

There have been several popular reports of various groups exploiting the deterministic nature of the game of roulette for profit. Moreover, through its history the inherent determinism in the game of roulette has attracted the attention of many luminaries of chaos theory. In this paper we provide a short review of that history and then set out to determine to what extent that determinism can really be exploited for profit. To do this, we provide a very simple model for the motion of a roulette wheel and ball and demonstrate that knowledge of initial position, velocity and acceleration is sufficient to predict the outcome with adequate certainty to achieve a positive expected return. We describe two physically realisable systems to obtain this knowledge both incognito and {\em in situ}. The first system relies only on a mechanical count of rotation of the ball and the wheel to measure the relevant parameters. By applying this techniques to a standard casino-grade European roulette wheel we demonstrate an expected return of at least 18%, well above the -2.7% expected of a random bet. With a more sophisticated, albeit more intrusive, system (mounting a digital camera above the wheel) we demonstrate a range of systematic and statistically significant biases which can be exploited to provide an improved guess of the outcome. Finally, our analysis demonstrates that even a very slight slant in the roulette table leads to a very pronounced bias which could be further exploited to substantially enhance returns.

preprint2012arXiv

The Impacts of Subsidy Policies on Vaccination Decisions in Contact Networks

Often, vaccination programs are carried out based on self-interest rather than being mandatory. Owing to the perceptions about risks associated with vaccines and the `herd immunity' effect, it may provide suboptimal vaccination coverage for the population as a whole. In this case, some subsidy policies may be offered by the government to promote vaccination coverage. But, not all subsidy policies are effective in controlling the transmission of infectious diseases. We address the question of which subsidy policy is best, and how to appropriately distribute the limited subsidies to maximize vaccine coverage. To answer these questions, we establish a model based on evolutionary game theory, where individuals try to maximize their personal payoffs when considering the voluntary vaccination mechanism. Our model shows that voluntary vaccination alone is insufficient to control an epidemic. Hence, two subsidy policies are systematically studied: (1) in the free subsidy policy the total amount of subsidies is distributed to some individuals and all the donees may vaccinate at no cost, and (2) in the part-offset subsidy policy each vaccinated person is offset by a certain proportion of the vaccination cost. Simulations suggest that, since the part-offset subsidy policy can encourage more individuals to be vaccinated, the performance of this policy is significantly better than that of the free subsidy policy.

preprint2012arXiv

The instability of followers and emergent vorticity in flocking behaviour for an experimental interaction rule

Computational models of collective behavior in birds has allowed us to infer interaction rules directly from experimental data. Using a generic form of these rules we explore the collective behavior and emergent dynamics of a simulated swarm. For a wide range of flock size and interaction extent (the fixed number of neighbors with which an individual will interact) we find that the computational collective is inherently stable --- individuals are attracted to one another and will position themselves a preferred distance from their fixed neighbors within a rigid lattice. Nonetheless, the irregular overall shape of the flock, coupled with the need for individuals on the boundary to move towards their neighbors creates a torque which leads the flock to rotate and then meander. We argue that this "rolling meander" is a very good proxy for real collective behavior in animal species and yet arises from a simple homogeneous and deterministic rule for interaction. Rather than then introduce leaders --- which has already been shown, quite straightforwardly, to drive collective swarms such as this --- we introduce a small number of "followers". Each follower is bound to consider a random fixed individual to be among their neighbors, irrespective of actual metric distance between them. We find that the introduction of a small number of such followers causes a phase transition that quickly leads to instability in the flock structure (as no stable configuration arises) and the previously rigid crystalline interaction among neighbors now becomes fluid: the distance between neighbors decreases, the flock ceases to rotate and meanders less.

preprint2012arXiv

The stability of networks --- towards a structural dynamical systems theory

The need to build a link between the structure of a complex network and the dynamical properties of the corresponding complex system (comprised of multiple low dimensional systems) has recently become apparent. Several attempts to tackle this problem have been made and all focus on either the controllability or synchronisability of the network --- usually analyzed by way of the master stability function, or the graph Laplacian. We take a different approach. Using the basic tools from dynamical systems theory we show that the dynamical stability of a network can easily be defined in terms of the eigenvalues of an homologue of the network adjacency matrix. This allows us to compute the stability of a network (a quantity derived from the eigenspectrum of the adjacency matrix). Numerical experiments show that this quantity is very closely related too, and can even be predicted from, the standard structural network properties. Following from this we show that the stability of large network systems can be understood via an analytic study of the eigenvalues of their fixed points --- even for a very large number of fixed points.

preprint2011arXiv

Dynamical modeling of collective behavior from pigeon flight data: flock cohesion and dispersion

Several models of flocking have been promoted based on simulations with qualitatively naturalistic behavior. In this paper we provide the first direct application of computational modeling methods to infer flocking behavior from experimental field data. We show that this approach is able to infer general rules for interaction, or lack of interaction, among members of a flock or, more generally, any community. Using experimental field measurements of homing pigeons in flight we demonstrate the existence of a basic distance dependent attraction/repulsion relationship and show that this rule is sufficient to explain collective behavior observed in nature. Positional data of individuals over time are used as input data to a computational algorithm capable of building complex nonlinear functions that can represent the system behavior. Topological nearest neighbor interactions are considered to characterize the components within this model. The efficacy of this method is demonstrated with simulated noisy data generated from the classical (two dimensional) Vicsek model. When applied to experimental data from homing pigeon flights we show that the more complex three dimensional models are capable of predicting and simulating trajectories, as well as exhibiting realistic collective dynamics. The simulations of the reconstructed models are used to extract properties of the collective behavior in pigeons, and how it is affected by changing the initial conditions of the system. Our results demonstrate that this approach may be applied to construct models capable of simulating trajectories and collective dynamics using experimental field measurements of herd movement. From these models, the behavior of the individual agents (animals) may be inferred.

preprint2011arXiv

Optimizing and controlling functions of complex networks by manipulating rich-club connections

Traditionally, there is no evidence suggesting that there are strong ties between the rich-club property and the function of complex networks. In this study, we find that whether a very small portion of rich nodes connected to each other or not can strongly affect the frequency of occurrence of basic building blocks (motif) within networks, and therefore the function, of a heterogeneous network. Conversely whether a homogeneous network has a rich-club property or not generally has no significant effect on its structure and function. These findings open the possibility to optimize and control the function of complex networks by manipulating rich-club connections. Furthermore, based on the subgraph ratio profile, we develop a more rigorous approach to judge whether a network has a rich-club or not. The new method does not calculate how many links there are among rich nodes but depends on how the links among rich nodes can affect the overall structure as well as function of a given network. These results can also help us to understand the evolution of dynamical networks and design new models for characterizing real-world networks.

preprint2011arXiv

Reciprocal relationships in collective flights of homing pigeons

Collective motion of bird flocks can be explained via the hypothesis of many wrongs, and/or, a structured leadership mechanism. In pigeons, previous studies have shown that there is a well-defined hierarchical structure and certain specific individuals occupy more dominant positions --- suggesting that leadership by the few individuals drives the behavior of the collective. Conversely, by analyzing the same data-sets, we uncover a more egalitarian mechanism. We show that both reciprocal relationships and a stratified hierarchical leadership are important and necessary in the collective movements of pigeon flocks. Rather than birds adopting either exclusive averaging or leadership strategies, our experimental results show that it is an integrated combination of both compromise and leadership which drives the group's movement decisions.

preprint2010arXiv

Fitness-driven deactivation in network evolution

Individual nodes in evolving real-world networks typically experience growth and decay --- that is, the popularity and influence of individuals peaks and then fades. In this paper, we study this phenomenon via an intrinsic nodal fitness function and an intuitive aging mechanism. Each node of the network is endowed with a fitness which represents its activity. All the nodes have two discrete stages: active and inactive. The evolution of the network combines the addition of new active nodes randomly connected to existing active ones and the deactivation of old active nodes with possibility inversely proportional to their fitnesses. We obtain a structured exponential network when the fitness distribution of the individuals is homogeneous and a structured scale-free network with heterogeneous fitness distributions. Furthermore, we recover two universal scaling laws of the clustering coefficient for both cases, $C(k) \sim k^{-1}$ and $C \sim n^{-1}$, where $k$ and $n$ refer to the node degree and the number of active individuals, respectively. These results offer a new simple description of the growth and aging of networks where intrinsic features of individual nodes drive their popularity, and hence degree.

preprint2010arXiv

Recurrence-based time series analysis by means of complex network methods

Complex networks are an important paradigm of modern complex systems sciences which allows quantitatively assessing the structural properties of systems composed of different interacting entities. During the last years, intensive efforts have been spent on applying network-based concepts also for the analysis of dynamically relevant higher-order statistical properties of time series. Notably, many corresponding approaches are closely related with the concept of recurrence in phase space. In this paper, we review recent methodological advances in time series analysis based on complex networks, with a special emphasis on methods founded on recurrence plots. The potentials and limitations of the individual methods are discussed and illustrated for paradigmatic examples of dynamical systems as well as for real-world time series. Complex network measures are shown to provide information about structural features of dynamical systems that are complementary to those characterized by other methods of time series analysis and, hence, substantially enrich the knowledge gathered from other existing (linear as well as nonlinear) approaches.

preprint2010arXiv

Rich-club connectivity dominates assortativity and transitivity of complex networks

Rich-club, assortativity and clustering coefficients are frequently-used measures to estimate topological properties of complex networks. Here we find that the connectivity among a very small portion of the richest nodes can dominate the assortativity and clustering coefficients of a large network, which reveals that the rich-club connectivity is leveraged throughout the network. Our study suggests that more attention should be payed to the organization pattern of rich nodes, for the structure of a complex system as a whole is determined by the associations between the most influential individuals. Moreover, by manipulating the connectivity pattern in a very small rich-club, it is sufficient to produce a network with desired assortativity or transitivity. Conversely, our findings offer a simple explanation for the observed assortativity and transitivity in many real world networks --- such biases can be explained by the connectivities among the richest nodes.

preprint2009arXiv

Mapping from Architecture to Dynamics: A Unified View of Dynamical Processes on Networks

Although it is unambiguously agreed that structure plays a fundamental role in shaping the dynamics of complex systems, this intricate relationship still remains unclear. We investigate a general computational transformation by which we can map the network topology directly to the dynamical patterns emergent on it -- independent of the nature of the dynamical process. We find that many seemingly diverse dynamical processes such as coupled oscillators and diffusion phenomena can all be understood and unified through this same procedure. Using the multiscale complexity measure derived form the structure-dynamics transformation, we find that the topological features like hierarchy, heterogeneity and modularity all result in higher complexity. This result suggests a universal principle: it is the desire for functional diversity that drives the evolution of network architecture.

preprint2009arXiv

Revising the simple measures of assortativity in complex networks

We find that traditional statistics for measuring degree mixing are strongly affected by superrich nodes. To counteract and measure the effect of superrich nodes, we propose a paradigm to quantify the mixing pattern of a real network in which different mixing patterns may appear among low-degree nodes and among high-degree nodes. The new paradigm and the simple revised measure uncover the true complex degree mixing patterns of complex networks with superrich nodes. The new method indicates that some networks show a false disassortative mixing induced by superrich nodes, and have no tendency to be genuinely disassortative. Our results also show that the previously observed fragility of scale-free networks is actually greatly exacerbated by the presence of even a very small number of superrich nodes.