BacArena (version 1.0.1)

constrain: Function for constraining the models based on metabolite concentrations

Description

The generic function constrain changes the constraints of the model representation of an organism.

Usage

constrain(object, reacts, lb, dryweight, time, scale)
"constrain"(object, reacts, lb, dryweight, time, scale)

Arguments

object

An object of class Organisms.

reacts

A character vector giving the names of reactions which should be constrained.

lb

A numeric vector giving the constraint values of lower bounds (e.g. avaible metabolite concentrations

dryweight

A number giving the current dryweight of the organism.

time

A number giving the time intervals for each simulation step.

scale

A numeric defining the scaling (units for linear programming has to be in certain range)

Value

Returns the lower bounds, which carry the constraints and names of relevant reactions.

Details

The constraints are calculated according to the flux definition as mmol/(gDW*hr) with the parameters dryweight and time.

See Also

Organism-class

Examples

data(Ec_core, envir = environment()) #get Escherichia coli core metabolic model
org <- Organism(Ec_core,deathrate=0.05,
           growthlimit=0.05,growtype="exponential") #initialize an organism
lobnds <- constrain(org,org@medium,org@lbnd[org@medium],1,1)