From c3c0007bd76292d976bd7624bf72aa75e168128c Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 8 Nov 2022 15:44:39 +0100 Subject: Same SFORB simplification for saem.mmkin fits --- R/saem.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/R/saem.R b/R/saem.R index 09a4832f..87718603 100644 --- a/R/saem.R +++ b/R/saem.R @@ -413,7 +413,7 @@ saemix_model <- function(object, solution_type = "auto", k_1output <- exp(psi[id, 2]) t <- xidep[, "time"] - sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 + k_12 * k_21 - (k_12 + k_1output) * k_21) + sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 - k_1output * k_21) b1 = 0.5 * (k_12 + k_21 + k_1output) + sqrt_exp b2 = 0.5 * (k_12 + k_21 + k_1output) - sqrt_exp @@ -427,7 +427,7 @@ saemix_model <- function(object, solution_type = "auto", k_1output <- psi[id, 2] t <- xidep[, "time"] - sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 + k_12 * k_21 - (k_12 + k_1output) * k_21) + sqrt_exp = sqrt(1/4 * (k_12 + k_21 + k_1output)^2 - k_1output * k_21) b1 = 0.5 * (k_12 + k_21 + k_1output) + sqrt_exp b2 = 0.5 * (k_12 + k_21 + k_1output) - sqrt_exp -- cgit v1.2.1