SimEvolEnzCons (version 2.0.0)

RNV.ranking.order.factor: Name and value of RNV-ranking-order factor

Description

Gives the name and values of the RNV-ranking-order factor

Usage

RNV.ranking.order.factor(A_fun,correl_fun,E_ini_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

E_ini_fun

Numeric vector of initial concentrations.

B_fun

Numeric vector of global co-regulation coefficients

Value

Invisible list of 2 elements:

• $value: numeric vector (length n) of the values of the RNV-ranking-order factor for each enzymes

• $name: character sting indicating the name of the RNV-ranking-order variable

Details

Factors governing ranking order of RNV

RNV-ranking-order factor depends on constraint. It would be:

• in case of independence (correl_fun="SC"): activities A power 1/3;

• in case of competition only (correl_fun="Comp"): activities A power 1/4;

• in case of regulation only (correl_fun="RegPos" or "RegNeg"): absolute value of inverse of global co-regulation coefficients |1/B|;

• in case of competition and regulation (correl_fun="CRPos" or "CRNeg"): absolute value of inverse of global co-regulation coefficients minus initial relative concentrations |1/B - e0|.

Examples

# NOT RUN {
A <- c(1,20,30)
E0 <- c(30,30,30)

#Independence
rank_var <- RNV.ranking.order.factor(A,"SC",E0)
all.equal(A^(1/3),rank_var$value) #TRUE

#Positive regulation
B <- 1/c(0.2,0.5,0.3)
rank_var <- RNV.ranking.order.factor(A,"RegPos",E0,B)
all.equal(1/B,rank_var$value) #TRUE


# }