From 8bdb4cd437a9d4663e542f95869e8692aa38dadb Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Thu, 7 May 2020 08:59:29 +0200 Subject: Static documentation rebuilt by pkgdown --- docs/reference/mkinfit.html | 171 ++++++++++++-------------------------------- 1 file changed, 47 insertions(+), 124 deletions(-) (limited to 'docs/reference/mkinfit.html') diff --git a/docs/reference/mkinfit.html b/docs/reference/mkinfit.html index 09837505..62823c58 100644 --- a/docs/reference/mkinfit.html +++ b/docs/reference/mkinfit.html @@ -10,23 +10,27 @@ - + - + - + + + + + - - + + - + - - + + @@ -43,7 +47,6 @@ kinetic model is solved using the function mkinpredict. The parameters of the selected error model are fitted simultaneously with the degradation model parameters, as both of them are arguments of the likelihood function." /> - @@ -61,7 +64,7 @@ likelihood function." /> - +
@@ -134,7 +142,7 @@ likelihood function." />

Fit a kinetic model to data with one or more state variables

- + Source: R/mkinfit.R
@@ -168,8 +176,8 @@ likelihood function.

rtol = 1e-10, n.outtimes = 100, error_model = c("const", "obs", "tc"), - error_model_algorithm = c("auto", "d_3", "direct", "twostep", "threestep", "fourstep", - "IRLS", "OLS"), + error_model_algorithm = c("auto", "d_3", "direct", "twostep", "threestep", + "fourstep", "IRLS", "OLS"), reweight.tol = 1e-08, reweight.max.iter = 10, trace_parms = FALSE, @@ -414,82 +422,7 @@ estimators.

Examples


 # Use shorthand notation for parent only degradation
-fit <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE)
-summary(fit)
#> mkin version used for fitting:    0.9.49.11 
-#> R version used for fitting:       3.6.3 
-#> Date of fit:     Mon Apr 20 18:52:22 2020 
-#> Date of summary: Mon Apr 20 18:52:22 2020 
-#> 
-#> Equations:
-#> d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent
-#> 
-#> Model predictions using solution type analytical 
-#> 
-#> Fitted using 222 model solutions performed in 0.473 s
-#> 
-#> Error model: Constant variance 
-#> 
-#> Error model algorithm: OLS 
-#> 
-#> Starting values for parameters to be optimised:
-#>          value   type
-#> parent_0  85.1  state
-#> alpha      1.0 deparm
-#> beta      10.0 deparm
-#> 
-#> Starting values for the transformed parameters actually optimised:
-#>               value lower upper
-#> parent_0  85.100000  -Inf   Inf
-#> log_alpha  0.000000  -Inf   Inf
-#> log_beta   2.302585  -Inf   Inf
-#> 
-#> Fixed parameter values:
-#> None
-#> 
-#> Optimised, transformed parameters with symmetric confidence intervals:
-#>           Estimate Std. Error    Lower   Upper
-#> parent_0  85.87000     1.8070 81.23000 90.5200
-#> log_alpha  0.05192     0.1353 -0.29580  0.3996
-#> log_beta   0.65100     0.2287  0.06315  1.2390
-#> sigma      1.85700     0.4378  0.73200  2.9830
-#> 
-#> Parameter correlation:
-#>             parent_0  log_alpha   log_beta     sigma
-#> parent_0   1.000e+00 -1.565e-01 -3.142e-01 4.758e-08
-#> log_alpha -1.565e-01  1.000e+00  9.564e-01 1.007e-07
-#> log_beta  -3.142e-01  9.564e-01  1.000e+00 8.568e-08
-#> sigma      4.758e-08  1.007e-07  8.568e-08 1.000e+00
-#> 
-#> Backtransformed parameters:
-#> Confidence intervals for internally transformed parameters are asymmetric.
-#> t-test (unrealistically) based on the assumption of normal distribution
-#> for estimators of untransformed parameters.
-#>          Estimate t value    Pr(>t)   Lower  Upper
-#> parent_0   85.870  47.530 3.893e-08 81.2300 90.520
-#> alpha       1.053   7.393 3.562e-04  0.7439  1.491
-#> beta        1.917   4.373 3.601e-03  1.0650  3.451
-#> sigma       1.857   4.243 4.074e-03  0.7320  2.983
-#> 
-#> FOCUS Chi2 error levels in percent:
-#>          err.min n.optim df
-#> All data   6.657       3  6
-#> parent     6.657       3  6
-#> 
-#> Estimated disappearance times:
-#>         DT50  DT90 DT50back
-#> parent 1.785 15.15     4.56
-#> 
-#> Data:
-#>  time variable observed predicted residual
-#>     0   parent     85.1    85.875  -0.7749
-#>     1   parent     57.9    55.191   2.7091
-#>     3   parent     29.9    31.845  -1.9452
-#>     7   parent     14.6    17.012  -2.4124
-#>    14   parent      9.7     9.241   0.4590
-#>    28   parent      6.6     4.754   1.8460
-#>    63   parent      4.0     2.102   1.8977
-#>    91   parent      3.9     1.441   2.4590
-#>   119   parent      0.6     1.092  -0.4919

