Dear Tiago, all,
I'm happy to have discovered the new partition overlap functionality. It works
very nicely. But is it possible to obtain the description length of a nested
consensus state?
Getting the description length of the consensus state works fine for an
unnested model:
from graph_tool.all import *
h = collection.data['polbooks']
num_models = 2
states = [minimize_blockmodel_dl(h) for i in range(num_models)]
bs = [states[i].get_blocks().a for i in range(num_models)]
c, r = partition_overlap_center(bs)
h_blocks = h.new_vp('int')
h_blocks.a = c
BlockState(h, h_blocks).entropy()
For nested partitions, getting the consensus state is a bit more complicated:
states_nested = [minimize_nested_blockmodel_dl(h) for i in range(num_models)]
bs_nested = []
for j in range(num_models):
l = []
for k in range(len(states_nested[j].get_levels())):
l.append(states_nested[j].get_bstack()[k].vp.b.a)
bs_nested.append(l)
c_nested, r_nested = nested_partition_overlap_center(bs_nested)
But how can I get its description length?
Many thanks for your attention
Haiko
GESIS, Cologne
_______________________________________________
graph-tool mailing list
[email protected]
https://lists.skewed.de/mailman/listinfo/graph-tool
