BASiCS_VarianceDecomp: Decomposition of gene expression variability according to...
In BASiCS: Bayesian Analysis of Single-Cell Sequencing data

Description Usage Arguments Details Value Author(s) References See Also Examples

Function to decompose total variability of gene expression into biological and technical components.

BASiCS_VarianceDecomp(
  Chain,
  OrderVariable = c("BioVarGlobal", "GeneName", "TechVarGlobal", "ShotNoiseGlobal"),
  Plot = TRUE,
  main = "Overall variance decomposition",
  ylab = "% of variance",
  beside = FALSE,
  palette = "Set1",
  legend = c("Biological", "Technical", "Shot noise"),
  names.arg = if (nBatch == 1) "Overall" else c("Overall", paste("Batch",
    seq_len(nBatch)))
)

`Chain`	an object of class `BASiCS_Chain`
`OrderVariable`	Ordering variable for output. Possible values: `'GeneName'`, `'BioVarGlobal'`, `'TechVarGlobal'` and `'ShotNoiseGlobal'`. Default: `OrderVariable = "BioVarGlobal"`.
`Plot`	If `TRUE`, a barplot of the variance decomposition (global and by batches, if any) is generated. Default: `Plot = TRUE`.
`main`	Plot title.
`ylab`	y axis label.
`beside`	If `TRUE`, bars are placed beside each other. If `FALSE`, bars are stacked.
`palette`	Palette to be passed to `scale_fill_brewer` to create a discrete colour mapping.
`legend`	Labels for variance components.
`names.arg`	X axis labels.

See vignette

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.

Catalina A. Vallejos cnvallej@uc.cl

Vallejos, Marioni and Richardson (2015). PLoS Computational Biology.

BASiCS_Chain

# For illustration purposes we load a built-in 'BASiCS_Chain' object
# (obtained using the 'BASiCS_MCMC' function)
data(ChainSC)

VD <- BASiCS_VarianceDecomp(ChainSC)