range_delta: Bounds of _delta_i_
In SimEvolEnzCons: Simulation of Evolution of Enzyme Under Constraints

Description Usage Arguments Details Value See Also Examples

Computes the bounds of the actual mutation effect δ_i such as all mutant concentrations are between 0 and total concentration, for a mutation targeting enzyme i

1	range_delta(E_res,alpha_fun,i_fun,tol_fun=0.0001)

`E_res`	Numeric vector of resident enzyme concentrations
`alpha_fun`	Numeric matrix of redistribution coefficients
`i_fun`	Integer number indicating the enzyme targeted by the mutation
`tol_fun`	Numeric and positive value. Accuracy for delta bounds. Default is `0.0001`

This function range.delta computes the bounds of δ_i such as all mutant concentrations are between 0 and total concentration Etot, for a mutation targeting enzyme i_fun. Mutant concentrations are equal to resident concentrations plus α_ij * δ_i (see function \code{\link{mut.E.indirect}}). For any enzyme j, mutant value is E_j^r + α_ij * δ_i.

The inferior (resp. superior) bound of δ_i corresponds to minimal (resp. maximal) value of δ_i such as all mutant concentrations are superior or equal to 0 and inferior or equal to Etot, with at least one mutant concentration equal to 0 or Etot.

tol_fun is the accuracy (or allowed tolerance) for δ bounds. It allows to avoid asymptote problem when computing the RNV.

Numeric vector of the inferior and the superior bounds of actual mutation effect δ_i

See function alpha_ij to compute matrix of redistribution coefficients alpha_fun.

beta <- matrix(c(1,10,5,0.1,1,0.5,0.2,2,1),nrow=3)
Er <- c(30,30,30)
correl <- "CRPos"
alpha <- alpha_ij(Er,correl,beta)

#mutant enzyme
i <- 1

range_delta(Er,alpha,i)