Efficient Community Detection in Boolean Composed Multiplex Networks

Networks (or graphs) are used to model the dyadic relations between entities in a complex system. In cases where there exists multiple relations between the entities, the complex system can be represented as a multilayer network, where the network in each layer represents one particular relation (or feature). The analysis of multilayer networks involves combining edges from specific layers and then computing a network property. Different subsets of the layers can be combined. For any Boolean combination operation (e.g. AND, OR), the number of possible subsets is exponential to the number of layers. Thus recomputing for each subset from scratch is an expensive process. In this paper, we propose to efficiently analyze multilayer networks using a method that we term network decomposition. Network decomposition is based on analyzing each network layer individually and then aggregating the analysis results. We demonstrate the effectiveness of using network decomposition for detecting communities on different combinations of network layers. Our results on multilayer networks obtained from real-world and synthetic datasets show that our proposed network decomposition method requires significantly lower computation time while producing results of high accuracy.



There are no comments yet.



Core Decomposition in Multilayer Networks: Theory, Algorithms, and Applications

Multilayer networks are a powerful paradigm to model complex systems, wh...

A Multilayer Structure Facilitates the Production of Antifragile Systems in Boolean Network Models

Antifragility is a property to not only resist stress and but also to be...

A New Community Definition For MultiLayer Networks And A Novel Approach For Its Efficient Computation

As the use of MultiLayer Networks (or MLNs) for modeling and analysis is...

Global and Individualized Community Detection in Inhomogeneous Multilayer Networks

In network applications, it has become increasingly common to obtain dat...

Consensus Community Detection in Multilayer Networks using Parameter-free Graph Pruning

The clustering ensemble paradigm has emerged as an effective tool for co...

ALMA: Alternating Minimization Algorithm for Clustering Mixture Multilayer Network

The paper considers a Mixture Multilayer Stochastic Block Model (MMLSBM)...

Multilayer Network of Language: a Unified Framework for Structural Analysis of Linguistic Subsystems

Recently, the focus of complex networks research has shifted from the an...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

Networks (or graphs) are used to represent the pair-wise relationship between entities in a system. The entities may be connected by multiple relations (or features). For example, traffic accidents can be related if they occur in the same location, or under the same light condition; actors can be related if they acted in movies of the same genre; authors can be related if they publish in the same conferences. The relationships pertaining to each feature can be represented as a network. This group of networks are together termed as a multiplex or homogeneous multilayer network, where each network (layer) denotes a distinct relationship based on a feature among the same set of entities.

Motivation. Often different sets of layers need to be combined and analyzed. For example, if a city with limited resources was trying to reduce accidents, they might be interested in finding the top two conditions that lead to accidents. In this case, the combination of all pairs of layers have to be tested. Another example would be to identify which subset of social networking platforms form the most strongly connected groups. In this case, all possible combinations of the layers representing the platforms have to analyzed.

In the current approach, multiplexes are analyzed by combining the required layers into a single network, which we term the composed network. Network analysis operations such as community detection or finding high centrality vertices is then applied to this composed network. It is easy to see that the analysis has to be recomputed from scratch for different combinations of layers, rendering the analysis process very expensive. As a solution, we propose an elegant method, which we term network decomposition, for finding communities in different compositions of the network layers.

Our Contribution. Our main contribution is to apply network decomposition for efficiently finding communities in different Boolean compositions of network layers. The principle is to first analyze the property, such as communities, in each individual layer and then aggregate the results to obtain the final results on the composed network (see Figure 1). Thus, we only need to analyze each network once, and then combine the results as per the aggregation method.

Figure 1: Illustration of network decomposition. Each layer is analyzed separately and then the results are aggregated.

Our goal is to apply network decomposition to efficiently find communities (i.e. strongly connected groups on nodes) in layers of multiplex that are combined using the Boolean AND and OR operations. For the AND (OR) operation, the combined network will contain an edge, if there exists an edge in all (any one) of the individual layers. An AND-Composition represents how multiple features together affect an analysis. For example, in identifying regions that become accident prone due to poor lighting conditions as well as bad roads. An OR-Composition represents how any one of the features affects a property. For example, in finding the group of people who are friends via least one of the social networking platforms among Facebook, Linkedin and Twitter.

