[BioC] flowClust results in subset filters

Bastian Angermann AngerB at gmx.de
Fri Apr 24 19:30:37 CEST 2009


Hi,

I am using R2.9.0 and the 2.4 version of Bioconductor. When I try to run
a tmixfilter subsetted on a rectangular gate the result is of class logicalFilterResult thus I am unable to access the individual clusters:
#################################################
library(flowCore)
library(flowClust)
data(rituximab)

r_filter <-  rectangleGate("FSC.H"=c(100,500),"SSC.H"=c(100,500)) 
# gate boundaries are not meant to be meaningful. 
t_filter <- tmixFilter("tmix", c("FSC.H", "SSC.H"), K=3)
result <- filter(rituximab,t_filter %subset% r_filter)

class(result)
################################################

Is that a bug, missing feature or am I missing something? 
I'd like to be able to see the result of the 
tmixFilter in the context of the ancestor populations,
thus the %subset%. I'am working around that by adding an index vector
to the original flowFrame, to be able to do manual subsetting, but
it would be nice to to that without this hack.

Thanks a lot,
Bastian Angermann
--



More information about the Bioconductor mailing list