tests/testthat/test_genefreq.R

In shashidhar22/LymphoSeq2: Analyze high-throughput sequencing of T and B cell receptors

context("Get VDJ gene frequency")
library(LymphoSeq2)
library(tidyverse)

test_that("Gene frequency charcterized by gene name sums to 1", {
  stable <- LymphoSeq2::readImmunoSeq("test_data", threads = 1)
  ntable <- LymphoSeq2::productiveSeq(stable, aggregate = "junction")
  gtable <- LymphoSeq2::geneFreq(ntable) %>%
    dplyr::group_by(repertoire_id, gene_type) %>%
    dplyr::summarise(freq_tot = sum(gene_frequency)) %>%
    dplyr::pull(freq_tot)
  expect_true(base::all(gtable == 1))
})

test_that("Gene frequency charcterized by gene name sums to 1", {
  stable <- LymphoSeq2::readImmunoSeq("test_data", threads = 1)
  ntable <- LymphoSeq2::productiveSeq(stable, aggregate = "junction")
  gtable <- LymphoSeq2::geneFreq(ntable, family = TRUE) %>%
    dplyr::group_by(repertoire_id, gene_type) %>%
    dplyr::summarise(freq_tot = sum(gene_frequency)) %>%
    dplyr::pull(freq_tot)
  expect_true(base::all(gtable == 1))
})