RKHSMetaMod-package: Set of Rcpp and R functions to produce a sequence of meta...

In RKHSMetaMod: Ridge Group Sparse Optimization Problem for Estimation of a Meta Model Based on Reproducing Kernel Hilbert Spaces

Description

Fits a meta model to an unknown model m by solving the ridge group sparse (or group lasso) optimization problem based on the reproducing kernel Hilbert spaces (RKHS), for the Gaussian regression model :

Y=m(X)+σ\varepsilon,

where variables X_{1},...,X_{d} are independent and uniformly distributed on [0,1] and are independent of \varepsilon's.

We define the ridge group sparse criteria by :

C(f_{0},θ)=\Vert Y-f_{0}I_{n}-∑_{v\in\mathcal{P}}K_{v}θ_{v}\Vert ^2+√{n}γ∑_{v\in\mathcal{P}}\Vert K_vθ_{v}\Vert +nμ∑_{v\in\mathcal{P}}\Vert K_{v}^{1/2}θ_{v}\Vert,

and the group lasso criteria is obtained by setting γ=0 in the criteria above. We set μ_g=√{n}μ to be the group lasso penalty parameter.

For each pair of the penalty parameters (μ,γ) in the ridge group sparse criteria, one meta model, called RKHS meta model, is calculated. The RKHS meta model is an additive model with at most vMax groups. It satisfies the properties of the Hoeffding decomposition, and its terms estimate the terms in the Hoeffding decomposition of the function m.

These estimators are evaluated using a testing dataset. That is, the prediction error is calculated for each RKHS meta model and the one with the minimum prediction error is the "best" estimator for the true model m. It provides a function that estimates the empirical sensitivity indices of the "best" RKHS meta model as an approximation of the true sensitivity indices.

Details

Details.

Author(s)

Halaleh Kamari.

Maintainer: [e-mail]

References

Kamari, H., Huet, S. and Taupin, M.-L. (2019) RKHSMetaMod : An R package to estimate the Hoeffding decomposition of an unknown function by solving RKHS Ridge Group Sparse optimization problem. <arXiv:1905.13695>

See Also

RKHSMetaMod

Examples

d <- 3
n <- 50
library(lhs)
X <- maximinLHS(n, d)
c <- c(0.2,0.6,0.8)
F <- 1;for (a in 1:d) F <- F*(abs(4*X[,a]-2)+c[a])/(1+c[a])
epsilon <- rnorm(n,0,1);sigma <- 0.2
Y <- F + sigma*epsilon
Dmax <- 3
kernel <- "matern"
frc <- c(10,100)
gamma <- c(.5,.01,.001,0)
result <- RKHSMetMod(Y,X,kernel,Dmax,gamma,frc,FALSE)  

Example output

RKHSMetaMod documentation built on July 7, 2019, 1:07 a.m.