# HG changeset patch # User crs4 # Date 1489681321 14400 # Node ID e7616348a1be12fe735732c4a36520dfe4dd7a12 # Parent ffa4ab20d08956f032e969c8d1668fae882164b4 Uploaded diff -r ffa4ab20d089 -r e7616348a1be exomedepth.R --- 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