Question

我想生成随机网络，并希望将网络与原始网络进行比较，该网络具有16809个节点和173393个边缘。因此，为了将其与不同的网络模型进行比较，我将不得不生成具有相同边数的网络模型。在erdos.renyi模型中，我可以生成指定边数的随机图。如何使用r。

中的igraph库生成具有相同边数的无标度和小世界网络

我的示例脚本如下。

library(igraph)
g_erdos_renyi <- erdos.renyi.game(16809, 173393 , type = "gnm" , directed = F , loops = F)
g_scale <- barabasi.game(16809 , m = 10)
g_small <- watts.strogatz.game(1, 16809, 10, 0.05)

如何生成随机网络g_scale和g_small，边数为173393 .. ??

Answer 1

简短的回答是，它有点繁琐，您需要决定使用什么策略来舍入或缩小边数。

小世界

对于小世界而言，因为它是一个高度有序的结构，所以很难准确指定每个节点以相同的度数开始所需的边数，并随机重新连接。我能想到的最好的方法是建立下一个最大的网络并随机删除边缘：

n <- 16809
m <- 173393

# Work out how theye divide into each other
rem <- m %% n
div <- m %/% n

set.seed(123)
if(rem != 0) {
  g <- sample_smallworld(1, n, div+1, p = 0.001)
# Randomly delete the unwanted edges. Should be quite homegenous
  g <- delete_edges(g, sample(1:gsize(g), size = gsize(g) - m))
} else {
  g <- sample_smallworld(1, n, div, p = 0.001)
}

优惠附件

对于BA网络，它还需要有一定数量的边缘。您可以使用out.seq参数指定每一步添加的边数：

# Barabasi - Albert --------------------------------------------------------

genOutSeq <- function(n, m) {
  n <- n-1 # Shift it along
  rem <- m %% n
  c(0, rep(m%/%n + 1, rem), rep(m%/%n, n - rem))
}

n <- 16809
m <- 173393

# This creates the right number of edges but some are multiple
set.seed(11)
g <- sample_pa(n, power = 0.5, out.seq = genOutSeq(n, m),
               algorithm = "psumtree-multiple", directed = FALSE)
gsize(g)

set.seed(11)
g <- sample_pa(n, power = 0.5, out.seq = genOutSeq(n, m),
               algorithm = "psumtree", directed = FALSE)

# Randomly add the remainder

nMulti <- m - gsize(g) # Number of multiple edges that were removed

for (i in 1:nMulti) {
  vPair <- sample(1:n, size = 2)
  while (get.edge.ids(g, vPair) > 0) {
    add_edges(g, vPair)
    vPair <- sample(1:n, size = 2)
  }
}

g

正如您所看到的，第一次运行使用了一个生成多边的算法。我通过随机添加它来解决这个问题，但这取决于你使用的策略。

生成r中具有指定边数的随机网络模型

1 个答案:

小世界

优惠附件