commonSeqsVenn: Common sequences Venn diagram

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

Common sequences Venn diagram

Description

Creates a Venn diagram comparing the number of common sequences in two or three repertoire_ids.

Usage

commonSeqsVenn(repertoire_ids, productive_aa)

Arguments

repertoire_ids	A character vector of two or three names of repertoire_ids in productiveSeq table to compare.
productive_aa	A tibble of amino acid sequences generated by the LymphoSeq function productiveSeq.

Value

Returns a a Venn diagram of the number of common sequences between two or three repertoire_ids.

See Also

commonSeqs

Examples

file_path <- system.file("extdata", "TCRB_sequencing", package = "LymphoSeq2")

stable <- readImmunoSeq(path = file_path)

atable <- productiveSeq(study_table = stable, aggregate = "junction_aa")

# Plot a triple Venn diagram
commonSeqsVenn(repertoire_ids = c("TRB_Unsorted_0", 
   "TRB_Unsorted_32", "TRB_Unsorted_83"), 
   productive_aa = atable)

# Plot a double Venn diagram
commonSeqsVenn(repertoire_ids = c("TRB_Unsorted_0", 
   "TRB_Unsorted_32"), productive_aa = atable)

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