Example evaluation of FOCUS Example Dataset D

+ + + + +

This is just a very simple vignette showing how to fit a degradation model for a parent compound with one transformation product using mkin. After loading the library we look a the data. We have observed concentrations in the column named value at the times specified in column time for the two observed variables named parent and m1.

library("mkin")

## Loading required package: minpack.lm

## Loading required package: rootSolve

## Loading required package: inline

## Loading required package: methods

## Loading required package: parallel

print(FOCUS_2006_D)

library("mkin")

## Lade nötiges Paket: minpack.lm

## Lade nötiges Paket: rootSolve

## Lade nötiges Paket: inline

## Lade nötiges Paket: methods

## Lade nötiges Paket: parallel

print(FOCUS_2006_D)

##      name time  value
 ## 1  parent    0  99.46
 ## 2  parent    0 102.04
@@ -179,27 +145,27 @@ $(document).ready(function () {
 ## 44     m1  120  33.31

Next we specify the degradation model: The parent compound degrades with simple first-order kinetics (SFO) to one metabolite named m1, which also degrades with SFO kinetics.

The call to mkinmod returns a degradation model. The differential equations represented in R code can be found in the character vector $diffs of the mkinmod object. If a C compiler (gcc) is installed and functional, the differential equation model will be compiled from auto-generated C code.

SFO_SFO <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))

SFO_SFO <- mkinmod(parent = mkinsub("SFO", "m1"), m1 = mkinsub("SFO"))

## Successfully compiled differential equation model from auto-generated C code.

print(SFO_SFO$diffs)

print(SFO_SFO$diffs)

##                                                       parent 
-## "d_parent = - k_parent_sink * parent - k_parent_m1 * parent" 
+## "d_parent = - k_parent_sink * parent - k_parent_m1 * parent" 
 ##                                                           m1 
-##             "d_m1 = + k_parent_m1 * parent - k_m1_sink * m1"

+## "d_m1 = + k_parent_m1 * parent - k_m1_sink * m1"

We do the fitting without progress report (quiet = TRUE).

fit <- mkinfit(SFO_SFO, FOCUS_2006_D, quiet = TRUE)

fit <- mkinfit(SFO_SFO, FOCUS_2006_D, quiet = TRUE)

A plot of the fit including a residual plot for both observed variables is obtained using the plot_sep method for mkinfit objects, which shows separate graphs for all compounds and their residuals.

plot_sep(fit, lpos = c("topright", "bottomright"))

plot_sep(fit, lpos = c("topright", "bottomright"))

Confidence intervals for the parameter estimates are obtained using the mkinparplot function.

mkinparplot(fit)

mkinparplot(fit)

A comprehensive report of the results is obtained using the summary method for mkinfit objects.

summary(fit)

## mkin version:    0.9.45 
-## R version:       3.4.0 
-## Date of fit:     Fri May  5 12:14:00 2017 
-## Date of summary: Fri May  5 12:14:00 2017 
+summary(fit)
+## mkin version:    0.9.46.1 
+## R version:       3.4.1 
+## Date of fit:     Thu Sep 14 12:15:01 2017 
+## Date of summary: Thu Sep 14 12:15:02 2017 
 ## 
 ## Equations:
 ## d_parent/dt = - k_parent_sink * parent - k_parent_m1 * parent
@@ -207,7 +173,7 @@ $(document).ready(function () {
 ## 
 ## Model predictions using solution type deSolve 
 ## 
-## Fitted with method Port using 153 model solutions performed in 1.033 s
+## Fitted with method Port using 153 model solutions performed in 1.14 s
 ## 
 ## Weighting: none
 ## 
@@ -318,25 +284,25 @@ $(document).ready(function () {
 ##   100       m1    33.13 3.198e+01  1.148e+00
 ##   120       m1    25.15 2.879e+01 -3.640e+00
 ##   120       m1    33.31 2.879e+01  4.520e+00
+


 
-
-
+  
+

Example evaluation of FOCUS Example Dataset D

Johannes Ranke

2017-05-05

Example evaluation of FOCUS Example Dataset D

Johannes Ranke

2017-09-14