# Efficient exact design using the RC heuristic

### Description

Computes an efficient exact design under linear resource constraints using the RC heuristic.

### Usage

1 2 |

### Arguments

`F` |
The |

`b, A` |
The vector of length |

`w0` |
The |

`crit` |
The optimality criterion. Possible values are |

`R` |
The region of summation for the IV-optimality criterion. The argument |

`w1` |
The |

`kappa` |
A small non-negative perturbation parameter. |

`tab` |
A vector determining the regressor components to be printed with the resulting design.
This argument should be a subvector of |

`graph` |
A vector determining the regressor components to be plotted with the resulting design.
This argument should be a subvector of |

`t.max` |
The time limit for the computation. |

### Details

This is an implementation of the algorithm proposed by Harman et al. employing the tabu search principle, and related to
the DETMAX procedure; see References. The inequalities `A%*%w<=b`

,
`w0<=w`

with the specific properties mentioned above, form the so-called resource constraints which encompass many
practical restrictions on the design, always permit a feasible solution, and lead to a bounded set of feasible solutions.

The information matrix of `w1`

should preferably have the reciprocal condition number of at least `1e-5`

. Note that the floor
of an optimal approximate design (computed using od.SOCP) is sometimes a good initial design. Alternatively,
the initial design can be the result of a different optimal design procedure, such as od.IQP. Even if no initial design
is provided, the model should be non-singular in the sense that there *exists* an exact design `w`

with a well
conditioned information matrix, satisfying all constraints. If this requirement is not satisfied, the computation may fail,
or it may produce a deficient design.

If the criterion of IV-optimality is selected, the region `R`

should be chosen such that the associated matrix `L`

(see the help page of the function od.crit) is non-singular, preferably with a reciprocal condition number of at least `1e-5`

.
If this requirement is not satisfied, the computation may fail, or produce a deficient design.

The perturbation parameter `kappa`

can be used to add `n*m`

iid random numbers from the uniform distribution
in `[-kappa,kappa]`

to the elements of `F`

before the optimization is executed. This can be helpful for
generating a random design from the potentially large set of optimal or nearly-optimal designs. However,
the RC heuristic uses a tabu principle based on the criterion values of designs, therefore in some problems a nonzero
`kappa`

can be detrimental to the optimization process.

The procedure always returns a permissible design, but in some cases, especially if `t.max`

is too small,
the resulting design can be inefficient. The performance depends on the problem and on the hardware used, but in most
cases the function can compute a nearly-optimal exact design for a problem with a hundred design points
within minutes of computing time. Because this is a heuristic method, we advise the user to verify the quality of the
resulting design by comparing it to the result of an alternative method (such as od.IQP and od.MISOCP)
and/or by computing its efficiency relative to the corresponding optimal approximate design (computed using od.SOCP).
In the special case of the single constraint on the size, it is generally more efficient to use the functions od.KL, or
od.RCs.

### Value

A list with the following components:

`method` |
The method used for computing the design |

`w.best` |
The best design found. |

`Phi.best` |
The value of the criterion of |

`t.act` |
The actual time taken by the computation. |

### Author(s)

Radoslav Harman, Lenka Filova

### References

Harman R, Bachrata A, Filova L (2016): Heuristic construction of exact experimental designs under multiple resource constraints, Applied Stochastic Models in Business and Industry, Volume 32, pp. 3-17

### See Also

`od.IQP, od.MISOCP, od.SOCP, od.KL, od.RCs`

### Examples

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | ```
# Consider the spring balance weighing model with 6 items of unknown
# weight. Suppose that we have already performed one weighing of each
# item separately. We will compute an A-efficient augmentation designs
# with 20 additional weighings. Then we will compute A-efficient designs
# under the additional restriction that no item can be used more
# than 8 times.
# Create the matrix of regressors of the model and the design
# to be augmented.
F.sbw <- F.cube(~x1 + x2 + x3 + x4 + x5 + x6 - 1, rep(0, 6),
rep(1, 6), rep(2, 6))
w0 <- rep(0, 64); w0[apply(F.sbw, 1, sum)==1] <- 1
# Compute an exact A-efficient augmentation design with 26
# total weighings.
b.sbw <- 26; A.sbw <- matrix(1, nrow=1, ncol=64)
res <- od.RC(F.sbw, b.sbw, A.sbw, w0=w0, crit="A", tab=1:6, t.max = 30)
# There are many A-optimal designs for this problem, which is possible
# to see by running the line above several times with a very small
# non-zero kappa. Note that each of the A-optimal experiments uses each
# item exactly 11 times. Suppose that we can use each item at
# most 8 times.
# Create the constraints A.eight * w <= b.eight representing
# the restriction that no item can be used more than eight times
# in the entire experiment.
b.eight <- rep(8, 6); A.eight <- t(F.sbw)
# Compute an exact A-efficient design with 26 total weighings under
# all constraints defined above.
b.sbw <- c(b.eight, 26); A.sbw <- rbind(A.eight, rep(1,64))
res <- od.RC(F.sbw, b.sbw, A.sbw, w0=w0, crit="A", tab=1:6, t.max = 60)
# To find a lower bound on the true efficiency of the resulting design,
# let us compare it to the A-optimal approximate design.
res.approx <- od.SOCP(F.sbw, b.sbw, A.sbw, w0=w0, crit="A",
tab=1:6, t.max = 20)
res$Phi.best / res.approx$Phi.best
``` |