reframed icon indicating copy to clipboard operation
reframed copied to clipboard

Infeasible or Unbounded in reframed while solvable by cobrapy

Open ccshao opened this issue 3 years ago • 0 comments

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!

ccshao avatar May 12 '22 09:05 ccshao