2014
DOI: 10.1016/j.csda.2013.11.017
|View full text |Cite
|
Sign up to set email alerts
|

Efficient optimization of the likelihood function in Gaussian process modelling

Abstract: Gaussian Process (GP) models are popular statistical surrogates used for emulating computationally expensive computer simulators. The quality of a GP model fit can be assessed by a goodness of fit measure based on optimized likelihood. Finding the global maximum of the likelihood function for a GP model is typically very challenging as the likelihood surface often has multiple local optima, and an explicit expression for the gradient of the likelihood function is typically unavailable. Previous methods for opt… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1
1

Citation Types

0
8
0

Year Published

2015
2015
2020
2020

Publication Types

Select...
4
2
1
1

Relationship

0
8

Authors

Journals

citations
Cited by 17 publications
(8 citation statements)
references
References 31 publications
0
8
0
Order By: Relevance
“…For instance, we should use a variety of test functions, repeat the simulations to average out the effect of initial design choice, find optimal (n 0 , n new ) combination in the sequential procedure, and so on. (2) The scalarization process should be further strengthened by using more informative discrepancy measure as compared to Euclidean distance. (3) Does this scalarization procedure affect the likeliness of finding the inverse solution?…”
Section: Discussionmentioning
confidence: 99%
See 2 more Smart Citations
“…For instance, we should use a variety of test functions, repeat the simulations to average out the effect of initial design choice, find optimal (n 0 , n new ) combination in the sequential procedure, and so on. (2) The scalarization process should be further strengthened by using more informative discrepancy measure as compared to Euclidean distance. (3) Does this scalarization procedure affect the likeliness of finding the inverse solution?…”
Section: Discussionmentioning
confidence: 99%
“…See [10] for more details. Popular implementations of the GP model like mlegp [11], GPfit ( [12]), GPmfit [13], and DiceKriging ( [14]) use some sort of numerical fix to overcome the computational instability issue. We used GPfit in R for all implementations of the GP model.…”
Section: Gaussian Process Modelsmentioning
confidence: 99%
See 1 more Smart Citation
“…MRST model: Finding the optimal drilling locations for production and injection wells in an oil reservoir is of utmost importance (see Onwunalu and Durlofsky, 2010). Butler et al (2014) used a Matlab Reservoir Simulator (MRST) (Lie et al, 2012) to generate the anticipated net present value (NPV) of the produced oil for a well to be drilled at a particular location. The goal here was to determine the configuration of wells that yields the best NPV.…”
Section: Real-life Computer Modelsmentioning
confidence: 99%
“…Chen et al (2015) discuss the importance of having extensive simulations and multiple data sets to validate a new statistical procedure. This working example and an incomplete list of other papers (Loeppky et al, 2009(Loeppky et al, , 2010Williams et al, 2011;Gramacy and Lee, 2012;Gramacy and Apley, 2014;Butler et al, 2014;Lam and Notz, 2008) all have a similar structure of a reasonably comprehensive set of simulations and number of methods that need to be assessed. The papers above typically show the results as a series of sideby-side boxplots or side-by-side dotcharts for each method, with one plot for each test function and sample size.…”
Section: Working Examplementioning
confidence: 99%