我已经开始使用R中的网络。我知道R中有一个名为igraph
的软件包,它具有许多用于网络分析和可视化的功能。我想计算一个小网络的程度。 e.g。
mydegree=c(2,3,5)
g=degree.sequence.game(mydegree)
现在要获得每个节点的程度,我可以简单地使用degree
中的igraph
命令并说出degree(g)
。但是,如果我不想使用这个度数命令,我可以使用什么公式来计算学位?
答案 0 :(得分:1)
你真的不想使用任何igraph函数吗?我假设您至少想要使用igraph加载或生成图形,否则这个问题没有多大意义。
因此,如果您想将生成的图形转换为不需要igraph处理的图形,那么您基本上有三个选项。
get.adjlist()
,正如您在另一个答案中看到的那样,会给您一个邻接列表,然后程度由列表中每个向量的长度给出:
g <- degree.sequence.game(rep(2,10))
al.g <- get.adjlist(g)
sapply(al.g, length)
# [1] 2 2 2 2 2 2 2 2 2 2
get.adjacency()
将为您提供一个邻接矩阵,然后可以通过简单求和行(或列,如果您有有向图)来计算度:
adj.g <- get.adjacency(g)
rowSums(adj.g)
# [1] 2 2 2 2 2 2 2 2 2 2
get.edgelist()
返回两列矩阵中的简单边列表。您需要计算每个顶点在此矩阵中出现的次数以获得它们的度数:
el.g <- get.edgelist(g)
table(el.g)
# el.g
# 1 2 3 4 5 6 7 8 9 10
# 2 2 2 2 2 2 2 2 2 2
答案 1 :(得分:0)
如果您阅读igraph文档,您会看到有很多方法可以做到这一点。这是一个。获取相关列表,计算每个元素的长度:
mydegree=c(2,3,5)
g=degree.sequence.game(mydegree)
unlist(lapply(get.adjlist(g),length))
[1] 2 3 5
另一个:
> mydegree = c(3,5,2,2,4,2,5,6,7)
> g=degree.sequence.game(mydegree)
> unlist(lapply(get.adjlist(g),length))
[1] 3 5 2 2 4 2 5 6 7
另一个:
> degree(graph.famous("Zachary"))
[1] 16 9 10 6 3 4 4 4 5 2 3 1 2 5 2 2 2 2 2 3 2 2 2 5 3
[26] 3 2 4 3 4 4 6 12 17
> unlist(lapply(get.adjlist(graph.famous("Zachary")),length))
[1] 16 9 10 6 3 4 4 4 5 2 3 1 2 5 2 2 2 2 2 3 2 2 2 5 3
[26] 3 2 4 3 4 4 6 12 17