# NOT RUN {
E <- c(30,30,30)
beta <- matrix(c(1,10,5,0.1,1,0.5,0.2,2,1),nrow=3)
correl <- "RegPos"
mu <- 1 #canonical size of mutation
i <- 3 #enzyme directly targeted by mutation
mut.E.direct(E,i,mu,correl,beta)
# }
Run the code above in your browser using DataLab