Creates a heat map from a Bhattacharyya, Similarity, Sorensen, or PSI matrix.
Arguments
- matrix
A Bhattacharyya, Similarity, Sorensen, or PSI matrix produced by the LymphoSeq2
scoringMatrix()function.
Details
The plot is made using the package ggplot2 and can be reformatted using ggplot2 functions. See examples below.
See also
An excellent resource for examples on how to reformat a ggplot can
be found in the R Graphics Cookbook online (http://www.cookbook-r.com/Graphs/).
The functions to create the similarity matrices can be found
here: scoringMatrix()
Examples
file_path <- system.file("extdata", "TCRB_sequencing",
package = "LymphoSeq2")
study_table <- LymphoSeq2::readImmunoSeq(path = file_path, threads = 1)
#> Dataset Analysis:
#> Files: 10, Total: 0.00 GB, Largest: 0.0 MB
#> Available memory: 12.0 GB
study_table <- LymphoSeq2::topSeqs(study_table, top = 100)
amino_table <- LymphoSeq2::productiveSeq(study_table,
aggregate = "junction_aa")
# Plot similarity using Similarity score
similarity_matrix <- LymphoSeq2::scoringMatrix(
productive_table = amino_table,
mode = "Similarity"
)
LymphoSeq2::pairwisePlot(matrix = similarity_matrix)
#> Warning: Ignoring unknown parameters: `linewidth`
# Plot similarity using Bhattacharyya score
bhattacharyya_matrix <- LymphoSeq2::scoringMatrix(
productive_table = amino_table,
mode = "Bhattacharyya"
)
LymphoSeq2::pairwisePlot(matrix = bhattacharyya_matrix)
#> Warning: Ignoring unknown parameters: `linewidth`
# Change plot color, title legend, and add title
LymphoSeq2::pairwisePlot(matrix = similarity_matrix) +
ggplot2::scale_fill_gradient(low = "#deebf7", high = "#3182bd") +
ggplot2::labs(fill = "Similarity score") +
ggplot2::ggtitle("Pairwise similarity score")
#> Warning: Ignoring unknown parameters: `linewidth`
#> Scale for fill is already present.
#> Adding another scale for fill, which will replace the existing scale.
