mygen = Generation()
mygen.set_JDD( myest.get_JDD('full') )
mygen.set_KTRI( myest.get_KTRI('full') )
mygen.construct_simple_2K()
mygen.mcmc_random_2_5_K(nmae_threshold=0.03)
mygen.save_graphs('%s_2Ksimple+MCMC_Full' % fname)
#######################################################
###### 30% sample - 2K with triangles + Improved MCMC
elif run_case == 3:
p_sample = 0.3
myest = Estimation()
myest.load_graph(fname)
myest.sample('rw', p_sample)
myest.estimate_JDD()
myest.estimate_CCK()
mygen = Generation()
mygen.set_JDD( myest.get_JDD('realizable') )
mygen.set_KTRI( myest.get_KTRI('estimate') )
mygen.construct_triangles_2K()
mygen.mcmc_improved_2_5_K(nmae_threshold=0.05)
mygen.save_graphs('%s_2KT+ImpMCMC_%.2fsample' % (fname, p_sample))
#######################################################
###### 30% sample - 2K simple + MCMC
elif run_case == 4:
p_sample = 0.3
myest = Estimation()
myest.load_graph(fname)
myest.sample('rw', p_sample)
mygen.set_JDD(myest.get_JDD('full'))
mygen.set_KTRI(myest.get_KTRI('full'))
mygen.construct_simple_2K()
mygen.mcmc_random_2_5_K(error_threshold=error_threshold)
mygen.save_graphs('%s_2Ksimple+MCMC_Full' % fname)
#######################################################
###### 30% sample - 2K with triangles + Improved MCMC
elif run_case == 3:
p_sample = 0.4
myest = Estimation()
myest.load_graph(fname)
myest.sample('rw', p_sample)
myest.estimate_JDD()
myest.estimate_CCK()
myest.estimation_summary()
mygen = Generation()
mygen.set_JDD(myest.get_JDD('realizable'))
mygen.set_KTRI(myest.get_KTRI('estimate'))
mygen.construct_triangles_2K()
mygen.mcmc_improved_2_5_K(error_threshold=error_threshold)
mygen.save_graphs('%s_2KT+ImpMCMC_%.2fsample' % (fname, p_sample))
#######################################################
###### 30% sample - 2K simple + MCMC
elif run_case == 4:
p_sample = 0.4
myest = Estimation()
myest.load_graph(fname)
