From e05656d57668688b971c28e32b4cfd4d3eac4662 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 18 Sep 2019 13:03:50 +0200 Subject: Correct and rebuild docs - Reconcile docs and code for max_twa_parent - Correct links to docs in twa vignette - Static documentation rebuilt by pkgdown --- docs/reference/logistic.solution.html | 42 +++++++++++++++++++++-------------- 1 file changed, 25 insertions(+), 17 deletions(-) (limited to 'docs/reference/logistic.solution.html') diff --git a/docs/reference/logistic.solution.html b/docs/reference/logistic.solution.html index eb4b0de0..91fe035e 100644 --- a/docs/reference/logistic.solution.html +++ b/docs/reference/logistic.solution.html @@ -16,13 +16,15 @@ - + + - - + + + @@ -48,6 +50,7 @@ + @@ -117,6 +120,7 @@ +
@@ -182,30 +186,30 @@

Examples

  # Reproduce the plot on page 57 of FOCUS (2014)
-  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.2),
-       from = 0, to = 100, ylim = c(0, 100),
-       xlab = "Time", ylab = "Residue")
  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.4),
-       from = 0, to = 100, add = TRUE, lty = 2, col = 2)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.8),
-       from = 0, to = 100, add = TRUE, lty = 3, col = 3)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.001, 0.2),
-       from = 0, to = 100, add = TRUE, lty = 4, col = 4)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.08, 0.2),
-       from = 0, to = 100, add = TRUE, lty = 5, col = 5)
  legend("topright", inset = 0.05,
-         legend = paste0("k0 = ", c(0.0001, 0.0001, 0.0001, 0.001, 0.08),
-                         ", r = ", c(0.2, 0.4, 0.8, 0.2, 0.2)),
+  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.2),
+       from = 0, to = 100, ylim = c(0, 100),
+       xlab = "Time", ylab = "Residue")
  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.4),
+       from = 0, to = 100, add = TRUE, lty = 2, col = 2)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.0001, 0.8),
+       from = 0, to = 100, add = TRUE, lty = 3, col = 3)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.001, 0.2),
+       from = 0, to = 100, add = TRUE, lty = 4, col = 4)
  plot(function(x) logistic.solution(x, 100, 0.08, 0.08, 0.2),
+       from = 0, to = 100, add = TRUE, lty = 5, col = 5)
  legend("topright", inset = 0.05,
+         legend = paste0("k0 = ", c(0.0001, 0.0001, 0.0001, 0.001, 0.08),
+                         ", r = ", c(0.2, 0.4, 0.8, 0.2, 0.2)),
          lty = 1:5, col = 1:5)

   # Fit with synthetic data
   logistic <- mkinmod(parent = mkinsub("logistic"))
 
-  sampling_times = c(0, 1, 3, 7, 14, 28, 60, 90, 120)
-  parms_logistic <- c(kmax = 0.08, k0 = 0.0001, r = 0.2)
+  sampling_times = c(0, 1, 3, 7, 14, 28, 60, 90, 120)
+  parms_logistic <- c(kmax = 0.08, k0 = 0.0001, r = 0.2)
   d_logistic <- mkinpredict(logistic,
-    parms_logistic, c(parent = 100),
+    parms_logistic, c(parent = 100),
     sampling_times)
   d_2_1 <- add_err(d_logistic,
     sdfunc = function(x) sigma_twocomp(x, 0.5, 0.07),
     n = 1, reps = 2, digits = 5, LOD = 0.1, seed = 123456)[[1]]
 
   m <- mkinfit("logistic", d_2_1, quiet = TRUE)
-  plot_sep(m)
  summary(m)$bpar
#>              Estimate   se_notrans   t value       Pr(>t)        Lower
+  plot_sep(m)
  summary(m)$bpar
#>              Estimate   se_notrans   t value       Pr(>t)        Lower
 #> parent_0 1.057896e+02 1.9023449649 55.610120 3.768361e-16 1.016451e+02
 #> kmax     6.398190e-02 0.0143201029  4.467978 3.841828e-04 3.929235e-02
 #> k0       1.612775e-04 0.0005866813  0.274898 3.940351e-01 5.846685e-08
@@ -236,19 +240,23 @@
+ + + -- cgit v1.2.1