https://github.com/cran/fields
Tip revision: 1e85525d7f4e727af643137ca8803bae3be1dd6d authored by Douglas Nychka on 16 May 2014, 16:10:31 UTC
version 7.1
version 7.1
Tip revision: 1e85525
qr.q2ty.R
# fields, Tools for spatial data
# Copyright 2004-2013, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
"qr.q2ty" <- function(qr, y) {
if (!is.matrix(y)) {
y <- as.matrix(y)
}
dy <- dim(y)
dq <- dim(qr$qr)
rank <- qr$rank
if (dy[1] != dq[1])
stop("y and qr$qr should have same number of rows")
qr.qty(qr, y)[(rank + 1):dy[1], ]
}