Using network decomposition, we generate communities for each of the individual layers and then combine the resultant communities based on the Boolean operators. As our results in Section IV demonstrate, our approach significantly reduces the computational costs while providing results of high accuracy.

Problem Statement: Given a set of layers , that are combined using a Boolean operation and a community detection algorithm , develop an aggregation algorithm , such that ).

Our goal is to find an aggregation algorithm , such that the results of finding the communities in the individual layers and then aggregating them via , should be the nearly the same as the communities obtained from the composed network where the layers are combined using the Boolean operator  111We use ”nearly equal” rather than ”equal” because community detection is non-unique and changes to the order in which vertices are processed can slightly change the results..

The remainder of the paper is organized as follows. In Section II, we provide a brief description of community detection in multiplexes. In Section III, we present our contribution of community detection using network decomposition. In Section IV we present the experimental results. In Section V, we discuss related research and conclude in Section VI with a discussion of our future plans.

Ii Community Detection in Multiplexes

We describe how multiplex networks are created from datasets, and how communities are obtained in AND and OR composition of layers.

Ii-a Creating Multiplex Networks

In a multirelational dataset, entities are related through multiple features. For example, people can be connected through different social networks. Authors of academic papers can be connected based on the common journals where they publish.

In a multiplex or multilayer network, each feature is represented as a separate network or layer. The set of entities (or nodes) remain the same in each layer. The edges, based on interactions between the entities change across the layers with respect to the corresponding feature. Table I lists the notations of the terms of multilayer networks discussed in this paper.

Number of layers
Set of entities
Set of features/layers
or The layer
Represntative node for entity in the layer
Set of nodes in the layer
An edge in the layer
Set of edges in the layer
or The community in the layer
Set of nodes in
Set of edges in
Table I: List of notations used for defining the concepts.

Ii-A1 Modeling the IMDb Dataset as a Multiplex

We use the Internet Movie Database (IMDb) to illustrate how a multiplex is constructed. The IMDb is an online database that contains information on television programs and movies including actors, directors, genre, and year of release [data/type2/IMDb].

We create a multiplex where the entities represent actors and two actors are connected to each other if they have acted in the same movie. Each layer in the multiplex represents a movie genre, such as comedy, drama, action, etc. As per the notations in Table I, the IMDb multiplex is defined follows:

  • : The set of actors form the entities, = .

  • : The set of features corresponds to the genres, = . The number of layers, , is equal to .

  • or : The layer represents the relations among the set of actors, , with respect to the genre .

  • : In the layer, , the actor, , is represented by a vertex, .

  • : If the and the actors have worked together in at least one movie that belongs to the genre, , , contains the edge .

Figure 2: Example of the IMDb Multiplex for co-actors with 16 actors and two genres, comedy and drama.

In Figure 2 we have selected two genres, comedy () and drama () to form the two layers, and , respectively. This multiplex shows the co-actor relationship among 16 actors (denoted by nodes numbered from 1 to 16) with respect to these genres. The same 16 actors are present in both layers. Note that each co-actor network has a distinct structure. By taking the information from the two networks together we can gain interesting insights to the data, as follows.

For example, actors and have never worked together in a drama, but have worked together in a comedy. Thus this pair of actors may together be more likely to be considered in a comedy, rather than a drama. Also observe that the actor is the actor with most connections in the drama genre, while in the case of comedies, actor is one of the nodes with the most connections, i.e. worked with most number of actors.

Ii-B Community Detection in Multiplexes

Community detection involves finding items with similar properties by identifying tightly connected groups of vertices. We consider non-overlapping communities, that is, there are no common vertices or edges between two communities. Figure 3 shows the communities for the composed layers, and for IMDB multiplex in Figure 2.

