optionally use rfunctions::geninv instead of MASS::ginv for increased…

… performance
daroczig · Jun 14, 2014 · ac185e6 · ac185e6
1 parent d0d472d
commit ac185e6
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -14,3 +14,5 @@ License: AGPL-3
 Depends:
     Matrix,
     MASS
+Suggests:
+    rfunctions
diff --git a/R/unbalanced.R b/R/unbalanced.R
@@ -49,7 +49,10 @@ WithinTransformation1c <- function(i = m[, 1], j = m[, 2], t = m[, 3], value = m
     Qa <- t(D2) %*% Da
 
     ## geninv of Qa
-    Qai <- ginv(as.matrix(Qa))
+    if (require('rfunctions'))
+        Qai <- geninv(as.matrix(Qa))
+    else
+        Qai <- ginv(as.matrix(Qa))
     Qai <- Matrix(Qai)
 
     ## full value (added NA)
@@ -58,7 +61,6 @@ WithinTransformation1c <- function(i = m[, 1], j = m[, 2], t = m[, 3], value = m
     ## this should work on smaller matrices (formula 22)
     ## Pa <- (It - (D1 %*% dN2i %*% t(D1))) - (Da %*% Qai %*% t(Da))
 
-
     fia <- Qai %*% t(Da) %*% value
 
     t <- as.numeric(factor(t))