+fit <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE)
#> Error in (function (t, parent_0, alpha, beta) {    parent = parent_0/(t/beta + 1)^alpha})(t = c(0, 1, 1.2020202020202, 2.4040404040404, 3, 3.60606060606061, 4.80808080808081, 6.01010101010101, 7, 7.21212121212121, 8.41414141414141, 9.61616161616162, 10.8181818181818, 12.020202020202, 13.2222222222222, 14, 14.4242424242424, 15.6262626262626, 16.8282828282828, 18.030303030303, 19.2323232323232, 20.4343434343434, 21.6363636363636, 22.8383838383838, 24.040404040404, 25.2424242424242, 26.4444444444444, 27.6464646464646, 28, 28.8484848484848, 30.050505050505, 31.2525252525253, 32.4545454545455, 33.6565656565657, 34.8585858585859, 36.0606060606061, 37.2626262626263, 38.4646464646465, 39.6666666666667, 40.8686868686869, 42.0707070707071, 43.2727272727273, 44.4747474747475, 45.6767676767677, 46.8787878787879, 48.0808080808081, 49.2828282828283, 50.4848484848485, 51.6868686868687, 52.8888888888889, 54.0909090909091, 55.2929292929293, 56.4949494949495, 57.6969696969697, 58.8989898989899, 60.1010101010101, 61.3030303030303, 62.5050505050505, 63, 63.7070707070707, 64.9090909090909, 66.1111111111111, 67.3131313131313, 68.5151515151515, 69.7171717171717, 70.9191919191919, 72.1212121212121, 73.3232323232323, 74.5252525252525, 75.7272727272727, 76.9292929292929, 78.1313131313131, 79.3333333333333, 80.5353535353535, 81.7373737373737, 82.9393939393939, 84.1414141414141, 85.3434343434343, 86.5454545454545, 87.7474747474747, 88.9494949494949, 90.1515151515152, 91, 91.3535353535353, 92.5555555555556, 93.7575757575758, 94.959595959596, 96.1616161616162, 97.3636363636364, 98.5656565656566, 99.7676767676768, 100.969696969697, 102.171717171717, 103.373737373737, 104.575757575758, 105.777777777778, 106.979797979798, 108.181818181818, 109.383838383838, 110.585858585859, 111.787878787879, 112.989898989899, 114.191919191919, 115.393939393939, 116.59595959596, 117.79797979798, 119), parent.0 = c(parent = 85.1),     alpha = 1, beta = 10): unbenutztes Argument (parent.0 = 85.1)
#> Timing stopped at: 0 0 0.001
summary(fit)
#> Error in summary(fit): Objekt 'fit' nicht gefunden

 # One parent compound, one metabolite, both single first order.
 # Use mkinsub for convenience in model formulation. Pathway to sink included per default.
 SFO_SFO <- mkinmod(
@@ -497,7 +430,8 @@ estimators.

   m1 = mkinsub("SFO"))
#> Successfully compiled differential equation model from auto-generated C code.
# Fit the model to the FOCUS example dataset D using defaults
 print(system.time(fit <- mkinfit(SFO_SFO, FOCUS_2006_D,
                            solution_type = "eigen", quiet = TRUE)))
#> Warning: Observations with value of zero were removed from the data
#>        User      System verstrichen 
-#>       1.465       0.003       1.474 
coef(fit)
#> NULL
endpoints(fit)
#> $ff
+#>       1.526       0.000       1.526 
parms(fit)
#>      parent_0 k_parent_sink   k_parent_m1     k_m1_sink         sigma 
+#>  99.598483069   0.047920122   0.050777612   0.005260651   3.125503875 
endpoints(fit)
#> $ff
 #> parent_sink   parent_m1     m1_sink 
 #>    0.485524    0.514476    1.000000 
 #> 