Bridge Edges. We term the external edges that connect two communities as bridge edges. Formally, if there exists an edge, , such that and , where , then this edge is a bridge edge. Bridge edges form links between two distinct communities. In the AND and OR composed-layer of Figure 3, the actors and belong to different communities, but are connected by a bridge edge. The actors and have a bridge edge in the OR composed network, but not in the AND-composed network.

Ii-B1 Communities in AND-Composed Layers.

AND composition of layers in a multiplex allows users to find communities that are related across multiple features. Examples of some queries that can be addressed by the AND composition are;

  • Groups of actors who have expertise in working together in both comedies and dramas (IMDb multiplex).

  • Author groups who publish in all of these conferences; ICDM, SIGMOD and VLDB (DBLP multiplex).

  • Groups of accidents that have similar conditions for all these features; light conditions, weather conditions, road conditions, and speed limit (Accident Multiplex).

The standard practice is to combine the layers using the AND operation, i.e. only edges that occur in all the layers are included. Then a community detection algorithm, such as Infomap, is executed on the combined network. This single graph approach, termed C-SG-AND, is given in Algorithm 1.

0:  Layers
0:  return - a list of communities
1:    AND …AND } { contains edges that are in all the networks , , …, .}
2:    = COMM() {Find communities in .}
Algorithm 1 Algorithm for C-SG-AND

Ii-B2 Communities in OR-Composed Layers

OR-composition forms a composed network that includes an edge if it appears in any of the layers. Algorithm 2 shows the steps of this single network based community detection using the OR operation, termed as C-SG-OR. Examples of queries that can be addressed by the OR composition are;

  • Groups of actors who have acted together in either a comedy or drama (IMDb multiplex).

  • Groups of authors who have published in at least one of these conferences, ICDM, VLDB, SIGMOD (DBLP multiplex).

  • Groups of accidents that have at least one condition in common (Accident Multiplex).

0:  Layers
0:  return - a list of communities
1:    OR …OR } { contains edges that are in at least one of the networks , , …, .}
2:    = COMM() {Find communities in .}
Algorithm 2 Algorithm for C-SG-OR
Figure 3: Composed Layer Communities of the IMDb Multiplex shown in Figure 2
Figure 4: Communities in each layer of the IMDb Multiplex
0:  Communities from layers and : COMM() = {, , …, }, COMM() = {, , …, }
0:  return - a list of communities
1:    {Initialize the set of communities to NULL.}
2:  for each community pair say, and  do
3:       =( ) {Create new combined community by taking the common vertices of every pair of communities.}
4:       = {Add new community to the set of communities.}
5:  end for
Algorithm 3 Algorithm for CV-AND

Iii Our Contribution: Network Decomposition for Efficient Community Detection on Multiplex

The Boolean composition of the layers of a multiplex provides in-depth analysis of the database. However, for any single Boolean operation, say AND, different combinations are possible. Thus the cost of finding the communities on each of them separately is very expensive. Moreover, if the networks do not change, several computations are rendered redundant. For example, consider finding the communities in the composed layer and . In this case, the composed layer related to remains unchanged, but has to be recomputed.

As a solution, we propose network decomposition for efficient community detection on multiplex networks. In network decomposition, the communities in each layer are identified separately and the results are then aggregated to obtain the results with respect to the composed network. Note that the storage required is only of the order of , where is the number of vertices in each layer and is the number of features/layers. Figure 4 shows the communities in each of the layers of the example IMDB network.

The challenge is to develop aggregation algorithms, , that can correctly aggregate the communities from each of the layers to obtain the communities over the composed network. We present the aggregation methods for AND and OR composition. For ease of understanding we will discuss the algorithms with respect to two layers. Note, however, that any binary operations can be easily extended to multiple layers.

Iii-a Vertex based Community Detection of AND Composed Layers (CV-AND)

An earlier paper on obtaining communities in AND-composed layers was presented in [ICCS/SantraBC17]. We discuss this work, termed CV-AND, here for completeness.

