caMonoBiExpFromIG {CalciOMatic} | R Documentation |
The caMonoBiExpFromIG
function returns a vector of
intracellular calcium concentration ([Ca^2+]). A dCa
jump occurs at
tOn
, followed by either a monoexponential or a biexponential
return to the baseline value Ca0
, depending on the fields of
the input list ig
. The caMonoBiExpFromIG
function is a
low-level function of the ratioFitFromCa
and directFit
functions.
caMonoBiExpFromIG(t = 1, tOn = 1, ig = NULL)
t |
a vector of time values at which [Ca^2+] is computed (in s) |
tOn |
the time of the [Ca^2+] jump (in s) |
ig |
an object of class "initial_guess" , giving the
parameters of the decay. This object is a list with the following
numerical fields: "log_Ca0" , "log_dCa" ,
"log_tau" in the case of a monoexponential decay. In the case
of a biexponential decay, the two following fields are also
included: "mu" and "log_dtau" |
A vector containing the [Ca^2+] values. The vector has the two following attributes:
Time |
a copy of argument t |
tOn |
a copy of argument tOn |
Sebastien Joucla sebastien.joucla@parisdescartes.fr
## Parameters of the calcium transient Ca0 <- 0.10 dCa <- 0.25 tau <- 1 ## List of initial guesses ig_mono <- list(log_Ca0 = log(Ca0), log_dCa = log(dCa), log_tau = log(tau)) class(ig_mono) <- "initial_guess" ## Build the calcium transient Ca <- caMonoBiExpFromIG(t = seq(0,10,0.1), tOn = 2, ig = ig_mono) ## Plot the calcium transient vs. time plot(attr(Ca,"Time"), Ca, type="l") ## Add a vertical dashed line at tOn abline(v = attr(Ca,"tOn"), lty = 2)