Mercurial > repos > crs4 > exdepth

changeset 3:e7616348a1be draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author crs4
date Thu, 16 Mar 2017 12:22:01 -0400
parents ffa4ab20d089
children bb79be61ff3e
files exomedepth.R
diffstat 1 files changed, 3 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/exomedepth.R	Thu Mar 16 12:17:11 2017 -0400
+++ b/exomedepth.R	Thu Mar 16 12:22:01 2017 -0400
@@ -1,5 +1,5 @@
 # Load ExomeDepth library (without warnings)
-suppressMessages(library(ExomeDepth))
+library(ExomeDepth)
 
 # Import parameters
 args  <- commandArgs(trailingOnly=TRUE)
@@ -19,8 +19,8 @@
 bam_label   <- gsub(" ", "_", bam_label)
 
 for(i in 1:length(bam)){
-#  stopifnot(file.symlink(bam[i], paste(bam_label[i], "bam", sep=".")))
-#  stopifnot(file.symlink(bam_bai[i], paste(bam_label[i], "bam.bai", sep=".")))
+  stopifnot(file.symlink(bam[i], paste(bam_label[i], "bam", sep=".")))
+  stopifnot(file.symlink(bam_bai[i], paste(bam_label[i], "bam.bai", sep=".")))
 }
 
 # Generate read count data