CV-AND (see Algorithm 3) heavily depends on the presence of self-preserving communities. A community is self preserving if the vertices in it are so tightly connected such that any connected subset of the vertices will form a smaller community rather than joining an existing larger community. Formally, consider a graph , with a community whose vertices are given by the set . Now consider a subset of vertices . If the vertices in form a community by themselves, for any subset of , where , then community is self preserving. The main result of  [ICCS/SantraBC17] was that if the communities from the layers are self preserving, then the communities of the AND-composed graph can be obtained by taking the vertex based intersection of the communities from the individual layers.

Drawbacks The main drawback of CV-AND is that for most networks, there is no guarantee that the communities will be self-preserving. If this algorithm is applied without testing for self-preserving communities, the results may not be accurate.

As an example, consider the community in the comedy layer of the network (Figure 4). This community is not self preserving, and when combined with community in the drama layer, which has the same vertices, it gives one large community, { , , ,,,}. In reality, as seen in Figures 3, two separate communities are formed, { ,,} and { , ,}.

This is a subtle but important difference because the community id determines whether two entities are similar. If two disconnected groups of vertices are placed in the same community (as is possible when using CV-AND), then, two dissimilar groups are marked to be similar, which is incorrect.

0:  Communities from layers and : COMM() = {, , …, }, COMM() = {, , …, }
0:  return - a list of communities
1:    {Initialize the set of communities to NULL.}
2:  for each community pair say, and  do
3:      {} = ( ) {Create list of k new communities by taking the common edges of every pair of communities.}
4:       = } {Add new communities to the set of communities.}
5:  end for
Algorithm 4 Algorithm for CE-AND

Iii-B Edge based Community Detection of AND Composed Layers (CE-AND)

We address these limitations by developing a community detection method, CE-AND (see Algorithm  4), that is based on the intersection of edges rather than vertices as follows.

For every pair of communities, from layer and from layer , the edge-based community intersection, , will produce k disconnected edge-sets, . These edge sets will form the AND-composed communities, .

Correctness and Limitations. Figure 5 shows how the communities are obtained for the example network using CE-AND. Comparing this result to that in Figure 3, we see that most of the communities are obtained with the exception of the singleton node 8 and the common bride edge (1, 5).

Lemma 1.

Let be a set of vertices that are assigned to the same community in all the layers of the multiplex. Let the set of common edges in all the layers whose both endpoints are in , be . If , then the community composed of vertices will be detected by CE-AND.

Proof. The proof is based on the observation that if a subgraph has vertices and edges, then the subgraph is connected. Since this subgraph is a community and is present in all the layers, therefore the community will be detected. Lemma 1 provides a lower bound () on the density of common edges of the communities that is needed for them to be detected by the CE-AND.

Lemma 2.

The CE-AND method returns all the communities of size that were returned by the CV-AND method.

Proof. The proof is derived from Lemma 1. All self-preserving communities have to be densely connected, and thus satisfying the lower bound. Thus CE-AND will detect all the common self preserving communities.

Lemma 3.

Common edges that form a bridge edge in at least one layer will not be detected by CE-AND.

Proof. This proof follows from Algorithm 4 in that both endpoints of the common edge have to be in the same community to be considered.

Lemma 3 highlights the limitations of finding communities in cases where the common communities do not share a dense substructure. Note that these limitations only occur for communities that are not dense. Since communities are informally defined as tightly connected groups of vertices, these sparse subpgraphs in most cases may not be relevant communities at all. We illustrate the issues in Figure 6.

The left-hand panels of Figure 6 shows two layers. The top right panel shows the communities obtained by the standard single network approach (C-SG-AND). The bottom right panel shows the communities obtained by CE-AND.

Figure 5: AND-Composition Communities of the IMDb Multiplex shown in Figure 2, using CE-AND method
Figure 6: Effect of Bridge Edges on AND Composition

Note that the community produced by C-SG-AND contains the edges (h, o) and (l, s) that act as bridges in both Layer L1 and L2. Thus CE-AND is not able to detect this community, and instead produces two communities, and , which should be merged into one by taking the bridge edges into account.

