MATLAB Answers

Invalid assignment syntax or parse error for event function (Simbiology)

48 views (last 30 days)
BK V
BK V on 28 Sep 2021 at 15:04
Commented: BK V on 5 Oct 2021 at 19:17
I am adding an event to a QSP model to impose a function to be a small positive number (say epsilon) if its value goes below epsilon. Snippet of the code is:
f_HPD1 = '((syn_CT.PD1_PDL1+syn_CT.PD1_PDL2)/PD1_50)';
addevent(model, [f_HPD1 '< 1e-10'], [f_HPD1 '= 1e-10']);
Where,
syn_CT is compartment name,
PD1_PDL1, PD1_PDL2 are species in the compartment syn_CT
PD1_50 is model parameter
But I get this error message:
Error using SimBiology.internal.simulate
--> Error reported from Expression Validation:
Invalid assignment syntax or parse error for event function '((syn_CT.PD1_PDL1+syn_CT.PD1_PDL2)/PD1_50) = 1e-10' in event with
trigger '((syn_CT.PD1_PDL1+syn_CT.PD1_PDL2)/PD1_50)< 1e-10'. Event functions must be valid MATLAB expressions and cannot end in
semicolons, commas, comments ('%' and optional text), or line continuations ('...' and optional text).
Error in sbiosimulate (line 140)
simResultsCell = SimBiology.internal.simulate(mobj, cs, variants, doses);
Error in initial_conditions (line 84)
simData = sbiosimulate(model,variant);
Please help me resolve this, any suggestion is really appreciated.
Thank you

Accepted Answer

Arthur Goldsipe
Arthur Goldsipe on 28 Sep 2021 at 16:04
Hi,
Events can only assign a value to a component, not to an expression. So your assignment ((syn_CT.PD1_PDL1+syn_CT.PD1_PDL2)/PD1_50) = 1e-10 is invalid because the left-hand side is not a single component name. A valid assignment would look more like syn_CT.PD1_PDL1 = 1e-10. Note that the component on the left-hand side also has to have its Constant property set to false.
Looking at the bigger picture of what you're trying to do, I'm guessing what you might want to do instead is to create a parameter (let's call it HPD1), set its Constant property set to false, and update its value with a repeated assignment rule of the form HPD1 = max(1e-10, ((syn_CT.PD1_PDL1+syn_CT.PD1_PDL2)/PD1_50)).
-Arthur
  9 Comments
BK V
BK V on 5 Oct 2021 at 19:17
Thank you, I did check the units are consistent so was not sure why do these warnings appear. I cannot turn off the unit conversion since some of the data that we input are in different units and so unit conversion is needed.

Sign in to comment.

More Answers (0)

Communities

More Answers in the  SimBiology Community

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!