Violin Plot of Variance Fractions

Violin plot of variance fraction for each gene and each variable

Usage

plotVarPart(
  obj,
  col = c(ggColorHue(base::ncol(obj) - 4), "grey85", "grey65"),
  label.angle = 45,
  ncol = 3,
  ...
)

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

Arguments

obj

object returned by fitVarPart()

col

vector of colors

label.angle

angle of labels on x-axis

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)

plotVarPart(sortCols(vp))