Also consider the community which consists of the edges (a, i) and (e, m) that are bridges in Layer L1, but are part of the community in Layer L2. As only those edges that are within community in all layers are considered, CE-AND produces two communities, and .

Iii-C Edge based Community Detection of OR Composed Layers

We now consider how to obtain communities in OR-composed networks. Note that the number of edges in the OR-composed network will at least as much, generally more, than the number of edges in each layer. Thus for any two layers and , the total number of edges is .

The computational complexity of community detection algorithms are at least proportional to the size of the graph. Thus the denser the graph, the more time will be required to find the communities. Thus for the OR-composed case, the goal is not only to lower the time by reducing the need to recompute different compositions of layers, but also to reduce the size of the graph to be analyzed.

To recreate the communities of OR Composed Layers, we propose the CE-OR algorithm (given in Algorithm 5 and illustrated in Figure 7). The CE-OR method reduces the size of the graph to be analyzed by leveraging the fact that the common communities across layers can be processed as a single node. The steps of the CE-OR algorithm are as follows;

0:  Communities from layers and : COMM() = {, , …, }, COMM() = {, , …, }
0:  return - a list of communities { Find common communities using CE-AND}
1:  Apply CE-AND on and to get Construct OR-MG {Assign nodes of each common community as a meta node}
2:  for each community  do
3:      =
4:  end for{ Assign the vertices not in any common community as a meta node}
5:  for each vertex  do
6:      {Create null set}
7:      {Add to the set}
8:      =
9:  end for{Add Edges in the metagraph. Two metanodes, are connected if there is an intra-community edge from one constituent node of to a constituent node of in any one of the layers.}
10:  for all  all metanode pairs  do
11:     if  : or , and  then
12:         =
13:     end if
14:  end for
15:  Insert weights on the edges of OR-MG
16:  L = COMM(OR-MG)
17:  Expand the community representative nodes in each community from L to get
Algorithm 5 Algorithm for CE-OR

Overview of CE-OR. Find the common communities in all the network layers (Line 1) by using CE-AND. Then construct a metagraph (OR-MG), as follows. Each metanode represents a set of vertices. Combine each of the vertices in a common community into a metanode (Line 2-4). Assign all remainder vertices, that are not part of any common community into singleton sets. Each of these sets is also a metanode (Line 5-9). Connect two metanodes, via a metaedge, if there exists an intra-community (within community) edge, in any one of the layers between an element (node) of and an element (node) of (Line 10-14). Apply appropriate weights to these edges (Line 15). Apply community detection on the metagraph (Line 16). The communities in the OR-composed network are obtained by expanding the metanodes in the communities obtained by the CE-OR algorithm.

Assigning Weights to Metaedges. Note that the metanodes represent vertex sets of varying sizes, and the number of edges between them represent the degree of similarity. Therefore although the original graph is composed on unweighted edges, the edges in the metagraph have to be weighted to quantify the extent of this similarity. A critical component of the CE-OR algorithm is based on correctly assigning these weights.

We propose two different weight metrics to quantity the similarity between two meta nodes. For any meta edge , let and be the set of entity nodes in the AND-composed communities, respectively. Further, let the set of intra-community edges between and be . Then the weight to the metaedge can be computed as follows;

  • Aggregation: The weight is the number of edges between the two communities;

  • Fractional: The weight is the fraction of connected nodes between the two communities; .

Correctness and Limitations: Figure 7 illustrates how the CE-OR algorithm is applied to identify communities in the OR-composed layers of the example IMDb graph. First the CE-AND communities obtained in Figure 5 and the remaining vertex are used to form the metanodes (Figure 7 (a)). Then these nodes are connected based on the intra-community edges. These edges are weighted in the meta graph using (Figure 7 (b)). A community detection algorithm on the metagraph produces the communities of the OR-composed layers (Figure 7 (c)). Comparing with the communities obtained by the C-SG-OR method in Figure 3, to those obtained by expanding the communities in the metanodes (Figure 7 (d)), we see that all the communities have been obtained. However, the bridge edges between the communities are missing.

