https://github.com/cran/fields
Tip revision: 6769ffc81115fbf0bf7d9c566cf7ac81be0049dc authored by Doug Nychka on 25 July 2005, 00:00:00 UTC
version 3.04
version 3.04
Tip revision: 6769ffc
summary.gcv.sreg.r
"summary.gcv.sreg" <-
function(object, lambda, cost = 1, nstep.cv = 20, offset = 0, verbose =
TRUE,...)
{
out<- object # hack S3
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
}