Fit one or more kinetic models with one or more state variables to one or more datasets

This function calls mkinfit on all combinations of models and datasets specified in its first two arguments.

mmkin(models, datasets,
      cores = round(detectCores()/2), cluster = NULL, ...)

Arguments

models	Either a character vector of shorthand names ("SFO", "FOMC", "DFOP", "HS", "SFORB"), or an optionally named list of mkinmod objects.
datasets	An optionally named list of datasets suitable as observed data for mkinfit.
cores	The number of cores to be used for multicore processing. This is only used when the cluster argument is NULL. On Windows machines, cores > 1 is not supported, you need to use the cluster argument to use multiple logical processors.
cluster	A cluster as returned by makeCluster to be used for parallel execution.
…	Further arguments that will be passed to mkinfit.

Value

A matrix of mkinfit objects that can be indexed using the model and dataset names as row and column indices.

See also

[.mmkin for subsetting, plot.mmkin for plotting.

Examples

m_synth_SFO_lin <- mkinmod(parent = mkinsub("SFO", "M1"),
                           M1 = mkinsub("SFO", "M2"),
                           M2 = mkinsub("SFO"), use_of_ff = "max")
#> Error in mkinmod(parent = mkinsub("SFO", "M1"), M1 = mkinsub("SFO", "M2"),     M2 = mkinsub("SFO"), use_of_ff = "max"): konnte Funktion "mkinmod" nicht finden

m_synth_FOMC_lin <- mkinmod(parent = mkinsub("FOMC", "M1"),
                            M1 = mkinsub("SFO", "M2"),
                            M2 = mkinsub("SFO"), use_of_ff = "max")
#> Error in mkinmod(parent = mkinsub("FOMC", "M1"), M1 = mkinsub("SFO", "M2"),     M2 = mkinsub("SFO"), use_of_ff = "max"): konnte Funktion "mkinmod" nicht finden

models <- list(SFO_lin = m_synth_SFO_lin, FOMC_lin = m_synth_FOMC_lin)
#> Error in eval(expr, envir, enclos): Objekt 'm_synth_SFO_lin' nicht gefunden
datasets <- lapply(synthetic_data_for_UBA_2014[1:3], function(x) x$data)
names(datasets) <- paste("Dataset", 1:3)

time_default <- system.time(fits.0 <- mmkin(models, datasets, quiet = TRUE))
#> Error in mmkin(models, datasets, quiet = TRUE): konnte Funktion "mmkin" nicht finden
#> Timing stopped at: 0 0 0.001
time_1 <- system.time(fits.4 <- mmkin(models, datasets, cores = 1, quiet = TRUE))
#> Error in mmkin(models, datasets, cores = 1, quiet = TRUE): konnte Funktion "mmkin" nicht finden
#> Timing stopped at: 0 0 0

time_default
#> Error in eval(expr, envir, enclos): Objekt 'time_default' nicht gefunden
time_1
#> Error in eval(expr, envir, enclos): Objekt 'time_1' nicht gefunden

endpoints(fits.0[["SFO_lin", 2]])
#> Error in endpoints(fits.0[["SFO_lin", 2]]): konnte Funktion "endpoints" nicht finden

# plot.mkinfit handles rows or columns of mmkin result objects
plot(fits.0[1, ])
#> Error in plot(fits.0[1, ]): Objekt 'fits.0' nicht gefunden
plot(fits.0[1, ], obs_var = c("M1", "M2"))
#> Error in plot(fits.0[1, ], obs_var = c("M1", "M2")): Objekt 'fits.0' nicht gefunden
plot(fits.0[, 1])
#> Error in plot(fits.0[, 1]): Objekt 'fits.0' nicht gefunden
# Use double brackets to extract a single mkinfit object, which will be plotted
# by plot.mkinfit and can be plotted using plot_sep
plot(fits.0[[1, 1]], sep_obs = TRUE, show_residuals = TRUE, show_errmin = TRUE)
#> Error in plot(fits.0[[1, 1]], sep_obs = TRUE, show_residuals = TRUE, show_errmin = TRUE): Objekt 'fits.0' nicht gefunden
plot_sep(fits.0[[1, 1]])
#> Error in plot_sep(fits.0[[1, 1]]): konnte Funktion "plot_sep" nicht finden
# Plotting with mmkin (single brackets, extracting an mmkin object) does not
# allow to plot the observed variables separately
plot(fits.0[1, 1])
#> Error in plot(fits.0[1, 1]): Objekt 'fits.0' nicht gefunden