Figure 7: Illustration of CE-OR Algorithm on the IMDB Example Graph
Lemma 4.

Let , , be a set of vertices. Let be the set of edges whose both endpoints are in , i.e. , and . If all , form an intra-community edge in at least one of the network layers and the vertices in form a community in the OR-composed network, then this community will be detected by the CE-OR algorithm.


As per Algorithm 5, any edge that is an intra-community edge in at least one of the layers will be included in the metagraph. If a community in the OR-composed layer is formed only of intra-edges, then in the metagraph, given appropriate weighting function, these set of edges would also form a dense subgraph and hence a community. ∎

Lemma 5.

If a community in the OR-composed layer contains significant number of edges that are bridge (inter-community) edges in all the layers, then the community cannot be detected in its entirety by the CE-OR algorithm.


Since the CE-OR algorithm analyzes only the intra-community edges, the bridge edges will not be part of the metagraph. Thus any community that is composed of mainly bridge edges will not be part of the metagraph and hence will not be detected by CE-OR. ∎

Figure 8: Effect of bridge edges on OR Composition

Figure 8, illustrates the effect of bridge nodes on the accuracy of the communities found by CE-OR. The left-hand of panels show two layers of the network. The top right panel shows the communities obtained by the standard single network approach (C-SG-OR). The bottom right panel shows the communities obtained by our proposed CE-OR method.

Consider the community generated by C-SG-OR approach that has edges (i, l), (h, m), (j, o) and (l, t) which are not intra-community edges in any of the layers, and are present as bridge edges in only one of the layers. These edges will not be part of the metagraph and thus CE-OR does not know that they exist. CE-OR, thus, generates three communities , and , instead of merging them into one community, as per the C-SG-OR method.

However in the community generated by C-SG-OR, the edges (a, b) and (d, f) are bridge edges in one layer but are intra-community edges in another layer. Therefore these edges will be part of the metagraph. Thus CE-OR can use these edges and correctly generate the community .

The primary limitations of our CE-AND and CE-OR algorithms is due to the non-inclusion of bridge edges. In the AND-composed network, we rationalize this non-inclusion by positing that communities formed solely of bridge edges cannot be dense, and hence are not strong communities. In the OR-composed network, note that we only exclude an edge if it is a bridge edge in all the layers. This is an infrequent case where bridge nodes from all layers come together to form communities. Our experiments in Section IV, justify this policy of not including bridge edges by demonstrating that the normalized mutual information (NMI) values between the communities returned by CE-AND and C-SG-AND are in general high.

Iv Empirical Results

In this section we compare the performance and accuracy of our proposed algorithms with the ground truth results obtained by the standard methods, C-SG-AND and C-SG-OR.

Iv-a Experimental Setup:

Since the results of community detection depend heavily on the type of algorithm used [Abrahao:2012], to control this parameter in the experiment we use the popular community detection algorithm Infomap [InfoMap2014], both to find the communities in the single network approach and the network decomposition approach. Our algorithms were implemented in C++ and were executed on a Linux machine with 8 GB RAM and installed with UBUNTU 16.10.

