runINext: Run iNEXT on repertoire_ids

In elulu3/LymphoSeqTest: Analyze high-throughput sequencing of T and B cell receptors

Run iNEXT on repertoire_ids

Description

Given a repertoire_id table, for each generate rarefaction curves to estimate repertoire diversity. The method used to generate the rarefaction curve is derived from Chao et al., (2014) using the iNEXT library

Usage

runINext(sample_table, color = "repertoire_id")

Arguments

sample_table	A tibble consisting antigen receptor sequencing data imported by the LymphoSeq function readImmunoSeq. "junction_aa", "duplicate_count", and "duplicate_frequency" are required columns.
color

Examples

file_path <- system.file("extdata", "TCRB_sequencing", package = "LymphoSeq2")
stable <- readImmunoSeq(path = file_path)
atable <- productiveSeq(stable, 
                        aggregate = "junction_aa", 
                        prevalence = TRUE)
atable <- atable %>% dplyr::filter(repertoire_id == "TRB_Unsorted_1320")
rtable <- runINext(atable)

elulu3/LymphoSeqTest documentation built on Aug. 27, 2022, 5:47 a.m.