## Not run:
# #########################
# ## 1 response variable ##
# #########################
#
# in.name <- c("x1","x2")
# nlev <- c(20, 20)
# lower <- c(-2.048, -2.048)
# upper <- c(2.048, 2.048)
# out.name <- "y"
# weight <- 1
# C <- 10
# pr.mut <- c(0.1, 0.07, 0.04, rep(0.01, C-3))
#
# ## simulated problem (with known objective function)
# tn <- emmat0(in.name, nlev, lower, upper, out.name, nd = 10, fn1 = ackley)
#
# for(t in 1:(C-1))
# {
# tn <- emmatn(t, tn, na = 5, opt = "mn", weight, pr.mut = pr.mut,
# graph = "yes", fn1 = ackley)
# tn <- emmacheck(tn, graph = "no", fn1 = ackley)
# }
#
#
# ## applicative problem (with unknown objective function)
# tn <- emmat0(in.name, nlev, lower, upper, out.name, nd = 10)
# ## use the measured response values to manually fill 'tn$ypop'
# tn$ypop<-ackley(tn$xpop)
#
# for(t in 1:(C-1))
# {
# tn <- emmatn(t, tn, na = 5, opt = "mn", weight, pr.mut = pr.mut,
# graph = "yes")
# tn$ypop <- ackley(tn$xpop)
# tn <- emmacheck(tn, graph = "no")
# if(tn$add == 1) tn$ypop<-ackley(tn$xpop)
# }
# ## End(Not run)
##########################
## 2 response variables ##
##########################
in.name <- c("x1", "x2")
nlev <- c(20, 20)
lower <- c(-3, -3)
upper <- c(3, 3)
out.name <- c("y1", "y2")
weight <- c(0.2, 0.8)
C <- 10
pr.mut <- c(0.1, 0.07, 0.04, rep(0.01, C-3))
tn <- emmat0(in.name, nlev, lower, upper, out.name, nd = 10, fn1 = ackley,
fn2 = peaks)
for(t in 1:(C-1))
{
tn <- emmatn(t, tn, na = 5, opt = c("mn", "mx"), weight,
pr.mut = pr.mut, graph = "yes", fn1 = ackley, fn2 = peaks)
tn <- emmacheck(tn, graph = "no", fn1 = ackley, fn2 = peaks)
}
Run the code above in your browser using DataLab