From 4a5f879782ea89aa9e9f087ec8047cba594536d8 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 2 May 2014 18:30:22 +0200 Subject: Update of the FOCUS Z vignette, release candidate --- vignettes/FOCUS_Z.Rnw | 57 +++++++++++++++++++------ vignettes/FOCUS_Z.pdf | Bin 204581 -> 189030 bytes vignettes/figure/FOCUS_2006_Z_fits_1.pdf | Bin 5670 -> 5670 bytes vignettes/figure/FOCUS_2006_Z_fits_10.pdf | Bin 7587 -> 7587 bytes vignettes/figure/FOCUS_2006_Z_fits_11.pdf | Bin 7605 -> 7598 bytes vignettes/figure/FOCUS_2006_Z_fits_11b.pdf | Bin 0 -> 5008 bytes vignettes/figure/FOCUS_2006_Z_fits_2.pdf | Bin 5670 -> 5670 bytes vignettes/figure/FOCUS_2006_Z_fits_3.pdf | Bin 5670 -> 5670 bytes vignettes/figure/FOCUS_2006_Z_fits_5.pdf | Bin 6155 -> 6155 bytes vignettes/figure/FOCUS_2006_Z_fits_6.pdf | Bin 6966 -> 6966 bytes vignettes/figure/FOCUS_2006_Z_fits_7.pdf | Bin 6967 -> 6967 bytes vignettes/figure/FOCUS_2006_Z_fits_8.pdf | Bin 6158 -> 6158 bytes vignettes/figure/FOCUS_2006_Z_fits_9.pdf | Bin 6785 -> 6785 bytes vignettes/figure/FOCUS_2006_Z_residuals_11.pdf | Bin 5943 -> 5940 bytes vignettes/figure/FOCUS_2006_Z_residuals_6.pdf | Bin 5959 -> 5959 bytes 15 files changed, 44 insertions(+), 13 deletions(-) create mode 100644 vignettes/figure/FOCUS_2006_Z_fits_11b.pdf (limited to 'vignettes') diff --git a/vignettes/FOCUS_Z.Rnw b/vignettes/FOCUS_Z.Rnw index 44cfa46..694c5ee 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: <>= 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}. <>= 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. + +<>= +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. <>= 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. <>= 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. +<>= +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. +<>= +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. <>= par(mfrow = c(2, 2)) diff --git a/vignettes/FOCUS_Z.pdf b/vignettes/FOCUS_Z.pdf index ca67191..e29de6c 100644 Binary files a/vignettes/FOCUS_Z.pdf and b/vignettes/FOCUS_Z.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_1.pdf b/vignettes/figure/FOCUS_2006_Z_fits_1.pdf index 9ee4546..346b9f2 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_1.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_1.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_10.pdf b/vignettes/figure/FOCUS_2006_Z_fits_10.pdf index c2776ca..26d2ec0 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_10.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_10.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_11.pdf b/vignettes/figure/FOCUS_2006_Z_fits_11.pdf index 2b1c973..577677b 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_11.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_11.pdf 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 0000000..57880f1 Binary files /dev/null and b/vignettes/figure/FOCUS_2006_Z_fits_11b.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_2.pdf b/vignettes/figure/FOCUS_2006_Z_fits_2.pdf index 511fe81..8b6fd71 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_2.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_2.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_3.pdf b/vignettes/figure/FOCUS_2006_Z_fits_3.pdf index e5472f9..9d04fed 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_3.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_3.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_5.pdf b/vignettes/figure/FOCUS_2006_Z_fits_5.pdf index 831b7f2..cf9e2ff 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_5.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_5.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_6.pdf b/vignettes/figure/FOCUS_2006_Z_fits_6.pdf index 6c4bcd3..3c157dd 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_6.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_6.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_7.pdf b/vignettes/figure/FOCUS_2006_Z_fits_7.pdf index e4ac1a8..b50959b 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_7.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_7.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_8.pdf b/vignettes/figure/FOCUS_2006_Z_fits_8.pdf index 1130397..21985bd 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_8.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_8.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_fits_9.pdf b/vignettes/figure/FOCUS_2006_Z_fits_9.pdf index 920353a..970aabd 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_fits_9.pdf and b/vignettes/figure/FOCUS_2006_Z_fits_9.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf b/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf index 40a6afb..f69ffd9 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf and b/vignettes/figure/FOCUS_2006_Z_residuals_11.pdf differ diff --git a/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf b/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf index a3006d8..22f62bc 100644 Binary files a/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf and b/vignettes/figure/FOCUS_2006_Z_residuals_6.pdf differ -- cgit v1.2.1