Thanks to DAVID SCHOCH, I leveraged his function and slightly modified to create the following. It works perfectly fine.
Now, all I need is to calculate these variables i.e. run functions on the graph object H in parallel, maybe using purrr or furrr, to run this faster. My real data is big with 30 more functions. The output should be same as H_indices shown in the last line of the code bolow.
library(igraph); library(sna); library(centiserve); library(tidygraph); library(tibble); library(expm)
H <- play_islands(5, 10, 0.8, 3)
all_indices <- function(g) {
tibble(
degree = igraph::degree(g),
flowbet = sna::flowbet(get.adjacency(g,sparse=F)),
communibet = centiserve::communibet(g))
}
H_indices <- all_indices(H)