inst/doc/PythonFunction.R

In caRamel: Automatic Calibration by Evolutionary Multi Objective Algorithm

## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(echo = TRUE)
knitr::opts_chunk$set(python.reticulate = FALSE)

## ----caRa, eval=F, echo=T-----------------------------------------------------
#  library(caRamel)

## ----reti, eval=F, echo=T-----------------------------------------------------
#  library(reticulate)

## ----load, eval=F, echo=T-----------------------------------------------------
#  use_python("/usr/local/bin/python3")
#  source_python("kursawe.py")

## ----wrap, eval=F, echo=T-----------------------------------------------------
#  wrapperFunction <- function(i) {
#    # load the package
#    library(reticulate)
#    # python path
#    use_python("/usr/local/bin/python3")
#    # source the Python function
#    source_python("kursawe.py")
#    # call the Python function and return the results
#    return(kursawe(x[i,]))
#  }

## ----kursawe_variable, eval=F, echo=T-----------------------------------------
#  nvar <- 3 # number of variables
#  bounds <- matrix(data = 1, nrow = nvar, ncol = 2) # upper and lower bounds
#  bounds[, 1] <- -5 * bounds[, 1]
#  bounds[, 2] <- 5 * bounds[, 2]

## ----kursawe_objectives, eval=F, echo=T---------------------------------------
#  nobj <- 2 # number of objectives
#  minmax <- c(FALSE, FALSE) # min and min

## ----kursawe_param, , eval=F, echo=T, fig.show="hide", results="hide"---------
#  popsize <- 100 # size of the genetic population
#  archsize <- 100 # size of the archive for the Pareto front
#  maxrun <- 1000 # maximum number of calls
#  prec <- matrix(1.e-3, nrow = 1, ncol = nobj) # accuracy for the convergence phase
#  
#  results <-
#    caRamel(nobj,
#            nvar,
#            minmax,
#            bounds,
#            wrapperFunction, # It's the wrapper function that will be called
#            popsize,
#            archsize,
#            maxrun,
#            prec)

## ----kursawe_OK_plot, eval=F, echo=T------------------------------------------
#  print(results$success==TRUE)
#  
#  plot(results$objectives[,1], results$objectives[,2], main="Kursawe Pareto front", xlab="Objective #1", ylab="Objective #2")

## ----kursawe_plot_conv, eval=F, echo=T----------------------------------------
#  matplot(results$save_crit[,1],cbind(results$save_crit[,2],results$save_crit[,3]),type="l",col=c("blue","red"), main="Convergence", xlab="Number of calls", ylab="Objectives values")