diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2022-11-08 15:44:39 +0100 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2022-11-08 15:44:39 +0100 |
| commit | c3c0007bd76292d976bd7624bf72aa75e168128c (patch) | |
| tree | f21df450f3ffad1cd644f04becd41a21e7fb4825 /R | |
| parent | 56415f3c0638725aaa77db2c02eeaf9d4983452b (diff) | |
Same SFORB simplification for saem.mmkin fits
Diffstat (limited to 'R')
| -rw-r--r-- | R/saem.R | 4 |
1 files changed, 2 insertions, 2 deletions
@@ -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 |