https://github.com/cran/fields
Revision a968cb727c3b72a67d1ddc6588c5edf9c0e3723a authored by Doug Nychka on 21 April 2013, 07:56:06 UTC, committed by cran-robot on 21 April 2013, 07:56:06 UTC
1 parent 759a1e4
Tip revision: a968cb727c3b72a67d1ddc6588c5edf9c0e3723a authored by Doug Nychka on 21 April 2013, 07:56:06 UTC
version 6.7.6
version 6.7.6
Tip revision: a968cb7
fields.derivative.poly.R
# fields, Tools for spatial data
# Copyright 2004-2011, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
fields.derivative.poly <- function(x, m, dcoef) {
# dimension of x locations
# goal is find partial derivative matrix
d <- ncol(x)
out <- fields.mkpoly(rbind(x[1, ]), m)
ptab <- attr(out, "ptab")
if (nrow(ptab) != length(dcoef)) {
stop(" rows of ptab not equal to length of dcoef")
}
hold <- matrix(NA, ncol = d, nrow = nrow(x))
for (k in 1:d) {
nonzero <- ptab[, k] != 0
ptemp <- matrix(ptab[nonzero, ], ncol = d)
dtemp <- dcoef[nonzero]
dtemp <- dtemp * ptemp[, k]
ptemp[, k] <- ptemp[, k] - 1
hold[, k] <- fields.evlpoly2(x, dtemp, ptemp)
}
return(hold)
}
Computing file changes ...