diff --git a/R/correlationMatrix.R b/R/correlationMatrix.R index 4f3db962ec732402d6fcdbc1bf498496c7cf4b30..3c21825f9f17d123ab30af455fe34e1142ce8142 100644 --- a/R/correlationMatrix.R +++ b/R/correlationMatrix.R @@ -138,6 +138,7 @@ 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) + testMat <- testMat[-which(apply(testMat, 1, sd) == 0, ] to_remove <- NULL for (i in 1:nrow(testMat)) { if (!(rownames(testMat)[i] %in% lrsid)) {