predLogProb: Predictive logarithmic probability of Kriging model

Description Usage Arguments Details Value Note Author(s) References Examples

View source: R/KDSNlossFunc.R

Description

Calculates the predictive logarithmic probability of Kriging model given the model parameters (likelihood function). It is used to choose better meta models in model-based optimization using Kriging. The observation used for prediction is left out in the training data and then estimated.

Usage

1
predLogProb(predMean, predSigma, y, X)

Arguments

predMean

Predicted leave one out mean of the Kriging model of all observations (numeric vector).

predSigma

Predicted leave one out variance of the Kriging model of all observations.

y

Numeric response vector of the outcome. Should be formated as a one column matrix.

X

Numeric design matrix of the covariates. All factors have to be prior encoded.

Details

Gaussian processes are conditionaly normal distributed and therefore the normal likelihood is used. The leave out parameters can be efficiently computed without performing leave one out cross-validation.

Value

Numeric Value of log-likelihood with leave-one-out parameters.

Note

This function is not intended to be called directly by the user. Should only be used by experienced users, who want to customize the model. It is called in the model selection process of the kernel deep stacking network, e.g. tuneMboLevelCvKDSN.

Author(s)

Thomas Welchowski welchow@imbie.meb.uni-bonn.de

References

Carl Edward Rasmussen and Christopher K. I. Williams, (2006), Gaussian Processes for Machine Learning Massachusetts Institute of Technology

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
library(DiceKriging)
# Generate design of experiments
design.fact <- expand.grid(x1=seq(0,1,length=4), x2=seq(0,1,length=4))
y <- apply(design.fact, 1, branin) 

# Estimate Kriging model
km_fit <- km(design=design.fact, response=y, control=list(trace=FALSE), 
nugget.estim=TRUE, iso=FALSE)

# Calculate leave one out parameters and performance measure
loo <- leaveOneOut.km(model=km_fit, type="UK", trend.reestim=TRUE)
predLogProbs <- predLogProb(predMean=loo$mean, predSigma=loo$sd^2, y=y, X=X)

kernDeepStackNet documentation built on May 2, 2019, 8:16 a.m.