MOLHD (version 0.2)

pfMp: 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.

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

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.

Examples

# NOT RUN {
#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
# }