aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--ChangeLog5
-rw-r--r--R/mkinfit.R2
2 files changed, 4 insertions, 3 deletions
diff --git a/ChangeLog b/ChangeLog
index d87fa0d9..a57894d3 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,8 +1,9 @@
-2013-10-17 Johannes Ranke <jranke@stiller> for mkin (0.9-22)
+2013-10-17 Johannes Ranke <jranke@uni-bremen.de> for mkin (0.9-22)
* Fix plot.mkinfit as it passed graphical arguments like main to the solver
* Do not use plot=TRUE in mkinfit example
* The first successful fits in the not so simple GUI
+ * Fix iteratively reweighted least squares for the case of many metabolites
2013-10-16 Johannes Ranke<jranke@uni-bremen.de> for mkin (0.9-22)
@@ -16,7 +17,7 @@
* Implement IRLS assuming different variances for observed variables
-2013-10-09 Johannes Ranke <jranke@uni-bremen.de> for mkin (0.9-22)
+2013-10-09 Johannes Ranke for version (0.9-22)
* Do not use 0 values at time zero for chi2 error level calculations.
This is the way it is done in KinGUII and it makes sense
diff --git a/R/mkinfit.R b/R/mkinfit.R
index 70ca9c0e..83ac0873 100644
--- a/R/mkinfit.R
+++ b/R/mkinfit.R
@@ -215,7 +215,7 @@ mkinfit <- function(mkinmod, observed,
while (reweight.diff > reweight.tol & n.iter < reweight.max.iter) {
n.iter <- n.iter + 1
sigma.old <- sqrt(fit$var_ms_unweighted)
- observed[err] <- sqrt(fit$var_ms_unweighted)[observed$name]
+ observed[err] <- sqrt(fit$var_ms_unweighted)[as.character(observed$name)]
fit <- modFit(cost, fit$par, method = method.modFit,
control = control.modFit, ...)
reweight.diff = sum((sqrt(fit$var_ms_unweighted) - sigma.old)^2)

Contact - Imprint