sample_hierarchical_sbm: Sample the hierarchical stochastic block model

In igraph: Network Analysis and Visualization

Sample the hierarchical stochastic block model

Description

Sampling from a hierarchical stochastic block model of networks.

Usage

sample_hierarchical_sbm(n, m, rho, C, p)

hierarchical_sbm(...)

Arguments

n	Integer scalar, the number of vertices.
m	Integer scalar, the number of vertices per block. n / m must be integer. Alternatively, an integer vector of block sizes, if not all the blocks have equal sizes.
rho	Numeric vector, the fraction of vertices per cluster, within a block. Must sum up to 1, and rho * m must be integer for all elements of rho. Alternatively a list of rho vectors, one for each block, if they are not the same for all blocks.
C	A square, symmetric numeric matrix, the Bernoulli rates for the clusters within a block. Its size must mach the size of the rho vector. Alternatively, a list of square matrices, if the Bernoulli rates differ in different blocks.
p	Numeric scalar, the Bernoulli rate of connections between vertices in different blocks.
...	Passed to sample_hierarchical_sbm().

Details

The function generates a random graph according to the hierarchical stochastic block model.

Value

An igraph graph.

Related documentation in the C library

igraph_hsbm_game(), igraph_hsbm_list_game().

Author(s)

Gabor Csardi [e-mail]

See Also

Random graph models (games) erdos.renyi.game(), sample_(), sample_bipartite(), sample_chung_lu(), sample_correlated_gnp(), sample_correlated_gnp_pair(), sample_degseq(), sample_dot_product(), sample_fitness(), sample_fitness_pl(), sample_forestfire(), sample_gnm(), sample_gnp(), sample_grg(), sample_growing(), sample_islands(), sample_k_regular(), sample_last_cit(), sample_pa(), sample_pa_age(), sample_pref(), sample_sbm(), sample_smallworld(), sample_traits_callaway(), sample_tree()

Examples

## 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[])
}

igraph documentation built on Oct. 20, 2024, 1:06 a.m.