https://github.com/cran/fields
Tip revision: 5c24a1f871f2c5d4b6da602c080c2b28df583875 authored by Douglas Nychka on 16 December 2016, 10:13:26 UTC
version 8.8
version 8.8
Tip revision: 5c24a1f
mKrigCheckXY.R
mKrigCheckXY <- function(x, y, weights, Z, na.rm)
{
#
# check for missing values in y or X.
#
# save logical indicating where there are NA's
# and check for NA's
#
ind <- is.na(y)
if (any(ind) & !na.rm) {
stop("Need to remove missing values or use: na.rm=TRUE in the call")
}
#
# coerce x to be a matrix
x <- as.matrix(x)
#
# coerce y to be a vector
#
y <- as.matrix(y)
#
#default weights ( reciprocal variance of errors).
#
if (is.null(weights))
weights <- rep(1, nrow(y))
#
# check that dimensions agree
#
if (nrow(y) != nrow(x)) {
stop(" length of y and number of rows of x differ")
}
if (nrow(y) != length(weights)) {
stop(" length of y and weights differ")
}
# if Z is not NULL coerce to be a matrix
# and check # of rows
if (!is.null(Z)) {
if (!is.matrix(Z)) {
Z <- as.matrix(Z)
}
if (length(y) != nrow(Z)) {
stop(" length of y and number of rows of Z differ")
}
}
# if NAs can be removed then remove them and warn the user
if (na.rm) {
ind <- is.na(y)
if(all(ind)){
stop("Oops! All y values are missing!")
}
if (any(ind)) {
y <- y[!ind]
x <- as.matrix(x[!ind, ])
if (!is.null(Z)) {
Z <- as.matrix(Z[!ind, ])
}
weights <- weights[!ind]
}
}
#
# check for NA's in x matrix -- there should not be any !
if (any(c(is.na(x)))) {
stop(" NA's in x matrix")
}
#
# check for NA's in Z matrix
if (!is.null(Z)) {
if (any(c(is.na(Z)))) {
stop(" NA's in Z matrix")
}
}
# save x, weights and y w/o NAs
N <- length(y)
return(list(N = N, y = y, x = x, weights = weights, Z = Z,
NA.ind = ind) )
}