Description Usage Arguments Value Author(s) References See Also Examples
finds a cycle free flux distribution given a flux distribution and a network model.
1 2 3 4 5 6 7 8 9 10 11 12 | cfFBA(model, wtflux, objVal = NA #min objval
,fixExchRxn=TRUE
,excReactPos=NULL
,lpdir = SYBIL_SETTINGS("OPT_DIRECTION")
,solver = SYBIL_SETTINGS("SOLVER")
,method = SYBIL_SETTINGS("METHOD")
,tol=SYBIL_SETTINGS("TOLERANCE")
,solverParm=NA
,verboseMode = 2
,safeBounds=FALSE #
######### ADDED BY GABRIEL ###############
,retOptSol = TRUE)
|
model |
An object of class |
wtflux |
initial flux distribution that may contain loops. |
objVal |
value of the objective function. |
fixExchRxn |
a logical value default is true, which indicates that the exchange reactions should be fixed or not. |
excReactPos |
list of the exchange reactions that can be sent to avoid recalculation |
lpdir |
Character value, direction of optimisation. Can be set to |
solver |
Single character value. The solver to use. See
|
method |
Single character value. The optimization algorithm to use. Possible values
depend on the setting in |
solverParm |
additional parameters to the solver |
tol |
tolerance of resulting solution |
safeBounds |
can be used to deal with false infeasibility state returned by solvers |
verboseMode |
level of displaying messages. |
retOptSol |
If set to |
return a list that contains the status of the solution , the objective value and the fluxes that are free of cycles.
Abdelmoneim Amer Desouki
Desouki, Abdelmoneim Amer, et al. "CycleFreeFlux: efficient removal of thermodynamically infeasible loops from flux distributions." Bioinformatics 31.13 (2015): 2159-2165.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
data(iAF1260)
model=iAF1260
fba=optimizeProb(model)
cfopt=cfFBA(model,wtflux=getFluxDist(fba),objVal=lp_obj(fba), retOptSol=FALSE)
llflx=cfopt$fluxes
flx=getFluxDist(fba)
## Not run:
layout(matrix(c(1,2,3,1,2,3), 2, 3, byrow = TRUE))
hist(log10(abs(flx[abs(llflx-flx)>1e-3])),main="loop flux",col="lightblue")
hist(log10(abs(llflx[abs(llflx-flx)>1e-3])),main="after removing loops",col="orange")
hist(log10(abs(flx[abs(llflx-flx)<1e-3])),main="fluxes not in loops",col="lightgreen")
## End(Not run)
## The function is currently defined as
"cfFBA"
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.