From 70591022c07f0e8fb4dd67789b7c8d78af8ebc18 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Thu, 2 May 2019 13:17:05 +0200 Subject: Better initials for error model parameters - Also make it possible to specify initial values for error model parameters. - Run tests - Rebuild docs --- docs/reference/mkinpredict.html | 56 ++++++++++++++++++++--------------------- 1 file changed, 28 insertions(+), 28 deletions(-) (limited to 'docs/reference/mkinpredict.html') diff --git a/docs/reference/mkinpredict.html b/docs/reference/mkinpredict.html index 56c1ed27..4b75603f 100644 --- a/docs/reference/mkinpredict.html +++ b/docs/reference/mkinpredict.html @@ -328,7 +328,7 @@ c(parent = 100, m1 = 0), seq(0, 20, by = 0.1), solution_type = "eigen")[201,]))
#> time parent m1 #> 201 20 4.978707 27.46227
#> User System verstrichen -#> 0.004 0.000 0.003
system.time( +#> 0.003 0.000 0.003
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")[201,]))
#> time parent m1 @@ -338,34 +338,34 @@ 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.04 0.00 0.04
+#> 0.021 0.000 0.021
# Predict from a fitted model - f <- mkinfit(SFO_SFO, FOCUS_2006_C)
#> Negative log-likelihood at call 1: 552.5739 -#> Negative log-likelihood at call 3: 552.5739 -#> Negative log-likelihood at call 4: 552.5739 -#> Negative log-likelihood at call 6: 279.9345 -#> Negative log-likelihood at call 7: 279.9344 -#> Negative log-likelihood at call 8: 279.9294 -#> Negative log-likelihood at call 9: 279.9294 -#> Negative log-likelihood at call 12: 200.3629 -#> Negative log-likelihood at call 13: 200.3629 -#> Negative log-likelihood at call 18: 197.9039 -#> Negative log-likelihood at call 25: 196.6754 -#> Negative log-likelihood at call 27: 196.6754 -#> Negative log-likelihood at call 32: 196.5742 -#> Negative log-likelihood at call 33: 196.5742 -#> Negative log-likelihood at call 34: 196.5742 -#> Negative log-likelihood at call 38: 196.5361 -#> Negative log-likelihood at call 40: 196.5361 -#> Negative log-likelihood at call 44: 196.5336 -#> Negative log-likelihood at call 45: 196.5336 -#> Negative log-likelihood at call 50: 196.5334 -#> Negative log-likelihood at call 51: 196.5334 -#> Negative log-likelihood at call 52: 196.5334 -#> Negative log-likelihood at call 56: 196.5334 -#> Negative log-likelihood at call 58: 196.5334 -#> Negative log-likelihood at call 59: 196.5334 -#> Negative log-likelihood at call 65: 196.5334 + f <- mkinfit(SFO_SFO, FOCUS_2006_C)
#> 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 +#> Sum of squared residuals at call 7: 279.9344 +#> Sum of squared residuals at call 8: 279.9294 +#> Sum of squared residuals at call 9: 279.9294 +#> 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 25: 196.6754 +#> Sum of squared residuals at call 27: 196.6754 +#> Sum of squared residuals at call 32: 196.5742 +#> Sum of squared residuals at call 33: 196.5742 +#> Sum of squared residuals at call 34: 196.5742 +#> Sum of squared residuals at call 38: 196.5361 +#> Sum of squared residuals at call 40: 196.5361 +#> Sum of squared residuals at call 44: 196.5336 +#> Sum of squared residuals at call 45: 196.5336 +#> Sum of squared residuals at call 50: 196.5334 +#> Sum of squared residuals at call 51: 196.5334 +#> Sum of squared residuals at call 52: 196.5334 +#> Sum of squared residuals at call 56: 196.5334 +#> 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 #> 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 -- cgit v1.2.1