R/methods_plotTenXVars.R

In yavorska/TenXAnalysisPackage: Analysis of the 10x Genomics 1.3 Million Cells Data

setMethod(f = "plotTenXVars",
          signature = signature(geneName = "character"),

          definition = function(geneName){

            index <- which(meanDF$geneName == geneName)

            ggplot(data = NULL,
                   aes(x = 1:133, y = as.numeric(varDF[index, -(1:4)]),
                       color = c(rep("A", 69), rep("B", 64)))) +
              geom_point() +
              theme_classic() +
              xlab("Library") +
              ylab(paste0("Gene expression variance of ", geneName)) +
              theme(legend.position = "none")

          })