Browse code

Bug fix for the case stratifyBy="mm" in the readCdfUnits function. It used to return the pm probes due to some misplaced negation. Option "pm" and "pmmm" is unaffected by this fix.

git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/affxparser@20979 bc3139a8-67e5-0310-9ffc-ced21a209358

Kasper D. Hansen authored on 03/11/2006 18:44:40
Showing 1 changed files

... ...
@@ -225,7 +225,7 @@ readCdfUnits <- function(filename, units=NULL, readXY=TRUE, readBases=TRUE, read
225 225
         if (ngroup == 0)
226 226
           next;
227 227
 
228
-        pm <- !.subset2(.subset2(isPm, uu), gg);
228
+        pm <- .subset2(.subset2(isPm, uu), gg);
229 229
         mm <- (1:length(pm))[!pm]; # Note: which(!pm) is about 60% slower!
230 230
         for (kk in 1:ngroup) {
231 231
           group[[kk]] <- .subset(.subset2(group, kk), mm);