| mkFunction4RatioFit {CalciOMatic} | R Documentation |
The function mkFunction4RatioFit returns a function
prediciting the exponential time course of an intracellular
concentration. The function, the arguments of which depend on the
value of type, can be used within the RSS function of a
ratiometric fit.
mkFunction4RatioFit(type = "mono")
type |
a character string (either "mono" or "bi"),
specifying which type of exponential decay should be considered |
A function that has five or seven arguments, depending on the
value of type. If type is set to "mono", the five
arguments are the following: t, tOn, log_Ca0,
log_dCa and log_tau. If type is set to
"bi", the output function has two more arguments: mu and
log_dtau. For details about the meaning of each argument, see
the caMonoExp or caBiExp functions
Sebastien Joucla sebastien.joucla@parisdescartes.fr
caMonoExp, caBiExp,
ratioFitFromCa, ratioFitFromDf
## Time parameters
tOn <- 1
Time <- seq(0,30,0.1)
Ca0 <- 0.10
dCa <- 0.25
tau <- 1.5
mu <- 2
dtau <- 10
## Define a calcium biexponential decay with the
## mkCa_MonoBiExp_4_RatioFit function
Ca_biexp_fct <- mkFunction4RatioFit(type = "bi")
Ca_biexp_decay1 <- Ca_biexp_fct(t = Time, tOn = tOn,
log_Ca0 = log(Ca0),
log_dCa = log(dCa),
log_tau = log(tau),
mu = mu,
log_dtau = log(dtau))
## Define a calcium biexponential decay with the
## caBiExp function
Ca_biexp_decay2 <- caBiExp(t = Time, tOn = tOn, Ca0 = Ca0,
dCa = dCa, tau = tau,
fact = 1/(1+exp(-mu)), dtau = dtau)
## Check that both decays are similar
plot(Time, Ca_biexp_decay1, type="l", lwd=2)
lines(Time, Ca_biexp_decay2, col="red", lty=2, lwd=2)