# NOT RUN {
## One group
beta <- matrix(c(1,10,5,0.1,1,0.5,0.2,2,1),nrow=3)
L_Phi <- class_group(beta)
#enzyme 2 is in group 1
search_group(2,L_Phi)
## Two groups
n <- 3
beta <- diag(1,n)
beta[1,2] <- -0.32
beta[2,1] <- 1/beta[1,2]
L_Phi <- class_group(beta)
search_group(2,L_Phi) #gives 1
search_group(3,L_Phi) #gives 2
# }
Run the code above in your browser using DataLab