tests/testthat/test-get_variants_in_genes.R

In frhl/our: genetic data analysis

context('get_variants_in_genes')

test_that("standard functionality", {

  data("hsapiens_gene_ensembl_grch38_chr1_24")
  df = data.frame(chr=c(21,21,21),bp=c(5000000,5022493,5050000),id=c('rs1','rs2','rs3'))
  expect_true(nrow(get_variants_in_genes(df , hsapiens_gene_ensembl_grch38_chr1_24)) == 1)

  data("hsapiens_gene_ensembl_grch37_chr1_24")
  df = data.frame(chr=c(21,21,21),bp=c(5000000,5022493,5050000),id=c('rs1','rs2','rs3'))
  expect_true(nrow(get_variants_in_genes(df , hsapiens_gene_ensembl_grch37_chr1_24)) == 0)

})