# HG changeset patch
# User kyost
# Date 1525281601 14400
# Node ID fbfe7b785ea757bb9b1986007926c8c1de501195
# Parent  6f5ce4e865a590257005b72fb178c157933dbdde
Add test for 0 standard deviation.

diff -r 6f5ce4e865a5 -r fbfe7b785ea7 find_qPCR_regions.R
--- a/find_qPCR_regions.R	Sat Apr 28 19:45:57 2018 -0400
+++ b/find_qPCR_regions.R	Wed May 02 13:20:01 2018 -0400
@@ -57,6 +57,10 @@
 	
 	a_split <- asplit[[i]]
 	b_split <- bsplit[[i]]
+	
+	if (sd(t(a_split[1,5:m])) == 0 || sd(t(b_split[-1,5:m])) == 0 ){
+		next
+	}
 
 	#calculate correlation between total peak reads and spanning reads in each window
     corb <- data.frame(cor(t(a_split[1,5:m])