Commit d41414fe authored by Lena Hersemann's avatar Lena Hersemann

skip beta visualization until we find a less memory intensive way

parent 7efc2bf5
......@@ -307,7 +307,7 @@ stopifnot(t(ncol(beta)) == ncol(design[,-c(1, condColIdx), drop=FALSE]))
# visualize effect of batch correction
load_pack(pheatmap)
# beta %>% pheatmap( show_rownames=F)
beta %*% t(design[,-c(1, condColIdx), drop=FALSE]) %>% { colnames(.)<-exDesign$replicate; . } %>% pheatmap( show_rownames=F, main="batch correction offsets")
# beta %*% t(design[,-c(1, condColIdx), drop=FALSE]) %>% { colnames(.)<-exDesign$replicate; . } %>% pheatmap( show_rownames=F, main="batch correction offsets")
bcMatrix <- countMatrix - beta %*% t(design[,-c(1, condColIdx), drop=FALSE])
......
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