Em comentário à minha resposta a esta pergunta, é pedido para calcular valores bootstrap da modularidade de um grafo agrupado com o algoritmo cluster_louvain
. Mas a dificuldade é que por vezes o código tentado dá mais que um valor para a modularidade dos grafos aleatórios. Isso pode ser visto com o código abaixo.
suppressPackageStartupMessages(
library(igraph)
)
# primeiro calcular o valor da modularidade
louvain <- cluster_louvain(matrixnetwork, weights = NULL, resolution = 0.88)$modularity
louvain
#> [1] 0.3911249
# bootstrap da modularidade louvain
Nperm <- 1e3
set.seed(123)
boot_mod <- replicate(Nperm, {
randomnet <- rewire(matrixnetwork, with = each_edge(0.5))
E(randomnet)$weight <- sample(E(matrixnetwork)$weight)
cluster_louvain(randomnet, weights = NULL, resolution = 0.88)$modularity
})
class(boot_mod)
#> [1] "list"
mean(lengths(boot_mod) > 1L)
#> [1] 0.309
plot(density(boot_mod), main = main)
#> Error in density.default(boot_mod): argument 'x' must be numeric
Created on 2023-12-09
O objeto boot_mod
é de classe "list"
porque em cerca de 30% dos casos a modularidade é mais que um valor. E a densidade dá erro.
Como não se pode escolher arbitrariamente um qualquer valor da modularidade, como resolver o problema?
Dados
Os dados são os mesmos dos da pergunta do link acima, com a modificação na minha resposta.
suppressPackageStartupMessages(
library(igraph)
)
### DADOS EXEMPLO ###
valores <- c(1, 5, 3, 8, 2, 9, 3, 2, 3 ,2, 5, 3, 6, 5, 1, 5, 3, 4, 2, 2, 5, 6, 5, 7, 1, 2,
8, 12, 5, 1, 5, 3, 6, 5, 9, 3,3,4,5,7,8,2,7,8,4,3,1,4,4 )
cor_pequena <- 0.1
matrixnet <- matrix(valores, 7, 7)
matcor <- cor(matrixnet, method = "spearman")
matcor[] <- ifelse(matcor < cor_pequena, 0, matcor)
matrixnetwork <- graph.adjacency(matcor, mode = "undirected", weighted = TRUE,
add.colnames = NULL, diag = FALSE)
Created on 2023-12-09