######################################################## # # # UoM iGEM 2013: E.c(oil)i Fatty Acid Biosynthesis Uncertainty Model # # # ######################################################## import pysces import random import numpy as np import scipy, scipy.stats # To convert filetype SBML to .psc : #pysces.PyscesInterfaces.Core2interfaces().convertSBML2PSC('test.xml',sbmldir=None, pscfile=None, pscdir=None) # Functions to determine the parameters for the lognormal distribution: def getNormalParamFromLogNormParam_median(median, s): mu = np.log(median) x = np.exp(2*mu) sigma_2 = np.log( (x + np.sqrt(x*(4*(s**2)+x))) / (2*x) ) sigma = np.sqrt(sigma_2) return mu, sigma def getNormalParamFromLogNormParam_maxci_median(median, maxci): z=scipy.stats.norm.ppf(1 - 0.05/2) mu = np.log(median) sigma=(np.log(maxci)-mu)/z return mu, sigma # Loading the model into PySCeS : mod=pysces.model('MODELfileBalanced - SINK', dir='C:\\Pysces\\psc') # Reptitive loop for generating n models r=open("ParametresModele.txt","w") r1=open("SteadyStateConcentrations.txt","w") r2=open("SteadyStateFluxes.txt","w") r3=open("SteadyStateCC.txt","w") ListeParam="Model n\tBKR_OxoC4_KNADPH\tBKR_OxoC6_KNADPH\tBKR_OxoC8_KNADPH\tBKR_OxoC10_KNADPH\tBKR_OxoC12_KNADPH\tBKR_OxoC14_KNADPH\tBKR_OxoC16_KNADPH\tBKR_OxoC18_KNADPH\tBKR_OxoC4_KOxoC4\tBKR_OxoC6_KOxoC6\tBKR_OxoC8_KOxoC8\tBKR_OxoC10_KOxoC10\tBKR_OxoC12_KOxoC12\tBKR_OxoC14_KOxoC14\tBKR_OxoC16_KOxoC16\tBKR_OxoC18_KOxoC18\tBKR_OxoC4_KOHC4\tBKR_OxoC6_KOHC6\tBKR_OxoC8_KOHC8\tBKR_OxoC10_KOHC10\tBKR_OxoC12_KOHC12\tBKR_OxoC14_KOHC14\tBKR_OxoC16_KOHC16\tBKR_OxoC18_KOHC18\tBKR_OxoC4_KNADP\tBKR_OxoC6_KNADP\tBKR_OxoC8_KNADP\tBKR_OxoC10_KNADP\tBKR_OxoC12_KNADP\tBKR_OxoC14_KNADP\tBKR_OxoC16_KNADP\tBKR_OxoC18_KNADP\tDH_OHC4_KOHC4\tDH_OHC6_KOHC6\tDH_OHC8_KOHC8\tDH_OHC10_KOHC10\tDH_OHC12_KOHC12\tDH_OHC14_KOHC14\tDH_OHC16_KOHC16\tDH_OHC18_KOHC18\tDH_OHC4_KEC4\tDH_OHC6_KEC6\tDH_OHC8_KEC8\tDH_OHC10_KEC10\tDH_OHC12_KEC12\tDH_OHC14_KEC14\tDH_OHC16_KEC16\tDH_OHC18_KEC18\tENR_EC4_KEC4\tENR_EC6_KEC6\tENR_EC8_KEC8\tENR_EC10_KEC10\tENR_EC12_KEC12\tENR_EC14_KEC14\tENR_EC16_KEC16\tENR_EC18_KEC18\tENR_EC4_KNADH\tENR_EC6_KNADH\tENR_EC8_KNADH\tENR_EC10_KNADH\tENR_EC12_KNADH\tENR_EC14_KNADH\tENR_EC16_KNADH\tENR_EC18_KNADH\tENR_EC4_KC4ACP\tENR_EC6_KC6ACP\tENR_EC8_KC8ACP\tENR_EC10_KC10ACP\tENR_EC12_KC12ACP\tENR_EC14_KC14ACP\tENR_EC16_KC16ACP\tENR_EC18_KC18ACP\tENR_EC4_KNAD\tENR_EC6_KNAD\tENR_EC8_KNAD\tENR_EC10_KNAD\tENR_EC12_KNAD\tENR_EC14_KNAD\tENR_EC16_KNAD\tENR_EC18_KNAD\tKASI_C2ACP_KC3ACP\tKASI_C4ACP_KC3ACP\tKASI_C6ACP_KC3ACP\tKASI_C8ACP_KC3ACP\tKASI_C10ACP_KC3ACP\tKASI_C12ACP_KC3ACP\tKASI_C14ACP_KC3ACP\tKASI_C16ACP_KC3ACP\tKASI_C2ACP_KC2ACP\tKASI_C4ACP_KC4ACP\tKASI_C6ACP_KC6ACP\tKASI_C8ACP_KC8ACP\tKASI_C10ACP_KC10ACP\tKASI_C12ACP_KC12ACP\tKASI_C14ACP_KC14ACP\tKASI_C16ACP_KC16ACP\tKASI_C16ACP_KOxoC18\tKASI_C16ACP_KACP\tKASI_C16ACP_KCO2\tAAT_C16ACP_KC16ACP\tAAT_C18ACP_KC18ACP\tAAT_C16ACP_KC16\tAAT_C18ACP_KC18\tAAT_C16ACP_KACP\tAAT_C18ACP_KACP\tFACS_C16_KC16\tFACS_C18_KC18\tFACS_C16_KCoA\tFACS_C18_KCoA\tFACS_C16_KATP\tFACS_C18_KATP\tFACS_C16_KAMP\tFACS_C18_KAMP\tFACS_C16_KPP\tFACS_C18_KPP\tFACS_C16_KC16CoA\tFACS_C18_KC18CoA\tBKR_OxoC4_Vmax\tBKR_OxoC6_Vmax\tBKR_OxoC8_Vmax\tBKR_OxoC10_Vmax\tBKR_OxoC12_Vmax\tBKR_OxoC14_Vmax\tBKR_OxoC16_Vmax\tBKR_OxoC18_Vmax\tDH_OHC4_Vmax\tDH_OHC6_Vmax\tDH_OHC8_Vmax\tDH_OHC10_Vmax\tDH_OHC12_Vmax\tDH_OHC14_Vmax\tDH_OHC16_Vmax\tDH_OHC18_Vmax\tENR_EC4_Vmax\tENR_EC6_Vmax\tENR_EC8_Vmax\tENR_EC10_Vmax\tENR_EC12_Vmax\tENR_EC14_Vmax\tENR_EC16_Vmax\tENR_EC18_Vmax\tKASI_C2ACP_Vmax\tKASI_C4ACP_Vmax\tKASI_C6ACP_Vmax\tKASI_C8ACP_Vmax\tKASI_C10ACP_Vmax\tKASI_C12ACP_Vmax\tKASI_C14ACP_Vmax\tKASI_C16ACP_Vmax\tAAT_C16ACP_Vmax\tAAT_C18ACP_Vmax\tFACS_C16_Vmax\tFACS_C18_Vmax\tBKR_OxoC4_Keq\tBKR_OxoC6_Keq\tBKR_OxoC8_Keq\tBKR_OxoC10_Keq\tBKR_OxoC12_Keq\tBKR_OxoC14_Keq\tBKR_OxoC16_Keq\tBKR_OxoC18_Keq\tDH_OHC4_Keq\tDH_OHC6_Keq\tDH_OHC8_Keq\tDH_OHC10_Keq\tDH_OHC12_Keq\tDH_OHC14_Keq\tDH_OHC16_Keq\tDH_OHC18_Keq\tENR_EC4_Keq\tENR_EC6_Keq\tENR_EC8_Keq\tENR_EC10_Keq\tENR_EC12_Keq\tENR_EC14_Keq\tENR_EC16_Keq\tENR_EC18_Keq\tKASI_C2ACP_Keq\tKASI_C4ACP_Keq\tKASI_C6ACP_Keq\tKASI_C8ACP_Keq\tKASI_C10ACP_Keq\tKASI_C12ACP_Keq\tKASI_C14ACP_Keq\tKASI_C16ACP_Keq\tAAT_C16ACP_Keq\tAAT_C18ACP_Keq\tFACS_C16_Keq\tFACS_C18_Keq\tKASIII_C2CoA_KACP\tKASIII_C2CoA_KC2CoA\tKASIII_C2CoA_Vmax\tBCCP_Holo_KATP\tBCCP_Holo_KHCO3\tBCCP_Holo_KHolo\tBCCP_Holo_Vmax\tACC_COOHbiotinCCP_KC2CoA\tACC_COOHbiotinCCP_KCOOHbiotinCCP\tACC_COOHbiotinCCP_KHolo\tACC_COOHbiotinCCP_KC3CoA\tACC_COOHbiotinCCP_Keq\tACC_COOHbiotinCCP_Vmax\tMCT_C3CoA_KACP\tMCT_C3CoA_KCoA\tMCT_C3CoA_KC3CoA\tMCT_C3CoA_KC3ACP\tMCT_C3CoA_Keq\tMCT_C3CoA_Vmax\tSINK18_KC18\tSINK18_Vmax\n" headerSpecies="Model n\tNADH@default\tNAD@default\tNADPH@default\tC3CoA@default\tFabH@default\tC2CoA@default\tH2O@default\tNADP@default\tH@default\tCOOHbiotinCCP@default\tC14@default\tC8@default\tACP@default\tATP@default\tC14ACP@default\tHCO3@default\tC4ACP@default\tC12ACP@default\tOHC18@default\tHolo@default\tOxoC4@default\tOHC14@default\tspecies_8@default\tspecies_9@default\tC18@default\tADP@default\tOHC4@default\tspecies_2@default\tspecies_3@default\tspecies_1@default\tspecies_6@default\tspecies_7@default\tspecies_4@default\tspecies_5@default\tP@default\tC8ACP@default\tOxoC10@default\tC2ACP@default\tOxoC6@default\tCO2@default\tC16@default\tEC18@default\tEC16\tOxoC14@default\tspecies_14@default\tspecies_15@default\tspecies_11@default\tOxoC12@default\tOxoC16@default\tC16ACP@default\tC18ACP@default\tC18@default\tEC10@default\tOHC8@default\tEC14@default\tOHC10@default\tOxoC18@default\tC3ACP@default\tC12@default\tC6ACP@default\tOHC16@default\tC10ACP@default\tOxoC8@default\tOHC6@default\tOHC12@default\tEC12@default\tEC4@default\tEC8@default\tEC18@default\tPalmitoleic_acid@default\tEC6@default\tspecies_10@default\tCoA@default\tEC16@default\tspecies_12@default\tC10@default\tspecies_13@default\n" headerReactions="Model n\tBCCP_Holo\tACC_COOHbiotinCCP\tKASI_C16ACP\tBKR_OxoC18\tDH_OHC18\tENR_EC18\tDH_OHC14\tBKR_OxoC14\tBKR_OxoC16\tDH_OHC16\tDH_OHC4\tENR_EC4\tDH_OHC12\tKASI_C12ACP\tENR_EC14\tENR_EC10\tKASI_C8ACP\tKASI_C10ACP\tKASIII_C2CoA\tMCT_C3CoA\tENR_EC16\tKASI_C14ACP\tFACS_C18\tFACS_C12\tBKR_OxoC4\tDH_OHC8\tBKR_OxoC10\tDH_OHC10\tBKR_OxoC12\tAAT_C8ACP\tENR_EC12\tFACS_C10\tAAT_C18ACP\tKASI_C2ACP\tAAT_C10ACP\tAAT_C12ACP\tAAT_C14ACP\tAAT_C16ACP\tENR_EC8\tDH_OHC6\tENR_EC6\tKASI_C6ACP\tKASI_C4ACP\tBKR_OxoC6\n" headerCC="Model n\tJAAT_C16ACP_AAT_C16ACP\tJAAT_C16ACP_BCCP_Holo\tJAAT_C16ACP_ACC_COOHbiotinCCP\tJAAT_C16ACP_KASI_C16ACP\tJAAT_C16ACP_BKR_OxoC18\tJAAT_C16ACP_DH_OHC18\tJAAT_C16ACP_ENR_EC18\tJAAT_C16ACP_DH_OHC14\tJAAT_C16ACP_BKR_OxoC14\tJAAT_C16ACP_BKR_OxoC16\tJAAT_C16ACP_DH_OHC16\tJAAT_C16ACP_DH_OHC4\tJAAT_C16ACP_ENR_EC4\tJAAT_C16ACP_DH_OHC12\tJAAT_C16ACP_KASI_C12ACP\tJAAT_C16ACP_ENR_EC14\tJAAT_C16ACP_ENR_EC10\tJAAT_C16ACP_KASI_C8ACP\tJAAT_C16ACP_KASI_C10ACP\tJAAT_C16ACP_KASIII_C2CoA\tJSHMT_PGAM\tJAAT_C16ACP_MCT_C3CoA\tJAAT_C16ACP_ENR_EC16\tJAAT_C16ACP_KASI_C14ACP\tJAAT_C16ACP_FACS_C18\tJAAT_C16ACP_FACS_C12\tJAAT_C16ACP_BKR_OxoC4\tJAAT_C16ACP_FACS_C16\tJAAT_C16ACP_FACS_C14\tJAAT_C16ACP_BKR_OxoC8\tJAAT_C16ACP_DH_OHC8\tJAAT_C16ACP_BKR_OxoC10\tJAAT_C16ACP_DH_OHC10\tJAAT_C16ACP_BKR_OxoC12\tJAAT_C16ACP_FACS_C8\tJAAT_C16ACP_AAT_C8ACP\tJAAT_C16ACP_ENR_EC12\tJAAT_C16ACP_FACS_C10\tJAAT_C16ACP_AAT_C18ACP\tJAAT_C16ACP_KASI_C2ACP\tJAAT_C16ACP_AAT_C10ACP\tJAAT_C16ACP_AAT_C12ACP\tJAAT_C16ACP_AAT_C14ACP\tJAAT_C16ACP_ENR_EC8\tJAAT_C16ACP_DH_OHC6\tJAAT_C16ACP_ENR_EC6\tJAAT_C16ACP_KASI_C6ACP\tJAAT_C16ACP_KASI_C4ACP\tJAAT_C16ACP_BKR_OxoC6\tJAAT_C18ACP_AAT_C18ACP\tJAAT_C18ACP_BCCP_Holo\tJAAT_C18ACP_ACC_COOHbiotinCCP\tJAAT_C18ACP_KASI_C16ACP\tJAAT_C18ACP_BKR_OxoC18\tJAAT_C18ACP_DH_OHC18\tJAAT_C18ACP_ENR_EC18\tJAAT_C18ACP_DH_OHC14\tJAAT_C18ACP_BKR_OxoC14\tJAAT_C18ACP_BKR_OxoC16\tJAAT_C18ACP_DH_OHC16\tJAAT_C18ACP_DH_OHC4\tJAAT_C18ACP_ENR_EC4\tJAAT_C18ACP_DH_OHC12\tJAAT_C18ACP_KASI_C12ACP\tJAAT_C18ACP_ENR_EC14\tJAAT_C18ACP_ENR_EC10\tJAAT_C18ACP_KASI_C8ACP\tJAAT_C18ACP_KASI_C10ACP\tJAAT_C18ACP_KASIII_C2CoA\tJAAT_C18ACP_MCT_C3CoA\tJAAT_C18ACP_ENR_EC16\tJAAT_C18ACP_KASI_C14ACP\tJAAT_C18ACP_FACS_C18\tJAAT_C18ACP_FACS_C12\tJAAT_C18ACP_BKR_OxoC4\tJAAT_C18ACP_FACS_C16\tJAAT_C18ACP_FACS_C14\tJAAT_C18ACP_BKR_OxoC8\tJAAT_C18ACP_DH_OHC8\tJAAT_C18ACP_BKR_OxoC10\tJAAT_C18ACP_DH_OHC10\tJAAT_C18ACP_BKR_OxoC12\tJAAT_C18ACP_FACS_C8\tJAAT_C18ACP_AAT_C8ACP\tJAAT_C18ACP_ENR_EC12\tJAAT_C18ACP_FACS_C10\tJAAT_C18ACP_AAT_C18ACP\tJAAT_C18ACP_KASI_C2ACP\tJAAT_C18ACP_AAT_C10ACP\tJAAT_C18ACP_AAT_C12ACP\tJAAT_C18ACP_AAT_C14ACP\tJAAT_C18ACP_AAT_C16ACP\tJAAT_C18ACP_ENR_EC8\tJAAT_C18ACP_DH_OHC6\tJAAT_C18ACP_KASI_C6ACP\tJAAT_C18ACP_KASI_C4ACP\tJAAT_C18ACP_BKR_OxoC6\n" r.write(ListeParam) r1.write(headerSpecies) r2.write(headerReactions) r3.write(headerCC) #Open r4=open("BBBSpecies.txt","w") r5=open("CCCRates.txt","w") r6=open("DDDPlots.txt","w") #Naming SpeciesNames=['time', 'C18', 'CoA', 'Holo', 'COOHbiotinCCP', 'C3CoA', 'C3ACP', 'C16ACP', 'OxoC18', 'OHC18', 'EC18', 'C18ACP', 'OHC14', 'EC14', 'OxoC14', 'OxoC16', 'OHC16', 'EC16', 'OHC4', 'EC4', 'C4ACP', 'OHC12', 'EC12', 'C12ACP', 'C14ACP', 'EC10', 'C10ACP', 'C8ACP', 'OxoC10', 'OxoC12', 'C2ACP', 'C12', 'C12CoA', 'OxoC4', 'C16', 'C16CoA', 'C14', 'C14CoA', 'OxoC8', 'OHC8', 'EC8', 'OHC10', 'C10', 'C10CoA', 'OHC6', 'EC6', 'C6ACP', 'OxoC6'] RateNames=['time', 'SINKC18', 'CoABalance', 'BCCP_Holo', 'ACC_COOHbiotinCCP', 'KASI_C16ACP', 'BKR_OxoC18', 'DH_OHC18', 'ENR_EC18', 'DH_OHC14', 'BKR_OxoC14', 'BKR_OxoC16', 'DH_OHC16', 'DH_OHC4', 'ENR_EC4', 'DH_OHC12', 'KASI_C12ACP', 'ENR_EC14', 'ENR_EC10', 'KASI_C8ACP', 'KASI_C10ACP', 'KASIII_C2CoA', 'MCT_C3CoA', 'ENR_EC16', 'KASI_C14ACP', 'FACS_C12', 'BKR_OxoC4', 'FACS_C16', 'FACS_C14', 'BKR_OxoC8', 'DH_OHC8', 'BKR_OxoC10', 'DH_OHC10', 'BKR_OxoC12', 'ENR_EC12', 'FACS_C10', 'AAT_C18ACP', 'KASI_C2ACP', 'AAT_C10ACP', 'AAT_C12ACP', 'AAT_C14ACP', 'AAT_C16ACP', 'ENR_EC8', 'DH_OHC6', 'ENR_EC6', 'KASI_C6ACP', 'KASI_C4ACP', 'BKR_OxoC6'] n=1000 i=1 a= 1.1 aa=1.12 b=1.2 bb=1.22 while i<=n: # Random sampling of the parameters from the defined plausible distributions r.write(str(i)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC4_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC6_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC8_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC10_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC12_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC14_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC16_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(0.179,0.37) mod.BKR_OxoC18_KNADPH=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_KNADPH)+"\t") t=getNormalParamFromLogNormParam_median(2.6,0.5) mod.BKR_OxoC4_KOxoC4=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_KOxoC4)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC6_KOxoC6=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_KOxoC6)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC8_KOxoC8=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_KOxoC8)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC10_KOxoC10=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_KOxoC10)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC12_KOxoC12=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_KOxoC12)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC14_KOxoC14=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_KOxoC14)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC16_KOxoC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_KOxoC16)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC18_KOxoC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_KOxoC18)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC4_KOHC4=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_KOHC4)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC6_KOHC6=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_KOHC6)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC8_KOHC8=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_KOHC8)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC10_KOHC10=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_KOHC10)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC12_KOHC12=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_KOHC12)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC14_KOHC14=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_KOHC14)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC16_KOHC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_KOHC16)+"\t") t=getNormalParamFromLogNormParam_median(4.20,22.30) mod.BKR_OxoC18_KOHC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_KOHC18)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC4_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC6_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC8_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC10_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC12_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC14_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC16_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,2.69) mod.BKR_OxoC18_KNADP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_KNADP)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC4_KOHC4=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC4_KOHC4)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC6_KOHC6=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC6_KOHC6)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC8_KOHC8=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC8_KOHC8)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC10_KOHC10=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC10_KOHC10)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC12_KOHC12=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC12_KOHC12)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC14_KOHC14=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC14_KOHC14)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC16_KOHC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC16_KOHC16)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC18_KOHC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC18_KOHC18)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC4_KEC4=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC4_KEC4)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC6_KEC6=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC6_KEC6)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC8_KEC8=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC8_KEC8)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC10_KEC10=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC10_KEC10)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC12_KEC12=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC12_KEC12)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC14_KEC14=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC14_KEC14)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC16_KEC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC16_KEC16)+"\t") t=getNormalParamFromLogNormParam_median(0.88,3.69) mod.DH_OHC18_KEC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC18_KEC18)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC4_KEC4=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_KEC4)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC6_KEC6=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_KEC6)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC8_KEC8=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_KEC8)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC10_KEC10=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_KEC10)+"\t") t=getNormalParamFromLogNormParam_median(0.03,0.06) mod.ENR_EC12_KEC12=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_KEC12)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC14_KEC14=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_KEC14)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC16_KEC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_KEC16)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC18_KEC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_KEC18)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC4_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC6_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC8_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC10_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC12_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC14_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC16_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.05) mod.ENR_EC18_KNADH=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_KNADH)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC4_KC4ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_KC4ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC6_KC6ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_KC6ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC8_KC8ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_KC8ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC10_KC10ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_KC10ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC12_KC12ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_KC12ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC14_KC14ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_KC14ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC16_KC16ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_KC16ACP)+"\t") t=getNormalParamFromLogNormParam_median(4.35,21.67) mod.ENR_EC18_KC18ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_KC18ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC4_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC6_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC8_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC10_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC12_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC14_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC16_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.96,5.12) mod.ENR_EC18_KNAD=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_KNAD)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C2ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C2ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C4ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C4ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C6ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C6ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C8ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C8ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C10ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C10ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C12ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C12ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C14ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C14ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0115,0.022) mod.KASI_C16ACP_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.031,0.055) mod.KASI_C2ACP_KC2ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C2ACP_KC2ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.051,0.055) mod.KASI_C4ACP_KC4ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C4ACP_KC4ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.051,0.055) mod.KASI_C6ACP_KC6ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C6ACP_KC6ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.051,0.055) mod.KASI_C8ACP_KC8ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C8ACP_KC8ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.012,0.055) mod.KASI_C10ACP_KC10ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C10ACP_KC10ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0032,0.005) mod.KASI_C12ACP_KC12ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C12ACP_KC12ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.071,0.055) mod.KASI_C14ACP_KC14ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C14ACP_KC14ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.138,0.055) mod.KASI_C16ACP_KC16ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_KC16ACP)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.KASI_C16ACP_KOxoC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_KOxoC18)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.KASI_C16ACP_KACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_KACP)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.KASI_C16ACP_KCO2=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_KCO2)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C16ACP_KC16ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C16ACP_KC16ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C18ACP_KC18ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C18ACP_KC18ACP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C16ACP_KC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C16ACP_KC16)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C18ACP_KC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C18ACP_KC18)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C16ACP_KACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C16ACP_KACP)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C18ACP_KACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C18ACP_KACP)+"\t") t=getNormalParamFromLogNormParam_median(0.0024,0.015) mod.FACS_C16_KC16=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KC16)+"\t") t=getNormalParamFromLogNormParam_median(0.0089,0.015) mod.FACS_C18_KC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KC18)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C16_KCoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KCoA)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C18_KCoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KCoA)+"\t") t=getNormalParamFromLogNormParam_median(2.11,25.49) mod.FACS_C16_KATP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KATP)+"\t") t=getNormalParamFromLogNormParam_median(2.11,25.49) mod.FACS_C18_KATP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KATP)+"\t") t=getNormalParamFromLogNormParam_median(5.05,7.70) mod.FACS_C16_KAMP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KAMP)+"\t") t=getNormalParamFromLogNormParam_median(5.05,7.70) mod.FACS_C18_KAMP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KAMP)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C16_KPP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KPP)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C18_KPP=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KPP)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C16_KC16CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_KC16CoA)+"\t") t=getNormalParamFromLogNormParam_median(1.00,4.16) mod.FACS_C18_KC18CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_KC18CoA)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC4_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC6_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC8_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC10_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC12_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC14_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC16_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.BKR_OxoC18_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC4_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC4_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC6_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC6_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC8_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC8_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC10_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC10_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC12_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC12_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC14_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC14_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC16_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC16_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.DH_OHC18_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC18_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC4_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC6_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC8_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC10_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC12_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC14_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ENR_EC16_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.0256,0.14) mod.ENR_EC18_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.0256,0.14) mod.KASI_C2ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C2ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASI_C4ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C4ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASI_C6ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C6ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASI_C8ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C8ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.02,0.14) mod.KASI_C10ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C10ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASI_C12ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C12ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASI_C14ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C14ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.0033,0.14) mod.KASI_C16ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C16ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C16ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.81,3.69) mod.AAT_C18ACP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C18ACP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.000685,0.14) mod.FACS_C16_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.000284,0.14) mod.FACS_C18_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC4_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC4_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC6_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC6_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC8_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC8_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC10_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC10_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC12_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC12_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC14_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC14_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC16_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC16_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.BKR_OxoC18_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.BKR_OxoC18_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC4_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC4_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC6_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC6_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC8_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC8_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC10_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC10_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC12_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC12_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC14_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC14_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC16_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC16_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.DH_OHC18_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.DH_OHC18_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC4_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC4_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC6_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC6_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC8_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC8_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC10_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC10_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC12_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC12_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC14_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC14_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC16_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC16_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ENR_EC18_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ENR_EC18_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C2ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C2ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C4ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C4ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C6ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C6ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C8ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C8ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C10ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C10ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C12ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C12ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C14ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C14ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.KASI_C16ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASI_C16ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.AAT_C16ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C16ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.AAT_C18ACP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.AAT_C18ACP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.FACS_C16_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C16_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.FACS_C18_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.FACS_C18_Keq)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.KASIII_C2CoA_KACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASIII_C2CoA_KACP)+"\t") t=getNormalParamFromLogNormParam_median(0.04,89.95) mod.KASIII_C2CoA_KC2CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASIII_C2CoA_KC2CoA)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.KASIII_C2CoA_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.KASIII_C2CoA_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.081,0.03) mod.BCCP_Holo_KATP=random.lognormvariate(t[0],t[1]) r.write(str(mod.BCCP_Holo_KATP)+"\t") t=getNormalParamFromLogNormParam_median(16.2,2.5) mod.BCCP_Holo_KHCO3=random.lognormvariate(t[0],t[1]) r.write(str(mod.BCCP_Holo_KHCO3)+"\t") t=getNormalParamFromLogNormParam_median(20.17,82.34) mod.BCCP_Holo_KHolo=random.lognormvariate(t[0],t[1]) r.write(str(mod.BCCP_Holo_KHolo)+"\t") t=getNormalParamFromLogNormParam_median(78.6,3.2) mod.BCCP_Holo_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.BCCP_Holo_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ACC_COOHbiotinCCP_KC2CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_KC2CoA)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ACC_COOHbiotinCCP_KCOOHbiotinCCP=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_KCOOHbiotinCCP)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ACC_COOHbiotinCCP_KHolo=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_KHolo)+"\t") t=getNormalParamFromLogNormParam_median(0.1,4.29) mod.ACC_COOHbiotinCCP_KC3CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_KC3CoA)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.ACC_COOHbiotinCCP_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.ACC_COOHbiotinCCP_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.ACC_COOHbiotinCCP_Vmax)+"\t") t=getNormalParamFromLogNormParam_median(0.351,89.95) mod.MCT_C3CoA_KACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_KACP)+"\t") t=getNormalParamFromLogNormParam_median(0.25,89.95) mod.MCT_C3CoA_KCoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_KCoA)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.MCT_C3CoA_KC3CoA=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_KC3CoA)+"\t") t=getNormalParamFromLogNormParam_median(34.73,225.64) mod.MCT_C3CoA_KC3ACP=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_KC3ACP)+"\t") t=getNormalParamFromLogNormParam_maxci_median(5,500) mod.MCT_C3CoA_Keq=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_Keq)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.MCT_C3CoA_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.MCT_C3CoA_Vmax)+"\t") t=getNormalParamFromLogNormParam_maxci_median(8.0,800) mod.SINK18_KC18=random.lognormvariate(t[0],t[1]) r.write(str(mod.SINK18_KC18)+"\t") t=getNormalParamFromLogNormParam_maxci_median(0.047,4.7) mod.SINK18_Vmax=random.lognormvariate(t[0],t[1]) r.write(str(mod.SINK18_Vmax)+"\t") # Outputs #Verify the existence of a steady state : condition=True try: mod.doEigen() except: condition=False if condition: r.write("1\n") else: r.write("0\n") #Steady-state analysis #mod.mode_solver(1) #tell it to use fintslv #mod.mode_solver_fallback(0) #dont let it fallback to a different solver mod.__settings__['fintslv_tol'] = 0.1 #set non-default tolerance mod.__settings__['fintslv_step'] = 2 #set non-default step mod.doState() #this should go after setting the tol and step otherwise it tries to find the steady state then change these parameters r1.write(str(i)+"\t") r2.write(str(i)+"\t") # Equilibrium concentrations j=0 while j