La méthode permet d'effectuer le partitionnement d'un réseau en optimisant la modularité. La modularité est une valeur comprise entre -1/2 et 1 qui mesure la densité d'arêtes à l'intérieur des communautés comparée à celle des arêtes reliant les communautés entre elles. L'optimisation de la modularité conduit théoriquement au meilleur partitionnement possible, cependant son calcul effectif pour chacun des groupements de nœuds possibles est trop long en pratique, d'où l'utilisation d'algorithmes heuristiques. La première phase de la méthode de Louvain consiste à trouver de petites communautés par optimisation locale de la modularité sur chacun des nœuds. Dans un second temps, les nœuds d'une même communauté sont regroupés en un unique nœud, et la première phase est répétée sur le réseau nouvellement obtenu. La méthode est similaire à la méthode antérieure proposée par Clauset, Newman et Moore[2], qui amalgame les communautés dont la fusion produit la plus grande augmentation de la modularité.
La valeur à optimiser est la modularité, définie entre -1/2 et 1, et mesurant la densité des arêtes à l'intérieur des communautés par rapport à la densité des arêtes entre les communautés[3]. Pour un graphe pondéré, la modularité Q est définie comme :
où :
représente le poids de l'arête entre les nœuds et ;
et sont la somme des poids des arêtes liées aux nœuds et respectivement ;
est la somme de l'ensemble des poids des arêtes du graphe ;
Afin de maximiser cette valeur de manière efficace, la méthode de Louvain a deux phases qui se répètent de manière itérative.
Tout d'abord, chaque nœud du réseau est affecté à sa propre communauté. Ensuite, pour chaque nœud , on calcule le changement de modularité occasionné par la suppression de de sa propre communauté et son déplacement dans la communauté de chacun des voisins de . Cette valeur est calculée en deux étapes : (1) la suppression de de sa communauté d'origine, et (2) l'insertion de communauté de . Les deux équations sont assez similaires, et l'équation pour l'étape (2) est :
Ici, est la somme des poids des arêtes à l'intérieur de la communauté de destination de , est la somme des poids de toutes les arêtes liées aux nœuds de la communauté de destination de , est le degré pondéré de , est la somme des poids des arêtes entre et les autres nœuds de la communauté de destination de , et est la somme des poids de l'ensemble des arêtes du graphe. Ensuite, une fois cette valeur calculée pour toutes les communautés auxquelles est connecté, est placé dans la communauté qui entraine la plus grande augmentation de la modularité. Si aucune augmentation n'est possible, reste dans sa communauté d'origine. Ce processus est appliqué de manière répétée et séquentielle sur tous les nœuds jusqu'à ce qu'aucune augmentation de la modularité ne se produise. Une fois que ce maximum local est atteint, la première phase se termine.
Dans la seconde phase de l'algorithme, tous les nœuds d'une même communauté sont regroupés, et un nouveau réseau est construit où les nœuds sont les communautés de la phase précédente. Les arêtes entre les nœuds d'une même communauté sont représentés par des boucles sur le nouveau nœud, et les arêtes d'une communauté à une autre sont représentées par des arêtes pondérées. Une fois le nouveau réseau créé, la seconde phase se termine et la première phase peut être appliquée de nouveau sur le réseau résultant.
Visualisation d'un sous-réseau d'interactions sociales sur Twitter, usant de la méthode louvain pour détecter des "communautés"Réseau socialTwitter (2,4 millions de nœuds, 38 millions de liens) par Josep Pujol, Vijay Erramilli, et Pablo Rodriguez[4].
Réseau de téléphonie mobile (4 millions de nœuds, de 100 millions de liens) par Derek Greene, Donal Doyle, et Padraig Cunningham[5].
Lors de la comparaison des méthodes d'optimisation de la modularité, les deux mesures d'importance sont la vitesse de calcul et la valeur de la modularité. Une vitesse élevée permet à la méthode d'être adaptée pour traiter de larges volumes de données. Une modularité élevée indique des communautés mieux définies. Les méthodes comparées sont l'algorithme de Clauset, Newman, et Moore[6], Pons et Latapy[7], et Watika et Tsurumi[8].
Comparaison de méthodes d'optimisation de la modularité[9]
-/- dans le tableau signifie que la méthode a pris plus de 24 heures pour s'exécuter. Ce tableau[10] montre que la méthode de Louvain surpasse de nombreuses méthodes d'optimisation de la modularité similaires, tant en termes de modularité qu'en temps de calcul.
Dans le cas de graphe avec des arêtes orientées , Leicht and Newman [11] ont proposé une définition
avec resp. la somme des arcs orientés entrants en resp. sortants de .
Dugué et Perez adaptent l'algorithme de Louvain pour optimiser la modularité , démontrant les performances de l'algorithme et sa pertinence dans le cadre orienté [12],[13].
↑Fast unfolding of communities in large networks, Vincent D Blondel, Jean-Loup Guillaume, Renaud Lambiotte, Étienne Lefebvre, Journal of Statistical Mechanics: Theory and Experiment 2008 (10), P10008 (12pp) doi: 10.1088/1742-5468/2008/10/P10008. ArXiv: https://arxiv.org/abs/0803.0476
↑Multilevel local optimization of modularity, T. Aynaud, V.D. Blondel, J.-L. Guillaume, R. Lambiotte - in Graph Partitioning, 315-345, Publisher John Wiley & Sons, 2011.
↑(en) Nicolas Dugué et A. Perez Directed Louvain: maximizing modularity in directed networks. (rapport), (lire en ligne)
↑(en) Nicolas Dugué et Anthony Perez, « Direction matters in complex networks: A theoretical and applied study for greedy modularity optimization », Physica A: Statistical Mechanics and its Applications, vol. 603, , p. 127798 (DOI10.1016/j.physa.2022.127798, lire en ligne)
![{\displaystyle Q={\frac {1}{2m}}\sum \limits _{ij}{\bigg [}A_{ij}-{\frac {k_{i}k_{j}}{2m}}{\bigg ]}\delta (c_{i},c_{j}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bd7498b437de5c6bc34f7cc44fba306122350202)
![{\displaystyle \Delta Q={\bigg [}{\frac {\Sigma _{in}+k_{i,in}}{2m}}-{\bigg (}{\frac {\Sigma _{tot}+k_{i}}{2m}}{\bigg )}^{2}{\bigg ]}-{\bigg [}{\frac {\Sigma _{in}}{2m}}-{\bigg (}{\frac {\Sigma _{tot}}{2m}}{\bigg )}^{2}-{\bigg (}{\frac {k_{i}}{2m}}{\bigg )}^{2}{\bigg ]}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f1152ebb31a1a547269ad696e708282f25c2c5c8)
![{\displaystyle Q_{d}={\frac {1}{m}}\sum \limits _{ij}{\bigg [}A_{ij}-{\frac {k_{i}^{in}k_{j}^{out}}{m}}{\bigg ]}\delta (c_{i},c_{j}),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e0acd81eca81ec0178797483d73c15bfbcb597b0)
