readTSVmod
. This function returns
an object of the class modelorg
. Read csv files (example files included):
mpath <- system.file(package = "sybil", "extdata")
model <- readTSVmod(prefix = "Ec_core",
fpath = mpath, quote = "\"")
Perform flux balance analysis (FBA):
ec_f <- simpleFBA(model)
Perform single gene deletion analysis:
ec_g <- oneGeneDel(model)
Plot the values of the objective function after optimization in a
histogramm:
hist(ec_g)
Perform flux variability analysis:
ec_v <- fluxVar(model)
Plot the result:
plot(ec_v)
Schellenberger, J., Park, J. O., Conrad, T. C., and Palsson, B. Ø., (2010) BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions. BMC Bioinformatics 11, 213.
The openCOBRA project
Becker, S. A., Feist, A. M., Mo, M. L., Hannum, G., Palsson, B. Ø. and Herrgard, M. J. (2007) Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox. Nat Protoc 2, 727--738.
Schellenberger, J., Que, R., Fleming, R. M. T., Thiele, I., Orth, J. D., Feist, A. M., Zielinski, D. C., Bordbar, A., Lewis, N. E., Rahmanian, S., Kang, J., Hyduke, D. R. and Palsson, B. Ø. (2011) Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc 6, 1290--1307.
readSBMLmod
to read
metabolic models written in SBML language.data(Ec_core)
Ec_ofd <- oneGeneDel(Ec_core)
hist(Ec_ofd)
Run the code above in your browser using DataLab