From 316f1360d1e12eefe491f86b0bbcd6dcf091c736 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Sat, 12 Jul 2014 08:32:08 +0200 Subject: First working fits with IORE model --- R/transform_odeparms.R | 14 +++++++++++++- 1 file changed, 13 insertions(+), 1 deletion(-) (limited to 'R/transform_odeparms.R') diff --git a/R/transform_odeparms.R b/R/transform_odeparms.R index 912a5c0a..a36b7eae 100644 --- a/R/transform_odeparms.R +++ b/R/transform_odeparms.R @@ -33,6 +33,8 @@ transform_odeparms <- function(parms, mkinmod, # Log transformation for rate constants if requested k <- parms[grep("^k_", names(parms))] + k.iore <- parms[grep("^k.iore_", names(parms))] + k <- c(k, k.iore) if (length(k) > 0) { if(transform_rates) { transparms[paste0("log_", names(k))] <- log(k) @@ -40,6 +42,10 @@ transform_odeparms <- function(parms, mkinmod, else transparms[names(k)] <- k } + # Do not transform exponents in IORE models + N <- parms[grep("^N", names(parms))] + transparms[names(N)] <- N + # Go through state variables and apply isometric logratio transformation to # formation fractions if requested mod_vars = names(spec) @@ -98,8 +104,10 @@ backtransform_odeparms <- function(transparms, mkinmod, # Exponential transformation for rate constants if(transform_rates) { trans_k <- transparms[grep("^log_k_", names(transparms))] + trans_k.iore <- transparms[grep("^log_k.iore_", names(transparms))] + trans_k = c(trans_k, trans_k.iore) if (length(trans_k) > 0) { - k_names <- gsub("^log_k_", "k_", names(trans_k)) + k_names <- gsub("^log_k", "k", names(trans_k)) parms[k_names] <- exp(trans_k) } } else { @@ -107,6 +115,10 @@ backtransform_odeparms <- function(transparms, mkinmod, parms[names(trans_k)] <- trans_k } + # Do not transform exponents in IORE models + N <- transparms[grep("^N", names(transparms))] + parms[names(N)] <- N + # Go through state variables and apply inverse isometric logratio transformation mod_vars = names(spec) for (box in mod_vars) { -- cgit v1.2.3