[Qpcr-commits] r123 - pkg/NormqPCR/R

[noreply at r-forge.r-project.org]

Author: jperkins
Date: 2010-09-09 19:35:35 +0200 (Thu, 09 Sep 2010)
New Revision: 123

Modified:
   pkg/NormqPCR/R/combineTechReps.R
Log:
Ridiculously difficult to find bug in combineTechReps, was grepping very lazily previously

Modified: pkg/NormqPCR/R/combineTechReps.R
===================================================================
--- pkg/NormqPCR/R/combineTechReps.R	2010-09-09 07:56:42 UTC (rev 122)
+++ pkg/NormqPCR/R/combineTechReps.R	2010-09-09 17:35:35 UTC (rev 123)
@@ -10,7 +10,10 @@
     newDetectors <- unique(gsub("_TechReps.\\d","", origDetectors))
     NewExpM <- matrix(nrow = length(newDetectors), ncol = dim(expM)[2], dimnames = list(newDetectors,colnames(expM)))
     for(detector in newDetectors){
+	cat(detector,"\t")
+        cat(expM[gsub("_TechReps.\\d","", origDetectors) %in% origDetectors],"\t")
       dValues <- colMeans(expM[grepl(detector, origDetectors),],na.rm=TRUE)
+	cat(dValues,"\n")
       NewExpM[detector,] <- dValues
     }
     NewExpM[is.na(NewExpM)] <- NA