PlotValues: Plot percent spliced-in (PSI) or gene expression values

In MARVEL: Revealing Splicing Dynamics at Single-Cell Resolution

Plot percent spliced-in (PSI) or gene expression values

Description

Plots percent spliced-in (PSI) or gene expression values across different groups of cells. This is a wrapper function for PlotValues.Exp and PlotValues.PSI.

Usage

PlotValues(
  MarvelObject,
  cell.group.list,
  feature,
  maintitle = "gene_short_name",
  xlabels.size = 8,
  level,
  min.cells = NULL,
  sigma.sq = 0.001,
  bimodal.adjust = NULL,
  seed = NULL,
  modality.column = "modality.bimodal.adj",
  scale.y.log = FALSE,
  max.cells.jitter = 10000,
  max.cells.jitter.seed = 1,
  cell.group.colors = NULL,
  point.alpha = 0.2
)

Arguments

MarvelObject	Marvel object. S3 object generated from TransformExpValues function.
cell.group.list	List of character strings. Each element of the list is a vector of cell IDs corresponding to a cell group. The name of the element will be the cell group label.
feature	Character string. tran_id or gene_id for plotting. Should match tran_id or gene_id column of MarvelObject$ValidatedSpliceFeature or MarvelObject$GeneFeature slot when level set to "splicing" or "gene", respectively.
maintitle	Character string. Column to use as plot main title as per MarvelObject$ValidatedSpliceFeature or MarvelObject$GeneFeature when level set to "splicing" or "gene", respectively. Default is "gene_short_name" column.
xlabels.size	Numeric value. Size of x-axis labels as per ggplot2 function. Default is 8.
level	Character string. Indicate "splicing" or "gene" for PSI or gene expression value plotting, respectively.
min.cells	Numeric value. Only applicable when level set to "splicing". The minimum no. of cells expressing the splicing event to be included for analysis.
sigma.sq	Numeric value. Only applicable when level set to "splicing". The variance threshold below which the included/excluded modality will be defined as primary sub-modality, and above which it will be defined as dispersed sub-modality. Please refer to AssignModality function help page for more details. Default is 0.001.
bimodal.adjust	Logical. Only applicable when level set to "splicing". When set to TRUE, MARVEL will identify false bimodal modalities and reassign them as included/excluded modality. Please refer to AssignModality function help page for more details.
seed	Numeric value. Only applicable when level set to "splicing". Ensure the fitdist function returns the same values for alpha and beta paramters each time this function is executed using the same random number generator. Please refer to AssignModality function help page for more details.
modality.column	Character string. Only applicable when level set to "splicing". Can take the value "modality", "modality.var" or "modality.bimodal.adj". Please refer to AssignModality function help page for more details. Default is "modality.bimodal.adj".
scale.y.log	Logical value. Only applicable when level set to "splicing". If set to TRUE, the y-axis of will log10-scaled. Useful when most PSI values are extremely small (< 0.02) or big (> 0.98). Default is FALSE.
max.cells.jitter	Numeric value. Only applicable when level set to "splicing". Maximum number of cells for jitter points. Cells are randomly downsampled to show on jitter plot. Useful when there are large number of cells so that individual jitter points do not overcrowd the violin plot. Specified together with max.cells.jitter.seed. To disable this option, specify a value large than the number of cells in each cell group.
max.cells.jitter.seed	Numeric value. Only applicable when level set to "splicing". Cells downsampled are reproducible. Specified together with max.cells.jitter.
cell.group.colors	Character string. Vector of colors for the cell groups specified for PCA analysis using cell.type.columns, cell.type.variable, and cell.type.labels. If not specified, default ggplot2 colors will be used.
point.alpha	Numeric value. Transparency of the data points. Takes any values between 0-1. Default value is 0.2.

Value

An object of class S3 with new slot $adhocPlot$PSI or MarvelObject$adhocPlot$Exp when level set to "splicing" or "gene", respectively.

Examples

marvel.demo <- readRDS(system.file("extdata/data", "marvel.demo.rds", package="MARVEL"))

# Define cell groups to plot
df.pheno <- marvel.demo$SplicePheno
cell.group.g1 <- df.pheno[which(df.pheno$cell.type=="iPSC"), "sample.id"]
cell.group.g2 <- df.pheno[which(df.pheno$cell.type=="Endoderm"), "sample.id"]
cell.group.list <- list(cell.group.g1, cell.group.g2)
names(cell.group.list) <- c("iPSC", "Endoderm")

# Plot
marvel.demo <- PlotValues(MarvelObject=marvel.demo,
                          cell.group.list=cell.group.list,
                          feature="chr17:8383254:8382781|8383157:-@chr17:8382143:8382315",
                          level="splicing",
                          min.cells=5,
                          xlabels.size=5
                          )

# Check output
marvel.demo$adhocPlot$PSI

MARVEL documentation built on Oct. 31, 2022, 5:07 p.m.