aboutsummaryrefslogtreecommitdiff
path: root/man/transform_odeparms.Rd
diff options
context:
space:
mode:
Diffstat (limited to 'man/transform_odeparms.Rd')
-rw-r--r--man/transform_odeparms.Rd51
1 files changed, 51 insertions, 0 deletions
diff --git a/man/transform_odeparms.Rd b/man/transform_odeparms.Rd
new file mode 100644
index 0000000..cc1c8f7
--- /dev/null
+++ b/man/transform_odeparms.Rd
@@ -0,0 +1,51 @@
+\name{transform_odeparms}
+\alias{transform_odeparms}
+\alias{backtransform_odeparms}
+\title{
+ Functions to transform and backtransform kinetic parameters for fitting
+}
+\description{
+ The transformations are intended to map parameters that should only take
+ on restricted values to the full scale of real numbers. For kinetic rate
+ constants and other paramters that can only take on positive values, a
+ simple log transformation is used. For compositional parameters, such as
+ the formations fractions that should always sum up to 1 and can not be
+ negative, the \code{\link{ilr}} transformation is used.
+}
+\usage{
+transform_odeparms(parms, mod_vars)
+backtransform_odeparms(transparms, mod_vars)
+}
+\arguments{
+ \item{parms}{
+ Parameters of kinetic models as used in the differential equations.
+}
+ \item{transparms}{
+ Transformed parameters of kinetic models as used in the fitting procedure.
+}
+ \item{mod_vars}{
+ Names of the state variables in the kinetic model. These are used for
+ grouping the formation fractions before \code{\link{ilr}} transformation.
+}
+}
+\value{
+ A vector of transformed or backtransformed parameters with the same names
+ as the original parameters.
+}
+\author{
+ Johannes Ranke
+}
+\examples{
+SFO_SFO <- mkinmod(
+ parent = list(type = "SFO", to = "m1", sink = TRUE),
+ m1 = list(type = "SFO"))
+# Fit the model to the FOCUS example dataset D using defaults
+fit <- mkinfit(SFO_SFO, FOCUS_2006_D)
+summary(fit, data=FALSE) # See transformed and backtransformed parameters
+initials <- fit$start$initial
+transformed <- fit$start$transformed
+names(initials) <- names(transformed) <- rownames(fit$start)
+transform_odeparms(initials, c("parent", "m1"))
+backtransform_odeparms(transformed, c("parent", "m1"))
+}
+\keyword{ manip }

Contact - Imprint