Result selections

The timecourseSelections on the roadrunner object define what is part of the simulation.
By default only time and concentrations are part of the timecourseSelections. To add additional selections just add the required sbml ids to the selection.

To get sets of ids you can use the getXXXIds() which you add to the timecourseSelections. An example is given below.

Same principle applies for the steadystateSelections.



In [1]:

    
%matplotlib inline
from __future__ import print_function
import tellurium as te



In [2]:

    
r = te.loada("""
r1: S1 -> S2; k1*S1-k2*S2;
S1 = 10.0; S2 = 0.0;
k1 = 0.2; k2 = 0.1;
""")



In [3]:

    
# default selections in results
print(r.timeCourseSelections)
s = r.simulate(start=0, end=10, steps=10)
print(s)









    



['time', '[S1]', '[S2]']
    time,    [S1],    [S2]
 [[    0,      10,       0],
  [    1, 8.27215, 1.72785],
  [    2, 6.99211, 3.00789],
  [    3,  6.0438,  3.9562],
  [    4, 5.34126, 4.65874],
  [    5, 4.82082, 5.17918],
  [    6, 4.43527, 5.56473],
  [    7, 4.14962, 5.85038],
  [    8,   3.938,   6.062],
  [    9, 3.78127, 6.21873],
  [   10, 3.66519, 6.33481]]



In [4]:

    
# add additional ids to the timecourse selections
r.timeCourseSelections += r.getGlobalParameterIds() + r.getCompartmentIds() + r.getReactionIds()
print(r.timeCourseSelections)









    



['time', '[S1]', '[S2]', 'k1', 'k2', 'default_compartment', 'r1']



In [5]:

    
# reset to initial concentrations
r.reset()
s2 = r.simulate(start=0, end=10, steps=10)
print(s2)









    



    time,    [S1],    [S2],  k1,  k2, default_compartment,        r1
 [[    0,      10,       0, 0.2, 0.1,                   1,         2],
  [    1, 8.27215, 1.72785, 0.2, 0.1,                   1,   1.48165],
  [    2, 6.99211, 3.00789, 0.2, 0.1,                   1,   1.09763],
  [    3,  6.0438,  3.9562, 0.2, 0.1,                   1,  0.813141],
  [    4, 5.34126, 4.65874, 0.2, 0.1,                   1,  0.602379],
  [    5, 4.82082, 5.17918, 0.2, 0.1,                   1,  0.446245],
  [    6, 4.43527, 5.56473, 0.2, 0.1,                   1,  0.330581],
  [    7, 4.14962, 5.85038, 0.2, 0.1,                   1,  0.244885],
  [    8,   3.938,   6.062, 0.2, 0.1,                   1,  0.181399],
  [    9, 3.78127, 6.21873, 0.2, 0.1,                   1,  0.134382],
  [   10, 3.66519, 6.33481, 0.2, 0.1,                   1, 0.0995582]]



In [ ]: