cfFBA: cycle free flux

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/cfFBA.R

Description

finds a cycle free flux distribution given a flux distribution and a network model.

Usage

 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)

Arguments

model

An object of class modelorg.

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 "min" or "max".
Default: SYBIL_SETTINGS("OPT_DIRECTION").

solver

Single character value. The solver to use. See SYBIL_SETTINGS for possible values.
Default: SYBIL_SETTINGS("SOLVER").

method

Single character value. The optimization algorithm to use. Possible values depend on the setting in solver. See SYBIL_SETTINGS for possible values.
Default: LP_METHOD(SYBIL_SETTINGS).

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 TRUE, the function returns an object of class optsol_optimizeProb, otherwise a list.

Value

return a list that contains the status of the solution , the objective value and the fluxes that are free of cycles.

Author(s)

Abdelmoneim Amer Desouki

References

Desouki, Abdelmoneim Amer, et al. "CycleFreeFlux: efficient removal of thermodynamically infeasible loops from flux distributions." Bioinformatics 31.13 (2015): 2159-2165.

See Also

modelorg

Examples

 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"

sybilcycleFreeFlux documentation built on Dec. 16, 2019, 1:19 a.m.