Encoding conditional assignments in SimBiology

Is it possible to define conditional assignments in SimBiology, i.e.

if a > 0

param1 = a + b

else

param1 = a + c

end

I am aware of events but it is insuffient as it lacks the 'else' expression.

Answer 1

[Ayush Singh]

Hey emjey,

In SimBiology, conditional assignments like the one you described are not directly supported in the same way as typical programming languages. However, you can achieve similar functionality using repeated assignment rules. Since events do not support an `else` clause directly, you will need to set up your model in a way that effectively mimics the conditional logic.

You can achieve this using repeated assignment rules with logical conditions. SimBiology allows you to define rules that continuously update the value of a parameter based on other parameters or species in the model.

You can follow below steps to achieve similar functionality

1. Define Parameters:

- Define the parameters 'a', 'b', 'c', and 'param1' in your model.

2. Create a Rule:

- Add a repeated assignment rule to your model that uses logical conditions.

Example:

Suppose you have parameters 'a', 'b', 'c', and 'param1'. You can define a repeated assignment rule for 'param1' using logical conditions as follows:

% Create a SimBiology model
model = sbiomodel('conditionalAssignmentModel');

% Add parameters
a = addparameter(model, 'a', 1);
b = addparameter(model, 'b', 2);
c = addparameter(model, 'c', 3);
param1 = addparameter(model, 'param1', 0);

% Add a repeated assignment rule
% The rule will use logical conditions to mimic the if-else structure
addrule(model, 'param1 = (a > 0) * (a + b) + (a <= 0) * (a + c)', 'repeatedAssignment');

% Simulate the model (optional)
cs = getconfigset(model);
cs.StopTime = 10;
simData = sbiosimulate(model);

% Plot the results (optional)
sbioplot(simData);

Below is the explanation for above

The rule `param1 = (a > 0) * (a + b) + (a <= 0) * (a + c)` works as follows:

• `(a > 0)` evaluates to `1` if `a` is greater than `0`, and `0` otherwise.
• `(a <= 0)` evaluates to `1` if `a` is less than or equal to `0`, and `0` otherwise.
• When `a > 0`, `(a > 0) * (a + b)` will be `a + b` and `(a <= 0) * (a + c)` will be `0`.
• When `a <= 0`, `(a > 0) * (a + b)` will be `0` and `(a <= 0) * (a + c)` will be `a + c`.

Thus, `param1` will be assigned the value of `a + b` if `a > 0`, and `a + c` if `a <= 0`.

Hope it helps!

Answer 2

[Arthur Goldsipe]

Hi emjey,

I see this problem solved in one of two ways:

1) Take advantage of true/false getting converted to 1/0 in mathematical operations. This is the approach suggested in the other answer. Another altenative for the rule expression is "param1 = a + b*(a>0) + c*(a<=0)". One disadvantage of this approach is that the rule expressions can easily get complicated and hard to debug, especially for more complicated "if, else if, else" expressions.

2) You can put any complicated logic, include if/else in a helper function that you call from the repeated assignment. That way, you can use if/else to write the code the most natural way. You ca find an example of that approach here. The main downside to this is that you must share this function file with anyone else who wants to use the model, since it isn't automatically bundled with the model or project file.