Directed Network Comparison Using Motifs

TL;DR

This paper introduces a new method to compare directed networks by analyzing their motif structures, which captures local and global differences more effectively than existing methods.

Contribution

The novel contribution is a motif-based approach for comparing directed networks using node-specific motif distributions and Jensen–Shannon divergence.

Findings

The proposed method outperforms existing baselines in comparing real directed networks and their perturbed versions.

The method is robust to different parameter settings and captures higher-order network differences effectively.

Abstract

Analyzing and characterizing the differences between networks is a fundamental and challenging problem in network science. Most previous network comparison methods that rely on topological properties have been restricted to measuring differences between two undirected networks. However, many networks, such as biological networks, social networks, and transportation networks, exhibit inherent directionality and higher-order attributes that should not be ignored when comparing networks. Therefore, we propose a motif-based directed network comparison method that captures local, global, and higher-order differences between two directed networks. Specifically, we first construct a motif distribution vector for each node, which captures the information of a node’s involvement in different directed motifs. Then, the dissimilarity between two directed networks is defined on the basis of a matrix, which is composed of the motif distribution vector of every node and the Jensen–Shannon divergence. The performance of our method is evaluated via the comparison of six real directed networks with their null models, as well as their perturbed networks based on edge perturbation. Our method is superior to the state-of-the-art baselines and is robust with different parameter settings.

Figures

Full text

1. Introduction

Many systems in various domains featuring intricate interaction relationships can be effectively represented in the form of complex networks [1], including social platforms [2,3], biological systems [4], and economic systems [5]. Due to the diversity of network forms [6,7] and the high-order features of networks [8,9], the precise measurement of similarity between different networks, namely, the design of an effective network comparison method, has emerged as a central focus in the field of network science. Network comparison aims to quantify the differences between two networks based on the network topological structure, allowing the effective handling of different types of tasks [10,11]. For example, in the field of pattern recognition, network comparison can be applied to classify contents such as images, documents, and videos [12]. In the biological domain, network comparison can be used to analyze which protein interactions may have equivalent functions [13]. In neuroscience, the comparison of brain networks contributes to understanding the functional differences between normal and pathological brains [14].

The original term used to compare networks was the graph isomorphism problem [15], which has been proven to fall within the NP complexity class [16]. In recent years, researchers have proposed various methodologies from different perspectives and technologies to measure the similarity between networks [17,18,19,20,21,22]. The majority of these methods have primarily concentrated on the comparison of undirected networks. However, interactions among distinct entities in the real world commonly exhibit asymmetry. In social networks, an instance of user i trusting user j does not necessarily imply reciprocal trust from j to i. The directionality of the interactions between nodes in a network, which can not be captured by an undirected network, has boosted the research of directed network comparison. For example, Bagrow and Bollt [23] utilized portrait divergence, a metric based on the distribution of the shortest path lengths, to evaluate the structural similarities between networks. Koutra et al. [24] proposed DeltaCon by calculating the Matusita distance of similarity matrices between two networks. Sarajlic et al. [25] extended network distance measures to directed networks using directed graphlets and demonstrated their efficacy in distinguishing various directed networks. Centrality-based methods, such as the degree [26], closeness [27], and the clustering coefficient [28], compare networks based on the centrality values of each node. Although these methods are capable of comparing networks effectively to some extent, most of them do not consider the higher-order structure, i.e., interactions among more than two nodes, of a network, which has been shown to be ubiquitous in various complex systems [9]. The interactions among multiple nodes have been modeled by simplexes, hypergraphs, and subgraphs such as motifs in work from different domains [29,30,31]. To capture the higher-order interactions between nodes in directed networks, we propose using direct motifs to quantify the dissimilarity between two networks. Motifs refer to recurring subgraphs in a network, where these subgraphs exhibit specific interaction patterns that facilitate the understanding of the functionality of networks [32]. Motifs have been widely used in different network tasks, i.e., community detection [33], link prediction [34], and node-ranking problems [35]. In contrast to traditional conventional methods, motif-based approaches consistently exhibit superior performance in tackling these problems.

