X-Git-Url: https://git.auder.net/?p=morpheus.git;a=blobdiff_plain;f=pkg%2FR%2FcomputeMu.R;fp=pkg%2FR%2FcomputeMu.R;h=f961fdebe40bb81fca5409ae1c4cc8aeb9099385;hp=bc52bb31df714ab1146c279eb1f22d239cf40105;hb=4b2f17bb108bab0f263619cfe00eabfb1e9b8860;hpb=35ffd7108a30ae1777dc9e5fc9d33442ad73548e

diff --git a/pkg/R/computeMu.R b/pkg/R/computeMu.R
index bc52bb3..f961fde 100644
--- a/pkg/R/computeMu.R
+++ b/pkg/R/computeMu.R
@@ -46,6 +46,8 @@ computeMu = function(X, Y, optargs=list())
     large_ratio <- ( abs(Σ[-d] / Σ[-1]) > 3 )
     K <- if (any(large_ratio)) max(2, which.min(large_ratio)) else d
   }
+  else if (K > d)
+    stop("K: integer >= 2, <= d")
 
   # Step 1: generate a family of d matrices to joint-diagonalize to increase robustness
   d = ncol(X)