Generalization of Effective Conductance Centrality for Egonetworks

Total Page:16

File Type:pdf, Size:1020Kb

Generalization of Effective Conductance Centrality for Egonetworks Generalization of Effective Conductance Centrality for Egonetworks∗ Heman Shakeri1, Behnaz Moradi-Jamei2, Pietro Poggi-Corradini3, Nathan Albin3, and Caterina Scoglio1 1Electrical and Computer Engineering Department, Kansas State University, Manhattan, Kansas, USA 2Department of Statistics, Kansas State University, Manhattan, Kansas, USA and 3Department of Mathematics, Kansas State University, Manhattan, Kansas, USA (Dated: October 9, 2018) We study the popular centrality measure known as effective conductance or in some circles as information centrality. This is an important notion of centrality for undirected networks, with many applications, e.g., for random walks, electrical resistor networks, epidemic spreading, etc. In this paper, we first reinterpret this measure in terms of modulus (energy) of families of walks on the network. This modulus centrality measure coincides with the effective conductance measure on simple undirected networks, and extends it to much more general situations, e.g., directed networks as well. Secondly, we study a variation of this modulus approach in the egocentric network paradigm. Egonetworks are networks formed around a focal node (ego) with a specific order of neighborhoods. We propose efficient analytical and approximate methods for computing these measures on both undirected and directed networks. Finally, we describe a simple method inspired by the modulus point-of-view, called shell degree, which proved to be a useful tool for network science. The concept of information centrality was first intro- where σ(e) > 0, is the conductance of the edge e. Thus duced in [1] and was later reinterpreted in terms of elec- modulus is a constrained convex optimization problem trical conductance in [2]. Given a network G = (V; E) that has a unique extremal density ρ∗ when 1 < p < . and a node a V , the effective conductance centrality of This point of view allows for much more flexibility, be-1 a is defined as2 cause it can be applied to a variety of different families of X 1 objects: walks, cycles, trees, etc, and also works when the eff(a) := : (1) underlying network is directed or weighted. Moreover, C eff(a; b) b V a 2 n R modulus has very useful properties of Γ-monotonicity and countable subadditivity. where (a; b) is effective resistance distance between a eff For undirected networks the effective conductance be- and b.R Note that this measure considers every possible tween a and b is connected to Mod (Γ(a; b)) as follow path that electrical current flow might take from a to an 2 [6, 7] arbitrary sink b. The situation can be clarified by introducing the no- 1 = Mod (Γ(a; b)): (5) tion of modulus of families of walks. This is a way of (a; b) 2 measuring the richness of certain families of walks on a Reff network (and beyond, see [3{5]). Given two nodes a and In the following, we reproduce a proof for this connec- b we may consider the connecting family Γ(a; b) of all tion and how to calculate Mod2(Γ(a; b)) in symmetric walks γ from a to b. Then, given edge density ρ : E R networks using the pseudoinverse of the Laplacian. ! for p [1; ], we define `ρ (Γ) := minγ Γ `ρ (γ) where Let F be the set of all unit flows f : E R that satisfy 2 1 2 ! `ρ(γ) is the ρ-length of a walk γ: Kirchoffs node law and pass through a network G from X a to b. Namely for v V `ρ (γ) := ρ (e) : (2) 2 8 e γ 1 v = a 2 <> The p-modulus of Γ is defined as ( :f)(v) = 1 v = b r − :>0 otherwise Modp (Γ) := min Energyp (ρ) (3) `ρ(Γ) 1 ≥ corresponds to the injected currents at each node. The Namely, we minimize the energy of candidate edge- arXiv:1705.02703v2 [physics.data-an] 26 Jul 2018 energy of f is densities ρ subject to the ρ-length of every walk in Γ X being greater than or equal one, i.e., `ρ(Γ) 1. These Energy(f) := (e)f(e)2 ≥ R densities can be interpreted as costs of using the given e E edge. The energy we consider is 2 X where (e) = 1 is the resistance of edge e. A unit Energy (ρ) = σ(e)ρ (e)p ; (4) R w(e) p current flow i F is a unit flow that also satisfies Ohm's e E 2 2 law, i.e., there is a function V : V R (called a poten- tial) such that for every edge (a; b):! ∗ Corresponding author: [email protected] (a; b)i(a; b) = V(b) V(a): R − 2 Let U : V R be a vertex potential function. We I. EGOCENTRIC EFFECTIVE CONDUCTANCE ! can define a density ρU as the gradient of U, i.e., for the CENTRALITY edge e = v; w f g U U As mentioned above, eff(a) is sociocentric in the sense ρU(e) := u w C j − j that it considers all walks from a to an arbitrary node in Then, ρU is admissible for walks from a to b, whenever G. However, in practice, it can be prohibitive to scale U(a) = 0; U(b) = 1. sociocentric methods to very large networks. Moreover, Conversely, if ρ is an admissible density, then we can in real-world situations it is not feasible to have access define a potential U(x) as the infimum of `ρ(γ) over all to the entire network. Rather, one can at best know walks from a to x. With this definition, ρU = ρ, see [7]. local information up to a few neighborhood levels. For In particular, assuming each edge has a unit resistance, instance, when data is anonymized to protect privacy of network entities, identifying the sociocentric picture is X 2 X Energy(ρU) = ρU(e) = U(u) U(w) : impossible, e.g., sexual networks may be limited to the j − j e E e= u;w E number of contacts of individuals. 2 f g2 An alternative approach is to consider measures that Hence, if we substitute U with V + C, where V eff(a;b) are adapted to egonetworks (also known as neighborhood is the electric potential when a unitR current flow i F networks). An ego network Ga(r) around a node a is 2 is passing through the network with source a and sink b constructed by collecting data (nodes and edges) start- and the effective resistance between a and b is eff, then, ing from the ego a and searching G out to a predefined R order of neighborhood r 1; (a) ; where (a) is the T 1 2 f ··· g Mod (a; b) = min ρUρU = : eccentricity of node a or the maximum distance from a 2 U (6) a=0 (a; b) U eff b=1 R to nodes in G. Egonetworks are often preferred because they support By Kirchhoff's law of current conservation: more flexible data collection methods [9] and often in- X volve less expensive computation costs. Egocentric mea- a (V V ) = ( :i)(i) i;j i − j r sures are more stable [10] against network sampling and j reliable (less sensitivity) with measurement errors [11]. N N where A = [aij] R × is the adjacency matrix of G, We concentrate on unweighted (binary) networks to sim- 2 with aij = 1 if and only if i; j E. In matrix form: plify the algebra, although, all of our methods and dis- 2 cussions can be easily generalized for weighted networks. LV = I (7) Thus, we let d(a; b) denote the shortest-path distance be- tween two nodes (smallest number of hops). The neigh- where L is the Laplacian matrix of G and I = :i. Be- borhood structure around an ego a is described by the cause V is defined up to an additive and the nullspacer of shells of order k: L is along the constant vector, we ground an arbitrary node k and thus reduce L by removing kth row and col- S(a; k) := y V : d(a; y) = k ; umn denoted by kL [8]. Now we can find solve (7): f 2 g and the corresponding families of walks Γ(v; S(a; k)), k k 1 k V = ( L)− I: consisting of simple walks that begin at ego v V and reach S(a; k) for the first time. Modulus allows2 a quan- k 1 we denote ( L)− by (reduced conductance matrix) tification of the richness of the family of walks, i.e., a G and obtain effective resistance between nodes a and b is family with many short walks has a larger modulus than k k a family with fewer and longer walks. Here we consider eff(a; b) = Va Vb R − (8) shell modulus Mod2(v; S(a; k)) which quantifies the ca- = + 2 Ga;a Gb;b − Ga;b pacity of walks emanating from the ego up to the shell and from (6): S(a; k) [5] without having to account the data outside Ga(k). 1 Mod (a; b) = ( + 2 )− (9) 2 Ga;a Gb;b − Ga;b Theorem 1. For undirected networks, we can calcu- Therefore, using (5), we can rewrite the effective con- late 2-modulus of Γ(v; S(a; k)) analytically without going ductance centrality in (1) in the Modulus language through the optimization problem in (3): PSs 1 1 X 1 + xs − (a) = Mod (Γ(a; b)): (10) j=S1 xi Ceff 2 Mod2(a; S(a; r)) = (11) b V a xs 2 n PSs 1 For the rest of this paper, we consider p = 2 due to its where xi = − ij .
Recommended publications
  • Learning on Hypergraphs: Spectral Theory and Clustering
    Learning on Hypergraphs: Spectral Theory and Clustering Pan Li, Olgica Milenkovic Coordinated Science Laboratory University of Illinois at Urbana-Champaign March 12, 2019 Learning on Graphs Graphs are indispensable mathematical data models capturing pairwise interactions: k-nn network social network publication network Important learning on graphs problems: clustering (community detection), semi-supervised/active learning, representation learning (graph embedding) etc. Beyond Pairwise Relations A graph models pairwise relations. Recent work has shown that high-order relations can be significantly more informative: Examples include: Understanding the organization of networks (Benson, Gleich and Leskovec'16) Determining the topological connectivity between data points (Zhou, Huang, Sch}olkopf'07). Graphs with high-order relations can be modeled as hypergraphs (formally defined later). Meta-graphs, meta-paths in heterogeneous information networks. Algorithmic methods for analyzing high-order relations and learning problems are still under development. Beyond Pairwise Relations Functional units in social and biological networks. High-order network motifs: Motif (Benson’16) Microfauna Pelagic fishes Crabs & Benthic fishes Macroinvertebrates Algorithmic methods for analyzing high-order relations and learning problems are still under development. Beyond Pairwise Relations Functional units in social and biological networks. Meta-graphs, meta-paths in heterogeneous information networks. (Zhou, Yu, Han'11) Beyond Pairwise Relations Functional units in social and biological networks. Meta-graphs, meta-paths in heterogeneous information networks. Algorithmic methods for analyzing high-order relations and learning problems are still under development. Review of Graph Clustering: Notation and Terminology Graph Clustering Task: Cluster the vertices that are \densely" connected by edges. Graph Partitioning and Conductance A (weighted) graph G = (V ; E; w): for e 2 E, we is the weight.
    [Show full text]
  • Multiplex Conductance and Gossip Based Information Spreading in Multiplex Networks
    1 Multiplex Conductance and Gossip Based Information Spreading in Multiplex Networks Yufan Huang, Student Member, IEEE, and Huaiyu Dai, Fellow, IEEE Abstract—In this network era, not only people are connected, different networks are also coupled through various interconnections. This kind of network of networks, or multilayer networks, has attracted research interest recently, and many beneficial features have been discovered. However, quantitative study of information spreading in such networks is essentially lacking. Despite some existing results in single networks, the layer heterogeneity and complicated interconnections among the layers make the study of information spreading in this type of networks challenging. In this work, we study the information spreading time in multiplex networks, adopting the gossip (random-walk) based information spreading model. A new metric called multiplex conductance is defined based on the multiplex network structure and used to quantify the information spreading time in a general multiplex network in the idealized setting. Multiplex conductance is then evaluated for some interesting multiplex networks to facilitate understanding in this new area. Finally, the tradeoff between the information spreading efficiency improvement and the layer cost is examined to explain the user’s social behavior and motivate effective multiplex network designs. Index Terms—Information spreading, multiplex networks, gossip algorithm, multiplex conductance F 1 INTRODUCTION N the election year, one of the most important tasks Arguably, this somewhat simplified version of multilayer I for presidential candidates is to disseminate their words networks already captures many interesting multi-scale and and opinions to voters in a fast and efficient manner. multi-component features, and serves as a good starting The underlying research problem on information spreading point for our intended study.
    [Show full text]
  • Rumour Spreading and Graph Conductance∗
    Rumour spreading and graph conductance∗ Flavio Chierichetti, Silvio Lattanzi, Alessandro Panconesi fchierichetti,lattanzi,[email protected] Dipartimento di Informatica Sapienza Universit`adi Roma October 12, 2009 Abstract if its completion time is poly-logarithmic in the size of We show that if a connected graph with n nodes has the network regardless of the source, and that it is slow conductance φ then rumour spreading, also known as otherwise. randomized broadcast, successfully broadcasts a mes- Randomized broadcast has been intensely investi- sage within O(log4 n/φ6) many steps, with high proba- gated (see the related-work section). Our long term bility, using the PUSH-PULL strategy. An interesting goal is to characterize a set of necessary and/or suffi- feature of our approach is that it draws a connection be- cient conditions for rumour spreading to be fast in a tween rumour spreading and the spectral sparsification given network. In this work, we provide a very general procedure of Spielman and Teng [23]. sufficient condition{ high conductance. Our main moti- vation comes from the study of social networks. Loosely 1 Introduction stated, we are looking after a theorem of the form \Ru- mour spreading is fast in social networks". Our result is Rumour spreading, also known as randomized broadcast a good step in this direction because there are reasons or randomized gossip (all terms that will be used as to believe that social networks have high conductance. synonyms throughout the paper), refers to the following This is certainly the case for preferential attachment distributed algorithm. Starting with one source node models such as that of [18].
    [Show full text]
  • Absorbing State 34, 36 Adaptive Boolean Networks 78 Adaptive Chemical Network 73 Adaptive Coupled Map Lattices 66 Adaptive Epide
    239 Index a b absorbing state 34, 36 b-value 107, 109, 112, 119 adaptive Boolean networks 78 Bak–Sneppen model 74, 75 adaptive chemical network 73 benchmark 216, 217, 234 adaptive coupled map lattices 66 Bethe–Peierls approach 225 adaptive epidemiological network betweenness of a node 68 95 bill-tracking website 6 adaptive network of coupled biodiversity 46 oscillators 83 biological evolution 73 adaptive networks 63–66, 70–72, bipartite networks 211 74, 83, 90, 98, 100, 102 bipartite structure 216 adaptive rewiring 91, 92 bipartition 218, 220 adaptive SIS model 96 bipartitioning problem 219 adjacency matrix 9, 66, 206, 207, birth-death processes 31, 37 209–211, 217, 223 bisection problem 232 aftershock 108, 113, 121, bistability 27 125–129, 131, 134, 135, 137 bivariate analysis 164 aftershock magnitudes 127 bivariate measures 162 aftershock sequence 121, 123 bivariate nonlinear approaches aftershock sequences 136 163 agent-based models 73 bivariate time-series analysis agent-based simulation 172, 173, 178 frameworks 3 block-structure 201, 203, 217, air transportation network 2, 10 222, 234, 235 anti-epileptic drugs 159, 163 Boolean networks 65, 76, 79 assortative 67 bounded rationality 25 attractor 164 box-counting technique 112, 113 autoregressive modeling 161, brain 161 172, 175 Brownian motion 6, 7 avalanche size 139 burst-firing 169, 170 average shortest path length 167 bursting 166, 167, 170 Reviews of Nonlinear Dynamics and Complexity. Edited by Heinz Georg Schuster Copyright c 2009 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim ISBN:
    [Show full text]
  • Random Boolean Networks As a Toy Model for the Brain
    UNIVERSITY OF GENEVA SCIENCE FACULTY VRIJE UNIVERSITEIT OF AMSTERDAM PHYSICS SECTION Random Boolean Networks as a toy model for the brain MASTER THESIS presented at the science faculty of the University of Geneva for obtaining the Master in Theoretical Physics by Chlo´eB´eguin Supervisor (VU): Pr. Greg J Stephens Co-Supervisor (UNIGE): Pr. J´er^ome Kasparian July 2017 Contents Introduction1 1 Biology, physics and the brain4 1.1 Biological description of the brain................4 1.2 Criticality in the brain......................8 1.2.1 Physics reminder..................... 10 1.2.2 Experimental evidences.................. 15 2 Models of neural networks 20 2.1 Classes of models......................... 21 2.1.1 Spiking models...................... 21 2.1.2 Rate-based models.................... 23 2.1.3 Attractor networks.................... 24 2.1.4 Links between the classes of models........... 25 2.2 Random Boolean Networks.................... 28 2.2.1 General definition..................... 28 2.2.2 Kauffman network.................... 30 2.2.3 Hopfield network..................... 31 2.2.4 Towards a RBN for the brain.............. 32 2.2.5 The model......................... 33 3 Characterisation of RBNs 34 3.1 Attractors............................. 34 3.2 Damage spreading........................ 36 3.3 Canonical specific heat...................... 37 4 Results 40 4.1 One population with Gaussian weights............. 40 4.2 Dale's principle and balance of inhibition - excitation..... 46 4.3 Lognormal distribution of the weights.............. 51 4.4 Discussion............................. 55 i 5 Conclusion 58 Bibliography 60 Acknowledgements 66 A Python Code 67 A.1 Dynamics............................. 67 A.2 Attractor search.......................... 69 A.3 Hamming Distance........................ 73 A.4 Canonical specific heat.....................
    [Show full text]
  • Boolean Analysis of MOS Circuits
    Bo olean Analysis of MOS Circuits Randal E Bryant Computer Science Department CarnegieMellon University Pittsburgh PA February Abstract The switchlevel mo del represents a digital metaloxide semiconductor MOS cir cuit as a network of charge storage no des connected by resistive transistor switches The functionality of such a network can b e expressed as a series of systems of Bo olean equations Solving these equations symb olically yields a set of Bo olean formulas that describ e the mapping from input and current state to the new network state This analysis supp orts the same class of networks as the switchlevel simulator MOSSIM I I and provides the same functionality including the handling of bidirectional eects and indeterminate X logic values In the worst case the analysis of an n no de network can yield a set of formulas containing a total of O n op erations However all but a limited set of dense passtransistor networks give formulas with O n total op era tions The analysis can serve as the basis of ecient programs for a varietyoflogic design tasks including logic simulation on b oth conventional and sp ecial purp ose computers fault simulation test generation and symb olic verication Keywords and phrases switchlevel networks symb olic analysis logic simulation fault simulation simulation accelerators Intro duction The switchlevel mo del has proved successful as an abstract representation of digital metaloxide semiconductor MOS circuits for a variety of applications This mo del repre sents a circuit in terms of its exact transistor
    [Show full text]
  • Percolation Theory Are Well-Suited To
    Models of Disordered Media and Predictions of Associated Hydraulic Conductivity A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science By L AARON BLANK B.S., Wright State University, 2004 2006 Wright State University WRIGHT STATE UNIVERSITY SCHOOL OF GRADUATE STUDIES Novermber 6, 2006 I HEREBY RECOMMEND THAT THE THESIS PREPARED UNDER MY SUPERVISION BY L Blank ENTITLED Models of Disordered Media and Predictions of Associated Hydraulic Conductivity BE ACCEPTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF Master of Science. _______________________ Allen Hunt, Ph.D. Thesis Advisor _______________________ Lok Lew Yan Voon, Ph.D. Department Chair _______________________ Joseph F. Thomas, Jr., Ph.D. Dean of the School of Graduate Studies Committee on Final Examination ____________________ Allen Hunt, Ph.D. ____________________ Brent D. Foy, Ph.D. ____________________ Gust Bambakidis, Ph.D. ____________________ Thomas Skinner, Ph.D. Abstract In the late 20th century there was a spill of Technetium in eastern Washington State at the US Department of Energy Hanford site. Resulting contamination of water supplies would raise serious health issues for local residents. Therefore, the ability to predict how these contaminants move through the soil is of great interest. The main contribution to contaminant transport arises from being carried along by flowing water. An important control on the movement of the water through the medium is the hydraulic conductivity, K, which defines the ease of water flow for a given pressure difference (analogous to the electrical conductivity). The overall goal of research in this area is to develop a technique which accurately predicts the hydraulic conductivity as well as its distribution, both in the horizontal and the vertical directions, for media representative of the Hanford subsurface.
    [Show full text]
  • Submodular Hypergraphs: P-Laplacians, Cheeger Inequalities and Spectral Clustering
    Submodular Hypergraphs: p-Laplacians, Cheeger Inequalities and Spectral Clustering Pan Li 1 Olgica Milenkovic 1 Abstract approximations have been based on the assumption that each hyperedge cut has the same weight, in which case the We introduce submodular hypergraphs, a family underlying hypergraph is termed homogeneous. of hypergraphs that have different submodular weights associated with different cuts of hyper- However, in image segmentation, MAP inference on edges. Submodular hypergraphs arise in cluster- Markov random fields (Arora et al., 2012; Shanu et al., ing applications in which higher-order structures 2016), network motif studies (Li & Milenkovic, 2017; Ben- carry relevant information. For such hypergraphs, son et al., 2016; Tsourakakis et al., 2017) and rank learn- we define the notion of p-Laplacians and derive ing (Li & Milenkovic, 2017), higher order relations between corresponding nodal domain theorems and k-way vertices captured by hypergraphs are typically associated Cheeger inequalities. We conclude with the de- with different cut weights. In (Li & Milenkovic, 2017), Li scription of algorithms for computing the spectra and Milenkovic generalized the notion of hyperedge cut of 1- and 2-Laplacians that constitute the basis of weights by assuming that different hyperedge cuts have new spectral hypergraph clustering methods. different weights, and that consequently, each hyperedge is associated with a vector of weights rather than a single scalar weight. If the weights of the hyperedge cuts are sub- 1. Introduction modular, then one can use a graph with nonnegative edge Spectral clustering algorithms are designed to solve a relax- weights to efficiently approximate the hypergraph, provided ation of the graph cut problem based on graph Laplacians that the largest size of a hyperedge is a relatively small con- that capture pairwise dependencies between vertices, and stant.
    [Show full text]
  • Artificial Neural Networks Based on Memristive Devices
    SCIENCE CHINA Information Sciences . REVIEW . June 2018, Vol. 61 060423:1–060423:14 https://doi.org/10.1007/s11432-018-9425-1 Special Focus on Memristive Devices and Neuromorphic Computing Artificial neural networks based on memristive devices Vignesh RAVICHANDRAN, Can LI, Ali BANAGOZAR, J. Joshua YANG & Qiangfei XIA* Department of Electrical and Computer Engineering, University of Massachusetts, Amherst MA 01003, USA Received 22 February 2018/Accepted 27 March 2018/Published online 15 May 2018 Abstract The advent of memristive devices and the continuing research and development in the field of neuromorphic computing show great potential as an alternative to traditional von Neumann computing and bring us ever closer to realizing a true “thinking machine”. Novel neural network architectures and algo- rithms inspired by the brain are becoming more and more attractive to keep up with computing needs, relying on intrinsic parallelism and reduced power consumption to outperform more conventional comput- ing methods. This article provides an overview of various neural networks with an emphasis on networks based on memristive emerging devices, with the advantages of memristor neural networks compared with pure complementary metal oxide semiconductor (CMOS) implementations. A general description of neural networks is presented, followed by a survey of prominent CMOS networks, and finally networks implemented using emerging memristive devices are discussed, along with the motivation for developing memristor based networks and the computational potential
    [Show full text]
  • Controllable Fabrication of Percolative Metal Nanoparticle Arrays Applied for Quantum Conductance-Based Strain Sensors
    materials Article Controllable Fabrication of Percolative Metal Nanoparticle Arrays Applied for Quantum Conductance-Based Strain Sensors Zhengyang Du, Ji’an Chen, Chang Liu, Chen Jin and Min Han * National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences and Collaborative Innovation Centre of Advanced Microstructures, Nanjing University, Nanjing 210093, China; [email protected] (Z.D.); [email protected] (J.C.); [email protected] (C.L.); [email protected] (C.J.) * Correspondence: [email protected]; Tel.: +86-25-83686248 Received: 13 October 2020; Accepted: 28 October 2020; Published: 29 October 2020 Abstract: We use gas phase deposition of well-defined nanoparticles (NPs) to fabricate closely-spaced Pd NP arrays on flexible membranes prepatterned with interdigital electrodes (IDEs). The evolution of the morphology and electron conductance of the NP arrays during deposition is analyzed. The growth of two-dimensional percolation clusters of interconnected NPs, which correlate with the percolation pathway for electron conduction in the NP deposits, is demonstrated. The percolative nature of the NP arrays permits us to finely control the percolation geometries and conductance of the NP film by controlling the NP deposition time so as to realize a precise and reproducible fabrication of sensing materials. Electron transport measurements reveal that the electrical conductance of the NP films is dominated by electron tunneling or hopping across the NP percolating networks. Based on the percolative and quantum tunneling nature, the closely-spaced Pd NP films on PET membranes are used as flexible strain sensors. The sensor demonstrates an excellent response ability to distinguish tiny deformations down to 5 10 4 strain and a high sensitivity with a large gauge factor of 200 up to × − 4% applied strain.
    [Show full text]
  • Approximating Graph Conductance: from Global to Local
    APPROXIMATING GRAPH CONDUCTANCE: FROM GLOBAL TO LOCAL YUEHENG ZHANG Abstract. In this expository paper, we introduce spectral algorithms and random walk based algorithms on graphs through the problem of finding sparse cuts in both global and local settings. First, we consider the global problem of finding a set with minimum conductance, and show how the Laplacian and normalized Laplacian matrices arise naturally when we relax this problem to an eigenvalue computation problem. Then we bound the relaxed quantity, completing a spectral approximation algorithm and a proof of Cheeger's in- equality. Next, we discuss limitations of this global algorithm and motivate the problem of finding local conductance, for which we introduce random walks on graphs. We discuss how to gain information about conductance from the con- vergence of random walks, and motivate convergence bounds related to local conductance. Finally, we briefly introduce algorithms utilizing such bounds. Throughout, we place an emphasis on providing motivation and explanations from an algorithmic point of view. Contents 1. Introduction 1 1.1. Notation and general definitions 3 1.2. Conductance of cuts and of a graph 4 2. Generalized Rayleigh quotient and variational characterization of eigenvalues 5 3. Spectral approximation of the conductance of a graph 7 3.1. Relaxing the problem of finding the conductance of a graph 7 3.1.1. Cut capacity and the Laplacian 7 3.1.2. Normalized Laplacian and the easy side of Cheeger 9 3.2. Lower bounding the relaxation and the hard side of Cheeger 11 4. Local low-conductance sets and random walks on a graph 15 4.1.
    [Show full text]
  • Graph Clustering Outline
    Graph Clustering Outline • Min s-t cut problem • Min cut problem • Multiway cut • Minimum k-cut • Other normalized cuts and spectral graph partitionings Min s-t cut • Weighted graph G(V,E) • An s-t cut C = (S,T) of a graph G = (V, E) is a cut partition of V into S and T such that s∈S and t∈T • Cost of a cut: Cost(C) = Σe(u,v) uЄS, v ЄT w(e) • Problem: Given G, s and t find the minimum cost s-t cut Max flow problem • Flow network – Abstraction for material flowing through the edges – G = (V,E) directed graph with no parallel edges – Two distinguished nodes: s = source, t= sink – c(e) = capacity of edge e Cuts • An s-t cut is a partition (S,T) of V with sЄS and tЄT • capacity of a cut (S,T) is cap(S,T) = Σe out of Sc(e) • Find s-t cut with the minimum capacity: this problem can be solved optimally in polynomial time by using flow techniques Flows • An s-t flow is a function that satisfies – For each eЄE 0≤f(e) ≤c(e) [capacity] – For each vЄV-{s,t}: Σe in to vf(e) = Σe out of vf(e) [conservation] • The value of a flow f is: v(f) = Σe out of s f(e) Max flow problem • Find s-t flow of maximum value Flows and cuts • Flow value lemma: Let f be any flow and let (S,T) be any s-t cut. Then, the net flow sent across the cut is equal to the amount leaving s Σe out of S f(e) – Σe in to S f(e) = v(f) Flows and cuts • Weak duality: Let f be any flow and let (S,T) be any s-t cut.
    [Show full text]