@@ -567,7 +501,8 @@ estimators.

 #> Sum of squared residuals at call 126: 371.2134
 #> Sum of squared residuals at call 135: 371.2134
 #> Negative log-likelihood at call 145: 97.22429
#> Optimisation successfully terminated.
#>        User      System verstrichen 
-#>       1.052       0.000       1.068 
coef(fit.deSolve)
#> NULL
endpoints(fit.deSolve)
#> $ff
+#>       1.084       0.000       1.085 
parms(fit.deSolve)
#>      parent_0 k_parent_sink   k_parent_m1     k_m1_sink         sigma 
+#>  99.598483072   0.047920122   0.050777612   0.005260651   3.125503874 
endpoints(fit.deSolve)
#> $ff
 #> parent_sink   parent_m1     m1_sink 
 #>    0.485524    0.514476    1.000000 
 #> 
@@ -583,25 +518,23 @@ estimators.

   parent = mkinsub("FOMC", "m1"),
   m1 = mkinsub("SFO"))
#> Successfully compiled differential equation model from auto-generated C code.
# Fit the model to the FOCUS example dataset D using defaults
 fit.FOMC_SFO <- mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
# Use starting parameters from parent only FOMC fit
-fit.FOMC = mkinfit("FOMC", FOCUS_2006_D, quiet = TRUE)
-fit.FOMC_SFO <- mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE,
-  parms.ini = fit.FOMC$bparms.ode)
#> Warning: Observations with value of zero were removed from the data

+fit.FOMC = mkinfit("FOMC", FOCUS_2006_D, quiet = TRUE)
#> Error in (function (t, parent_0, alpha, beta) {    parent = parent_0/(t/beta + 1)^alpha})(t = c(0, 0.757575757575758, 1, 1.51515151515152, 2.27272727272727, 3, 3.03030303030303, 3.78787878787879, 4.54545454545454, 5.3030303030303, 6.06060606060606, 6.81818181818182, 7, 7.57575757575758, 8.33333333333333, 9.09090909090909, 9.84848484848485, 10.6060606060606, 11.3636363636364, 12.1212121212121, 12.8787878787879, 13.6363636363636, 14, 14.3939393939394, 15.1515151515152, 15.9090909090909, 16.6666666666667, 17.4242424242424, 18.1818181818182, 18.9393939393939, 19.6969696969697, 20.4545454545455, 21, 21.2121212121212, 21.969696969697, 22.7272727272727, 23.4848484848485, 24.2424242424242, 25, 25.7575757575758, 26.5151515151515, 27.2727272727273, 28.030303030303, 28.7878787878788, 29.5454545454545, 30.3030303030303, 31.0606060606061, 31.8181818181818, 32.5757575757576, 33.3333333333333, 34.0909090909091, 34.8484848484849, 35, 35.6060606060606, 36.3636363636364, 37.1212121212121, 37.8787878787879, 38.6363636363636, 39.3939393939394, 40.1515151515151, 40.9090909090909, 41.6666666666667, 42.4242424242424, 43.1818181818182, 43.9393939393939, 44.6969696969697, 45.4545454545455, 46.2121212121212, 46.969696969697, 47.7272727272727, 48.4848484848485, 49.2424242424242, 50, 50.7575757575758, 51.5151515151515, 52.2727272727273, 53.030303030303, 53.7878787878788, 54.5454545454545, 55.3030303030303, 56.0606060606061, 56.8181818181818, 57.5757575757576, 58.3333333333333, 59.0909090909091, 59.8484848484849, 60.6060606060606, 61.3636363636364, 62.1212121212121, 62.8787878787879, 63.6363636363636, 64.3939393939394, 65.1515151515152, 65.9090909090909, 66.6666666666667, 67.4242424242424, 68.1818181818182, 68.9393939393939, 69.6969696969697, 70.4545454545455, 71.2121212121212, 71.969696969697, 72.7272727272727, 73.4848484848485, 74.2424242424242, 75), parent.0 = c(parent = 100.75), alpha = 1,     beta = 10): unbenutztes Argument (parent.0 = 100.75)
#> Timing stopped at: 0 0 0.001
fit.FOMC_SFO <- mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE,
+  parms.ini = fit.FOMC$bparms.ode)
#> Warning: Observations with value of zero were removed from the data
#> Error in mkinfit(FOMC_SFO, FOCUS_2006_D, quiet = TRUE, parms.ini = fit.FOMC$bparms.ode): Objekt 'fit.FOMC' nicht gefunden

 # Use stepwise fitting, using optimised parameters from parent only fit, SFORB
 SFORB_SFO <- mkinmod(
   parent = list(type = "SFORB", to = "m1", sink = TRUE),
   m1 = list(type = "SFO"))
