maxproLHD: Generate a MaxPro Latin-hypercube design

View source: R/maxpro.R

maxproLHDR Documentation

Generate a MaxPro Latin-hypercube design

Description

This function generates a MaxPro Latin-hypercube design.

Usage

maxproLHD(
  n,
  p,
  design = NULL,
  max.sa.iter = 1e+06,
  temp = 0,
  decay = 0.95,
  no.update.iter.max = 400,
  num.passes = 10,
  max.det.iter = 1e+06,
  method = "full",
  scaled = TRUE
)

Arguments

n

design size.

p

design dimension.

design

an initial LHD. If design=NULL, a random LHD is generated.

max.sa.iter

maximum number of swapping involved in the simulated annealing (SA) algorithm.

temp

initial temperature of the simulated annealing algorithm. If temp=0, it will be automatically determined.

decay

the temperature decay rate of simulated annealing.

no.update.iter.max

the maximum number of iterations where there is no update to the global optimum before SA stops.

num.passes

the maximum number of passes of the whole design matrix if deterministic swapping is used.

max.det.iter

maximum number of swapping involved in the deterministic swapping algorithm.

method

choice of "deterministic", "sa", or "full". If the method="full", the design is first optimized by SA and then deterministic swapping.

scaled

whether the design is scaled to unit hypercube. If scaled=FALSE, the design is represented by integer numbers from 1 to design size. Leave it as TRUE when no initial design is provided.

Details

maxproLHD generates a MaxPro Latin-hypercube design (Joseph, V. R., Gul, E., & Ba, S. 2015). The major difference with the MaxPro packages is that we have a deterministic swap algorithm, which can be enabled by setting method="deterministic" or method="full". For optimization details, see the detail section in customLHD.

Value

design

final design points.

total.iter

total number of swaps in the optimization.

criterion

final optimized criterion.

crit.hist

criterion history during the optimization process.

References

Joseph, V. R., Gul, E., & Ba, S. (2015). Maximum projection designs for computer experiments. Biometrika, 102(2), 371-380.

Examples

n = 20
p = 3
D = maxproLHD(n, p)


SFDesign documentation built on June 22, 2025, 1:06 a.m.