def main_exp(): L = 10 num_sites = 16 alphas = [0.5/16]*16 mr_systems = [mr_system(alphas,use_annealing=True,scale=1,iterations=500000, sse_epsilon=0.0001)[-1] for i in range(10)] mr_ics = map(motif_ic,mr_systems) mr_ginis = map(motif_gini,mr_systems) init_states = [random_motif_with_dirty_bits(L,num_sites)] ic = 12 M = 100 tau = 1 timesteps = 100 epsilon = 0.1 bins = [-10] + myrange(0,ic,epsilon) + [ic,ic + epsilon]+ [ic + 10] control_motifs = concat([weighted_ensemble(mutate_motif_with_dirty_bits, motif_ic_with_dirty_bits, init_states, bins, M, tau, timesteps, final_bin_index=None,verbose=1) for i in range(5)]) control_ics = map(lambda ((m,ics),p):motif_ic(m),concat(control_motifs)) control_ginis = map(lambda ((m,ics),p):motif_gini(m),concat(control_motifs)) binned_controls = [[gini for ic,gini in zip(control_ics,control_ginis) if abs(ic-mr_ic) < 0.1] for mr_ic in mr_ics]
def motif_prob_scan_exp(): L = 10 num_sites = 16 alphas = [0.5/16]*16 verbose=False mr_systems = [[mr_system(alphas,use_annealing=True,scale=1,iterations=5000000, sse_epsilon=0.0001,motif_prob=alpha,verbose=verbose)[-1] for i in range(5)] for alpha in myrange(0.3,1,0.1)]