Update correlationMatrix.R

R/correlationMatrix.R

@@ -138,9 +138,11 @@ getGenoCorMatrix <- function(lrsid, lchr, lpos, lrefall, pop, path = NULL, web =
         lind <- lind[lind$Paternal.ID == 0 & lind$Maternal.ID == 0, 2]
         
         testMat <- getAdditivelyCodedMatrix(genoMat,lind)
-        to_remove <- which(apply(testMat, 1, sd) == 0)
-        testMat <- testMat[-to_remove, ]
-        genoMap <- genoMap[-to_remove, ]
+        to_remove <- which(apply(testMat, 1, sd) == 0))
+        if (length(to_remove) > 0) {
+            testMat <- testMat[-to_remove, ]
+            genoMap <- genoMap[-to_remove, ]
+        }
         to_remove <- NULL
         for (i in 1:nrow(testMat)) {
             if (!(rownames(testMat)[i] %in% lrsid)) {