tests/testthat/test-constant.R

In rbentham/MCbiclust: Massive correlating biclusters for gene expression data and associated methods

data(CCLE_small)
data(Mitochondrial_genes)

mito.loc <- which(row.names(CCLE_small) %in% Mitochondrial_genes)
CCLE.mito <- CCLE_small[mito.loc,]

suppressWarnings(RNGversion("3.5.0")) 
set.seed(101)
CCLE.seed <- FindSeed(gem = CCLE.mito,seed.size = 10,
                      iterations = 100,messages = 50)

CCLE.sort <- SampleSort(gem = CCLE.mito,seed = CCLE.seed,sort.length = 11)

CCLE.samp.sort <- MCbiclust:::Vignette_sort[[1]]

CCLE.pc1 <- PC1VecFun(top.gem = CCLE.mito,
                      seed.sort = CCLE.samp.sort,
                      n = 10)

CCLE.hicor.genes <- as.numeric(HclustGenesHiCor(CCLE.mito,
                                                CCLE.seed,
                                                cuts = 8))

CCLE.cor.mat <- cor(t(CCLE.mito[CCLE.hicor.genes,CCLE.seed]))

gene.set1 <- labels(as.dendrogram(hclust(dist(CCLE.cor.mat)))[[1]])
gene.set2 <- labels(as.dendrogram(hclust(dist(CCLE.cor.mat)))[[2]])

gene.set1.loc <- which(row.names(CCLE.mito) %in% gene.set1)
gene.set2.loc <- which(row.names(CCLE.mito) %in% gene.set2)

CCLE.ps <- PointScoreCalc(CCLE.mito,gene.set1.loc,gene.set2.loc)

test_that('PointScoreCalc returns error and warnings in right situation',{
  expect_error(PointScoreCalc(CCLE.mito[,1],gene.set1.loc,gene.set2.loc))
  expect_error(PointScoreCalc(CCLE.mito[,1],gene.set1.loc,gene.set2.loc))
  expect_warning(PointScoreCalc(CCLE.mito[,seq_len(10)],gene.set1.loc,gene.set2.loc))
}
)

CCLE.mito[gene.set1.loc,1][10] <- NA 

test_that("HclustGenesHiCor and CVEval work with matrices",{
  expect_error(as.numeric(HclustGenesHiCor(as.matrix(CCLE.mito),
                                                CCLE.seed,
                                                cuts = 58)), NA)
  
  expect_error(CVEval(gem.part = as.matrix(CCLE.mito),
         gem.all = as.matrix(CCLE_small),
         seed = CCLE.seed,
         splits = 58), NA)
  
})

test_that("PC1VecFun returns as expected",{
  expect_equal(length(CCLE.pc1), 500)
})

CCLE.cor.vec <-  CVEval(gem.part = CCLE.mito,
                        gem.all = CCLE_small,
                        seed = CCLE.seed,
                        splits = 10)

CCLE.bic <- ThresholdBic(cor.vec = CCLE.cor.vec,sort.order = CCLE.samp.sort,
                         pc1 = as.numeric(CCLE.pc1))

CCLE.pc1 <- PC1Align(gem = CCLE_small, pc1 = CCLE.pc1,
                     cor.vec = CCLE.cor.vec ,
                     sort.order = CCLE.samp.sort,
                     bic =CCLE.bic)

CCLE.fork <- ForkClassifier(CCLE.pc1, samp.num = length(CCLE.bic[[2]]))

test_that("Basic functioning remains constant",{
  expect_equal_to_reference(CCLE.seed,"ccle_seed.rds")
  expect_equal_to_reference(CCLE.sort,"ccle_sort.rds")
  expect_equal_to_reference(CCLE.cor.vec,"ccle_cv.rds")
  expect_equal_to_reference(CCLE.bic,"ccle_bic.rds")
  #expect_equal_to_reference(CCLE.pc1,"ccle_pc1.rds")
  expect_equal_to_reference(CCLE.fork,"ccle_fork.rds")
  expect_equal_to_reference(CCLE.ps,"ccle_ps.rds")
})

CCLE.mito[gene.set1.loc,1][10] <- NA 
test_that("PointScoreCalc handles NA values",{
  expect_error(PointScoreCalc(CCLE.mito[,seq_len(11)],gene.set1.loc,gene.set2.loc),
               NA)
})

test_that("FindSeed returns as expected",{
  expect_equal(length(CCLE.seed), 10)
  expect_message(FindSeed(gem = CCLE.mito,seed.size = 10,
                          iterations = 100,messages = 50),"Iteration")
  expect_error(FindSeed(gem = CCLE.mito,seed.size = 501,
                        iterations = 100,messages = 50))
})

test_that("SampleSort returns as expected",{
  expect_equal(length(CCLE.sort), 11)
})

test_that("CVEval returns as expected",{
  expect_equal(length(CCLE.cor.vec), 1000)
})

test_that("ThresholdBic returns as expected",{
  expect_equal(length(CCLE.bic), 2)
  expect_equal(length(CCLE.bic[[1]]), 1000)
  expect_equal(length(CCLE.bic[[2]]), 58)
})

test_that("PC1Align returns as expected",{
  expect_equal(length(CCLE.pc1), 500)
})


CCLE.bic[[2]] <- CCLE.bic[[2]][1]

test_that("PC1Align works with bicluster with single sample",{
  expect_error(PC1Align(gem = CCLE_small, pc1 = CCLE.pc1,
                        cor.vec = CCLE.cor.vec ,
                        sort.order = CCLE.samp.sort,
                        bic =CCLE.bic), NA)
          }
          )