Local sensitivity for a MEMC model

Estimate the local effect of a parameter on a memc model output

Usage

memc_sensfunc(config, t, pars, ...)

Arguments

config

a memc model configuration object, either one of the pre-built configurations listed in model_configs or created using configure_model

t

vector of the time steps to run the model at

pars

vector of the parameters to test

...

additional arguments passed to FME::sensFun

Value

the results of the FME::sensFun

See also

Other sensitivity: format_sensout(), memc_sensrange()

Examples

if (FALSE) {
# Test the sensitivity of the MEND output for V.p, K.p, V.m
pars <- c("V.d" = 3.0e+00,"V.p" = 1.4e+01,"V.m" = 2.5e-01)
out <- memc_sensfunc(config = MEND_model, t = t, pars = pars)
pairs(out)
plot(out)
# Using the helper functions to make nice ggplots
to_plot <- format_sensout(out)
ggplot(data = to_plot) +
   geom_line(aes(time, value, color = parameter)) +
   facet_wrap("variable", scales = "free")
}