Il metodo raccomandato per calcolare il rango di una matrice in R sembra essere qr
:Il modo più veloce per calcolare il grado di matrice 2 * 2?
X <- matrix(c(1, 2, 3, 4), ncol = 2, byrow=T)
Y <- matrix(c(1.0, 1, 1, 1), ncol = 2, byrow=T)
qr(X)$rank
[1] 2
qr(Y)$rank
[1] 1
sono stato in grado di migliorare l'efficienza modificando questa funzione per il mio caso specifico:
qr2 <- function (x, tol = 1e-07) {
if (!is.double(x))
storage.mode(x) <- "double"
p <- as.integer(2)
n <- as.integer(2)
res <- .Fortran("dqrdc2", qr = x, n, n, p, as.double(tol),
rank = integer(1L), qraux = double(p), pivot = as.integer(1L:p),
double(2 * p), PACKAGE = "base")[c(1, 6, 7, 8)]
class(res) <- "qr"
res}
qr2(X)$rank
[1] 2
qr2(Y)$rank
[1] 1
library(microbenchmark)
microbenchmark(qr(X)$rank,qr2(X)$rank,times=1000)
Unit: microseconds
expr min lq median uq max
1 qr(X)$rank 41.577 44.041 45.580 46.812 1302.091
2 qr2(X)$rank 19.403 21.251 23.099 24.331 80.997
Utilizzando R , è possibile calcolare il rango di una matrice 2 * 2 ancora più veloce?
'fortune (98)' - beh, volte 4 Suppongo. – BenBarnes
@BenBarnes: Ho usato il tempo che ho risparmiato per guardare i lolcats sul web. –
Sto ottimizzando le prestazioni di una funzione che ho bisogno di eseguire alcune milioni di volte in una simulazione. 'qr' è usato all'interno di un ciclo while in questa funzione. Quindi, alla fine quei microsecondi alcuni fino a ore. – Roland