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 <[email protected]>
signature.asc
Description: OpenPGP digital signature
_______________________________________________ graph-tool mailing list [email protected] https://lists.skewed.de/mailman/listinfo/graph-tool
