tests/testthat/test-calculate_physicochemical_properties.R

In currocam/taxa2hmmer: Get more out of homologous sequence search with HMMER

test_that("calculate works", {
  data <- read_hmmer_from_xml(
    xml_file_paths = "phmmer_example.xml",
    fullseq_fasta_paths = "phmmer_example.fa")
  HMMER_tidy_tbl <- extract_from_HMMER_data_tbl(data)
  suppressWarnings({
    calculate_physicochemical_properties(HMMER_tidy_tbl$hits.fullseqfasta) %>%
      write.csv("properties.csv")
  }
  )
  testthat::expect_snapshot_file("properties.csv")
})