SimBiology

フォロー

To perform a sensitivity analysis I am running *.sbproj coded models using 'biosimulate' in a for-loop and see for some parameter combinations the following warning: Warning: The right-hand side of the system of SimBiology ODEs results in complex numbers. The imaginary part of the result will be ignored. > In sbioodeflux_298d036b_ba03_4469_9210_25b32c7bbb2b (line 44) In sbiosimulate (line 140) In GSA (line 107) Can I trust the solutions? If not, how can I exclude the solutiuons which resulted in the warning above, i.e. is there a way to find programatically if a warning was produced? Update I added assignment of NaN's in case of fatal errors but still don't know how to handle the warnings. try ... [t,y] = sbiosimulate(m1); f1(i,j) = trapz(t,y(:,varNo)); % f1 AUC f2(i,j) = y(end,varNo); % f2 endPoint catch me disp( getReport( me, 'extended', 'hyperlinks', 'on' ) ) f1(i,j) = NaN; f2(i,j) = NaN; disp('Found an error-----------------------') end

It says "The rule expression is invalid. Check rule expression for completeness." I am trying to make an initial assignment, and my expression is this: [Static Glucagon Secr] = -[p]*[[Static Glucagon Secr] - max(((o2*([Hypoglycemic Threshold] - [Plasma Glu]))/([Plasma Ins Conc]+[Basal Plasma Ins Conc]) +[Basal Glucagon Secr]),0)] Eventually, my goal is to create sort of an "if then" situation with the model. So that if G>Gb, this expression is true, otherwise it is not. I am using initial assignments, as well as two events for this. How can I fix this warning?