Merge branch 'qrX-pivot'

vegandevs · Mar 13, 2024 · a88584b · a88584b
2 parents ef9dd75 + 58236af
commit a88584b
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/R/intersetcor.R b/R/intersetcor.R
@@ -6,7 +6,7 @@
     if (is.null(object$CCA) || !object$CCA$rank)
         stop("no constrained ordination or rank of constraints is zero")
     wa <- object$CCA$wa
-    X <- qr.X(object$CCA$QR)
+    X <- qr.X(object$CCA$QR, ncol = length(object$CCA$QR$pivot))
     ## remove conditions (partial terms)
     if (!is.null(object$pCCA)) {
         X <- X[, -seq_along(object$pCCA$envcentre), drop = FALSE]

diff --git a/R/vif.cca.R b/R/vif.cca.R
@@ -1,5 +1,5 @@
 `vif.cca` <-
-    function(object) 
+    function(object)
 {
     if (is.null(object$CCA))
         stop("can be used only with constrained ordination")
@@ -8,7 +8,7 @@
     names(out)[Q$pivot] <- colnames(Q$qr)
     rank <- Q$rank
     V <- chol2inv(Q$qr, size = rank)
-    X <- qr.X(Q)[, Q$pivot[1:rank], drop=FALSE]
+    X <- qr.X(Q, ncol = length(Q$pivot))[, Q$pivot[1:rank], drop=FALSE]
     Vi <- crossprod(X)
     v1 <- diag(V)
     v2 <- diag(Vi)