From 51a1fb28416de58d1161924ec55cf8ac24e5fc75 Mon Sep 17 00:00:00 2001
From: Johannes Ranke
library("mkin")
-
-
-## Loading required package: minpack.lm
-## Loading required package: rootSolve
-
-
-FOCUS_2006_L1 = data.frame(
+FOCUS_2006_L1 = data.frame(
t = rep(c(0, 1, 2, 3, 5, 7, 14, 21, 30), each = 2),
parent = c(88.3, 91.4, 85.6, 84.5, 78.9, 77.6,
72.0, 71.9, 50.3, 59.4, 47.0, 45.1,
@@ -231,8 +225,8 @@ summary(m.L1.SFO)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:12 2014
-## Date of summary: Mon Jul 14 12:36:12 2014
+## Date of fit: Mon Jul 14 20:32:20 2014
+## Date of summary: Mon Jul 14 20:32:20 2014
##
## Equations:
## [1] d_parent = - k_parent_sink * parent
@@ -333,8 +327,8 @@ summary(m.L1.FOMC, data = FALSE)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:12 2014
-## Date of summary: Mon Jul 14 12:36:12 2014
+## Date of fit: Mon Jul 14 20:32:20 2014
+## Date of summary: Mon Jul 14 20:32:20 2014
##
## Equations:
## [1] d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent
@@ -425,8 +419,8 @@ summary(m.L2.SFO)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:12 2014
-## Date of summary: Mon Jul 14 12:36:13 2014
+## Date of fit: Mon Jul 14 20:32:20 2014
+## Date of summary: Mon Jul 14 20:32:20 2014
##
## Equations:
## [1] d_parent = - k_parent_sink * parent
@@ -534,8 +528,8 @@ mkinresplot(m.L2.FOMC)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:13 2014
-## Date of summary: Mon Jul 14 12:36:13 2014
+## Date of fit: Mon Jul 14 20:32:21 2014
+## Date of summary: Mon Jul 14 20:32:21 2014
##
## Equations:
## [1] d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent
@@ -619,8 +613,8 @@ plot(m.L2.DFOP)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:13 2014
-## Date of summary: Mon Jul 14 12:36:13 2014
+## Date of fit: Mon Jul 14 20:32:23 2014
+## Date of summary: Mon Jul 14 20:32:23 2014
##
## Equations:
## [1] d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * parent
@@ -705,8 +699,8 @@ plot(m.L3.SFO)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:14 2014
-## Date of summary: Mon Jul 14 12:36:14 2014
+## Date of fit: Mon Jul 14 20:32:23 2014
+## Date of summary: Mon Jul 14 20:32:23 2014
##
## Equations:
## [1] d_parent = - k_parent_sink * parent
@@ -790,8 +784,8 @@ plot(m.L3.FOMC)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:14 2014
-## Date of summary: Mon Jul 14 12:36:14 2014
+## Date of fit: Mon Jul 14 20:32:24 2014
+## Date of summary: Mon Jul 14 20:32:24 2014
##
## Equations:
## [1] d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent
@@ -862,8 +856,8 @@ plot(m.L3.DFOP)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:14 2014
-## Date of summary: Mon Jul 14 12:36:14 2014
+## Date of fit: Mon Jul 14 20:32:24 2014
+## Date of summary: Mon Jul 14 20:32:24 2014
##
## Equations:
## [1] d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * parent
@@ -952,8 +946,8 @@ plot(m.L4.SFO)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:15 2014
-## Date of summary: Mon Jul 14 12:36:15 2014
+## Date of fit: Mon Jul 14 20:32:27 2014
+## Date of summary: Mon Jul 14 20:32:27 2014
##
## Equations:
## [1] d_parent = - k_parent_sink * parent
@@ -1026,8 +1020,8 @@ plot(m.L4.FOMC)
## mkin version: 0.9.31
## R version: 3.1.1
-## Date of fit: Mon Jul 14 12:36:15 2014
-## Date of summary: Mon Jul 14 12:36:15 2014
+## Date of fit: Mon Jul 14 20:32:28 2014
+## Date of summary: Mon Jul 14 20:32:28 2014
##
## Equations:
## [1] d_parent = - (alpha/beta) * ((time/beta) + 1)^-1 * parent
diff --git a/vignettes/FOCUS_Z.pdf b/vignettes/FOCUS_Z.pdf
index 38490f09..559504cc 100644
Binary files a/vignettes/FOCUS_Z.pdf and b/vignettes/FOCUS_Z.pdf differ
diff --git a/vignettes/mkin.pdf b/vignettes/mkin.pdf
index f41a9d6a..f06d38f9 100644
Binary files a/vignettes/mkin.pdf and b/vignettes/mkin.pdf differ
--
cgit v1.2.1
From 34ace6b447062a492e526e5480562f495f096310 Mon Sep 17 00:00:00 2001
From: Johannes Ranke
Date: Mon, 14 Jul 2014 20:44:46 +0200
Subject: Small improvements to usage examples
---
README.md | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/README.md b/README.md
index 4d394808..e171180f 100644
--- a/README.md
+++ b/README.md
@@ -44,8 +44,8 @@ A very simple usage example would be
)
SFO <- mkinmod(parent = list(type = "SFO"))
SFO.fit <- mkinfit(SFO, example_data)
+ plot(SFO.fit, show_residuals = TRUE)
summary(SFO.fit)
- plot(SFO.fit)
A fairly complex usage example using a built-in dataset:
@@ -61,7 +61,8 @@ A fairly complex usage example using a built-in dataset:
fit <- mkinfit(model, data, method.modFit = "Port")
plot(fit, show_residuals = TRUE)
- s <- summary(fit)
+ summary(fit)
+ mkinparplot(fit)
For more examples and to see results, have a look at the examples provided in the
[`mkinfit`](http://kinfit.r-forge.r-project.org/mkin_static/mkinfit.html)
--
cgit v1.2.1
From 2abca7e5eadfdf7a9d1a95c448c9f62435e49745 Mon Sep 17 00:00:00 2001
From: Johannes Ranke
Date: Wed, 16 Jul 2014 09:40:40 +0200
Subject: Insert empty lines to have R make separate entries
---
NEWS.md | 30 ++++++++++++++++++++++++++++++
1 file changed, 30 insertions(+)
diff --git a/NEWS.md b/NEWS.md
index 94bc566a..3c0d72de 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -35,41 +35,61 @@
# CHANGES in mkin VERSION 0.9-29
- R/mkinresplot.R: Make it possible to specify `xlim`
+
- R/geometric_mean.R, man/geometric_mean.Rd: Add geometric mean function
+
- R/endpoints.R, man/endpoints.Rd: Calculate additional (pseudo)-DT50 values for FOMC, DFOP, HS and SFORB. Avoid calculation of formation fractions from rate constants when they are directly fitted
# CHANGES in mkin VERSION 0.9-28
- Do not backtransform confidence intervals for formation fractions if more than one compound is formed, as such parameters only define the pathways as a set
+
- Add historical remarks and some background to the main package vignette
+
- Correct 'isotropic' into 'isometric' for the ilr transformation
# CHANGES in mkin VERSION 0.9-27
- Fork the GUI into a separate package [gmkin](http://github.com/jranke/gmkin)
+
- DESCRIPTION, NAMESPACE, TODO: Adapt and add copyright information
+
- Remove files belonging to the GUI
+
- Possibility to fit without parameter transformations, using bounds as implemented in FME
+
- Add McCall 2,4,5-T dataset
+
- Enable selection of observed variables in plotting
+
- Add possibility to show residual plot in `plot.mkinfit`
+
- R/mkinparplot.R, man/mkinparplot.Rd: plot parameters with confidence intervals
# CHANGES in mkin VERSION 0.9-25
- Change vignette format from Sweave to knitr
+
- Split examples vignette to FOCUS_L and FOCUS_Z
+
- Remove warning about constant formation fractions in mkinmod as it was based on a misconception
+
- Restrict the unit test with the Schaefer data to parent and primary metabolites as formation fraction and DT50 for A2 are higly correlated and passing the test is platform dependent. For example, the test fails in 1 out of 14 platforms on CRAN as of today.
+
- Add Eurofins Regulatory AG copyright notices
+
- Import FME and deSolve instead of depending on them to have clean startup
+
- Add a starter function for the GUI: `gmkin()`
+
- Change the format of the workspace files of gmkin so they can be distributed and documented in the package
+
- Add gmkin workspace datasets FOCUS_2006_gmkin and FOCUS_2006_Z_gmkin
# CHANGES in mkin VERSION 0.9-24
- Bugfix re-enabling the fixing of any combination of initial values for state variables
+
- Default values for kinetic rate constants are not all 0.1 any more but are "salted" with a small increment to avoid numeric artefacts with the eigenvalue based solutions
# CHANGES in mkin VERSION 0.9-23
@@ -79,15 +99,25 @@
# CHANGES in mkin VERSION 0.9-22
- Get rid of the optimisation step in `mkinerrmin` - this was unnecessary. Thanks to KinGUII for the inspiration - actually this is equation 6-2 in FOCUS kinetics p. 91 that I had overlooked originally
+
- 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
+
- Unify naming of initial values of state variables
+
- Unify naming in dataframes of optimised and fixed parameters in the summary
+
- Show the weighting method for residuals in the summary
+
- Correct the output of the data in the case of manual weighting
+
- Implement IRLS assuming different variances for observed variables
+
- 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. It does impact the chi2 error levels in the output. Generally they seem to be lower for metabolites now, presumably because the mean of the observed values is higher
For a detailed list of changes to the mkin source please consult the commit history on http://github.com/jranke/mkin
--
cgit v1.2.1