predict_th: Prediction of theoretical equilibrium

Description Usage Arguments Details Value Special results See Also Examples

View source: R/predict_th.R

Description

Gives the theoretical equilibrium for relative concentrations

Usage

1
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

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.

Value

List of two elements:

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.

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
#### 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

SimEvolEnzCons documentation built on Oct. 29, 2021, 1:07 a.m.