To explore the similarity between different directed network structures, in this paper, we propose a motif-based directed network comparison method, $[eqn]$ , i.e., using motifs to examine smaller components of directed networks to assess the similarity between networks. We start by constructing a node motif distribution matrix, where the elements in the matrix are obtained by computing the distribution of nodes appearing in different directed motifs. Due to computational complexity, we consider the motifs composed of 2 to 4 nodes and thus obtain 35 different directed motifs. Later on, we use the Jensen–Shannon divergence to quantify the dissimilarity between two directed networks both locally and globally. We validate the effectiveness of $[eqn]$ in six real directed networks. Compared to the baseline methods, $[eqn]$ exhibits notable distinguishability and robustness in comparing networks.

The rest of this paper is organized as follows. Section 2 introduces the definition of motifs in a directed network and details the motif-based directed network comparison method. We provide a clear description of the baseline methods and directed network datasets in Section 3. All experimental results are presented in Section 4. Section 5 summarizes the full paper.

2. Method

2.1. The Definition of Motifs in a Directed Network

A directed unweighted network is represented by $[eqn]$ , where $[eqn]$ and $[eqn]$ are the sets of nodes and edges, respectively. The numbers of nodes and edges are given by N and M. The adjacent relationship between two nodes in G is given by the adjacency matrix A, with $[eqn]$ indicating that there is a directed edge between $[eqn]$ and $[eqn]$ and $[eqn]$ implying that there are no edges between them. We note that the directionality of G determines that A is an asymmetric matrix.

Motifs are the most common graphical patterns in complex networks, consisting of a group of closely connected nodes and edges. Due to the high complexity of computing motifs in a network, we normally consider motifs formed by 2 to 4 nodes. Motifs play a crucial role in the study of complex networks, acting as fundamental building blocks for large complex networks, analogous to genes in biology. In a directed network, the motifs are formed by nodes with directed edges. We show examples of directed motifs in Figure 1. There are 35 directed motifs, each comprising 2 to 4 nodes, individually represented by $[eqn]$ to $[eqn]$ , respectively. For instance, there are two kinds of motifs if we consider two nodes, which are given by $[eqn]$ and $[eqn]$ in the figure.

2.2. The Motif-Based Directed Network Comparison Method

Motifs contain important topological information of a network and thus are essential for network comparison. Based on the distinctive topological properties of directed motifs, we first compute the motif distribution in a directed network. As the time complexity of computing motifs is quite high, we will use the motifs listed in Figure 1 that are formed by $[eqn]$ , and 4 nodes for the computation of motif distribution. Specifically, we use $[eqn]$ to represent the motif distribution of node $[eqn]$ , where $[eqn]$ represents the fraction of motif j that contains $[eqn]$ . Consequently, an $[eqn]$ matrix $[eqn]$ can be constructed based on the motif distribution of every node. We further define the directed network node dispersion ( $[eqn]$ ) to measure the connectivity heterogeneity between nodes [22]. A larger $[eqn]$ indicates greater heterogeneity in the connectivity of nodes within the network, while a smaller $[eqn]$ suggests a more uniform distribution of node connections. And $[eqn]$ is given by the following formula:

$$[eqn]$$

where $[eqn]$ is the Jensen–Shannon divergence of the N motif distributions and is given by

$$[eqn]$$

where $[eqn]$ represents the average value of N motif distributions, and the specific calculation is as follows:

$$[eqn]$$

Given two directed networks $[eqn]$ and $[eqn]$ , the structural dissimilarity between them can be calculated based on their motif distribution matrices $[eqn]$ and $[eqn]$ . We use $[eqn]$ to represent the dissimilarity between $[eqn]$ and $[eqn]$ , and thus,

$$[eqn]$$

where

$$[eqn]$$

