From 1ef7008be2a72a0847064ad9c2ddcfa16b055482 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 3 May 2019 19:14:15 +0200 Subject: Improve error model fitting Now we have a three stage fitting process for nonconstant error models: - Unweighted least squares - Only optimize the error model - Optimize both Static documentation rebuilt by pkgdown --- docs/articles/FOCUS_D.html | 42 +++++++++++++++++++++--------------------- 1 file changed, 21 insertions(+), 21 deletions(-) (limited to 'docs/articles/FOCUS_D.html') diff --git a/docs/articles/FOCUS_D.html b/docs/articles/FOCUS_D.html index 3ce94de7..ae895cc3 100644 --- a/docs/articles/FOCUS_D.html +++ b/docs/articles/FOCUS_D.html @@ -88,7 +88,7 @@

Example evaluation of FOCUS Example Dataset D

Johannes Ranke

-

2019-05-02

+

2019-05-03

@@ -168,8 +168,8 @@
summary(fit)
## mkin version used for fitting:    0.9.49.4 
 ## R version used for fitting:       3.6.0 
-## Date of fit:     Thu May  2 18:49:13 2019 
-## Date of summary: Thu May  2 18:49:14 2019 
+## Date of fit:     Fri May  3 19:08:45 2019 
+## Date of summary: Fri May  3 19:08:46 2019 
 ## 
 ## Equations:
 ## d_parent/dt = - k_parent_sink * parent - k_parent_m1 * parent
@@ -177,7 +177,7 @@
 ## 
 ## Model predictions using solution type deSolve 
 ## 
-## Fitted using 396 model solutions performed in 0.999 s
+## Fitted using 389 model solutions performed in 0.983 s
 ## 
 ## Error model:
 ## Constant variance 
@@ -216,13 +216,13 @@
 ## log_k_parent_sink  6.067e-01         1.000e+00      -8.550e-02
 ## log_k_parent_m1   -6.372e-02        -8.550e-02       1.000e+00
 ## log_k_m1_sink     -1.688e-01        -6.252e-01       4.731e-01
-## sigma             -3.368e-07         3.365e-08       8.420e-07
+## sigma              1.164e-09        -8.908e-10       1.652e-08
 ##                   log_k_m1_sink      sigma
-## parent_0             -1.688e-01 -3.368e-07
-## log_k_parent_sink    -6.252e-01  3.365e-08
-## log_k_parent_m1       4.731e-01  8.420e-07
-## log_k_m1_sink         1.000e+00  1.958e-08
-## sigma                 1.958e-08  1.000e+00
+## parent_0             -1.688e-01  1.164e-09
+## log_k_parent_sink    -6.252e-01 -8.908e-10
+## log_k_parent_m1       4.731e-01  1.652e-08
+## log_k_m1_sink         1.000e+00 -1.340e-10
+## sigma                -1.340e-10  1.000e+00
 ## 
 ## Backtransformed parameters:
 ## Confidence intervals for internally transformed parameters are asymmetric.
@@ -258,10 +258,10 @@
 ##     0   parent   102.04  99.59848  2.442e+00
 ##     1   parent    93.50  90.23787  3.262e+00
 ##     1   parent    92.50  90.23787  2.262e+00
-##     3   parent    63.23  74.07319 -1.084e+01
-##     3   parent    68.99  74.07319 -5.083e+00
-##     7   parent    52.32  49.91206  2.408e+00
-##     7   parent    55.13  49.91206  5.218e+00
+##     3   parent    63.23  74.07320 -1.084e+01
+##     3   parent    68.99  74.07320 -5.083e+00
+##     7   parent    52.32  49.91207  2.408e+00
+##     7   parent    55.13  49.91207  5.218e+00
 ##    14   parent    27.27  25.01257  2.257e+00
 ##    14   parent    26.64  25.01257  1.627e+00
 ##    21   parent    11.50  12.53462 -1.035e+00
@@ -271,23 +271,23 @@
 ##    50   parent     0.69   0.71624 -2.624e-02
 ##    50   parent     0.63   0.71624 -8.624e-02
 ##    75   parent     0.05   0.06074 -1.074e-02
-##    75   parent     0.06   0.06074 -7.381e-04
+##    75   parent     0.06   0.06074 -7.382e-04
 ##     1       m1     4.84   4.80296  3.704e-02
 ##     1       m1     5.64   4.80296  8.370e-01
 ##     3       m1    12.91  13.02400 -1.140e-01
 ##     3       m1    12.96  13.02400 -6.400e-02
 ##     7       m1    22.97  25.04476 -2.075e+00
 ##     7       m1    24.47  25.04476 -5.748e-01
-##    14       m1    41.69  36.69003  5.000e+00
-##    14       m1    33.21  36.69003 -3.480e+00
+##    14       m1    41.69  36.69002  5.000e+00
+##    14       m1    33.21  36.69002 -3.480e+00
 ##    21       m1    44.37  41.65310  2.717e+00
 ##    21       m1    46.44  41.65310  4.787e+00
 ##    35       m1    41.22  43.31312 -2.093e+00
 ##    35       m1    37.95  43.31312 -5.363e+00
-##    50       m1    41.19  41.21832 -2.832e-02
-##    50       m1    40.01  41.21832 -1.208e+00
-##    75       m1    40.09  36.44704  3.643e+00
-##    75       m1    33.85  36.44704 -2.597e+00
+##    50       m1    41.19  41.21831 -2.831e-02
+##    50       m1    40.01  41.21831 -1.208e+00
+##    75       m1    40.09  36.44703  3.643e+00
+##    75       m1    33.85  36.44703 -2.597e+00
 ##   100       m1    31.04  31.98163 -9.416e-01
 ##   100       m1    33.13  31.98163  1.148e+00
 ##   120       m1    25.15  28.78984 -3.640e+00
-- 
cgit v1.2.1