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
COR.r
"COR" <-
function (dat)
{
m <- ncol(dat)
M <- nrow(dat)
ntemp <- !is.na(dat)
ntemp <- t(ntemp) %*% (ntemp)
hold.mean <- apply(dat, 2, "mean", na.rm = TRUE)
temp <- (dat - matrix(c(hold.mean), ncol = m, byrow = TRUE,
nrow = M))
temp[is.na(temp)] <- 0
temp <- ifelse(ntemp > 1, t(temp) %*% (temp)/(ntemp - 1),
NA)
sd <- sqrt(diag(temp))
list(cor = t(t(temp)/sd)/sd, N = ntemp, sd = sd, mean = hold.mean)
}