[BioC] Help with flowCore

[Finak, Greg]

I would add one more caveat.. double check that the $SPILL matrix is not just the identity matrix (sometimes this is the case). If so you'll need the single stained controls and the unstained control to calculate the spillover matrix with spillover(), and then proceed as Mike suggests.

Greg Finak
Staff Scientist
Fred Hutchinson Cancer Research Center


Mike <wjiang2 at fhcrc.org> wrote:


Andrea,

We should compensate the data outside of read.flowSet call. Here is what
you can do:

#read the FCS files into a flowSet
fs <- read.flowSet(...)

#loop through flowset and do the compensation for each flowFrame
#and return the compensated flowFrame,fsApply call will construct a new
flowSet for you

fs_comp <- fsApply(fs,function(frame){
     #extract compensation matrix from keywords
     comp <- keyword(frame)$`SPILL`
     new_frame <- compensate(frame,comp)
     new_frame
})

Let me know if you have more questions.

Mike
On 04/03/2013 03:43 AM, Valle Andrea wrote:
> Dear Dr. Jiang,
>
> I'm a PhD student in Milan and I'm trying flowCore package. First of all let me tell you that I think it has been a great step forward in flow cytometry!
> I'm writing you cause I have an issue about compensation and, although I tried a lots of possible solutions, I didn't succeed. I have to say that I'm not an expert in R langauge since this is my first approach with this subject.
> When I load an FCS file it's is not compensated, although it include $SPILL matrix. Of course I can use "compensate" to apply the $SPILL matrix to that file. However I can't do it upon flowset creation. Once I created a flowset, the $SPILL keyword is no more included in each file.
> How can I create a flowset in which each file is automatically compensated with its own compensation?
>
> I attached a representative file
>
> Many many thanks!
>
> Andrea Valle
> --------------------------------------------------------------------------
> LA TUA CURA E' SCRITTA NEL TUO DNA. AL SAN RAFFAELE LA STIAMO REALIZZANDO.
> AIUTA LA RICERCA, DAI IL TUO 5XMILLE - CF: 07636600962
> info:www.5xmille at hsr.it - www.5xmille.org
>
> Disclaimer added by CodeTwo Exchange Rules 2007
> http://www.codetwo.com
>

_______________________________________________
Bioconductor mailing list
Bioconductor at r-project.org
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor