1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
\name{mmkin}
\alias{mmkin}
\title{
Fit one or more kinetic models with one or more state variables to one or more datasets
}
\description{
This function calls \code{\link{mkinfit}} on all combinations of models and datasets
specified in its first two arguments.
}
\usage{
mmkin(models, datasets,
cores = round(detectCores()/2), cluster = NULL, ...)
}
\arguments{
\item{models}{
Either a character vector of shorthand names ("SFO", "FOMC", "DFOP",
"HS", "SFORB"), or an optionally named list of \code{\link{mkinmod}}
objects.
}
\item{datasets}{
An optionally named list of datasets suitable as observed data for
\code{\link{mkinfit}}.
}
\item{cores}{
The number of cores to be used for multicore processing. This is only
used when the \code{cluster} argument is \code{NULL}.
}
\item{cluster}{
A cluster as returned by \code{\link{makeCluster}} to be used for parallel
execution.
}
\item{\dots}{
Further arguments that will be passed to \code{\link{mkinfit}}.
}
}
\value{
A matrix of \code{\link{mkinfit}} objects that can be indexed using the model
and dataset names as row and column indices.
}
\seealso{
\code{\link{[.mmkin}} for subsetting, \code{\link{plot.mmkin}} for plotting.
}
\author{
Johannes Ranke
}
\examples{
\dontrun{
m_synth_SFO_lin <- mkinmod(parent = list(type = "SFO", to = "M1"),
M1 = list(type = "SFO", to = "M2"),
M2 = list(type = "SFO"), use_of_ff = "max")
m_synth_FOMC_lin <- mkinmod(parent = list(type = "FOMC", to = "M1"),
M1 = list(type = "SFO", to = "M2"),
M2 = list(type = "SFO"), use_of_ff = "max")
models <- list(SFO_lin = m_synth_SFO_lin, FOMC_lin = m_synth_FOMC_lin)
datasets <- lapply(synthetic_data_for_UBA_2014[1:3], function(x) x$data)
time_default <- system.time(fits.0 <- mmkin(models, datasets))
time_1 <- system.time(fits.1 <- mmkin(models, datasets, cores = 1))
time_default
time_1
endpoints(fits[["SFO_lin", 2]])
# Plot.mkinfit handles rows or columns of mmkin result objects
plot(fits.0[1, ])
# Double brackets to select a single mkinfit object, which will be
# plotted by plot.mkinfit
plot(fits.0[[1, 1]], sep_obs = TRUE, show_residuals = TRUE)
}
}
\keyword{ optimize }
|
