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)