Sample from a DCSBM

Sample an adjacency matrix from a degree-corrected block model (DCSBM)

sample_dcsbm(z, B, theta = 1)

Arguments

z

Node labels (\(n * 1\))

B

Connectivity matrix (\(K * K\))

theta

Node connectivity propensity vector (\(n * 1\))

Value

An adjacency matrix following DCSBM

Examples

B = pp_conn(n = 10^3, oir = 0.1, lambda = 7, pri = rep(1,3))$B
head(sample_dcsbm(sample(1:3, 10^3, replace = TRUE), B, theta = rexp(10^3)))
#> 6 x 1000 sparse Matrix of class "dgCMatrix"
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . 1 . . . . . . . . . 1 . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . 1 .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . .
#> [3,] . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . 1 . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . 1 . . . . . . 1 . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . 1 . . . . . . . . . . . . . . 1 . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>                                                                               
#> [1,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [2,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [3,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [4,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [5,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#> [6,] . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
#>       
#> [1,] .
#> [2,] .
#> [3,] .
#> [4,] .
#> [5,] .
#> [6,] .