https://github.com/cran/fields
Tip revision: c71fb7f6ffa323303affebf0e35a0070faa9c24d authored by Doug Nychka on 10 May 2004, 00:00:00 UTC
version 1.7.2
version 1.7.2
Tip revision: c71fb7f
Krig.find.gcvmin.r
"Krig.find.gcvmin" <-
function (info, lambda.grid, gcv, gcv.fun, tol, verbose = FALSE,
give.warnings=TRUE)
{
ind <- !is.na(gcv)
lambda.grid <- lambda.grid[ind]
gcv <- gcv[ind]
nstep.cv <- length(lambda.grid)
il <- order(gcv)[1]
lambda.gcv <- lambda.grid[il]
gcv.raw <- min(gcv)
if (verbose) {
cat("Results of coarse search lambda and GCV:",
lambda.grid[il],gcv.raw, fill=TRUE)
}
if ((il > 1) & (il < nstep.cv)) {
out <- golden.section.search(lambda.grid[il - 1], lambda.grid[il],
lambda.grid[il + 1], gcv.fun, f.extra = info, tol = tol *
gcv.raw)
return(out$x)
}
else {
if( give.warnings){
warning("GCV search gives a minumum at the endpoints of the grid search")}
return(lambda.gcv)
}
}