Controlling Epidemics Through Optimal Allocation of Test Kits and Vaccine Doses Across Networks Mingtao Xia, Lucas Bottcher,¨ Tom Chou

Total Page:16

File Type:pdf, Size:1020Kb

Controlling Epidemics Through Optimal Allocation of Test Kits and Vaccine Doses Across Networks Mingtao Xia, Lucas Bottcher,¨ Tom Chou 1 Controlling epidemics through optimal allocation of test kits and vaccine doses across networks Mingtao Xia, Lucas Bottcher,¨ Tom Chou Abstract—Efficient testing and vaccination protocols are crit- optimal vaccination strategies have been derived for a rapidly ical aspects of epidemic management. To study the optimal spreading disease in a highly mobile urban population using allocation of limited testing and vaccination resources in a PMP [12]. Complementing these control-theory-based inter- heterogeneous contact network of interacting susceptible, re- covered, and infected individuals, we present a degree-based ventions, a recent work [13] developed methods relying on testing and vaccination model for which we use control-theoretic reinforcement learning (RL) to identify infectious high-degree methods to derive optimal testing and vaccination policies using nodes (“superspreaders”) in temporal networks and reduce control-theoretic methods. Within our framework, we find that the overall infection rate with limited medical resources. optimal intervention policies first target high-degree nodes before The application of optimal control methods and PMP to shifting to lower-degree nodes in a time-dependent manner. Using such optimal policies, it is possible to delay outbreaks a heterogeneous node-based susceptible-infected-recovered- and reduce incidence rates to a greater extent than uniform susceptible (SIRS) model with applications to rumor spreading and reinforcement-learning-based interventions, particularly on was studied in [14]. certain scale-free networks. The machine-learning-based interventions of [13] showed Index Terms—Infection networks, epidemics, testing, vaccina- that RL is able to outperform intervention policies derived tion, optimal control, reinforcement learning from purely structural node characterizations that are, for instance, based on centrality measures. However, the methods I. INTRODUCTION of [13] were applied to rather small networks with a maximum Limiting the spread of novel pathogens such as SARS-CoV- number of nodes of about 400. Here, we focus on a com- 2 requires efficient testing [1], [2] and quarantine strategies [3], plementary approach by formulating optimal control and RL- especially when vaccines are not available or effective. Even if based target policies for a degree-based epidemic model [15] effective vaccines become available at scale, their population- that is constrained only by the maximum degree and not by the wide distribution is a complex and time-consuming endeavor, system size (i.e., number of nodes). Early work by May and influenced by, for example, age-structure [4]–[6], vaccine Anderson [16] employed such effective degree models to study hesitancy [7], and different objectives [8]. the population-level dynamics of human immunodeficiency Until a sufficient level of immunity within a population is virus (HIV) infections. These degree-based models and later reached, distancing and quarantine policies can also be used adaptations [17]–[19] do not account for degree correlations. to help slow the spread and evolutionary dynamics [9] of Effective degree models for susceptible-infected-susceptible infectious diseases. Epidemic modeling and control-theoretic (SIS) dynamics with degree correlations were derived in [20] approaches are useful for identifying both efficient testing and and applied to SIR dynamics in [21]. A further generalization vaccination policies. For an epidemic model of SARS-CoV- of these methods to model SIR dynamics with networked and 2 transmission, Pontryagin’s maximum principle (PMP) has well-mixed transmission pathways was presented in [22]. For a been used to derive optimal distancing and testing strategies detailed summary of degree-based epidemic models, see [23]. that minimize the number of COVID-19 cases and intervention In the next section, we propose and justify a degree-based costs [10]. Optimal control theory has also been applied epidemic, testing, and quarantining model. An optimal control to a multi-objective control problem that uses isolation and framework for this model is presented in Sec. III and, given vaccination to limit epidemic size and duration [11]. Both limited testing resources, an optimal testing strategy is calcu- of these recent investigations describe the underlying infec- lated. We extend the same underlying disease model to include tious disease dynamics through compartmental models without vaccination in Sec. IV and find optimal vaccination strategies underlying network structure, meaning that all interactions that minimize infection given a limited vaccination rate. We among different individuals are assumed to be homogeneous. summarize and discuss our results and how they depend on For a structured susceptible-infected-recovered (SIR) model, network and dynamical features of the model in Sec. V. For comparison, we also present in the Appendix a reinforcement- Mingtao Xia is in the Dept. of Mathematics at UCLA. learning-based algorithm that is able to approximate optimal E-mail: [email protected] testing strategies for the model introduced in Sec. II. Lucas Bottcher¨ is in the Dept. of Computational Medicine at UCLA and at the Frankfurt School of Finance and Management. E-mail: [email protected] Tom Chou is in the Depts. of Computational Medicine and Mathematics at II. DEGREE-BASED EPIDEMIC AND TESTING MODEL UCLA E-mail: [email protected] For the formulation of optimal testing policies that allocate Manuscript received July 18, 2021; revised August 31, 2021. testing resources to different individuals in a contact network, 2 we adopt an effective degree model of SIR dynamics with testing in a static network of N nodes. Nodes represent individuals, and edges between nodes represent corresponding contacts. Therefore, the degree of a node represents the number of its contacts. If K is the maximum degree across all nodes, we can divide the population into K distinct subpopu- lations, each of size Nk (k = 1; 2;:::;K) such that all nodes th PK in the k group have degree k. Therefore, N = k=1 Nk. In our epidemic model, we distinguish between untested u ∗ and tested infected individuals. Let Sk(t), Ik (t), Ik (t), and Rk(t) denote the numbers of susceptible, untested infected, tested infected, and recovered nodes with degree k at time t, respectively. Since these subpopulations together represent the entire population (the total number of nodes N), both N and Nk are constants in our model. Their values satisfy u ∗ the normalization condition Sk + Ik + Ik + Rk = Nk. The corresponding fractions are u u sk(t) = Sk(t)=N; ik(t) = Ik (t)=N; ∗ ∗ (1) ik(t) = Ik (t)=N; rk(t) = Rk(t)=N; P u ∗ Fig. 1. Degree distribution of a Barabasi–Albert´ network and a stochastic such that k(sk + ik + ik + rk) = 1. Using an effective- block model. (a) The degree distribution of a Barabasi–Albert´ network with degree approach [16], [22], we describe the evolution of the 99,939 nodes. Each new node is connected to m = 2 existing nodes (i.e., the above subpopulations by degree of each node is at least 2) using preferential attachment. Then nodes with degrees larger than 100 [26] are removed from the network. The grey K solid line is a guide-to-the-eye with slope -3 [27]. The inset shows a realization dsk(t) X P (`jk) = − ks (t) βuiu(t) + β∗i∗(t); (2) of a Barabasi–Albert´ network with 100 nodes. Node size scales with their dt k P (`) ` ` `=1 betweenness centrality. (b) The conditional probability P (`jk) associated with the Barabasi–Albert´ network generated in (a). (c) The degree distribution of a K diu(t) X P (`jk) stochastic block model with four blocks and 100,000 nodes. The inset shows k =ks (t) βuiu(t) + β∗i∗(t) (3) dt k P (`) ` ` a realization of a stochastic block model with 800 nodes, but using the same `=1 block probability matrix. (d) The conditional probability P (`jk) associated u u fk(t) u with the SBM. In both (b) and (d), all elements that are strictly zero are − γ ik(t) − ik(t); uncolored. Nk ∗ dik(t) ∗ ∗ fk(t) u = − γ ik(t) + ik(t); (4) dt Nk lead to quarantining and do not affect the disease dynamics. dr (t) u u ∗ u k u u ∗ ∗ However a fraction I =(Sk + I + I + Rk) ≡ I =Nk of =γ ik(t) + γ ik(t); (5) k k k k dt these fk(t)∆t tests will be administered to untested infecteds. where P (`) = N`=N is the degree distribution and P (`jk) is Once infected nodes have been identified by testing, they can the probability that a chosen node with degree k is connected be quarantined and removed from the disease transmission to a node with degree `. Our degree-based formulation of SIR dynamics. If infected individuals who already have been tested dynamics with testing, Eqs. (2)–(5), is an approximation of strictly avoid future testing, more tests will be available for the full node-based dynamics assuming that nodes of the same the other subpopulations, increasing the rate at which the degree are equally likely to be infected at any given time [15]. remaining untested infecteds will be tested. In this case, the Susceptible individuals become infected through contact fraction of tests administered to untested infecteds is modified: u u u u ∗ with untested and tested infected individuals at rates β Ik =(Sk +Ik +Rk) ≡ Ik =(Nk −Ik ). After normalizing by the ∗ u and β , respectively. Untested and tested infected individuals total population N, we arrive at the testing terms −fk(t)ik=Nk u ∗ u u recover at rates γ and γ , respectively. Differences in the (Eqs. (2) and (3)) or −fk(t)ik=[Nk(1−Ik =Nk)], respectively. recovery rates γu and γ∗ reflect differences in disease severity Biased testing can also be represented by using a testing u b u b ∗ of and treatment options for untested and tested infected fraction of the form Ik e =(Ik e +Sk +Ik +Rk), where b > 0 individuals.
Recommended publications
  • Arxiv:1705.10225V8 [Stat.ML] 6 Feb 2020
    Bayesian stochastic blockmodelinga Tiago P. Peixotoy Department of Mathematical Sciences and Centre for Networks and Collective Behaviour, University of Bath, United Kingdom and ISI Foundation, Turin, Italy This chapter provides a self-contained introduction to the use of Bayesian inference to ex- tract large-scale modular structures from network data, based on the stochastic blockmodel (SBM), as well as its degree-corrected and overlapping generalizations. We focus on non- parametric formulations that allow their inference in a manner that prevents overfitting, and enables model selection. We discuss aspects of the choice of priors, in particular how to avoid underfitting via increased Bayesian hierarchies, and we contrast the task of sampling network partitions from the posterior distribution with finding the single point estimate that maximizes it, while describing efficient algorithms to perform either one. We also show how inferring the SBM can be used to predict missing and spurious links, and shed light on the fundamental limitations of the detectability of modular structures in networks. arXiv:1705.10225v8 [stat.ML] 6 Feb 2020 aTo appear in “Advances in Network Clustering and Blockmodeling,” edited by P. Doreian, V. Batagelj, A. Ferligoj, (Wiley, New York, 2019 [forthcoming]). y [email protected] 2 CONTENTS I. Introduction 3 II. Structure versus randomness in networks 3 III. The stochastic blockmodel (SBM) 5 IV. Bayesian inference: the posterior probability of partitions 7 V. Microcanonical models and the minimum description length principle (MDL) 11 VI. The “resolution limit” underfitting problem, and the nested SBM 13 VII. Model variations 17 A. Model selection 18 B. Degree correction 18 C.
    [Show full text]
  • On Community Structure in Complex Networks: Challenges and Opportunities
    On community structure in complex networks: challenges and opportunities Hocine Cherifi · Gergely Palla · Boleslaw K. Szymanski · Xiaoyan Lu Received: November 5, 2019 Abstract Community structure is one of the most relevant features encoun- tered in numerous real-world applications of networked systems. Despite the tremendous effort of a large interdisciplinary community of scientists working on this subject over the past few decades to characterize, model, and analyze communities, more investigations are needed in order to better understand the impact of community structure and its dynamics on networked systems. Here, we first focus on generative models of communities in complex networks and their role in developing strong foundation for community detection algorithms. We discuss modularity and the use of modularity maximization as the basis for community detection. Then, we follow with an overview of the Stochastic Block Model and its different variants as well as inference of community structures from such models. Next, we focus on time evolving networks, where existing nodes and links can disappear, and in parallel new nodes and links may be introduced. The extraction of communities under such circumstances poses an interesting and non-trivial problem that has gained considerable interest over Hocine Cherifi LIB EA 7534 University of Burgundy, Esplanade Erasme, Dijon, France E-mail: hocine.cherifi@u-bourgogne.fr Gergely Palla MTA-ELTE Statistical and Biological Physics Research Group P´azm´any P. stny. 1/A, Budapest, H-1117, Hungary E-mail: [email protected] Boleslaw K. Szymanski Department of Computer Science & Network Science and Technology Center Rensselaer Polytechnic Institute 110 8th Street, Troy, NY 12180, USA E-mail: [email protected] Xiaoyan Lu Department of Computer Science & Network Science and Technology Center Rensselaer Polytechnic Institute 110 8th Street, Troy, NY 12180, USA E-mail: [email protected] arXiv:1908.04901v3 [physics.soc-ph] 6 Nov 2019 2 Cherifi, Palla, Szymanski, Lu the last decade.
    [Show full text]
  • Community Detection with Dependent Connectivity∗
    Submitted to the Annals of Statistics arXiv: arXiv:0000.0000 COMMUNITY DETECTION WITH DEPENDENT CONNECTIVITY∗ BY YUBAI YUANy AND ANNIE QUy University of Illinois at Urbana-Champaigny In network analysis, within-community members are more likely to be connected than between-community members, which is reflected in that the edges within a community are intercorrelated. However, existing probabilistic models for community detection such as the stochastic block model (SBM) are not designed to capture the dependence among edges. In this paper, we propose a new community detection approach to incorporate intra-community dependence of connectivities through the Bahadur representation. The pro- posed method does not require specifying the likelihood function, which could be intractable for correlated binary connectivities. In addition, the pro- posed method allows for heterogeneity among edges between different com- munities. In theory, we show that incorporating correlation information can achieve a faster convergence rate compared to the independent SBM, and the proposed algorithm has a lower estimation bias and accelerated convergence compared to the variational EM. Our simulation studies show that the pro- posed algorithm outperforms the popular variational EM algorithm assuming conditional independence among edges. We also demonstrate the application of the proposed method to agricultural product trading networks from differ- ent countries. 1. Introduction. Network data has arisen as one of the most common forms of information collection. This is due to the fact that the scope of study not only focuses on subjects alone, but also on the relationships among subjects. Networks consist of two components: (1) nodes or ver- tices corresponding to basic units of a system, and (2) edges representing connections between nodes.
    [Show full text]
  • Developing a Credit Scoring Model Using Social Network Analysis
    Developing a Credit Scoring Model Using Social Network Analysis Ahmad Abd Rabuh UP821483 This thesis is submitted in fulfilment of the requirements for the award of the degree of Doctor of Philosophy in the School of Business and Law at the University of Portsmouth November 2020 1 Author’s Declaration Whilst registered as a candidate for the above degree, I have not been registered for any other research award. The results and conclusions embodied in this thesis are the work of the named candidate and have not been submitted for any other academic award. Signature: Name: Ahmad Abd Rabuh Date: the 30th of April 2020 Word Count: 63,809 2 Acknowledgements As I witness this global pandemic, my greatest gratitude goes to my mother, Hanan, for her moral support and encouragement through the good, the bad and the ugly. During the period of my PhD revision, I received mentoring and help from my future lifetime partner, Kawthar. Also, I have been blessed with the help of wonderful people from Scholars at Risk, Rose Anderson and Sarina Rosenthal, who assisted me in the PhD application and scholarship processes. Additionally, I cannot be thankful enough to the Associate Dean of Research, Andy Thorpe, who supported me in every possible way during my time at the University of Portsmouth. Finally, many thanks to my supervisors, Mark Xu and Renatas Kizys for their continuous guidance. 3 Developing a Credit Scoring Model Using Social Network Analysis Table of Contents Abstract ......................................................................................................................................... 10 CHAPTER 1: INTRODUCTION ................................................................................................. 12 1.1. Research Questions ........................................................................................................ 16 1.2. Aims and Objectives .....................................................................................................
    [Show full text]
  • Exploratory Analysis of Networks
    Exploratory Analysis of Networks James D. Wilson Data Institute Conference 2017 James D. Wilson (USF) Exploratory Analysis of Networks 1 / 35 Summarizing Networks Goals 1 Provide one number summaries of network 2 Develop hypotheses about observed data 3 Motivate predictive models 4 Augment standard multivariate analyses James D. Wilson (USF) Exploratory Analysis of Networks 2 / 35 At first glance: Network Summary Statistics Measures of Connectivity Degree: number of edges incident to each node (popularity) For directed networks, we consider In- and Out-Degree Geodesic distance: shortest path between two nodes Diameter: longest geodesic distance Clustering coefficient: fraction of triangles to total triples Reciprocity: fraction of reciprocated ties (directed graphs only) James D. Wilson (USF) Exploratory Analysis of Networks 3 / 35 Examples in Social Networks James D. Wilson (USF) Exploratory Analysis of Networks 4 / 35 At first glance: Network Summary Statistics Measures of Nodal Influence / Centrality Centrality: how ”central” is a node in the observed network? Degree: popularity of a node Eigenvector: self- neighbor- importance Betweenness: extent+ to which a node lies on the shortest path between two nodes Closeness: mean distance from node to all other vertices Authorities: vertices with useful information on a topic Hubs: vertices that point to authorities James D. Wilson (USF) Exploratory Analysis of Networks 5 / 35 Examples Internet Search PageRank Algorithm Popular web-search scoring and search algorithm (Google) Score based on eigenvector centrality HITS Algorithm (Kleinberg, 1999) Hyperlink-induced topic search Uses hub and authority scores as basis for measuring importance of webpages James D. Wilson (USF) Exploratory Analysis of Networks 6 / 35 Example: Recommendation Systems Figure: eye2data.blogspot.com Want individual (node) with most influence James D.
    [Show full text]
  • On Spectral Embedding Performance and Elucidating Network Structure in Stochastic Block Model Graphs
    On spectral embedding performance and elucidating network structure in stochastic block model graphs Joshua Cape and Minh Tang and Carey E. Priebe Department of Applied Mathematics and Statistics Johns Hopkins University August 16, 2018 Abstract Statistical inference on graphs often proceeds via spectral methods involving low- dimensional embeddings of matrix-valued graph representations, such as the graph Laplacian or adjacency matrix. In this paper, we analyze the asymptotic information- theoretic relative performance of Laplacian spectral embedding and adjacency spectral embedding for block assignment recovery in stochastic block model graphs by way of Chernoff information. We investigate the relationship between spectral embedding per- formance and underlying network structure (e.g. homogeneity, affinity, core-periphery, (un)balancedness) via a comprehensive treatment of the two-block stochastic block model and the class of K-block models exhibiting homogeneous balanced affinity struc- ture. Our findings support the claim that, for a particular notion of sparsity, loosely speaking, “Laplacian spectral embedding favors relatively sparse graphs, whereas adja- cency spectral embedding favors not-too-sparse graphs.” We also provide evidence in support of the claim that “adjacency spectral embedding favors core-periphery network structure.” 2010 Mathematics Subject Classification: Primary 62H30; Secondary 62B10. arXiv:1808.04855v1 [math.ST] 14 Aug 2018 Keywords: Statistical network analysis; random graphs; stochastic block model; Laplacian spectral embed- ding; adjacency spectral embedding; Chernoff information; vertex clustering. 1 1 Preface The stochastic block model (SBM) (Holland et al., 1983) is a simple yet ubiquitous network model capable of capturing community structure that has been widely studied via spectral methods in the mathematics, statistics, physics, and engineering communities.
    [Show full text]
  • Link Prediction in Large-Scale Complex Networks (Application to Bibliographical Networks)
    N◦ attribué par la bibliothèque | _ | _ | _ | _ | _ | _ | _ | _ | _ | _ | Université Paris Nord Doctoral Thesis Link Prediction in Large-scale Complex Networks (Application to Bibliographical Networks) Author: Manisha Pujari A thesis submitted in fulfillment of the requirements for the degree of Doctor of Philosophy in Computer Science in the research team Apprentissage Artificiel et Applications LIPN CNRS UMR-7030 Jury: Reviewer: Céline Robardet Professor INSA Lyon Bénédicte le Grand Professor Université Paris 1 Panthéon Sorbonne Examiner: Aldo Gangemi Professor SPC, Université Paris 13 Christophe Prieur Associate Professor, HDR Université Paris Diderot Director: Céline Rouveirol Professor SPC, Université Paris 13 Supervisor: Rushed Kanawati Associate Professor SPC, Université Paris 13 “ The more I learn, the more I realize how much I don’t know. ” Albert Einstein Abstract Link Prediction in Large-scale Complex Networks (Application to Bibliographical Networks) In this work, we are interested to tackle the problem of link prediction in complex net- works. In particular, we explore topological dyadic approaches for link prediction. Dif- ferent topological proximity measures have been studied in the scientific literature for finding the probability of appearance of new links in a complex network. Supervised learning methods have also been used to combine the predictions made or information provided by different topological measures. They create predictive models using various topological measures. The problem of supervised learning for link prediction is a difficult problem especially due to the presence of heavy class imbalance. In this thesis, we search different alternative approaches to improve the performance of different dyadic approaches for link prediction.
    [Show full text]
  • Dynamic Stochastic Block Models, Clustering and Segmentations in Dynamic Graphs Marco Corneli
    Dynamic stochastic block models, clustering and segmentations in dynamic graphs Marco Corneli To cite this version: Marco Corneli. Dynamic stochastic block models, clustering and segmentations in dynamic graphs. Mathematics [math]. Université Paris 1 - Panthéon Sorbonne, 2017. English. tel-01644866 HAL Id: tel-01644866 https://hal.archives-ouvertes.fr/tel-01644866 Submitted on 22 Nov 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Université Paris 1 Pahntéon-Sorbonne Laboratoire SAMM Thèse présentée en première version en vu d’obtenir le grade de Docteur, spécialité "Mathématiques Appliquées" par Marco Corneli Dynamic stochastic block models, clustering and segmentation in dynamic graphs Thèse soutenue le 17 Novembre 2017 devant le jury composé de : M. Christophe Biernacki Université Lille 1 Sciences et Technologies (Rapporteur) M. Thomas Brendan Murphy University College Dublin (Rapporteur) M. Allou-Badara Samé IFSTTAR (Examinateur) Mme. Clémence Magnien CNRS (Examinatrice) Mme. Sophie Lèbre Université Paul Valéry Montpellier 3 (Examinatrice) M. Fabrice Rossi Université Paris 1 Panthéon-Sorbonne (Directeur) M. Pierre Latouche Université Paris 1 Panthéon-Sorbonne (Co-encadrant) Alla mia famiglia... Remerciements e remercie Fabrice Rossi et Pierre Latouche pour m’avoir fait confiance J et avoir encadré ce travail de thèse.
    [Show full text]
  • Efficient Inference in Stochastic Block Models with Vertex Labels
    1 Efficient inference in stochastic block models with vertex labels Clara Stegehuis and Laurent Massoulie´ Abstract—We study the stochastic block model with two connectivity matrix and d the average degree in the network. communities where vertices contain side information in the form Underneath the Kesten-Stigum threshold, no algorithm is able of a vertex label. These vertex labels may have arbitrary label to infer the community memberships better than a random distributions, depending on the community memberships. We analyze a linearized version of the popular belief propagation guess, even though a community structure may be present [14]. algorithm. We show that this algorithm achieves the highest In this setting, it is even impossible to discriminate between accuracy possible whenever a certain function of the network a graph generated by the stochastic block model and an parameters has a unique fixed point. Whenever this function has Erdos˝ Renyi´ random graph with the same average degree, multiple fixed points, the belief propagation algorithm may not even though a community structure is present [13]. Above the perform optimally. We show that increasing the information in the vertex labels may reduce the number of fixed points and Kesten-Stigum threshold, the communities can be efficiently hence lead to optimality of belief propagation. reconstructed [11], [15]. A popular algorithm for community detection is belief propagation [5] (BP). This algorithm starts with initial beliefs I. INTRODUCTION on the community memberships, and iterates until these beliefs converge to a fixed point. Above the Kesten-Stigum threshold, Many real-world networks contain community structures: a fixed point that is correlated with the true community mem- groups of densely connected nodes.
    [Show full text]
  • A Review of Stochastic Block Models and Extensions for Graph Clustering
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Newcastle University E-Prints Lee and Wilkinson Applied Network Science (2019) 4:122 Applied Network Science https://doi.org/10.1007/s41109-019-0232-2 REVIEW Open Access A review of stochastic block models and extensions for graph clustering Clement Lee1* and Darren J. Wilkinson2,3 *Correspondence: [email protected] Abstract 1Department of Mathematics and There have been rapid developments in model-based clustering of graphs, also known Statistics, Lancaster University, Lancaster, United Kingdom as block modelling, over the last ten years or so. We review different approaches and Full list of author information is extensions proposed for different aspects in this area, such as the type of the graph, the available at the end of the article clustering approach, the inference approach, and whether the number of groups is selected or estimated. We also review models that combine block modelling with topic modelling and/or longitudinal modelling, regarding how these models deal with multiple types of data. How different approaches cope with various issues will be summarised and compared, to facilitate the demand of practitioners for a concise overview of the current status of these areas of literature. Keywords: Model-based clustering, Stochastic block models, Mixed membership models, Topic modelling, Longitudinal modelling, Statistical inference Introduction Stochastic block models (SBMs) are an increasingly popular class of models in statistical analysis of graphs or networks. They can be used to discover or understand the (latent) structure of a network, as well as for clustering purposes.
    [Show full text]
  • Interpretable Stochastic Block Influence Model: Measuring Social Influence Among Homophilous Communities
    Interpretable Stochastic Block Influence Model: measuring social influence among homophilous communities YAN LENG, MIT Media Lab, McCombs School of Business at UT-Austin TARA SOWRIRAJAN, Harvard University, MIT Media Lab ALEX PENTLAND, MIT Media Lab Decision-making on networks can be explained by both homophily and social influence. While homophily drives the formation of communities with similar characteristics, social influence occurs both within and between communities. Social influence can be reasoned through role theory, which indicates that the influences among individuals depend on their roles and the behavior of interest. To operationalize these social science theories, we empirically identify the homophilous communities and use the community structures to capture the “roles”, which affect the particular decision-making processes. We propose a generative model named Stochastic Block Influence Model and jointly analyze both the network formation and the behavioral influence within and between different empirically-identified communities. To evaluate the performance and demonstrate the interpretability of our method, we study the adoption decisions of microfinance in an Indian village. We show that although individuals tend to form links within communities, there are strong positive and negative social influences between communities, supporting the weak tie theory. Moreover, we find that communities with shared characteristics are associated with positive influence. In contrast, the communities with a lack of overlap are associated with negative influence. Our framework facilitates the quantification of the influences underlying decision communities and is thus a useful tool for driving information diffusion, viral marketing, and technology adoptions. Key Words: Social influence; Homophily; Stochastic Block Model; Community structure; Generative model 1 INTRODUCTION We are living in an increasingly connected society [4, 26, 28, 37].
    [Show full text]
  • Barabasi-Albert Scale-Free Networks • Repeat Measurements Multiple Times and Plot Histograms of Assortativity
    Random Networks © Mario Freese Hiroki Sayama [email protected] Small-World Phenomenon “Bacon Number” • BN = 1 for those who co-starred with Kevin Bacon in a film • BN = 2 for those who co-starred with actors/actresses with BN=1 … • Mostly BN <= 3 !! • The largest finite BN = 8 !! 3 “Erdös Number” • EN = 1 for those who co-authored a paper with a Hungarian mathematician Paul Erdös (1913-1996) • EN = 2 for those who co-authored a paper with authors with EN=1 … • Mostly EN <= 7 !! • The largest finite EN = 13 !! FYI - Hiroki’s EN=4 (by Bing/Microsoft Academic Search) 4 3.5 degrees of separation in FB • https://research.facebook.com/blog/three-and-a- half-degrees-of-separation/ 5 “Small-world” phenomenon • Most real-world networks are remarkably “small” – Despite a huge number of nodes involved – Even if connections are relatively sparse • Why? 6 Random Networks Classical explanation: Erdös-Rényi random network model • A network made of N nodes • Each node pair is connected randomly and independently with probability p • A small characteristic path length is realized because of randomness – Number of nodes reachable from a single node within k steps increases exponentially with k 8 Exercise • Create and plot a few ER random networks using NetworkX • Measure their properties – Network density – Characteristic path length – Clustering coefficient – Degree distribution – etc. 9 Limitation of ER networks • ER random networks have very few loops or local clusters if connection probability is small • Real-world networks are often clustered with
    [Show full text]