Decomposition of gene expression variability according to BASiCS
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
Vector of characters containing gene names (must be in the same order as in the Counts matrix).
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