https://github.com/cran/fields
Tip revision: 8eab500c3dad2103092ff68706417414fe53e16b authored by Doug Nychka on 22 September 2009, 20:23:49 UTC
version 6.01
version 6.01
Tip revision: 8eab500
summary.gcv.sreg.R
# fields, Tools for spatial data
# Copyright 2004-2007, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
"summary.gcv.sreg" <- function(object, lambda, cost = 1,
nstep.cv = 20, offset = 0, verbose = TRUE, ...) {
out <- object
shat.pure.error <- out$shat.pure.error
pure.ss <- out$pure.ss
nt <- 2
np <- out$np
N <- out$N
out$cost <- cost
out$offset <- offset
lambda.est <- rep(NA, 6)
names(lambda.est) <- c("lambda", "trA", "GCV", "GCV.one",
"GCV.model", "shat")
#
# fill in stuff for this lambda
lambda.est[1] <- lambda
temp <- sreg.fit(lambda, out)
lambda.est[2] <- temp$trace
lambda.est[3] <- temp$gcv
lambda.est[4] <- temp$gcv.one
if (!is.na(shat.pure.error)) {
lambda.est[5] <- temp$gcv.model
}
lambda.est[6] <- temp$shat
lambda.est
}