Running Simulations

Once a model is successfully loaded we can run a time course simulation. To run a simulation we use the simulate() method:

result = rr.simulate()

The output will be in a Python numpy array. The first column will contain time points and the remaining columns will include all the floating species amounts/concentrations. In the simulate method we didn’t specify how long to run the simulation or how many points to generate. By default the starting time is set to zero, ending time to 10 and the number of points to 51. There are two ways to set these values. The easiest way is to change the positional arguments in simulate() in the following manner:

result = rr.simulate (0, 10, 100)

This will set the starting time to zero, the ending time to 10 and generate 100 points. This means that the result will be out in time intervals of 1/99.

The simulate method also accepts the steps keyword argument instead of points:

result = rr.simulate(0, 10, steps=99)

For more details of the simulate method see simulate(). The following table summarizes the various options.

Option

Description

start

Starting time for simulation

end

Ending time for simulation. Setting ‘end’ will automatically change ‘duration’

points

Number of rows to include in the output matrix

selections

(Optional) A list of variables to include in the output, e.g. ['time','A'] for a model with species A. More below.

steps

(Optional keyword argument) Number of steps at which the output is sampled where the samples are evenly spaced. Steps = points-1. Steps and points may not both be specified.

One important point to note about simulate(): When simulate() is run, the concentration of the floating species will naturally change. If simulate() is called a second time, the simulation will start the simulation from the previous simulated values. This can be used to easily follow on simulations. However there will be times when we wish to run the same simulation again but perhaps with slightly different parameters values. For this we must reset the model to its initial conditions. To do that we run the command reset():

rr.reset()

RoadRunner also has two other reset methods: resetAll() and resetToOrigin(). These are typically only used in advanced scenarios. This is because RoadRunner maintains its own copy of the “initial” value of every quantity, which can be set via rr_instance.setValue('init(quantity)',123). When resetAll() is called, it resets the quantity to RoadRunner’s internal copy, not the SBML-defined value. resetToOrigin() completely reverts everything in the model back to the SBML-specified values, whereas resetAll() preserves any changes you have made to initial values.