# HG changeset patch
# User davidecangelosi
# Date 1559040823 14400
# Node ID 9227f7b1c50b1c4346cbd23bb5e0f90e16ca935e
# Parent  8503457bd0c4bb6076205cafb4c0a6cc82b486bf
planemo upload for repository https://github.com/igg-molecular-biology-lab/pipe-t.git commit ff16cfc0a6c2b83c1e32d1204e91f3a20d6dc231

diff -r 8503457bd0c4 -r 9227f7b1c50b pipe-t.R
--- a/pipe-t.R	Tue May 28 04:35:18 2019 -0400
+++ b/pipe-t.R	Tue May 28 06:53:43 2019 -0400
@@ -694,7 +694,7 @@
 			#						xx <- log2(subset(x, x<Ct.max))
 			#						2^mean(xx)})
       geo.mean	<- apply(data, 2, function(x) {
-									xx <- subset(x, x<=Ct.max)
+									xx <- subset(x, x<Ct.max)
 									geometric.mean(xx)})
 			# Which sample to scale to
 			#if (missing(geo.mean.ref))
@@ -710,7 +710,7 @@
 		} # switch
     ,globalmean = {
         glo <- apply(data, 2, function(x) {
-            xx <- subset(x, x <= Ct.max)
+            xx <- subset(x, x <Ct.max)
             mean(xx)
         })
         data.norm <- t(t(data) - glo)