DiceEval-package | R Documentation |
Construction and evaluation of metamodels.
Package: | DiceEval |
Type: | Package |
Version: | 1.4 |
Date: | 2015-06-15 |
License: | GPL-3 |
This package is dedicated to the construction of metamodels. A validation procedure is also proposed using usual criteria (RMSE, MAE etc.) and cross-validation procedure. Moreover, graphical tools help to choose the best value for the penalty parameter of a stepwise or a PolyMARS model. Another routine is dedicated to the comparison of metamodels.
This work was conducted within the frame of the DICE (Deep Inside Computer Experiments) Consortium between ARMINES, Renault, EDF, IRSN, ONERA and TOTAL S.A. (http://emse.dice.fr/).
Functions gam
, mars
and polymars
are required for the construction of metamodels. km
provides Kriging models.
D. Dupuy & C. Helbert
Dupuy D., Helbert C., Franco J. (2015), DiceDesign and DiceEval: Two R-Packages for Design and Analysis of Computer Experiments, Journal of Statistical Software, 65(11), 1–38, https://www.jstatsoft.org/v65/i11/.
Friedman J. (1991), Multivariate Adaptative Regression Splines (invited paper), Annals of Statistics, 10/1, 1-141.
Hastie T. and Tibshirani R. (1990), Generalized Additive Models, Chapman and Hall, London.
Hastie T., Tibshirani R. and Friedman J. (2001), The Elements of Statistical Learning : Data Mining, Inference and Prediction, Springer.
Helbert C. and Dupuy D. (2007-09-26), Retour d'exp?riences sur m?tamod?les : partie th?orique, Livrable r?dig? dans le cadre du Consortium DICE.
Kooperberg C., Bose S. and Stone C.J. (1997), Polychotomous Regression, Journal of the American Statistical Association, 92 Issue 437, 117-127.
Rasmussen C.E. and Williams C.K.I. (2006), Gaussian Processes for Machine Learning, the MIT Press, www.GaussianProcess.org/gpml.
Stones C., Hansen M.H., Kooperberg C. and Truong Y.K. (1997), Polynomial Splines and their Tensor Products in Extended Linear Modeling, Annals of Statistics, 25/4, 1371-1470.
modelFit
, modelPredict
,
crossValidation
and modelComparison
Different space-filling designs can be found in the DiceDesign
package and we refer to the DiceKriging
package for the construction of kriging models. This package takes part of a toolbox inplemented during the Dice consortium.
## Not run: rm(list=ls()) # A 2D example Branin <- function(x1,x2) { x1 <- 1/2*(15*x1+5) x2 <- 15/2*(x2+1) (x2 - 5.1/(4*pi^2)*(x1^2) + 5/pi*x1 - 6)^2 + 10*(1 - 1/(8*pi))*cos(x1) + 10 } # A 2D uniform design with n points in [-1,1]^2 n <- 50 X <- matrix(runif(n*2,-1,1),ncol=2,nrow=n) Y <- Branin(X[,1],X[,2]) Z <- (Y-mean(Y))/sd(Y) # Construction of a PolyMARS model with a penalty parameter equal to 2 library(polspline) modPolyMARS <- modelFit(X,Z,type = "PolyMARS",gcv=2.2) # Prediction and comparison between the exact function and the predicted one xtest <- seq(-1, 1, length= 21) ytest <- seq(-1, 1, length= 21) Zreal <- outer(xtest, ytest, Branin) Zreal <- (Zreal-mean(Y))/sd(Y) Zpredict <- modelPredict(modPolyMARS,expand.grid(xtest,ytest)) m <- min(floor(Zreal),floor(Zpredict)) M <- max(ceiling(Zreal),ceiling(Zpredict)) persp(xtest, ytest, Zreal, theta = 30, phi = 30, expand = 0.5, col = "lightblue",main="Branin function",zlim=c(m,M), ticktype = "detailed") persp(xtest, ytest, matrix(Zpredict,nrow=length(xtest), ncol=length(ytest)), theta = 30, phi = 30, expand = 0.5, col = "lightblue",main="PolyMARS Model",zlab="Ypredict",zlim=c(m,M), ticktype = "detailed") # Comparison of models modelComparison(X,Y,type=c("Linear", "StepLinear","PolyMARS","Kriging"), formula=Y~X1+X2+X1:X2+I(X1^2)+I(X2^2),penalty=log(dim(X)[1]), gcv=4) # see also the demonstration example in dimension 5 (source: IRSN) demo(IRSN5D) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.