Description Usage Arguments Details Value Special results See Also Examples
Gives the theoretical equilibrium for relative concentrations
1 | predict_th(A_fun,correl_fun,B_fun=NULL)
|
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 |
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
.
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
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
.
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.
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
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.