diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2019-02-21 14:34:45 +0100 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2019-02-21 14:50:27 +0100 |
| commit | d89e3d22eb9dc383897b09e9c5aa1b57f65cdbf0 (patch) | |
| tree | e81237fcbd996390eadd439295f2fb4b3874b0ab /R/endpoints.R | |
| parent | f134599b4d2cd3558b887b7f06faf1dfb599196e (diff) | |
Add the logistic model
Diffstat (limited to 'R/endpoints.R')
| -rw-r--r-- | R/endpoints.R | 8 |
1 files changed, 8 insertions, 0 deletions
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)
|