commonSeqsVenn: Common sequences Venn diagram

In davidcoffey/LymphoSeq: Analyze high-throughput sequencing of T and B cell receptors

Description

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

Usage

commonSeqsVenn(samples, productive.seqs)

Arguments

samples	A character vector of two or three names of samples in productive.seqs to compare.
productive.seqs	A list of productive 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 samples.

See Also

commonSeqs

Examples

file.path <- system.file("extdata", "TCRB_sequencing", package = "LymphoSeq")

file.list <- readImmunoSeq(path = file.path)

productive.aa <- productiveSeq(file.list = file.list, aggregate = "aminoAcid")

# Plot a triple Venn diagram
commonSeqsVenn(samples = c("TRB_Unsorted_0", 
   "TRB_Unsorted_32", "TRB_Unsorted_83"), 
   productive.seqs = productive.aa)

# Plot a double Venn diagram
commonSeqsVenn(samples = c("TRB_Unsorted_0", 
   "TRB_Unsorted_32"), productive.seqs = productive.aa)

# Save Venn diagram as a .png file to working directory
png(filename = "Venn diagram.png", res = 300, units = "in", height = 5, width = 5)

commonSeqsVenn(samples = c("TRB_Unsorted_0", "TRB_Unsorted_32"), 
   productive.seqs = productive.aa)

dev.off()

davidcoffey/LymphoSeq documentation built on Dec. 31, 2019, 9:52 p.m.