pfMp {MOLHD} | R Documentation |
Generate the Pareto front for optimal Latin Hypercube Designs based on both the Maximin and the MaxPro criteria.
Description
Generate the Pareto front for optimal Latin Hypercube Designs based on both the Maximin and the MaxPro criteria.
Usage
pfMp(n, p, crlim, nstarts = 1, times = 300, maxiter = 1e+06, temp0 = 0,
wtset = cbind(c(1, 0), c(0.8, 0.2), c(0.6, 0.4), c(0.4, 0.6), c(0.2, 0.8),
c(0, 1)))
Arguments
n |
number of runs desired |
p |
number of design factors desired |
crlim |
a matrix saving the best and worst values for each criterion to be used for defining the scaling choices for converting the natural criteria values onto a desirability scale between 0 and 1. Each column corresponds to one criterion. The best and worst values are recommended based on the values from each single criterion search. It is recommended that slightly wider range is used for defining the scaling choice for the Pareto front search. |
nstarts |
Optional, default is "1". The number of random starts |
times |
Optional, default is "300". The maximum number of non-improving searches allowed before terminating the search. Lower this parameter if you expect the search to converge faster. |
maxiter |
Optional, default is "1e+06".The maximum total number of iterations. Lower this number if the design is prohibitively large and you want to terminate the search earlier to report the best design found. |
temp0 |
Initial temperature for simualted annealing |
wtset |
Optional, default is "cbind(c(1,0),c(0.8,0.2),c(0.6,0.4),c(0.4,0.6),c(0.2,0.8),c(0,1))". The set of weight combinations to guide the search in varied directions. Each column is a weight vector that guides the search in a certain direction. |
Details
This function is to search the Pareto front and the Pareto set of designs based on the Maximin and Maxpro criteria. Each design on Pareto front is not dominated by any other design.This function utilizes a version of simulated annealing algorithm and several computational shortcuts to efficiently generate the optimal Latin hypercube designs. Choose a lower maximum limit of the criteria but high enough for Pareto front designs will save the computing time.
Value
pfdes |
The column bind of Pareto front designs whose criteria values are on the Pareto front. |
pfvals |
The Pareto front of criteria values based on the Maximin and MaxPro criteria. Columns are the optimization criteria. |
time_rec |
Time to complete the search |
Examples
#Generate the Pareto designs of maximin distance and maximum projection for LHD(10,5)
## Not run:
D1=MmLHD(n=10,p=5,nstarts=30)
D2=mpLHD(n=10,p=5,nstarts=30)
Mmlim=c(D1$criterion-0.2,D1$criterion-0.2+2)
mplim=c(D2$criterion-2,D2$criterion-2+10)
crlim=cbind(Mmlim,mplim)
pf=pfMp(10,5,crlim,nstarts = 30)
pf$pfvals
pf$pfdes
pf$time_rec
## End(Not run)