Use this if you are using igraph from R
| sample_hierarchical_sbm {igraph} | R Documentation | 
Sampling from a hierarchical stochastic block model of networks.
sample_hierarchical_sbm(n, m, rho, C, p)
hierarchical_sbm(...)
n | 
 Integer scalar, the number of vertices.  | 
m | 
 Integer scalar, the number of vertices per block.   | 
rho | 
 Numeric vector, the fraction of vertices per cluster, within a
block. Must sum up to 1, and   | 
C | 
 A square, symmetric numeric matrix, the Bernoulli rates for the
clusters within a block. Its size must mach the size of the   | 
p | 
 Numeric scalar, the Bernoulli rate of connections between vertices in different blocks.  | 
... | 
 Passed to   | 
The function generates a random graph according to the hierarchical stochastic block model.
An igraph graph.
Gabor Csardi csardi.gabor@gmail.com
## Ten blocks with three clusters each
C <- matrix(c(1  , 3/4,   0,
              3/4,   0, 3/4,
              0  , 3/4, 3/4), nrow=3)
g <- sample_hierarchical_sbm(100, 10, rho=c(3, 3, 4)/10, C=C, p=1/20)
g
if (require(Matrix)) { image(g[]) }