Description Usage Details References See Also Examples
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
.
mtkIshigamiEvaluator()
mtkNativeEvaluator(model="Ishigami")
mtkEvaluator(protocol = "R", site = "mtk", service = "Ishigami")
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
.
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èque R pour l'exploration numérique des modèles.
In: Analyse de sensibilité et exploration de modèles : Application aux sciences de la nature et de l'environnement
(R. Faivre, B. Iooss, S. Mahévas, D. Makowski, H. Monod, Eds). Editions Quae, Versailles.
help(Ishigami.factors),help(ishigami.fun, sensitivity)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 | ### 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.