aboutsummaryrefslogtreecommitdiff
path: root/R/endpoints.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2020-09-01 08:05:43 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2020-09-01 08:05:43 +0200
commit6fbb2341f998b322f21fe4e7cd655e0ca61cdbe4 (patch)
tree50f25047bb108a9185ed65bac76b8392eee33c9b /R/endpoints.R
parenta7e209b553e2e99ffc9aa2b003f9b73b9ac4b985 (diff)
Back-calculated DT50 values for biphasic models
Diffstat (limited to 'R/endpoints.R')
-rw-r--r--R/endpoints.R11
1 files changed, 10 insertions, 1 deletions
diff --git a/R/endpoints.R b/R/endpoints.R
index cccd2ff7..e4813db9 100644
--- a/R/endpoints.R
+++ b/R/endpoints.R
@@ -22,8 +22,10 @@
#' fit <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE)
#' endpoints(fit)
#' \dontrun{
-#' fit_2 <- mkinfit("SFORB", FOCUS_2006_C, quiet = TRUE)
+#' fit_2 <- mkinfit("DFOP", FOCUS_2006_C, quiet = TRUE)
#' endpoints(fit_2)
+#' fit_3 <- mkinfit("SFORB", FOCUS_2006_C, quiet = TRUE)
+#' endpoints(fit_3)
#' }
#'
#' @export
@@ -135,7 +137,9 @@ endpoints <- function(fit) {
silent = TRUE)
if (inherits(DT50, "try-error")) DT50 = NA
if (inherits(DT90, "try-error")) DT90 = NA
+ DT50_back = DT90 / (log(10)/log(2)) # Backcalculated DT50 as recommended in FOCUS 2011
+ ep$distimes[obs_var, c("DT50back")] = DT50_back
ep$distimes[obs_var, c("DT50_k1")] = DT50_k1
ep$distimes[obs_var, c("DT50_k2")] = DT50_k2
}
@@ -152,8 +156,10 @@ endpoints <- function(fit) {
}
DT50 <- DTx(50)
DT90 <- DTx(90)
+ DT50_back = DT90 / (log(10)/log(2)) # Backcalculated DT50 as recommended in FOCUS 2011
DT50_k1 = log(2)/k1
DT50_k2 = log(2)/k2
+ ep$distimes[obs_var, c("DT50back")] = DT50_back
ep$distimes[obs_var, c("DT50_k1")] = DT50_k1
ep$distimes[obs_var, c("DT50_k2")] = DT50_k2
}
@@ -191,6 +197,8 @@ endpoints <- function(fit) {
DT90 = if (inherits(log_DT90, "try-error")) NA
else exp(log_DT90)
+ DT50_back = DT90 / (log(10)/log(2)) # Backcalculated DT50 as recommended in FOCUS 2011
+
for (k_out_name in k_out_names)
{
ep$ff[[sub("k_", "", k_out_name)]] = parms.all[[k_out_name]] / k_1output
@@ -200,6 +208,7 @@ endpoints <- function(fit) {
ep$SFORB[[paste(obs_var, "b1", sep="_")]] = b1
ep$SFORB[[paste(obs_var, "b2", sep="_")]] = b2
+ ep$distimes[obs_var, c("DT50back")] = DT50_back
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
}

Contact - Imprint