Datasets Used. We performed our experiments on multiplexes created from three real-world datasets and one synthetic dataset created using the RMAT [RMAT] graph generator. We selected the real-world datasets such that they were sufficiently large and contained communities. To test on larger networks with more vertices, we created the synthetic RMAT dataset. The details of the datasets are as follows (also see Table II);

  • IMDb: From the IMDB dataset [data/type2/IMDb], we created the following three layers in the multiplex, with the nodes representing the actors. In the first layer, (L1, co-acting) two nodes are connected if they co-acted in at least one movie. In the second layer, (L2, rating) two nodes are connected if the average ratings of the movies where they acted were similar. In the third layer, (L3, genre) two nodes are connected if they acted in movies of similar genres.

  • DBLP: From the DBLP dataset of academic publications [data/type2/DBLP], we selected all papers published from 2000-2018 in top three conferences VLDB (L1), SIGMOD(L2) and ICDM (L3). The nodes were the authors. Two authors in each layer were connected if they had published a paper in the conference corresponding to the layer.

  • Accident: From the dataset of road accidents that occurred in the United Kingdom in 2014 [UKDataset2014], we represented each accident as a node. Two nodes are connected in a layer if they occurred within 10 miles of each other and have similar Light (L1), Weather (L2) or Road Surface Conditions (L3).

  • RMAT: The RMAT generator creates networks based on the Kronecker product of a matrix. We set the number of vertices to

    and the edges to roughly eight times the number of vertices. We set the probabilities in each quadrant of the matrix as a=0.65, b=c=d=0.15 to create a scale-free graph.

    The first layer (L1) was the graph obtained by the generator. We applied cross perturbation to the other layers. That is we selected two edges (a, b) and (c, d), and replaced them with new edges (a, c) and (b, d). Thus the number of edges remain the same, but the degree distribution and the structure of the graph changes.

    In layer L2 we applied perturbation to 1% of the edges and in layer L3 to 5% of the edges. We limited the number of perturbations because if the network structure is significantly changed between the layers, the number of bridge edges also increase.

Ground Truth and Accuracy Metrics: Since our goal is to achieve the results obtained by the standard C-SG-AND and C-SG-OR methods, we use the communities obtained from these methods as the ground truth. We disregard communities of just one vertex, since these result due to an artifact of the algorithm rather than provide any meaningful analysis. We use two metrics to evaluate the accuracy of the communities - i) Normalized Mutual Information (NMI) that measures the quality with respect to the participating entity nodes only and ii) modified-NMI that also takes into account the topology of the communities. For both metrics higher is better, with maximum value of 1 and minimum of 0(definitions in  [Labatut13]).

Each multiplex has 3 layers. Thus, a total of 4 compositions are possible (3 for 2-layers and 1 3-layers). Thus we compare results for 8 (4 combinations X 2 Boolean operations) composed networks.

Name Vertices Edges Edges Edges
in L1 in L2 in L3
IMDB 9,485 45,581 13,945,912 996,527
DBLP 17,204 5,831 17,737 12,986
Accident 5000 193,860 235,175 216,397
RMAT 32,768 230,445 230,445 230,445
Table II: Summary of the sizes of the multiplexes.

Iv-B Accuracy of the Aggregation Algorithms.

For the AND-composed networks we show in Figure 9, the average NMI and m-NMI of all the four multiplexes with respect to the ground truth for the CV-AND and CE-AND methods. The results show that the accuracy obtained with CE-AND is higher than that from CV-AND.

Multiplex L1, L2 L1, L3 L2, L3 L1, L2, L3
Accuracy Values using CE-AND
IMDB .97 .93 .98 .97 .88 .86 .99 .99
DBLP .92 .84 .99 .96 .98 .96 .98 .95
Accident .96 .98 .94 .98 .91 .96 .88 .95
RMAT .92 .82 .90 .79 .90 .78 .90 .77
Accuracy Values using CE-OR using Fractional Weights
IMDB .01 .01 .97 .99 1 1 1 1
DBLP .83 .79 .87 .80 .75 .60 .71 .56
Accident .88 .93 .94 .98 .98 .99 .86 .93
RMAT .74 .64 .76 .59 .75 .55 .73 .54
Table III: Accuracy Values using CE-AND on the different compositions of the datasets.
Figure 9: Comparison of Accuracy of CE-AND and CV-AND based on NMI and m-NMI.

For the OR-composed networks we show in Figure 10, the average NMI and m-NMI of all the four multiplexes with respect to the ground truth for the two weighting metrics; Aggregation () and Fractional (). The results show that the accuracy obtained using both the metrics are similar.

Figure 10: Accuracy of CE-OR with Different Weighting schemes based on NMI and m-NMI.

In Table III we provide the accuracy values for all the different layer compositions with respect to CE-AND for the AND composition and CE-OR with Fractional Weights. As can be seen nearly all the values are high, .

