Simple Parmameter Scanning Functions

parameterScanning.simpleTimeCourseScan(r, parameter, variable, lowRange, highRange, numberOfScans, timeEnd=10, numberOfPoints=100, formatStr='{:10.6f}', legendLoc='upper right')[source]

Run a time course simulation at different parameter values, observe a single variable

Parameters:

  • r (reference) – Roadrunner instance

  • parameter (string) – The name of the parameter to change

  • variable (string) – The name of the variable to record during the scan

  • lowRange (float) – The starting value for the parameter

  • highRange (float) – The final value for the parameter

  • numberOfScans (integer) – The number of values of the parameter to try

  • timeEnd (float) – Optional: Simulate a time course up to this time

  • numberOfPoints – (integer): Optional: Generate this number of points for each time course

  • formatStr (string) – Optional: The format string for values listed in the plot legend

Returns:

numpy array – first column being time remaining columns correspond to each parameter scan.

Example

import tellurium as te
import teUtils as tu

r = te.loada('''
J1: $Xo -> S1; k10*Xo - k11*S1;
J2: S1 -> S2;  k20*S1 - k21*S2;
J3: S2 -> $X1; k30*S2 - k31*X1;

k10 = 1.21;  k11 = 0.69
k20 = 1.03;  k21 = 0.13
k30 = 1.89;  k31 = 0.10
Xo = 6.00
X1 = 0
S1 = 0; S2 = 0;
''')

tu.parameterScanning.simpleTimeCourseScan(r, 'k20', 'S1',
        3, 12, 7, timeEnd=6, numberOfPoints=200, formatStr='{:4.1f}')