SimEvolEnzCons (version 2.0.0)

predict_th: Prediction of theoretical equilibrium

Description

Gives the theoretical equilibrium for relative concentrations

Usage

predict_th(A_fun,correl_fun,B_fun=NULL)

Arguments

A_fun

Numeric vector of activities

correl_fun

Character string indicating the abbreviation of the constraint applied on the system

B_fun

Numeric vector of global co-regulation coefficients

Value

List of two elements:

• $pred_e: numeric vector of relative concentrations at theoretical equilibrium. Same length as A_fun

• $pred_r: numeric vector of response coefficients at theoretical equilibrium. Same length as A_fun

Special results

In case of negative regulation (correl_fun = "RegNeg" or "CRNeg"), relative concentrations would be negative.

In case of competition plus regulation (correl_fun = "CRPos" or "CRNeg"), response coefficients is not defined and $pred_r returns NaN.

Details

Gives values at theoretical equilibrium for relative concentrations and response coefficients. This equilibrium corresponds to null derivative for relative concentrations, without conditions on flux.

When there are regulation groups, preferably use predict_grp.

See Also

Use function activities to compute enzyme activities.

Use function is.correl.authorized to see allowed constraints for correl_fun.

Use function predict_grp to predict equilibria when there are regulation groups.

Examples

# NOT RUN {
#### For independancy "SC" or competition "Comp"
A <- c(1,10,30)

eq_th <- predict_th(A,"SC")

eq_th$pred_e
eq_th$pred_r

###### In presence of regulation
A <- c(1,10,30)
beta <- matrix(c(1,10,5,0.1,1,0.5,0.2,2,1),nrow=3)
B <- apply(beta,1,sumbis) 

eq_th <- predict_th(A,"CRPos",B)

eq_th$pred_e
eq_th$pred_r


# }