aboutsummaryrefslogtreecommitdiff
path: root/R/saem.R
diff options
context:
space:
mode:
Diffstat (limited to 'R/saem.R')
-rw-r--r--R/saem.R4
1 files 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

Contact - Imprint