implement changes to store signalValue

as implemented in #1493 but putting these commits on `release_1.0.2` instead of master
GreenleafLab · Jun 30, 2022 · 1c2affd · 1c2affd
1 parent 6a0ec0c
commit 1c2affd
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/ReproduciblePeakSet.R b/R/ReproduciblePeakSet.R
@@ -789,8 +789,8 @@ addReproduciblePeakSet <- function(
 	run <- system2(pathToMacs2, cmd, wait=TRUE, stdout=NULL, stderr=NULL)
 
 	#Read Summits!
-	out <- data.table::fread(summitsFile, select = c(1,2,3,5))
-	out <- GRanges(out$V1, IRanges(out$V2 + 1, out$V3), score = out$V5)
+	out <- data.table::fread(narrowPeaksFile, select = c(1,2,5,7,10))
+	out <- GRanges(out$V1, IRanges(out$V2 + out$V10 + 1, out$V2 + out$V10 + 1), score = out$V5/10, signalValue = out$V7)
 
 	#Remove Files
 	r2 <- suppressWarnings(file.remove(summitsFile, narrowPeaksFile, xlsFile))