tests/testthat/test-domain_annotation.R

In almeidasilvaf/tfhunter: Identification and classification of plant transcription factors

#----Load data------------------------------------------------------------------
data(gsu)
seq <- gsu[1:2]

#----Start tests----------------------------------------------------------------
test_that("annotate_pfam() returns a data frame of domain annotation", {
    
    annot_local <- data.frame(
        Gene = c("gsu04730.1", "gsu00340.1"),
        Domain = c("PF02042", "PF00076")
    )
    if(hmmer_is_installed()) {
        annot_local <- annotate_pfam(seq)
    }
    
    expect_equal(class(annot_local), "data.frame")
    expect_equal(ncol(annot_local), 2)
    expect_equal(names(annot_local), c("Gene", "Domain"))
    expect_true(nrow(annot_local) >= 2)
})