https://github.com/cran/fields
Raw File
Tip revision: 1e85525d7f4e727af643137ca8803bae3be1dd6d authored by Douglas Nychka on 16 May 2014, 16:10:31 UTC
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], ]
}
back to top