-
Notifications
You must be signed in to change notification settings - Fork 70
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
"CupSodaSimulator doesn't support model with expression" error. #515
Comments
cupSODA generates the ODEs from the the stoichiometry matrix from the reactions, so expressions, so it does not support the use of expressions. |
@JamesPino I changed the script to print the expressions imported with the model and I get that So I think that model has expressions. |
It seems that the way the sbml is written and how it is converted to pysb (through bngl), the only way to do it safely is to use expressions. For example, for a reaction A -> B with a rate constant k, the ratelaw in SBML is k*[A]. In BNG, we need to pull out the k from the ratelaw, so Atomizer tries to divide out the concentration [A] from it. Sometimes there’s a danger of divide-by-zero, so it makes a function like if(A>0, ratelaw/[A], 0). My guess is that’s the problem they’re running into. After converting the model is pysb you can export to a file, then edit out the expressions.
Sadly, it looks like this is the only way since cupSODA doesn't handle expressions, and SBML encodes the model with ratelaws. I can't think of another solution. The Tyson model you are using should be the same that is in pysb.examples.tyson_oscillator |
@JamesPino Okay Thanks. After exporting in a new file, like did you do, how can I run a cupSODA simulation? Because it gives me the same error... |
So you'll have to modify that file to remove expressions. I am not sure your experience building models, but basically you have to add reactions corresponding to those expressions. Is your project to start from SBML or are you just trying to use that model? I honestly don't know if it is easier to create the model from the information on biomodels from scratch or modify the exported file. Either way, the exported file should give you the monomer species and parameters, then its just creating the reactions from there. If your goal is just to use that model for analysis, then you might be able to make the PySB implementation of it work for you. There are still a lot of issues to be worked out to convert SBML models to rule based models in a robust way. |
@JamesPino Unfortunatelly I'm not exeperinced in building model from scratch, so I don't know how to fix this problem. In these days I tried to learn something abou, but still I don't understand how to build reactions from expressions. Do you have any idea ? |
Hi, I've a problem simulating this model BIOMD0000000005 using the cupsoda simulator. I wrote this simple python script
and the error is
Some general informations:
Where I'm wrong? I don't understand!
Thansk to everyone!!!
The text was updated successfully, but these errors were encountered: