# HG changeset patch # User Ido Tamir # Date 1376837816 -7200 # Node ID 7bdd29cdfed8a6066c7e5edb3cfb1d15d087a5b6 # Parent 1eb46757d8bb91a9e1f40d8c807cdd99d4ce8b3d more debug messages diff -r 1eb46757d8bb -r 7bdd29cdfed8 corr.R --- a/corr.R Sun Aug 18 09:42:48 2013 -0400 +++ b/corr.R Sun Aug 18 16:56:56 2013 +0200 @@ -30,7 +30,8 @@ ## adds mappable regions as 0 coverage to track covWith0 <- function(cov, mappable){ - c0 <- setdiff(mappable, cov) + print(paste("cov:",cov)) + c0 <- setdiff(mappable, cov) cus <- if(length(c0) > 0){ elementMetadata(c0)$score <- 0 sort(c(c0, cov)) @@ -49,7 +50,8 @@ ## correlation dist corrDist <- function(covRles, outnames){ - vl <- length(covRles) + print("calculating correlation") + vl <- length(covRles) v <- 1:vl o <- matrix(NA,vl,vl) colnames(o) <- outnames diff -r 1eb46757d8bb -r 7bdd29cdfed8 mkCorr.sh --- a/mkCorr.sh Sun Aug 18 09:42:48 2013 -0400 +++ b/mkCorr.sh Sun Aug 18 16:56:56 2013 +0200 @@ -4,7 +4,9 @@ #echo "running corr" #echo $PWD - +module () { + eval `/usr/bin/modulecmd bash $*` +} if [ -f /biosw/modules/init/modules-env.sh ]; then module load R/3.0.0