isothermal_OED_limit: OED of isothermal microbial inactivation with detection limit

In bioOED: Sensitivity Analysis and Optimum Experiment Design for Microbial Inactivation

Description

Calculates an Optimal Experiment for an isothermal microbial inactivation experiment considering the maximum duration of the experiment according to the detection limit.

Usage

isothermal_OED_limit(model, pars, limit, n_points, min_time, max_time,
  min_temp, max_temp, criterion = "D", opts = NULL, x_0 = NULL)

Arguments

model	character string defining the inactivation model to use.
pars	list defining the nominal model parameters.
limit	numerical value describing the maximum number of log-reductions that can be identified in the experiment limit = logDL - logN0, where DL is the detection limit.
n_points	numerical stating the number of data points.
min_time	numerical stating the lower limit for the time points.
max_time	numerical stating the upper limit for the time points.
min_temp	numerical stating the lower limit for the temperature.
max_temp	numerical stating the upper limit for the temperature.
criterion	character string defining the criterion to use.
opts	options for the MEIGO algorithm. By default, a maximum of 2000 function evaluations with local finish with the DHC algorithm (see help from MEIGO).
x_0	initial point for the MEIGO algorithm. By default, it is NULL.

Value

A MEIGO object

Examples

pars <- list(z = 4.2, D_R = 3.9, temp_ref = 55)
opts <- list(maxeval=2000,local_finish="DHC")
## Not run: 
OED <- isothermal_OED_limit("Bigelow", pars, n_points = 5, criterion = "E-mod", limit = 6,
                      min_time = 0, max_time = 100, min_temp = 52.5, max_temp = 60,
                      opts = opts)
plot(OED)

## End(Not run)

bioOED documentation built on Aug. 7, 2019, 5:03 p.m.