sybil (version 2.2.0)

changeBounds: Change Variable Bounds in a Metabolic Network

Description

The function changes the upper and/or lower bounds of a given metabolic network model to new values.

Usage

changeBounds(model, react, lb = NULL, ub = NULL)

Arguments

model

An object of class '>modelorg.

react

An object of class '>reactId, character or integer. Specifies the fluxes (variables) for which to change the upper and/or lower bounds.

lb

Numeric vector giving the lower bounds for the fluxes mentioned in react. If missing, lower bounds are set to zero. If lb has a length of 1, the value of lb will be used for all reactions in react.

ub

Numeric vector giving the upper bounds for the fluxes mentioned in react. If missing, upper bounds are set to zero. If ub has a length of 1, the value of ub will be used for all reactions in react.

Value

Returns the given model (an object of the same class as the argument lpmodel) containing the new objective function.

Details

The argument react will be evaluated by the function checkReactId.

See Also

checkReactId

Examples

# NOT RUN {
  ## change the E.coli core model to lactate input:
  data(Ec_core)
  Ec_new <- changeBounds(Ec_core,
                         c("EX_glc", "EX_lac"),
                         lb = c(0, -20), ub = 1000)
# }