aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2014-05-02 18:30:22 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2014-05-02 18:30:22 +0200
commit4a5f879782ea89aa9e9f087ec8047cba594536d8 (patch)
tree1f572cd22b8d13e72084d5fb3a1be2b3505faf61
parent3ac16665a0b543b128d993bc8a767eea7a2a84ad (diff)
Update of the FOCUS Z vignette, release candidate
-rw-r--r--DESCRIPTION2
-rw-r--r--vignettes/FOCUS_Z.Rnw57
-rw-r--r--vignettes/FOCUS_Z.pdfbin204581 -> 189030 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_1.pdfbin5670 -> 5670 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_10.pdfbin7587 -> 7587 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_11.pdfbin7605 -> 7598 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_11b.pdfbin0 -> 5008 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_2.pdfbin5670 -> 5670 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_3.pdfbin5670 -> 5670 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_5.pdfbin6155 -> 6155 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_6.pdfbin6966 -> 6966 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_7.pdfbin6967 -> 6967 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_8.pdfbin6158 -> 6158 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_fits_9.pdfbin6785 -> 6785 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_residuals_11.pdfbin5943 -> 5940 bytes
-rw-r--r--vignettes/figure/FOCUS_2006_Z_residuals_6.pdfbin5959 -> 5959 bytes
16 files changed, 45 insertions, 14 deletions
diff --git a/DESCRIPTION b/DESCRIPTION
index 618344a8..594d68c7 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -3,7 +3,7 @@ Type: Package
Title: Routines for fitting kinetic models with one or more state
variables to chemical degradation data
Version: 0.9-26
-Date: 2014-04-28
+Date: 2014-05-02
Authors@R: c(person("Johannes", "Ranke", role = c("aut", "cre", "cph"),
email = "jranke@uni-bremen.de"),
person("Katrin", "Lindenberger", role = "ctb"),
diff --git a/vignettes/FOCUS_Z.Rnw b/vignettes/FOCUS_Z.Rnw
index 44cfa468..694c5eea 100644
--- a/vignettes/FOCUS_Z.Rnw
+++ b/vignettes/FOCUS_Z.Rnw
@@ -85,7 +85,8 @@ is negligible. Accordingly, the exact magnitude of the fitted parameter
This suggests, in agreement with the analysis in the FOCUS kinetics report, to simplify
the model by removing the pathway to sink.
-A similar result can be obtained when formation fractions are used in the model formulation:
+A similar result can be obtained when formation fractions are used in the model
+formulation:
<<FOCUS_2006_Z_fits_2, echo=TRUE, fig.height=4>>=
Z.2a.ff <- mkinmod(Z0 = list(type = "SFO", to = "Z1"),
@@ -97,13 +98,13 @@ plot(m.Z.2a.ff)
summary(m.Z.2a.ff, data = FALSE)
@
-Here, the ilr transformed formation fraction fitted in the model takes a very large value,
-and the backtransformed formation fraction from parent Z to Z1 is practically unity. Again,
-the covariance matrix is not returned as the model is overparameterised.
+Here, the ilr transformed formation fraction fitted in the model takes a very
+large value, and the backtransformed formation fraction from parent Z to Z1 is
+practically unity. Again, the covariance matrix is not returned as the model is
+overparameterised.
The simplified model is obtained by setting the list component \texttt{sink} to
-\texttt{FALSE}. This model definition is not supported when formation fractions
-are used.
+\texttt{FALSE}.
<<FOCUS_2006_Z_fits_3, echo=TRUE, fig.height=4>>=
Z.3 <- mkinmod(Z0 = list(type = "SFO", to = "Z1", sink = FALSE),
@@ -113,6 +114,19 @@ plot(m.Z.3)
summary(m.Z.3, data = FALSE)
@
+This model definition is not supported when formation fractions
+are used, but the formation fraction can be fixed to unity.
+
+<<FOCUS_2006_Z_fits_3.ff, echo=TRUE, fig.height=4>>=
+Z.3.ff <- mkinmod(Z0 = list(type = "SFO", to = "Z1"),
+ Z1 = list(type = "SFO"), use_of_ff = "max")
+m.Z.3.ff <- mkinfit(Z.3.ff, FOCUS_2006_Z_mkin,
+ parms.ini = c(f_Z0_to_Z1 = 1),
+ fixed_parms = "f_Z0_to_Z1",
+ quiet = TRUE)
+summary(m.Z.3.ff, data = FALSE)
+@
+
\section{Including metabolites Z2 and Z3}
As suggested in the FOCUS report, the pathway to sink was removed for metabolite Z1 as
@@ -196,7 +210,7 @@ for the same reasons as in the original analysis.
<<FOCUS_2006_Z_fits_9, echo=TRUE, fig.height=4>>=
Z.mkin.3 <- mkinmod(Z0 = list(type = "SFORB", to = "Z1", sink = FALSE),
- Z1 = list(type = "SFO", to = "Z2"),
+ Z1 = list(type = "SFO", to = "Z2", sink = FALSE),
Z2 = list(type = "SFO"))
m.Z.mkin.3 <- mkinfit(Z.mkin.3, FOCUS_2006_Z_mkin, quiet = TRUE)
plot(m.Z.mkin.3)
@@ -206,8 +220,8 @@ summary(m.Z.mkin.3, data = FALSE)
This results in a much better representation of the behaviour of the parent
compound Z0.
-Finally, Z3 is added as well. This model appears overparameterised (no
-covariance matrix returned) if the sink for Z1 is left in the model.
+Finally, Z3 is added as well. These models appear overparameterised (no
+covariance matrix returned) if the sink for Z1 is left in the models.
<<FOCUS_2006_Z_fits_10, echo=TRUE, fig.height=4>>=
Z.mkin.4 <- mkinmod(Z0 = list(type = "SFORB", to = "Z1", sink = FALSE),
@@ -234,12 +248,29 @@ m.Z.mkin.5 <- mkinfit(Z.mkin.5, FOCUS_2006_Z_mkin,
parms.ini = m.Z.mkin.4$bparms.ode[1:5],
quiet = TRUE)
plot(m.Z.mkin.5)
-summary(m.Z.mkin.5, data = FALSE)
+summary(m.Z.mkin.5, data = FALSE)$bpar
+@
+
+The summary view of the backtransformed parameters shows that we get no
+confidence intervals due to overparameterisation. As the optimized
+\texttt{k\_Z3\_bound\_free} is excessively small, it is reasonable to fix it to
+zero.
+<<FOCUS_2006_Z_fits_11a, echo=TRUE, fig.height=4>>=
+m.Z.mkin.5a <- mkinfit(Z.mkin.5, FOCUS_2006_Z_mkin,
+ parms.ini = c(m.Z.mkin.4$bparms.ode[1:5],
+ k_Z3_bound_free = 0),
+ fixed_parms = "k_Z3_bound_free",
+ quiet = TRUE)
+summary(m.Z.mkin.5a, data = FALSE)$bpar
+@
+
+A graphical representation of the confidence intervals can finally be obtained.
+<<FOCUS_2006_Z_fits_11b, echo=TRUE, fig.height=4>>=
+mkinparplot(m.Z.mkin.5a)
@
-Looking at the confidence intervals of the SFORB model parameters of Z3, it is
-clear that nothing can be said about the degradation rate of Z3 towards the end
-of the experiment. However, this appears to be a feature of the data.
+It is clear that nothing can be said about the degradation rate of Z3 towards
+the end of the experiment. However, this appears to be a feature of the data.
<<FOCUS_2006_Z_residuals_11>>=
par(mfrow = c(2, 2))
diff --git a/vignettes/FOCUS_Z.pdf b/vignettes/FOCUS_Z.pdf
index ca67191e..e29de6ca 100644
--- a/vignettes/FOCUS_Z.pdf
+++ b/vignettes/FOCUS_Z.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_1.pdf b/vignettes/figure/FOCUS_2006_Z_fits_1.pdf
index 9ee45463..346b9f27 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_1.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_1.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_10.pdf b/vignettes/figure/FOCUS_2006_Z_fits_10.pdf
index c2776cab..26d2ec09 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_10.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_10.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_11.pdf b/vignettes/figure/FOCUS_2006_Z_fits_11.pdf
index 2b1c973c..577677bf 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_11.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_11.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_11b.pdf b/vignettes/figure/FOCUS_2006_Z_fits_11b.pdf
new file mode 100644
index 00000000..57880f13
--- /dev/null
+++ b/vignettes/figure/FOCUS_2006_Z_fits_11b.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_2.pdf b/vignettes/figure/FOCUS_2006_Z_fits_2.pdf
index 511fe81c..8b6fd715 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_2.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_2.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_3.pdf b/vignettes/figure/FOCUS_2006_Z_fits_3.pdf
index e5472f9a..9d04fed7 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_3.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_3.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_5.pdf b/vignettes/figure/FOCUS_2006_Z_fits_5.pdf
index 831b7f26..cf9e2ffb 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_5.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_5.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_6.pdf b/vignettes/figure/FOCUS_2006_Z_fits_6.pdf
index 6c4bcd33..3c157dd3 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_6.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_6.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_7.pdf b/vignettes/figure/FOCUS_2006_Z_fits_7.pdf
index e4ac1a83..b50959bd 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_7.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_7.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_8.pdf b/vignettes/figure/FOCUS_2006_Z_fits_8.pdf
index 11303976..21985bd3 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_8.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_8.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_fits_9.pdf b/vignettes/figure/FOCUS_2006_Z_fits_9.pdf
index 920353a7..970aabd3 100644
--- a/vignettes/figure/FOCUS_2006_Z_fits_9.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_fits_9.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf b/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf
index 40a6afb1..f69ffd9f 100644
--- a/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf
Binary files differ
diff --git a/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf b/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf
index a3006d8f..22f62bcb 100644
--- a/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf
+++ b/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf
Binary files differ

Contact - Imprint