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
print.summary.krig.image.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
"print.summary.krig.image" <- function(x, ...) {
digits <- x$digits
c1 <- "Number of Observations:"
c2 <- x$num.observation
c1 <- c(c1, "Degree of polynomial null space ( base model):")
c2 <- c(c2, x$m - 1)
c1 <- c(c1, "MLE sigma ")
c2 <- c(c2, format(signif(x$shat.MLE, digits)))
c1 <- c(c1, "MLE rho ")
c2 <- c(c2, format(signif(x$rhohat, digits)))
c1 <- c(c1, "Scale used for covariance (rho)")
c2 <- c(c2, signif(x$rho, digits))
c1 <- c(c1, "Scale used for nugget (sigma^2)")
c2 <- c(c2, signif(x$sigma2, digits))
c1 <- c(c1, "lambda (sigma2/rho)")
c2 <- c(c2, signif(x$lambda, digits))
sum <- cbind(c1, c2)
dimnames(sum) <- list(rep("", dim(sum)[1]), rep("", dim(sum)[2]))
res.quantile <- x$res.quantile
names(res.quantile) <- c("min", "1st Q", "median", "3rd Q",
"max")
cat("Call:\n")
dput(x$call)
print(sum, quote = FALSE)
if ((x$correlation.model)) {
cat("Y is standardized before spatial estimate is found",
fill = TRUE)
}
cat(" Residuals:", "\n")
print(signif(res.quantile, digits))
cat("Covariance function name:", x$cov.function, fill = TRUE)
invisible(x)
}