LymphoSeq (version 1.0.2)

commonSeqsVenn: Common sequences Venn diagram

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("TCRB_Day0_Unsorted", 
   "TCRB_Day32_Unsorted", "TCRB_Day83_Unsorted"), 
   productive.seqs = productive.aa)

# Plot a double Venn diagram
commonSeqsVenn(samples = c("TCRB_Day0_Unsorted", 
   "TCRB_Day32_Unsorted"), 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("TCRB_Day0_Unsorted", "TCRB_Day32_Unsorted"), 
   productive.seqs = productive.aa)

dev.off()