Function to decompose total variability of gene expression into biological and technical components.
BASiCS_VarianceDecomp(Data, object, OrderVariable = "BioVarGlobal",
Plot = TRUE, ...)
an object of class BASiCS_Data-class
an object of class BASiCS_Chain-class
Ordering variable for output. Must take values in c("GeneNames", "BioVarGlobal", "TechVarGlobal", "ShotNoiseGlobal")
.
If TRUE
, a barplot of the variance decomposition (global and by batches, if any) is generated
Other arguments to be passed to barplot
A data.frame
whose first 4 columns correspond to
GeneName
Gene name (as indicated by user)
BioVarGlobal
Percentage of variance explained by a biological cell-to-cell heterogeneity component (overall across all cells)
TechVarGlobal
Percentage of variance explained by the technical cell-to-cell heterogeneity component (overall across all cells)
ShotNoiseGlobal
Percentage of variance explained by the shot noise component (baseline, overall across all cells)
If more than 1 batch of cells are being analysed, the remaining columns contain the corresponding variance decomposition calculated within each batch.
See vignette
Vallejos, Marioni and Richardson (2015). Bayesian Analysis of Single-Cell Sequencing data.
# NOT RUN {
# See
help(BASiCS_MCMC)
# }
Run the code above in your browser using DataLab