The dissimilarity $[eqn]$ comprises two terms, and we use the parameter $[eqn]$ to adjust their weights. The first term illustrates the difference between the average motif distributions, that is, $[eqn]$ and $[eqn]$ , and predominantly signifies the global distinctions between the two networks. The second term mainly describes the difference between the DNNDs of the two networks, indicating the local difference between them. A lower value of $[eqn]$ indicates a higher network similarity and vice versa.

3. Baselines and Datasets

3.1. Baselines

Portrait-based directed network comparison method [23]: For a directed network G, we construct a portrait matrix B based on the distance between nodes. Each element $[eqn]$ represents the number of nodes that have k nodes at distance l, where $[eqn]$ , $[eqn]$ , and d represents the diameter of G. We note that we utilize the shortest directed path length to calculate the distance between nodes. In addition, B is independent of the ordering and labeling of the nodes. Based on $[eqn]$ , we can derive the probability that a randomly selected node has k nodes at a distance of l and is given by

$$[eqn]$$

For two directed networks, $[eqn]$ and $[eqn]$ , the probability distributions $[eqn]$ and $[eqn]$ are employed to interpret the rows of the network portraits for each of them. The similarity between $[eqn]$ and $[eqn]$ is represented by $[eqn]$ and is defined as follows:

$$[eqn]$$

where $[eqn]$ , and $[eqn]$ represents the Kullback–Liebler divergence between two distributions.

DeltaCon-based directed network comparison method [24]: DeltaCon considers the similarity between two networks by quantifying the difference between r-step paths other than the edges. Given a directed and unweighted network G and its adjacency matrix A, the r-step paths are encoded in the similarity matrix $[eqn]$ , where D and I are diagonal matrices with diagonal elements equal to the node degree and 1, respectively, and $[eqn]$ . We assume that the similarity matrices for two directed and unweighted networks $[eqn]$ and $[eqn]$ are denoted by S and $[eqn]$ , and the dissimilarity $[eqn]$ between them is given by the following equation:

$$[eqn]$$

Closeness-based directed network comparison method: Centrality measures, such as degree, betweenness, and closeness, were used to compare networks [27]. However, in this part of the experiments, we found that closeness centrality surpasses other centrality methods in network comparison. Therefore, we omit the other centrality measures and only use closeness for directed network comparison. Closeness centrality measures the importance of a node within a network by evaluating the proximity of its connections to other nodes. The closeness centrality of a node is defined by

$$[eqn]$$

where $[eqn]$ represents the directed shortest path length from node $[eqn]$ to node $[eqn]$ . For two directed networks $[eqn]$ and $[eqn]$ , we assume that the closeness centrality vectors for them are given by $[eqn]$ and $[eqn]$ . Therefore, the dissimilarity between $[eqn]$ and $[eqn]$ based on closeness centrality is given as follows:

$$[eqn]$$

3.2. Description of Directed Network Datasets

To evaluate the performance of our proposed methods and the state-of-the-art baselines, we selected six real-world directed networks from diverse domains, including biological networks, transportation networks, and social networks. The descriptions of each of the datasets are as follows:

Mac [36] describes the interactions between adult female Japanese macaques and is about the dominance behavior between them. Each node denotes a macaque, and a directed edge from node $[eqn]$ to $[eqn]$ indicates the dominance of $[eqn]$ over $[eqn]$ .

Caenorhabditis elegans (Elegans) [37] is a neural network of Caenorhabditis elegans. It uses directed edges to represent neural connections among neurons in the nervous system of Caenorhabditis elegans.

Physicians [38] is a directed network that describes the spread of innovation among physicians. A directed edge $[eqn]$ between two physicians $[eqn]$ and $[eqn]$ implies that $[eqn]$ would turn to $[eqn]$ if he or she needs suggestions or is interested in a discussion.

Email-Eu-core (Email) [39] is an email network that captures email interactions between institution members in a large European research institution. A directed edge between two staff members $[eqn]$ and $[eqn]$ means that staff member $[eqn]$ has sent an email to staff member $[eqn]$ .

US airport [40] illustrates the flight connections between US airports. A directed edge $[eqn]$ between two airports $[eqn]$ and $[eqn]$ illustrates that there is at least a flight from airport $[eqn]$ to $[eqn]$ .