#> Successfully compiled differential equation model from auto-generated C code.
# Fit the model to the FOCUS example dataset D using defaults
 fit.SFORB_SFO <- mkinfit(SFORB_SFO, FOCUS_2006_D, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
fit.SFORB_SFO.deSolve <- mkinfit(SFORB_SFO, FOCUS_2006_D, solution_type = "deSolve",
                                  quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
# Use starting parameters from parent only SFORB fit (not really needed in this case)
-fit.SFORB = mkinfit("SFORB", FOCUS_2006_D, quiet = TRUE)
-fit.SFORB_SFO <- mkinfit(SFORB_SFO, FOCUS_2006_D, parms.ini = fit.SFORB$bparms.ode, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
# }
+fit.SFORB = mkinfit("SFORB", FOCUS_2006_D, quiet = TRUE)
#> Error in (function (t, parent_0, k_12, k_21, k_1output) {    sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 + k_12 *         k_21 - (k_12 + k_1output) * k_21)    b1 = 0.5 * (k_12 + k_21 + k_1output) + sqrt_exp    b2 = 0.5 * (k_12 + k_21 + k_1output) - sqrt_exp    parent = parent_0 * (((k_12 + k_21 - b1)/(b2 - b1)) * exp(-b1 *         t) + ((k_12 + k_21 - b2)/(b1 - b2)) * exp(-b2 * t))})(t = c(0, 0.757575757575758, 1, 1.51515151515152, 2.27272727272727, 3, 3.03030303030303, 3.78787878787879, 4.54545454545454, 5.3030303030303, 6.06060606060606, 6.81818181818182, 7, 7.57575757575758, 8.33333333333333, 9.09090909090909, 9.84848484848485, 10.6060606060606, 11.3636363636364, 12.1212121212121, 12.8787878787879, 13.6363636363636, 14, 14.3939393939394, 15.1515151515152, 15.9090909090909, 16.6666666666667, 17.4242424242424, 18.1818181818182, 18.9393939393939, 19.6969696969697, 20.4545454545455, 21, 21.2121212121212, 21.969696969697, 22.7272727272727, 23.4848484848485, 24.2424242424242, 25, 25.7575757575758, 26.5151515151515, 27.2727272727273, 28.030303030303, 28.7878787878788, 29.5454545454545, 30.3030303030303, 31.0606060606061, 31.8181818181818, 32.5757575757576, 33.3333333333333, 34.0909090909091, 34.8484848484849, 35, 35.6060606060606, 36.3636363636364, 37.1212121212121, 37.8787878787879, 38.6363636363636, 39.3939393939394, 40.1515151515151, 40.9090909090909, 41.6666666666667, 42.4242424242424, 43.1818181818182, 43.9393939393939, 44.6969696969697, 45.4545454545455, 46.2121212121212, 46.969696969697, 47.7272727272727, 48.4848484848485, 49.2424242424242, 50, 50.7575757575758, 51.5151515151515, 52.2727272727273, 53.030303030303, 53.7878787878788, 54.5454545454545, 55.3030303030303, 56.0606060606061, 56.8181818181818, 57.5757575757576, 58.3333333333333, 59.0909090909091, 59.8484848484849, 60.6060606060606, 61.3636363636364, 62.1212121212121, 62.8787878787879, 63.6363636363636, 64.3939393939394, 65.1515151515152, 65.9090909090909, 66.6666666666667, 67.4242424242424, 68.1818181818182, 68.9393939393939, 69.6969696969697, 70.4545454545455, 71.2121212121212, 71.969696969697, 72.7272727272727, 73.4848484848485, 74.2424242424242, 75), parent.0 = c(parent_free = 100.75), k_1output = 0.1,     k_12 = 0.1, k_21 = 0.02): unbenutztes Argument (parent.0 = 100.75)
#> Timing stopped at: 0.001 0 0.001
fit.SFORB_SFO <- mkinfit(SFORB_SFO, FOCUS_2006_D, parms.ini = fit.SFORB$bparms.ode, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
#> Error in mkinfit(SFORB_SFO, FOCUS_2006_D, parms.ini = fit.SFORB$bparms.ode,     quiet = TRUE): Objekt 'fit.SFORB' nicht gefunden
# }
 
 # \dontrun{
 # Weighted fits, including IRLS
 SFO_SFO.ff <- mkinmod(parent = mkinsub("SFO", "m1"),
                       m1 = mkinsub("SFO"), use_of_ff = "max")
#> Successfully compiled differential equation model from auto-generated C code.
f.noweight <- mkinfit(SFO_SFO.ff, FOCUS_2006_D, quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
summary(f.noweight)
#> mkin version used for fitting:    0.9.49.11 
-#> R version used for fitting:       3.6.3 
-#> Date of fit:     Mon Apr 20 18:52:37 2020 
-#> Date of summary: Mon Apr 20 18:52:37 2020 
+#> R version used for fitting:       4.0.0 
+#> Date of fit:     Thu May  7 08:59:03 2020 
+#> Date of summary: Thu May  7 08:59:03 2020 
 #> 
 #> Equations:
 #> d_parent/dt = - k_parent * parent
@@ -609,7 +542,7 @@ estimators.

 #> 
 #> Model predictions using solution type deSolve 
 #> 
-#> Fitted using 422 model solutions performed in 1.087 s
+#> Fitted using 422 model solutions performed in 1.106 s
 #> 
 #> Error model: Constant variance 
 #> 
@@ -716,9 +649,9 @@ estimators.

 #>   100       m1    33.13  31.98163  1.148e+00
 #>   120       m1    25.15  28.78984 -3.640e+00
 #>   120       m1    33.31  28.78984  4.520e+00
f.obs <- mkinfit(SFO_SFO.ff, FOCUS_2006_D, error_model = "obs", quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
summary(f.obs)
#> mkin version used for fitting:    0.9.49.11 
-#> R version used for fitting:       3.6.3 
-#> Date of fit:     Mon Apr 20 18:52:40 2020 
-#> Date of summary: Mon Apr 20 18:52:40 2020 
+#> R version used for fitting:       4.0.0 
+#> Date of fit:     Thu May  7 08:59:06 2020 
+#> Date of summary: Thu May  7 08:59:06 2020 
 #> 
 #> Equations:
 #> d_parent/dt = - k_parent * parent
@@ -726,7 +659,7 @@ estimators.

 #> 
 #> Model predictions using solution type deSolve 
 #> 
-#> Fitted using 979 model solutions performed in 2.518 s
+#> Fitted using 979 model solutions performed in 2.604 s
 #> 
 #> Error model: Variance unique to each observed variable 
 #> 
@@ -848,9 +781,9 @@ estimators.

 #>   100       m1    33.13  31.98773  1.142e+00
 #>   120       m1    25.15  28.80429 -3.654e+00
 #>   120       m1    33.31  28.80429  4.506e+00
f.tc <- mkinfit(SFO_SFO.ff, FOCUS_2006_D, error_model = "tc", quiet = TRUE)
#> Warning: Observations with value of zero were removed from the data
summary(f.tc)
#> mkin version used for fitting:    0.9.49.11 
-#> R version used for fitting:       3.6.3 
-#> Date of fit:     Mon Apr 20 18:52:50 2020 
-#> Date of summary: Mon Apr 20 18:52:50 2020 
+#> R version used for fitting:       4.0.0 
+#> Date of fit:     Thu May  7 08:59:16 2020 
+#> Date of summary: Thu May  7 08:59:16 2020 
 #> 
 #> Equations:
 #> d_parent/dt = - k_parent * parent
@@ -858,7 +791,7 @@ estimators.

 #> 
 #> Model predictions using solution type deSolve 
 #> 
-#> Fitted using 2552 model solutions performed in 10.186 s
+#> Fitted using 2552 model solutions performed in 10.544 s
 #> 
 #> Error model: Two-component variance function 
 #> 
@@ -977,20 +910,10 @@ estimators.
- @@ -1001,7 +924,7 @@ estimators.

-

Site built with pkgdown 1.4.1.

+

Site built with pkgdown 1.5.1.

-- cgit v1.2.1