Check correlation between number of genic reads and genes and transcripts detected per cell

gistfile1.txt

## Read in summary file from Drop-seq pipe
infile = 

DGE_info <- read.table(infile,
                       sep="\t",
                       header=T)

## Correlation Number of genic reads vs num of genes detected per cell
ggplot(DGE_info,aes(NUM_GENIC_READS,NUM_GENES,col=NUM_TRANSCRIPTS)) +
  geom_point()

## Correlation Number of genic reads vs num of transcripts (UMIs) detected per cell
ggplot(DGE_info,aes(NUM_GENIC_READS,NUM_TRANSCRIPTS,col=NUM_GENES)) +
  geom_point()