The Ishigami model is an example evaluator implemented in the native mtk. It corresponds to the Ishigami function described in Saltelli et al., 2000. The behavior of the model is influenced by three factors x1, x2, x3.
Usage
mtkIshigamiEvaluator()
mtkNativeEvaluator(model="Ishigami")
mtkEvaluator(protocol = "R", site = "mtk", service = "Ishigami")
Details
The implementation of the Ishigami model includes the object Ishigami.factors on the input factors and the class mtkIshigamiEvaluator to run the simulations.
In mtk, there are a few ways to build an evaluator of the Ishigami model, but we usually recommend the following class constructors: mtkIshigamiEvaluator , mtkNativeEvaluator.
References
T. Ishigami and T. Homma (1990). An importance quantification technique in uncertainty analysis for computer models,
In: International Symposium on Uncertainity Modelling and Analysis (ISUMA'90) (1990).
A. Saltelli, K. Chan and E. M. Scott (2000). Sensitivity Analysis. Wiley, New York.
J. Wang, H. Richard, R. Faivre, H. Monod (2013). Le package mtk, une biblioth<c3><83><c2><a8>que R pour l'exploration num<c3><83><c2><a9>rique des mod<c3><83><c2><a8>les.
In: Analyse de sensibilit<c3><83><c2><a9> et exploration de mod<c3><83><c2><a8>les : Application aux sciences de la nature et de l'environnement
(R. Faivre, B. Iooss, S. Mah<c3><83><c2><a9>vas, D. Makowski, H. Monod, Eds). Editions Quae, Versailles.
### Run simulations of the "Ishigami" model
### for a random sample of input combinations
## Example I: by using the class constructor: mtkIshigamiEvaluator()
#
# Input the factors used in the "Ishigami" model
data(Ishigami.factors)
# Build the workflow:
# 1) specify the design process
exp1.designer <- mtkNativeDesigner(design = "BasicMonteCarlo",
information = list(size=20) )
# 2) specify the evaluation process;
exp1.evaluator <- mtkIshigamiEvaluator()
# 3) specify the workflow
exp1 <- mtkExpWorkflow(expFactors = Ishigami.factors,
processesVector = c(design=exp1.designer,
evaluate=exp1.evaluator) )
# Run the workflow and report the results.
run(exp1)
print(exp1)
## Example II: by using the class constructor: mtkNativeEvaluator()
# Generate the Ishigami input factors
data(Ishigami.factors)
# Build the workflow:
# 1) specify the design process
exp1.designer <- mtkNativeDesigner(design = "BasicMonteCarlo",
information = list(size=20) )
# 2) specify the evaluation process;
exp1.evaluator <- mtkNativeEvaluator(model="Ishigami")
# 3) specify the workflow
exp1 <- mtkExpWorkflow(expFactors = Ishigami.factors,
processesVector = c(design=exp1.designer, evaluate=exp1.evaluator) )
# Run the workflow and report the results.
run(exp1)
print(exp1)
## Example III: by using the generic class constructor: mtkEvaluator()
# Generate the Ishigami input factors
data(Ishigami.factors)
# Build the workflow:
# 1) specify the design process
exp1.designer <- mtkNativeDesigner(
design = "BasicMonteCarlo", information = list(size=20) )
# 2) specify the evaluation process;
exp1.evaluator <- mtkEvaluator(protocol = "R", site = "mtk", service = "Ishigami")
# 3) specify the workflow
exp1 <- mtkExpWorkflow(expFactors = Ishigami.factors,
processesVector = c(design=exp1.designer, evaluate=exp1.evaluator) )
# Run the workflow and report the results.
run(exp1)
print(exp1)