[BioC] flowClust results in subset filters
Florian Hahne
fhahne at fhcrc.org
Mon Apr 27 07:49:02 CEST 2009
Hi Bastian,
this is indeed a feature which wasn't implemented yet. I fixed the
code and everything seems to be working as intended now, but I want to
give it some more rigorous testing before submitting. I am planning to
push that back into the release branch since it is more of a bug fix
than feature addition.
If you want to be more independent of the types of comparisons between
gating steps, you might want to take a look at the workFlow stuff (as
described in the package vignette). There you can compute summary
statistics for a particular gating operation based on another
reference view in the workFlow (via the additional 'reference'
argument of the summary method)
bw,
Florian
On 24.04.2009, at 10:30, Bastian Angermann wrote:
> 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
> --
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list