tests/testthat/test-wilcox_rowscore.R

In montilab/CaDrA: Candidate Driver Analysis

test_that("wilcox_rowscore returns correct results", {
  
  mat <- matrix(c(1,0,1,0,0,0,0,0,1,0, 
                  0,0,1,0,1,0,1,0,0,0,
                  0,0,0,0,1,0,1,0,1,0), nrow=3)
  
  colnames(mat) <- 1:10
  row.names(mat) <- c("TP_1", "TP_2", "TP_3")
  
  set.seed(42)
  input_score = rnorm(n = ncol(mat))
  names(input_score) <- colnames(mat)
  
  result <- wilcox_rowscore(
    FS = mat, 
    input_score = input_score, 
    meta_feature = NULL,
    alternative = "less"
  )
  
  testthat::expect_length(result, 3L)
  testthat::expect_type(result, "double")
  testthat::expect_identical(result, c(TP_1=3, TP_2=13,TP_3=5))
  
})

## --------------------------------------------------- ##
test_that("wilcox_rowscore issues error messages when needed", {
  
  FS <-  data.frame(a = rnorm(10), b = rnorm(10))
  set.seed(42)
  input_score = rnorm(n = ncol(FS))
  names(input_score) <- colnames(FS)
  
  expect_error( 
    wilcox_rowscore(
      FS = FS,  
      input_score = input_score,
      meta_feature = NULL
    )
  )
  
})


## --------------------------------------------------- ##
test_that("wilcox_score issues correct errors", {
  
  row <- c(1, 0, 1, 0, 0, 0, 0, 0, 1, 0)
  ranks <- 1:length(row)
  x = ranks[which(row == 1)]
  y = ranks[which(row == 0)]
  
  mu = NA
  expect_error(wilcox_score(x, y, mu))
  
  mu = 0
  expect_error(wilcox_score(x=c(), y, mu))
  expect_error(wilcox_score(x=x, y=c(), mu))
  
})
  
## --------------------------------------------------- ##
test_that("wilcox_score returns correct results", {
  
  row <- c(1, 0, 1, 0, 0, 0, 0, 0, 1, 0)
  ranks <- 1:length(row)
  x = ranks[which(row == 1)]
  y = ranks[which(row == 0)]
  mu = 0
  
  result <- wilcox_score(x=x, y=y, mu=mu)
  expect_named(result, c("score.W", "p_value"))
  expect_equal(result[1], c("score.W"=7.0))
  
  expect_equal(result[2], 
               c("p_value"= pnorm( (7 - 21/2 + 0.5)/sqrt(21 * 11/12))))
  
})