From 4c6f29fe2a3ece5a85160b891c89ce0f55299c11 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Wed, 23 Jul 2014 08:34:59 +0200
Subject: Parallel metabolite formation with formation fractions in mkinerrmin

---
 R/mkinerrmin.R | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

(limited to 'R')

diff --git a/R/mkinerrmin.R b/R/mkinerrmin.R
index 9ebac6a4..09724730 100644
--- a/R/mkinerrmin.R
+++ b/R/mkinerrmin.R
@@ -65,11 +65,12 @@ mkinerrmin <- function(fit, alpha = 0.05)
     # Formation fractions are attributed to the target variable, so look
     # for source compartments with formation fractions
     for (source_var in fit$obs_vars) {
+      n.ff.source = length(grep(paste("^f", source_var, sep = "_"),
+                                 names(parms.optim)))
+      n.paths.source = length(fit$mkinmod$spec[[source_var]]$to)
       for (target_var in fit$mkinmod$spec[[source_var]]$to) {
         if (obs_var == target_var) {
-          n.ff.optim <- n.ff.optim +  
-                        length(grep(paste("^f", source_var, sep = "_"),
-                                    names(parms.optim))) 
+          n.ff.optim <- n.ff.optim + n.ff.source/n.paths.source
         }
       }
     }
-- 
cgit v1.2.1