Bar Plot of Variance Fractions

Bar plot of variance fraction for each gene and each variable

Usage

plotPercentBars(
  obj,
  col = c(ggColorHue(base::ncol(obj) - 4), "grey85", "grey65"),
  ncol = 3,
  ...
)

# S4 method for class 'data.frame'
plotPercentBars(
  obj,
  col = c(ggColorHue(base::ncol(obj) - 4), "grey85", "grey65"),
  ncol = 3,
  ...
)

Arguments

obj

object returned by fitVarPart()

col

vector of colors

ncol

number of columns in the plot

...

additional arguments

Value

Violin plot showing variance fractions

Examples

library(SingleCellExperiment)

# Load example data
data(example_sce, package="muscat")
sce <- example_sce

# Compute library size for each cell
sce$libSize <- colSums(counts(sce))

# Specify regression formula and cell annotation 
form <- ~ group_id + (1|sample_id)
fit <- lucida(sce, form, "cluster_id", verbose=FALSE)
#> B cells 
#> CD14+ Monocytes 
#> CD4 T cells 
#> CD8 T cells 
#> FCGR3A+ Monocytes 

# Model with only intercept and random effect
form <- ~ (1|sample_id)
fit.null <- lucida(sce, form, "cluster_id", verbose=FALSE)
#> B cells 
#> CD14+ Monocytes 
#> CD4 T cells 
#> CD8 T cells 
#> FCGR3A+ Monocytes 

# Variance partitioning analysis
vp <- fitVarPart(fit, fit.null)

# Bar plots of a subset of genes
library(dplyr)
vp %>%
  sortCols %>%
  filter(ID %in% c('ISG15', 'C19orf10')) %>%
  plotPercentBars