reframed
reframed copied to clipboard
Infeasible or Unbounded in reframed while solvable by cobrapy
When I tried to perform FBA with both reframed and cobrapy with bio_mass as the objective function, the former gives "Infeasible or Unbounded" status and the latter works fine.
I used the human GEM model, v1.11.0; reframed is installed via downloading the repo and manual installation; cobrapy, v0.23.0. Solver is gurobi.
Here are the codes in reframed.
#- Works with test data.
import reframed
model = reframed.load_cbmodel('../tutorial/reframed/reframed-master/tests/data/e_coli_core.xml.gz')
sol = reframed.FBA(model)
# Objective: 0.8739215069684304
# Status: Optimal
model = reframed.load_cbmodel('../tutorial/reframed/reframed-master/tests/data/iML1515.xml.gz')
sol = reframed.FBA(model)
# Objective: 0.876997214426969
# Status: Optimal
model = reframed.load_cbmodel("../01_model/model/Human-GEM.xml")
sol = reframed.FBA(model)
# Objective: None
# Status: Infeasible or Unbounded
In cobrapy.
import cobra
model = cobra.io.read_sbml_model('../tutorial/reframed/reframed-master/tests/data/e_coli_core.xml.gz')
sol = model.optimize()
# <Solution 0.874 at 0x7feb5cb28310>
model = cobra.io.read_sbml_model('../tutorial/reframed/reframed-master/tests/data/iML1515.xml.gz')
sol = model.optimize()
<Solution 0.877 at 0x7feb5ae3b880>
model = cobra.io.read_sbml_model("../01_model/model/Human-GEM.xml")
sol = model.optimize()
<Solution 187.354 at 0x7feb52c5f790>
As shown above the solutions for two test model are same, but reframed failed to solve with human-GEM.
I would like to use the GIMME function in reframed on Human-GEM, so any suggestions on the how to make reframed with human model are really appreciated. Thanks!