isothermal_OED_limit {bioOED}R Documentation

OED of isothermal microbial inactivation with detection limit

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)


[Package bioOED version 0.2.1 Index]