compute.beta.from.B: Co-regulation coefficient computation
In SimEvolEnzCons: Simulation of Evolution of Enzyme Under Constraints

Description Usage Arguments Details Value See Also Examples

Computes the matrix of co-regulation coefficients M_β from a vector of global co-regulation coefficients B

1	compute.beta.from.B(B_fun,L_Phi_fun=NULL)

`B_fun`	Numeric vector of global co-regulation coefficients
`L_Phi_fun`	List of regulation groups. Default is `NULL`.

Default value L_Phi_fun = NULL is appropriate only if enzymes are supposed to be all co-regulated, i.e. no β value is null. Return the same result as as a list of one element such as L_Phi_fun = list(1:n), where n is the number of enzymes, which is also the length of B_fun.

L_Phi_fun must be a list of p elements (the number of regulation groups) containing numbers between 1 and n (number of enzymes), where each list element contains the numbers of the enzyme in the group. Each enzyme only occurs in one group. See function class_group to have an idea of L_Phi_fun structure. For compute.beta.from.B, independent enzymes can be not contained in L_Phi_fun, and thus, L_Phi_fun may be smaller than B_fun.

Numeric matrix n*n of the co-regulation coefficients, where n is the length of B_fun.

If beta_fun is NULL, compute.beta.from.B returns NULL.

Use function is.B.accurate to verify B_fun conformity.

See function class_group to know format of L_Phi_fun.

B <- 1/c(0.5,0.2,0.3)
correl <- "RegPos"

is.B.accurate(B,3,correl)

beta <- compute.beta.from.B(B)

#Seven enzymes and three groups
n <- 7
p <- 3
B <- c(1.1824,  3.695, -8.593023,  1.3, 13,  6.5, 1)
L_Phi <- list(1:3) #firt three enzymes in first group
L_Phi[[2]] <- 4:6
L_Phi[[3]] <- 7 #last enzyme independent

beta <- compute.beta.from.B(B,L_Phi)