tests/testthat/test-smedian.sample.R

In KrasnitzLab/CNprep: Pre-process DNA Copy Number (CN) Data for Detection of CN Events

### Unit tests for smedian.sample function

library(CNprep)


### Tests smedian.sample() results

context("smedian.sample() results")

test_that("smedian.sample() must return expected results 01", {
    
    RNGkind(kind="Mersenne-Twister", normal.kind="Inversion", 
            sample.kind="Rejection")
    
    position <- c(1, 5)
    values <- c(0.072073840, 0.119913919, 0.154459489, 0.040994620, -0.082843732,
                0.111907384, 0.001913919, 0.032259489, 0.140994620, -0.000843732)
    
    set.seed(2111)
    
    results <- CNprep:::smedian.sample(position, values)
    
    expected <- 0.119913919000
    
    expect_equal(results, expected)
})


test_that("smedian.sample() must return expected results 02", {
    
    RNGkind(kind="Mersenne-Twister", normal.kind="Inversion", 
            sample.kind="Rejection")
    
    position <- c(1, 5)
    values <- c(0.072073840, 0.119913919, 0.154459489, 0.040994620, -0.082843732,
                0.111907384, 0.001913919, 0.032259489, 0.140994620, -0.000843732)
    
    set.seed(311)
    
    results <- CNprep:::smedian.sample(position, values)
    
    expected <- 0.0409946200
    
    expect_equal(results, expected)
})

test_that("smedian.sample() must return expected results when NA present", {
    
    RNGkind(kind="Mersenne-Twister", normal.kind="Inversion", 
            sample.kind="Rejection")
    
    position <- c(1, 5)
    values <- c(0.072073840, NA, NA, NA, NA,
                0.111907384, NA, NA, NA, -0.000843732)
    
    set.seed(3311)
    
    results <- CNprep:::smedian.sample(position, values)
    
    expected <- 0.07207384000
    
    expect_equal(results, expected)
})