# HG changeset patch
# User modencode-dcc
# Date 1360296396 18000
# Node ID 495a6d033ca1a5ca752f12682e6b0612331e9a54
# Parent  c573810591798e1774d003d63640171536d0f157
Uploaded

diff -r c57381059179 -r 495a6d033ca1 run_spp.R
--- a/run_spp.R	Thu Feb 07 23:06:21 2013 -0500
+++ b/run_spp.R	Thu Feb 07 23:06:36 2013 -0500
@@ -860,14 +860,14 @@
 	# Write to narrowPeak file
 	if (!is.na(iparams$output.npeak.file)) {
 		write.narrowpeak.binding(narrow.peaks,iparams$output.npeak.file,margin=round(crosscorr$whs/2),npeaks=iparams$npeak)
-		system(paste('gzip -f ',iparams$output.npeak.file))
+		#system(paste('gzip -f ',iparams$output.npeak.file))
 	}
 	
 	# Compute and write regionPeak file
 	if (!is.na(iparams$output.rpeak.file)) {
 		region.peaks <- add.broad.peak.regions(chip.data,control.data,narrow.peaks,window.size=max(50,round(crosscorr$whs/4)),z.thr=10)
 		write.narrowpeak.binding(region.peaks,iparams$output.rpeak.file,margin=round(crosscorr$whs/2),npeaks=iparams$npeak)
-		system(paste('gzip -f ',iparams$output.rpeak.file))
+		#system(paste('gzip -f ',iparams$output.rpeak.file))
 	}
 	
 	# Save Rdata file