aboutsummaryrefslogtreecommitdiff
path: root/R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-11-08 15:44:39 +0100
committerJohannes Ranke <jranke@uni-bremen.de>2022-11-08 15:44:39 +0100
commitc3c0007bd76292d976bd7624bf72aa75e168128c (patch)
treef21df450f3ffad1cd644f04becd41a21e7fb4825 /R
parent56415f3c0638725aaa77db2c02eeaf9d4983452b (diff)
Same SFORB simplification for saem.mmkin fits
Diffstat (limited to '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