Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

Table of Contents

...

SoftwareAdvantagesDisadvantagesCost
R

Fast

Opensource

HPC

Great community contribution

Scalable to large datasets

Most automated population clustering options

Learning curve

Support is from the community

FREE
RStudio

Same as R but with freiendlier user interface

Code can be deployed on HPC

Same as R

Learning curve (less than R due to GUI)

FREE
RToolsUsed to deploy newly published methodsComplicatedFREE
FlowJo

Commercial software

Powerful

Easy to use

Large number of advanced users

Works well

Cumbersome to implement new advanced analysis methods$$
SeqGeq

Useful for single cell workflows (RNA/DNA)

Commercial software

Powerful

Easy to use

Large number of advanced users

Many automated population clustering options


$$
Matlab

Powerful

Some published workflows use Matlab

Code is easier to understand compared to R (personal advantage, although this will likely change soon)

Flow cytometry analysis community is heavily developing in R compared to Matlab

...

In FlowJo, after importing your cytometry data, click the 'Check sample quality' button in the tools workspace. This will flag samples that should be checked. FlowJo plots the median values for each parameter over time and flags any that are outside of 2 standard deviations. Green is good. Any thing else should be reviewed. 

Active sample quality check - FlowAI in FlowJo

...

  1. Open RStudio (ensure you run as administrator, otherwise you may run into permission errors).
  2. Make a new project (a project can be considered a workspace where everything can be saved to)
  3. You want to download and load packages for the analysis that you will be doing (not efficient to load all the libraries, as RStudio is limited to 100, I learnt by experience, to see loaded libraries, type library())
  4. Some basic packages/libraries should be installed while working with flow data, this list isn't comprehensive but it contains some useful basicsfor the workflows listed, you can load the libraries that you will use

    Code Block
    #this will download the core Bioconductor packages
    source("https://bioconductor.org/biocLite.R")
    biocLite()
    
    #this will download flowCore (useful core for working with flow data)
    source("https://bioconductor.org/biocLite.R")
    biocLite("FlowSOMflowCore",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download flowUtils (another set of useful tools)
    source("https://bioconductor.org/biocLite.R")
    biocLite("flowCoreflowUtils",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download flowViz
    source("https://bioconductor.org/biocLite.R")
    biocLite("flowViz",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download flowSOM 
    source("https://bioconductor.org/biocLite.R")
    biocLite("flowUtilsFlowSOM",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download geneplotter 
    source("https://bioconductor.org/biocLite.R")
    biocLite("geneplotter",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download Seurat 
    source("https://bioconductor.org/biocLite.R")
    biocLite("Seurat",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download stringi
    source("https://bioconductor.org/biocLite.R")
    biocLite("stringi",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download yaml 
    source("https://bioconductor.org/biocLite.R")
    biocLite("yaml",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download dplyr
    source("https://bioconductor.org/biocLite.R")
    biocLite("dplyr",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download openCyto 
    source("https://bioconductor.org/biocLite.R")
    biocLite("openCyto",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download tsne 
    source("https://bioconductor.org/biocLite.R")
    biocLite("tsne",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download Rtsne 
    source("https://bioconductor.org/biocLite.R")
    biocLite("Rtsne",dependencies=TRUE,suppressUpdates=TRUE)
    
    #this will download edgeR 
    source("https://bioconductor.org/biocLite.R")
    biocLite("edgeR",dependencies=TRUE,suppressUpdates=TRUE)
    
    install.packages("dplyr",dependencies=TRUE,suppressUpdates=TRUE)
    install.packages("yaml",dependencies=TRUE,suppressUpdates=TRUE) 
    install.packages("devtools",dependencies=TRUE,suppressUpdates=TRUE)
    
    #use this to load downloaded packages
    library(flowCore)
    library(FlowSOM)
    library(flowViz)
    library(flowUtils)
    library(geneplotter)
    library(Seurat)
    library(dplyr)
    library(yaml)
    library(stringi)
    library(openCyto)
    library(tsne)
    library(Rtsne)


    Note: There are many Vignettes in the packages which are ever so helpful. Vignettes are help guides that can help to show you how to use different tools/functions.

  5. Some commonly used codes - good to get familiar with

    Code Block
    #gets to working directory
    getwd()
    #sets the working directory
    setwd('C:/Users/utopi/Desktop/testdata')
    #assign a file to a variable
    fileName <- "C:/Users/utopi/Desktop/testdata/sample.fcs"
    #assign a folder to a variable
    folderName <- "C:/Users/utopi/Desktop/testdata/"
    #read a FCS file to a variable using read.FCS from flowCore
    data1 = read.FCS('A1.fcs')
    #assign a variable an example data file i.e. from a vignette example
    fileName <- system.file("extdata","lymphocytes.fcs",package="FlowSOM")


  6. Examples to try...
    flowClust workflow in vignette
    flowSOM workflow in vignette 
    openCyto workflow in vignette
    Rtsne - https://github.com/lmweber/FlowSOM-Rtsne-example
    https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4980779/
    https://onlinelibrary.wiley.com/doi/full/10.1002/cyto.a.23030  - good comparison paper

    Lets analyse some data..
  7. Put your data in a folder on your computer
  8. In RStudio go to file and then new R script (we will be writing the script so that we can rerun it if needed)
  9. Lets assign the folder to a variable

...