26 Jan
2017
26 Jan
'17
12:49 p.m.
On 26.01.2017 12:33, P-M wrote:
Hi Tiago,
This should hopefully do it:
import graph_tool.all as gt
g = gt.load_graph('graph_no_multi_reac_type.gt') gt.remove_parallel_edges(g)
print 'minimising state' state = gt.minimize_nested_blockmodel_dl(g, deg_corr=True)
print 'equilibrating Markov chain' # We will first equilibrate the Markov chain gt.mcmc_equilibrate(state, wait=1000, mcmc_args=dict(niter=10))
As explained in the documentation, before doing the mcmc, you have to enable sampling with: state = state.copy(sampling=True) (The error message is very cryptic and needs to be improved.) Best, Tiago -- Tiago de Paula Peixoto <tiago@skewed.de>