Function to decompose total variability of gene expression into biological and technical components.
BASiCS_VarianceDecomp(Chain, OrderVariable = "BioVarGlobal", Plot = TRUE,
...)
Ordering variable for output.
Possible values: 'GeneName'
, 'BioVarGlobal'
,
'TechVarGlobal'
and 'ShotNoiseGlobal'
.
Default: OrderVariable = "BioVarGlobal"
.
If TRUE
, a barplot of the variance decomposition
(global and by batches, if any) is generated. Default: Plot = TRUE
.
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 component (overall across all cells)
TechVarGlobal
Percentage of variance explained by the technical component (overall across all cells)
ShotNoiseGlobal
Percentage of variance explained by the shot noise component (baseline Poisson noise, 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). PLoS Computational Biology.
# NOT RUN {
# See
help(BASiCS_MCMC)
# }
Run the code above in your browser using DataLab