def test_modification(): st1 = Phosphorylation(Agent('BRAF'), Agent('MAP2K1'), 'S', '222') sa = SifAssembler([st1]) sa.make_model(True, True, True) assert (len(sa.graph.nodes()) == 2) assert (len(sa.graph.edges()) == 1) sa.save_model('test_sif.sif', True) with open('test_sif.sif', 'rb') as fh: txt = fh.read().decode('utf-8') assert txt == 'BRAF 0 MAP2K1\n', txt
bn_str = boolean2.modify_states(bn_str, turnon=on, turnoff=off) model = boolean2.Model(text=bn_str, mode='async') for i in range(nsim): model.initialize() model.iterate(steps=nsteps) coll.collect(states=model.states, nodes=model.nodes) avgs = coll.get_averages(normalize=True) return avgs if __name__ == '__main__': # Build Boolean net for basic pathway st = ac.load_statements('ras_pathway.pkl') sa = SifAssembler(st) sa.make_model(use_name_as_key=True) sa.save_model('ras_pathway.sif') bn_str = sa.print_boolean_net('ras_pathway_bn.txt') # Build Boolean net for extended pathway st_ext = ac.load_statements('ras_pathway_extension.pkl') sa = SifAssembler(st + st_ext) sa.make_model(use_name_as_key=True) sa.save_model('ras_pathway_extension.sif') bn_str = sa.print_boolean_net('ras_pathway_extension_bn.txt') # Condition 1 off = [] on = ['GROWTH-FACTOR'] avgs = get_sim_avgs(bn_str, off=off, on=on) jun_basic_noinh = avgs['JUN']