r之R中的网络模块化计算

Network Modularity 的方程在其维基百科页面(和著名的书籍)上给出。我想看到它在一些代码中工作。我发现这可以使用 modularity library对于 igraph与 R (The R Foundation for Statistical Computing) 一起使用。

我想看看代码中用于计算模块性的下面的示例(或类似的示例)。该库给出了示例，但这并不是我真正想要的。

让我们有一组顶点 V = {1, 2, 3, 4, 5} 和边 E = {(1,5), (2,3), (2,4), (2,5) (3 ,5)} 形成一个无向图。

将这些顶点划分为两个社区:c1 = {2,3} 和 c2 = {1,4,5}。要计算的是这两个社区的模块化。

请您参考如下方法:

library(igraph) 
g <- graph(c(1,5,2,3,2,4,2,5,3,5)) 
membership <- c(1,2,2,1,1) 
modularity(g, membership)