Skip to content

Infeasible or Unbounded in reframed while solvable by cobrapy #13

@ccshao

Description

@ccshao

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!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions