plotGame: Plot equilibrium search result (2-objectives only)
In GPGame: Solving Complex Game Problems using Gaussian Processes

Description Usage Arguments Value Examples

View source: R/plotGame.R

Plot equilibrium search result (2-objectives only)

plotGame(
  res,
  equilibrium = "NE",
  add = FALSE,
  UQ_eq = TRUE,
  simus = NULL,
  integcontrol = NULL,
  simucontrol = NULL,
  Nadir = NULL,
  Shadow = NULL,
  ncores = 1,
  calibcontrol = NULL
)

`res`	list returned by `solve_game`
`equilibrium`	either "`NE`" for Nash, "`KSE`" for Kalai-Smoridinsky and "`NKSE`" for Nash-Kalai-Smoridinsky
`add`	logical; if `TRUE` adds the first graphical output to an already existing plot; if `FALSE`, (default) starts a new plot
`UQ_eq`	logical; should simulations of the equilibrium be displayed?
`simus`	optional matrix of conditional simulation if `UQ_Eq` is `TRUE`
`integcontrol`	list with `n.s` element (maybe n.s should be returned by solve_game). See `solve_game`.
`simucontrol`	optional list for handling conditional simulations. See `solve_game`.
`Nadir, Shadow`	optional vectors of size `nobj`. Replaces the nadir point for `KSE`. If only a subset of values needs to be defined, the other coordinates can be set to `Inf` (resp. `-Inf`).
`ncores`	number of CPU available (> 1 makes mean parallel `TRUE`)
`calibcontrol`	an optional list for calibration problems, containing `target` a vector of target values for the objectives and `log` a Boolean stating if a log transformation should be used or not.

No value returned, called for visualization.

library(GPareto)
library(parallel)

# Turn off on Windows
parallel <- FALSE # TRUE
ncores <- 1
if(parallel) ncores <- detectCores()
cov.reestim <- TRUE
n.sim <- 20
n.ynew <- 20
IS <- TRUE
set.seed(1)

pb <- "P1" # 'P1' 'PDE' 'Diff'
fun <- P1

equilibrium = "NE"

d <- 2
nobj <- 2
n.init <- 20
n.ite <- 4
model.trend <- ~1
n.s <- rep(31, 2) #31
x.to.obj   <- c(1,2)
gridtype <- 'cartesian'
nsimPoints <- 800
ncandPoints <- 200
sur_window_filter <- NULL
sur_pnash_filter  <- NULL
Pnash_only_filter <- NULL
res <- solve_game(fun, equilibrium = equilibrium, crit = "sur", model = NULL, n.init=n.init,
  n.ite = n.ite, nobj=nobj, x.to.obj = x.to.obj, integcontrol=list(n.s=n.s, gridtype=gridtype),
  simucontrol=list(n.ynew=n.ynew, n.sim=n.sim, IS=IS), ncores = ncores, d = d,
  filtercontrol=list(filter=sur_window_filter, nsimPoints=nsimPoints, ncandPoints=ncandPoints),
  kmcontrol=list(model.trend=model.trend), trace=3,
  seed=1)
plotGame(res, equilibrium = equilibrium)

dom <- matrix(c(0,0,1,1),2)
plotGameGrid("P1", graphs = "objective", domain = dom, n.grid = 51, equilibrium = equilibrium)
plotGame(res, equilibrium = equilibrium, add = TRUE)