chordDiagramVDJ {LymphoSeq} | R Documentation

## Chord diagram of VJ or DJ gene associations

### Description

Creates a chord diagram showing VJ or DJ gene associations from one or more samples.

### Usage

    chordDiagramVDJ(sample, association = "VJ", colors = c("red", "blue"))

### Arguments

`sample` | A data frame consisting of frequencies of antigen receptor sequences. "vFamilyName", "jFamilyName", and if applicable, "dFamilyName" are a required columns. Using output from the LymphoSeq function topSeqs is recommended.
`association` | A character vector of gene familes to associate. Options include "VJ" or "DJ".
`colors` | A character vector of 2 colors corresponding to the V/D and J gene colors respectively.

### Details

The size of the ribbons connecting VJ or DJ genes correspond to the number of samples or number of sequences that make up that recombination event. The thicker the ribbon, the higher the frequency of the recombination.

### Value

Returns a chord diagram showing VJ or DJ gene associations from one or more samples.

### See Also

`topSeqs`

### Examples

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

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

    productive.nt <- productiveSeq(file.list = file.list, aggregate = "nucleotide")

    top.seqs <- topSeqs(productive.seqs = productive.nt, top = 1)

    chordDiagramVDJ(sample = top.seqs, association = "VJ", colors = c("red", "blue"))

    # Remove "TCRB" from gene family name
    top.seqs <- as.data.frame(apply(top.seqs, 2, function(x) gsub("TCRB", "", x)))

    chordDiagramVDJ(sample = top.seqs, association = "VJ", colors = c("red", "blue"))

[Package LymphoSeq version 1.15.0 Index]