Some low values occur for the CE-OR method. An egregious example is IMDb (L1, L2) for which the accuracy results are less than 1%! In this case the metagraph had 193 nodes, and on running the community detection algorithm 56 communities were obtained. However, the ground truth communities obtained by C-SG-OR had only 2 communities. This happened because there existed many bridge edges in the layers that were not included in the metagraph. Moreover, because the communities represented in the metanodes were small in size, the weights were also lower and could not combine the communities.

Iv-C Performance of the Aggregation Algorithms

We now compare the time taken to obtain the communities using the aggregation methods (CV-AND, CE-AND and CE-OR) with respect to C-SG-AND and C-SG-OR.

Figure 11 shows that the time to compute the communities over all the 4-composed layers is significantly lower for both CV-AND and CE-AND methods than C-SG-AND. When the layers are sparse, CE-AND will be faster than CV-AND, as can be seen for DBLP multiplex. However if the network layers are dense, then the edge-based intersection approach of CE-AND has a higher cost as compared to the CV-AND.

Figure 11: Efficiency of CV-AND and CE-AND as compared to C-SG-AND
Figure 12: Efficiency of CE-OR as compared to C-SG-OR

Figure 12 gives the time for executing CE-OR. For CE-OR, CE-AND is used as a subroutine. One scan of community edges is required to generate the meta graph (OR-MG) on which we apply Infomap. If the layers are sparse and the multiplex contains many bridge nodes, then cost of generating the meta graph and applying Infomap will become an overhead as compared to simply applying Infomap on OR graph (C-SG-OR approach). This can be seen from the DBLP multiplex where sparse layers (density of densest layer (SIGMOD) = 0.0001!) make the CE-OR 67% less efficient as compared to C-SG-OR. However, for multiplexes with fewer bridge edges (IMDb, Accident), CE-OR is significantly faster.

V Related Work

Analysis of Multilayer Networks. Homogeneous multilayer networks [MultiLayerSurveyKivelaABGMP13, Boccaletti20141] are used to handle interactions among the same set of entities such as co-authorship [MinSubMulLayer2012], citation across different topics [ng2011multirank], and relationships across different social media platforms [magnani2013formation]

. Multiplexes have also been analyzed in form of adjacency tensors

[de2013mathematical, maruhashi2011multiaspectforensics].

Community detection in temporal homogeneous multilayer networks use spectral optimization of the modularity function [zhang2017modularity, bazzi2016community]. Techniques based on information theory have been proposed for reducing the number of layers in multilayer protein-protein interactions  [LayerAggDomenicoNAL14]. Software for creating or analysing multiplex networks include Muxviz [de2014muxviz, softwareMN/muxviz], MAMMULT [battiston2014structural, softwareMN/MAMMULT], MultiTensor [de2017community] and Pymnet [softwareMN/Pymnet]. However, they support only a few simple analysis functions, not community detection.

Community Detection is a well-studied problem in monoplex (single network) analysis. A large number of competitive approximation algorithms exist (see reviews in [Fortunato2009, Xie2013]) including algorithms for  [Yang10], weighted  [Liu2011], dynamic networks  [porterchaos13], as well as parallel algorithms  [StaudtSM14, BhowmickSrinivasan13].

Recent work has also extended community detection algorithms for multilayer networks based on matrix factorization [dong2012clustering], pattern mining [silva2012mining], cluster expansion philosophy [li2008scalable] and Bayesian probabilistic models [xu2012model], but all of these use the single layer approach. To the best of our knowledge, we are the first to infer the communities of the combined network from communities of individual layers.

Vi Conclusion and Future Work

In this paper, we presented algorithms for efficiently finding communities in Boolean composed layers of multiplex networks. The results show that for most cases our algorithms are significantly faster than the standard methods and produce results of similar quality. The only cases that our algorithm fails is when the layers have significantly more bridge edges.

In future plan, we will investigate how to include some percentage of bridge edges without reducing the computation time. We also plan to explore adaptive techniques that can select between the network decomposition and standard methods as suitable. Finally, we also aim to develop methods to efficiently find communities when the NOT operator is used.