Commit 91a0d85f authored by Vincent  LAVILLE's avatar Vincent LAVILLE

Update correlationMatrix.R

parent f80f60f3
......@@ -138,7 +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)
to_remove <- which(apply(testMat, 1, sd) == 0))
to_remove <- which(apply(testMat, 1, sd) == 0)
if (length(to_remove) > 0) {
testMat <- testMat[-to_remove, ]
genoMap <- genoMap[-to_remove, ]
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment