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/reference/mkinpredict.html | 18 ++++++++++-------- 1 file changed, 10 insertions(+), 8 deletions(-) (limited to 'docs/reference/mkinpredict.html') diff --git a/docs/reference/mkinpredict.html b/docs/reference/mkinpredict.html index 7b75b85d..074ac3ac 100644 --- a/docs/reference/mkinpredict.html +++ b/docs/reference/mkinpredict.html @@ -333,14 +333,14 @@ c(parent = 100, m1 = 0), seq(0, 20, by = 0.1), solution_type = "deSolve")[201,]))
#> time parent m1 #> 201 20 4.978707 27.46227
#> User System verstrichen -#> 0.002 0.000 0.002
system.time( +#> 0.002 0.000 0.001
system.time( print(mkinpredict(SFO_SFO, c(k_parent_m1 = 0.05, k_parent_sink = 0.1, k_m1_sink = 0.01), c(parent = 100, m1 = 0), seq(0, 20, by = 0.1), solution_type = "deSolve", use_compiled = FALSE)[201,]))
#> time parent m1 #> 201 20 4.978707 27.46227
#> User System verstrichen -#> 0.021 0.000 0.021
+#> 0.022 0.000 0.021
# Predict from a fitted model - f <- mkinfit(SFO_SFO, FOCUS_2006_C)
#> Sum of squared residuals at call 1: 552.5739 + f <- mkinfit(SFO_SFO, FOCUS_2006_C)
#> Ordinary least squares optimisation
#> Sum of squared residuals at call 1: 552.5739 #> Sum of squared residuals at call 3: 552.5739 #> Sum of squared residuals at call 4: 552.5739 #> Sum of squared residuals at call 6: 279.9345 @@ -350,6 +350,7 @@ #> Sum of squared residuals at call 12: 200.3629 #> Sum of squared residuals at call 13: 200.3629 #> Sum of squared residuals at call 18: 197.9039 +#> Sum of squared residuals at call 23: 197.9039 #> Sum of squared residuals at call 25: 196.6754 #> Sum of squared residuals at call 27: 196.6754 #> Sum of squared residuals at call 32: 196.5742 @@ -366,15 +367,16 @@ #> Sum of squared residuals at call 58: 196.5334 #> Sum of squared residuals at call 59: 196.5334 #> Sum of squared residuals at call 65: 196.5334 +#> Sum of squared residuals at call 73: 196.5334 #> Negative log-likelihood at call 75: 26.64668 #> Negative log-likelihood at call 103: 26.64668 #> Optimisation successfully terminated.
head(mkinpredict(f))
#> time parent m1 #> 1 0.0 82.49216 0.000000 -#> 2 0.1 80.00563 1.179956 -#> 3 0.2 77.59404 2.312583 -#> 4 0.3 75.25515 3.399424 -#> 5 0.4 72.98675 4.441974 -#> 6 0.5 70.78673 5.441685
+#> 2 0.1 80.00563 1.179955 +#> 3 0.2 77.59404 2.312580 +#> 4 0.3 75.25515 3.399419 +#> 5 0.4 72.98675 4.441969 +#> 6 0.5 70.78673 5.441679