R/get_sweep_pair_counts.R
get_sweep_pair_counts.RdRetrieves the pairwise connection propensity counts from latest MCMC sweep on an sbm_network object. Needs to have
mcmc_sweep(track_pairs = TRUE) run on the object before.
get_sweep_pair_counts(sbm)
| sbm | Object of class |
|---|
pairwise connection propensity counts from latest MCMC sweep
Other helpers:
build_score_fn(),
get_combination_indices(),
get_sweep_results(),
print.sbm_network(),
rolling_mean(),
verify_model()
# Start with a small random network net <- sim_basic_block_network(n_blocks = 3, n_nodes_per_block = 15) %>% initialize_blocks(num_blocks = 4) %>% mcmc_sweep(num_sweeps = 4, variable_num_blocks = FALSE) # Retrieve the sweep results from network get_sweep_results(net)#> $nodes_moved #> [1] "g1_8" "g2_3" "g3_9" "g3_13" "g3_11" "g1_13" "g1_9" "g3_1" "g1_2" #> [10] "g3_8" "g3_3" "g3_15" "g3_12" "g3_10" "g3_5" "g1_5" "g2_14" "g3_6" #> [19] "g2_7" "g1_2" "g3_11" "g3_10" "g3_1" "g2_2" "g2_7" "g3_4" "g2_8" #> [28] "g3_13" "g3_15" "g3_14" "g3_6" "g3_5" "g1_7" "g3_10" "g3_7" "g2_10" #> [37] "g1_2" "g3_4" "g2_1" "g1_3" "g1_13" "g3_1" "g3_15" "g1_9" "g2_14" #> [46] "g3_6" "g1_10" "g3_8" "g2_13" "g1_7" "g1_5" "g3_7" "g1_15" "g3_5" #> [55] "g3_4" "g3_13" #> #> $sweep_info #> entropy_delta num_nodes_moved #> 1 -11.9101128 19 #> 2 -23.1109183 14 #> 3 0.9113516 13 #> 4 1.7059270 10 #>