file.path <- system.file("extdata", "TCRB_sequencing", package = "LymphoSeq")
file.list <- readImmunoSeq(path = file.path)
productive.aa <- productiveSeq(file.list = file.list, aggregate = "aminoAcid")
similarity.matrix <- similarityMatrix(productive.seqs = productive.aa)
pairwisePlot(matrix = similarity.matrix)
bhattacharyya.matrix <- bhattacharyyaMatrix(productive.seqs = productive.aa)
pairwisePlot(matrix = bhattacharyya.matrix)
# Change plot color, title legend, and add title
pairwisePlot(matrix = similarity.matrix) +
ggplot2::scale_fill_gradient(low = "#deebf7", high = "#3182bd") +
ggplot2::labs(fill = "Similarity score") + ggplot2::ggtitle("Figure Title")
Run the code above in your browser using DataLab