Update SEACR_1.2.R

Added a d0 term to percentile threshold part of loop so that write.table on line 144 does not fail due to an unspecified d0
FredHutch · Nov 5, 2019 · 169c97f · 169c97f
1 parent d30a4c6
commit 169c97f
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/SEACR_1.2.R b/SEACR_1.2.R
@@ -137,6 +137,7 @@ if(is.na(numtest)){ ## If 2nd field is a bedgraph, calculate empirical threshold
 	ctrl<-as.vector(as.numeric(paste(0,argsL$ctrl,sep="")))
 	x0<-min(frame$values[frame$percentile <= ctrl[1]])
 	z0<-min(frame2$values[frame2$percentile <= ctrl[1]])
+	d0<-0
 	fdr<-ctrl[1] ## New for SEACR_1.1
 }
 invis <- gc(verbose=FALSE)