compute.beta.from.B: Co-regulation coefficient computation

Description

Computes the matrix of co-regulation coefficients \(M_\beta\) from a vector of global co-regulation coefficients B

Usage

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

Arguments

B_fun

Numeric vector of global co-regulation coefficients

L_Phi_fun

List of regulation groups. Default is NULL.

Value

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.

Details

Default value L_Phi_fun = NULL is appropriate only if enzymes are supposed to be all co-regulated, i.e. no \(\beta\) 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.

Examples

Run this code

# NOT RUN {
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)



# }