tests/testthat/testpIBjell.R

In ypriverol/pIR: pIR: Estimating Isoelectric Point (pI) of Peptide and Proteins from Amino Acid Sequence.

library(pIR)

context("Isoelectric point estimation: Bjellvist Method")

test_that("Test of the Bjellvist Methods.. ", {
    # Test for the salomon pkSet
    seq <- "AGAAPYVQAFDSLLAGPVAE"

    # Test for expasy pKSet
    expect_equal(pIBjell(sequence = seq, pkSetMethod = "expasy"), 3.666)

    # Test for skoog pKSet
    expect_equal(pIBjell(sequence = seq, pkSetMethod = "skoog"), 2.935)

    # Test for calibrated pKSet
    expect_equal(pIBjell(sequence = seq, pkSetMethod = "calibrated"), 4.085)
}
)