Used in the M-step (M_step.R). Note: predict(fit.krige, ...) has weird beheviour when making predictions for a single value in terms of standard error. We work around this issue in this implementation.
EI(theta, test.fun, fit.krige, theta.hash, dir)
The expected improvement.
Vector of coefficients.
Test function (cf. EstimationAlgorithmBei.R
).
Fitted Kriging model.
Tentative optimal value for theta, i.e., the largest or smallest feasible value for theta (if dir = 1 or dir = -1, respectively). A 'feasible value' is one that satisfies all moment restrictions.
Search direction. dir = 1
corresponds to looking for an
upper bound. dir = -1
corresponds to looking for a lower bound.