diff options
Diffstat (limited to 'man')
| -rw-r--r-- | man/saemix.Rd | 26 |
1 files changed, 25 insertions, 1 deletions
diff --git a/man/saemix.Rd b/man/saemix.Rd index eedbf537..5f96e84f 100644 --- a/man/saemix.Rd +++ b/man/saemix.Rd @@ -43,7 +43,7 @@ sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"), \dontrun{ f_mmkin <- mmkin(list("SFO-SFO" = sfo_sfo), ds, quiet = TRUE) library(saemix) -m_saemix <- saemix_model(f_mmkin) +m_saemix <- saemix_model(f_mmkin, cores = 1) d_saemix <- saemix_data(f_mmkin) saemix_options <- list(seed = 123456, save = FALSE, save.graphs = FALSE, displayProgress = FALSE, @@ -51,4 +51,28 @@ saemix_options <- list(seed = 123456, f_saemix <- saemix(m_saemix, d_saemix, saemix_options) plot(f_saemix, plot.type = "convergence") } +# Synthetic data with two-component error +sampling_times = c(0, 1, 3, 7, 14, 28, 60, 90, 120) +dt50_sfo_in <- c(800, 900, 1000, 1111.11, 1250) +k_in <- log(2) / dt50_sfo_in + +SFO <- mkinmod(parent = mkinsub("SFO")) + +pred_sfo <- function(k) { + mkinpredict(SFO, c(k_parent = k), + c(parent = 100), sampling_times) +} + +ds_sfo_mean <- lapply(k_in, pred_sfo) +set.seed(123456L) +ds_sfo_syn <- lapply(ds_sfo_mean, function(ds) { + add_err(ds, sdfunc = function(value) sqrt(1^2 + value^2 * 0.07^2), + n = 1)[[1]] + }) +f_mmkin_syn <- mmkin("SFO", ds_sfo_syn, error_model = "tc", quiet = TRUE) +m_saemix_tc <- saemix_model(f_mmkin_syn, cores = 1) +d_saemix_tc <- saemix_data(f_mmkin_syn) +f_saemix_tc <- saemix(m_saemix_tc, d_saemix_tc, saemix_options) +plot(f_saemix_tc, plot.type = "convergence") + } |