9.31556123e+00, 3.22439204e-08, 5.71925258e-05, 3.70698682e-01, 3.85329914e+04, 2.41529096e+04, 2.43090986e+04, 2.86685025e+04, 1.11990152e+03, 1.46303304e+03, 9.01419378e+05, 9.16429780e+04, 7.99801158e+04, 1.69026159e+04, 5.09448820e-01, 7.84471874e-05, 7.47324740e-07, 9.20472237e-03, 7.66006690e-01, 1.17392697e-08, 1.92114387e-02, 1.51333432e+00, 5.89673352e-09, 5.41091972e-04, 1.81968382e+00, 1.80402454e-06, 2.24990216e-03, 1.43036340e+00, 1.00064585e-07, 7.26550483e-03, 1.00000000e-02, 3.64962469e-06, 2.09630718e-01, 3.55723180e-07, 2.61192682e-02, 2.62023580e-05, 1.54476111e-03, 1.01089158e-06, 8.08297025e-03, 1.10796394e-07, 6.87439549e-04, 1.28003247e-05, 2.11267099e-02, 4.86731639e-06, 2.37896838e-02, 7.98431070e-08, 2.03144337e-01, 4.92250704e-06, 2.97646552e-02, 4.87999878e-07, 8.02754496e-03, 2.80093210e-06, 1.33217623e-04, 1.17453801e-04, 3.07145084e-03, 6.03418156e-04, 1.50334426e-02, 6.65686789e-05, 1.64949367e-02, 7.59172307e-04, 5.24390424e-03, 1.89608777e-04, 1.80147595e-03, 1.30296705e-04, 1.47706376e-03, 1.14915557e+00, 6.26994247e-06, 6.42049097e-05, 9.89376804e-01, 1.58266178e-05, 4.84818743e-03, 6.05651201e+00, 7.66416641e-06, 2.52347502e-04, 9.60186500e+00]) tspan = np.linspace(0, 20160, 1000) run_flux_visualization(model, tspan, parameters=parames, verbose=False) print('finished')
from corm import model from visualization.species_visualization_bidirectional_equilibration import run_flux_visualization import numpy as np # tipe 1 cluster: 8509 # type 2 cluster: 7848 # type 3 cluster: 5113 all_dream_log_parames = np.load("/home/oscar/PycharmProjects/CORM/results/2015_02_02_COX2_all_traces.npy")[1784] pysb_sampled_parameter_names = ['kr_AA_cat2', 'kcat_AA2', 'kr_AA_cat3', 'kcat_AA3', 'kr_AG_cat2', 'kr_AG_cat3', 'kcat_AG3', 'kr_AA_allo1', 'kr_AA_allo2', 'kr_AA_allo3', 'kr_AG_allo1', 'kr_AG_allo2'] generic_kf = np.log10(1.5e4) param_dict = {pname: pvalue for pname, pvalue in zip(pysb_sampled_parameter_names, all_dream_log_parames)} for pname, pvalue in param_dict.items(): # Sub in parameter values at current location in parameter space if 'kr' in pname: model.parameters[pname].value = 10 ** (pvalue + generic_kf) elif 'kcat' in pname: model.parameters[pname].value = 10 ** pvalue tspan = np.linspace(0, 10, num=100) run_flux_visualization(model, tspan) print('finished')