我有一个有向图(顶点= 20,边<100)。使用了edge.betweenness.community()函数,找到了5个集群。簇的大小是不同的:14,3,1,1,1 脚本在下面。
ga.data<- read.csv('q2_4.csv', header=TRUE)
g <- graph.data.frame(ga.data, directed=TRUE)
gnc <- edge.betweenness.community(g, directed=TRUE)
windows()
plot(gnc, g, vertex.frame.color="black", edge.color="black",
edge.arrow.size=0.5, edge.width=.1, vertex.label.font=NA)
# print the lables of vertices
for (i in 1: length(gnc)) {
print(V(g2)[clusters(g2)$membership == i])
}
输出结果为:
# + 14/20 vertices, named:
# [1] A19 A15 A09 A16 A20 A01 A06 A13 A10 A18 A08 A14 A02 A03
# + 1/20 vertex, named:
# [1] A17
# + 1/20 vertex, named:
# [1] A11
# + 3/20 vertices, named:
# [1] A12 A04 A07
# + 1/20 vertex, named:
# [1] A05
我想监督群集的数量并将所有顶点划分为4个群集。我在剧本中添加了一行:
print(gnh<- cutat(gnc, no = 4)) # set number of clusters = 4
#[1] 1 1 3 4 2 1 1 2 1 1 2 1 1 1 1 1 1 1 1 1
最后一行中的一定数量的社区。 如何为结果指定名称?