Chess [40] is a network that characterizes the interaction between players in an international chess game within a month. A directed edge is formed from a white player to a black player in this network.

Table 1 shows the basic properties of the directed networks mentioned above, including the number of nodes $[eqn]$ , the number of edges $[eqn]$ , the average degree $[eqn]$ , the average shortest path length $[eqn]$ , and the network diameter $[eqn]$ .

4. Experimental Results

4.1. The Dissimilarity between a Real Network and Its Null Models

The null model is widely used as a tool for the comparison of network topology [41] and retains specific network properties, such as the degree distribution or clustering coefficient, via the random reshuffling of network connections. In this section, we propose three null models for directed networks to gradually change the network topology and use our comparison method to compare each directed network and its null models.

We extend the $[eqn]$ -series null models that were originally proposed for undirected networks to directed networks [42], which retain the degree distributions, correlations, and clustering of a real directed network to some extent. Concretely, the models are illustrated as follows: $[eqn]$ preserves the out-degree and in-degree of a node by randomly rewiring each directed edge. Therefore, the degree sequence of the original network is preserved in the reshuffling process. $[eqn]$ reshuffles every edge in the network while maintaining the out-degree, in-degree, and joint degree distribution of the original network. $[eqn]$ rewires every edge by preserving the distribution of the degree-dependent clustering coefficient. We note that the newly formed directed edges should never have existed in the original network.

We show examples of how to generate the null models in Figure 2a–c, in which the blue dashed lines indicate the newly connected edges. The left panel shows the original network, and the right panel shows the network after the rewiring process in each of the figures. Figure 2a shows an instance for $[eqn]$ . Specifically, we disconnect the edges $[eqn]$ and $[eqn]$ and form new edges, i.e., $[eqn]$ , $[eqn]$ and $[eqn]$ , $[eqn]$ . Therefore, the in-degree and out-degree of each node are preserved in this process. Figure 2b demonstrates the generation of a random network through $[eqn]$ , which is more strict than $[eqn]$ . For example, if we disconnect the directed edge between $[eqn]$ and $[eqn]$ , that is, $[eqn]$ , we need to find a node that has the same in-degree and out-degree as $[eqn]$ , and the appropriate node is $[eqn]$ . Accordingly, we connect $[eqn]$ and $[eqn]$ and form a new directed edge $[eqn]$ . Therefore, $[eqn]$ maintains the degree sequence and the joint degree distribution of a network. In Figure 2c, the degree (sum of in-degree and out-degree) and clustering coefficient for each node are $[eqn]$ and $[eqn]$ , respectively. Therefore, the average clustering coefficients for nodes that have degrees of $[eqn]$ are $[eqn]$ , respectively, which are also called degree-dependent clustering coefficients. We disconnect the directed edges $[eqn]$ and $[eqn]$ and form new directed edges as $[eqn]$ and $[eqn]$ . In the rewired network, the degree-dependent clustering coefficient distribution is the same as that in the original network.

A lower value of k implies a greater disruption of the original network structure. In Figure 3, we use the four directed network comparison methods mentioned above to quantify the dissimilarity between each of the directed networks and its three null models. Figure 3a–d represent the comparison results via $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ , respectively. The experimental results for six networks suggest that as k increases, the similarity between the original network and its null models gradually increases for our proposed method. The dissimilarity observed in our approach aligns with the generation of null models, providing further confirmation of the effectiveness and stability of our model in comparing directed networks from different domains. For the baseline models, $[eqn]$ shows a similar trend to our method. However, $[eqn]$ and $[eqn]$ show bad performance in networks such as Email and Physicians, respectively.

4.2. The Comparison of the Directed Network and Its Perturbed Network

