tests/testthat/test-predict_ic50s.R

In tzina97/netMHCpanW: Interface to 'NetMHCpan'

test_that("use", {
  if (!is_netmhcpan_installed()) return()
  protein_sequence <- "FANTASTICALLY"
  peptide_length <- 13
  temp_fasta_filename <- tempfile()
  temp_xls_filename <- tempfile()

  t <- predict_ic50s(
    protein_sequence = protein_sequence,
    peptide_length = peptide_length,
    mhc_haplotype = "HLA-A02:329",
    temp_fasta_filename = temp_fasta_filename,
    temp_xls_filename = temp_xls_filename
  )

  # Cleanup
  expect_true(!file.exists(temp_fasta_filename))
  expect_true(!file.exists(temp_xls_filename))

  expect_equal(colnames(t), c("peptide", "ic50"))
  expect_true(is.character(t$peptide))
  expect_true(!is.factor(t$peptide))
  expect_true(is.numeric(t$ic50))
  expect_true(!is.factor(t$ic50))
  n_peptides <- nchar(protein_sequence) - peptide_length + 1
  expect_equal(n_peptides, nrow(t))
})