aboutsummaryrefslogtreecommitdiff
path: root/man/saemix.Rd
diff options
context:
space:
mode:
Diffstat (limited to 'man/saemix.Rd')
-rw-r--r--man/saemix.Rd23
1 files changed, 18 insertions, 5 deletions
diff --git a/man/saemix.Rd b/man/saemix.Rd
index 23b0a4ad..b41796ca 100644
--- a/man/saemix.Rd
+++ b/man/saemix.Rd
@@ -5,13 +5,16 @@
\alias{saemix_data}
\title{Create saemix models from mmkin row objects}
\usage{
-saemix_model(object)
+saemix_model(object, cores = parallel::detectCores())
saemix_data(object, ...)
}
\arguments{
\item{object}{An mmkin row object containing several fits of the same model to different datasets}
+\item{cores}{The number of cores to be used for multicore processing.
+On Windows machines, cores > 1 is currently not supported.}
+
\item{\dots}{Further parameters passed to \link[saemix:saemixData]{saemix::saemixData}}
}
\value{
@@ -25,17 +28,27 @@ object for use with the saemix package. An mmkin row object is essentially a
list of mkinfit objects that have been obtained by fitting the same model to
a list of datasets.
}
+\details{
+Starting values for the fixed effects (population mean parameters, argument psi0 of
+\code{\link[saemix:saemixModel]{saemix::saemixModel()}} are the mean values of the parameters found using
+mmkin. Starting variances of the random effects (argument omega.init) are the
+variances of the deviations of the parameters from these mean values.
+}
\examples{
ds <- lapply(experimental_data_for_UBA_2019[6:10],
function(x) subset(x$data[c("name", "time", "value")]))
names(ds) <- paste("Dataset", 6:10)
sfo_sfo <- mkinmod(parent = mkinsub("SFO", "A1"),
A1 = mkinsub("SFO"))
-f_mmkin <- mmkin(list("SFO-SFO" = sfo_sfo), ds, quiet = TRUE, cores = 5)
+\dontrun{
+f_mmkin <- mmkin(list("SFO-SFO" = sfo_sfo), ds, quiet = TRUE)
+library(saemix)
m_saemix <- saemix_model(f_mmkin)
d_saemix <- saemix_data(f_mmkin)
-saemix_options <- list(seed = 123456, save = FALSE, save.graphs = FALSE)
-\dontrun{
- saemix(m_saemix, d_saemix, saemix_options)
+saemix_options <- list(seed = 123456,
+ save = FALSE, save.graphs = FALSE, displayProgress = FALSE,
+ nbiter.saemix = c(200, 80))
+f_saemix <- saemix(m_saemix, d_saemix, saemix_options)
+plot(f_saemix, plot.type = "convergence")
}
}

Contact - Imprint