plotTrack: Plot clone

In elulu3/LymphoSeqTest: Analyze high-throughput sequencing of T and B cell receptors

Plot clone

Description

Creates alluvial tracking amino acid frequencies across multiple samples

Usage

plotTrack(
  clone_table,
  alist = NULL,
  apal = NULL,
  breaks = 2,
  alphas = 0,
  breaks_pal = c("#7fc97f", "#beaed4")
)

Arguments

clone_table	A tibble of productive amino acid sequences generated by LymphoSeq function cloneTrack
alist	An optional list of amino acid, if a list is provided only those sequences will be highlighted
apal	An optional list of palatte colors used for the amino acids to be highlighted
breaks	Add an additional band behind the alluvial plot to highligh group, specificies after which bar the break should appear
alphas	Set alpha for the band
breaks_pal	Specify palette for the bands

Details

The plot is made using the package ggplot2 and can be reformatted using ggplot2 functions. See examples below.

Value

An alluvial diagram tracking particular clone across samples.

Examples

file_path <- system.file("extdata", "TCRB_sequencing", package = "LymphoSeq2")
stable <- LymphoSeq2::readImmunoSeq(path = file_path)
atable <- LymphoSeq2::productiveSeq(stable, aggregate = "junction_aa")
ctable <- LymphoSeq2::cloneTrack(study_table = atable,
                                sample_list = c("TRB_CD8_949", "TRB_CD8_CMV_369"))
LymphoSeq2::plotTrack(ctable)

elulu3/LymphoSeqTest documentation built on Aug. 27, 2022, 5:47 a.m.