In this section, we report perturbation experiments performed on the edges of six real directed networks to further assess the stability and applicability of the motif-based comparison method. Specifically, for each given network, we randomly add or remove edges with a certain proportion, f, where the range of f is $[eqn]$ . The positive value of f indicates that we randomly add an $[eqn]$ fraction of directed edges to the network, and the negative value of f means that we randomly remove an $[eqn]$ fraction of directed edges. We compare the original network with the perturbed network by adding or removing edges using different network comparison methods, as shown in Figure 4. The four comparison methods ( $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ ) show similar trends; that is, the increase in $[eqn]$ will make the perturbed network have a greater difference from the original network, which is consistent with intuition. This conclusion is especially significant when f is negative, where our method, as well as the baselines, can significantly distinguish between the original network and the ones after perturbation. However, the motif-based comparison method is much better than the rest of the baselines for positive values of f. The curves of the other three baselines for $[eqn]$ are flatter than those of our method. Taking the Mac network as an example (Figure 4a), the values of $[eqn]$ range from $[eqn]$ to $[eqn]$ for $[eqn]$ , and the values of $[eqn]$ are the same for $[eqn]$ and $[eqn]$ , which is unreasonable. $[eqn]$ and $[eqn]$ also show insignificant dissimilarities between networks in Figure 4a–f. The baseline methods, such as $[eqn]$ and $[eqn]$ , are based on the distance between nodes, and $[eqn]$ considers the r-step paths of a network for network comparison. However, they do not consider the higher-order organization of a network, i.e., at the level of subgraphs, and thus may result in poor performance in network comparison. For example, the distance between nodes is obtained through pairwise interactions (i.e., edges) between nodes. In addition, the r-step paths are also constructed using edges.

4.3. Parameter Sensitivity Analysis

The motif-based directed network comparison method involves a parameter, denoted by $[eqn]$ , that determines how much importance is given to the global or local differences between two networks, with a larger value of $[eqn]$ indicting that we consider more of global difference and vice versa. Therefore, we performed a parameter analysis for $[eqn]$ in the six real-world directed networks via the comparison of the original network and its perturbed networks. The results are given in Figure 5, in which we use curves with different colors to indicate different chosen values of $[eqn]$ . The figure displays curves that exhibit a similar trend for different values of $[eqn]$ , and there is small deviation among the curves when $[eqn]$ . However, the network dissimilarity for different f is more significant for $[eqn]$ in most networks (except Physicians and Email), which means that we need to consider the global or local differences between networks for comparison. Therefore, we used $[eqn]$ in the above analysis.

5. Conclusions

In this paper, we introduce a comparison method $[eqn]$ that utilizes network motifs to assess similarities between directed networks. The method, which considers both local and global differences between two directed networks as well as higher-order information, is based on node motif distributions and employs the Jensen–Shannon divergence. In detail, we use motifs with sizes up to 4 that are listed in Figure 1 to compute the motif distribution of nodes in a directed network. Based on the Jensen–Shannon divergence and motif distributions of nodes, we define the dispersion of directed network nodes ( $[eqn]$ ) to quantify the heterogeneity of connectivity between nodes. Lastly, for two given directed networks, the similarity between them is further defined by the combination of the $[eqn]$ metrics and the average motif distributions. Our method aims to better understand the internal connection patterns of the network nodes by capturing essential subgraph structures. To show the effectiveness of our method, we compare a directed network with its null models, which gradually change the structure of the original network. In addition, we further compare our method with the baselines to characterize the similarity between an original network and its perturbed networks. The results show that our method outperforms these baseline methods across networks from different domains.

Motifs have been widely used to address a range of tasks. In our analysis, we take into account the directionality of edges by utilizing directed motifs to compare directed networks. We limit our analysis to motifs with sizes up to 4 due to the high computational expenses involved. It is worthwhile to investigate the impact on network comparison performance when analyzing motifs with varying numbers of nodes. Although considering larger motifs could potentially enhance the effectiveness of our approach, it may pose scalability challenges when dealing with large networks containing millions of nodes. Given the success of motifs in network comparison, we believe that developing efficient algorithms for computing motifs could be a promising avenue for research. This has the potential not only to enhance network comparison but also to improve other network tasks, such as community detection, node classification, influence maximization, and more.