# NOT RUN {
data(Ec_core)
# retrieve all exchange reactions
ex <- findExchReact(Ec_core)
# perform flux balance analysis
opt <- optimizeProb(Ec_core, algorithm = "fba")
# get flux distribution of all exchange reactions
fd <- getFluxDist(opt, ex)
# group exchange reactions
getNetFlux(fd)
# }
Run the code above in your browser using DataLab