MSSearch {multiDoE}R Documentation

Local search algorithm for high quality design generation

Description

The MSSearch function can be used to obtain an optimal multi-stratum experimental design considering one or more optimality criteria, up to a maximum of six criteria simultaneously.
This function implements the procedure MS-Opt proposed by Sambo, Borrotti, Mylona e Gilmour (2016) as an extension of the Coordinate-Exchange (CE) algorithm for constructing approximately optimal designs. This innovative procedure is able to handle all possible multi-stratum experimental structures and, instead of minimizing a single objective function as in the original CE algorithm, it seeks to minimize the following scalarization of the objective functions for all considered criteria:

f_W = \sum_{c \in C}{\alpha_c f_c(d; \eta)=\overline{\alpha} \cdot \overline{f}},

with

\sum_{c \in C} \alpha_c = 1,

where C is the set of criteria to be minimized, f_c is the objective function for the c criterion and \overline{\alpha} is the vector that controls the relative weights of the objective functions.

Usage

MSSearch(msopt, alpha, ...)

Arguments

msopt

A list as returned by the MSOpt function.

alpha

A vector of weights, whose elements must sum to one. length(alpha) must be equal to the number of criteria considered, that is it must be equal to the length of the criteria element of msopt.

...

optional arguments (see Details).

Details

MSSearch by default does not apply any normalization to the individual objective functions f_c before the calculation of f_w is performed. However, it is possible to subject the vector of objective functions \overline{f} to the following transformation:

\overline{f}_{norm} = \frac{\overline{f} - CritTR}{CritSC},

by specifying CritTR and CritSC vectors as additional parameters, as described below.

Additional arguments can be specified as follows:

Value

MSSearch returns a list, whose elements are:

References

M. Borrotti and F. Sambo and K. Mylona and S. Gilmour. A multi-objective coordinate-exchange two-phase local search algorithm for multi-stratum experiments. Statistics & Computing, 2016.

Examples

library(multiDoE)

## To check the number of digits to be printed.
backup_options <- options()
options(digits = 10)

## Definition of parameters for experimental setup
facts <- list(1, 2:5)
units <- list(21, 2)
level <- 3
etas <- list(1)
model2 <- "quadratic"

## Single-objective optimization
criteria_S <- c('I')
msopt_S <- MSOpt(facts, units, level, etas, criteria_S, model2)

mssearch_S <- MSSearch(msopt_S, alpha = 1, "Restarts", 100)


## Multi-objective optimization
criteria_M <- c('Id', 'Ds', 'As')
msopt_M <- MSOpt(facts, units, level, etas, criteria_M, model2)

mssearch_M <- MSSearch(msopt_M, alpha = c(1/2, 1/4, 1/4), "Restarts", 100)


options(backup_options)

## To reduce the computational cost of MSSearch function, you may reduce the number of restarts.


[Package multiDoE version 0.9.4 Index]