battlefield/minmaxPareto_battle.R

In samuel-rosa/spsann: Optimization of Spatial Samples via Simulated Annealing

# Default example ##############################################################
rm(list = ls())
gc()
devtools::load_all()
# Example
require(sp)
data(meuse.grid)
candi <- meuse.grid[, 1:2]
covars <- meuse.grid[, c(1, 2)]

# General (greedy) cooling schedule
schedule <- scheduleSPSANN(
  initial.acceptance = c(0.01, 0.99), chains = 100,
  x.max = 1540, y.max = 2060, x.min = 0,
  y.min = 0, cellsize = 40)

# CORR
schedule$initial.temperature <- 0.1
set.seed(2001)
osc_corr <- optimCORR(
  points = 10, candi = candi, covars = covars,
  schedule = schedule)

# DIST
schedule$initial.temperature <- 0.1
set.seed(2001)
osc_dist <- optimDIST(
  points = 10, candi = candi, covars = covars,
  schedule = schedule)

# PPL
schedule$initial.temperature <- 0.1
set.seed(2001)
osc_ppl <- optimPPL(
  points = 10, candi = candi,
  schedule = schedule)

# MSSD
schedule$initial.temperature <- 0.1
set.seed(2001)
osc_mssd <- optimMSSD(
  points = 10, candi = candi,
  schedule = schedule)

# Pareto
pareto <- minmaxPareto(
  osc = list(
    DIST = osc_dist,
    CORR = osc_corr,
    PPL = osc_ppl,
    MSSD = osc_mssd),
  candi = candi,
  covars = covars)
round(pareto, 4)