From d89e3d22eb9dc383897b09e9c5aa1b57f65cdbf0 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Thu, 21 Feb 2019 14:34:45 +0100
Subject: Add the logistic model

---
 R/endpoints.R | 8 ++++++++
 1 file changed, 8 insertions(+)

(limited to 'R/endpoints.R')

diff --git a/R/endpoints.R b/R/endpoints.R
index ac1e3e7c..80450185 100644
--- a/R/endpoints.R
+++ b/R/endpoints.R
@@ -159,6 +159,14 @@ endpoints <- function(fit) {
       ep$distimes[obs_var, c(paste("DT50", obs_var, "b1", sep = "_"))] = DT50_b1
       ep$distimes[obs_var, c(paste("DT50", obs_var, "b2", sep = "_"))] = DT50_b2
     }
+    if (type == "logistic") {
+      # FOCUS kinetics (2014) p. 67
+      kmax = parms.all["kmax"]
+      k0 = parms.all["k0"]
+      r = parms.all["r"]
+      DT50 = (1/r) * log(1 - ((kmax/k0) * (1 - 2^(r/kmax))))
+      DT90 = (1/r) * log(1 - ((kmax/k0) * (1 - 10^(r/kmax))))
+    }
     ep$distimes[obs_var, c("DT50", "DT90")] = c(DT50, DT90)
   }
   return(ep)
-- 
cgit v1.2.1