tests/testthat/test-galgo.R

In harpomaxx/GSgalgoR: An Evolutionary Framework for the Identification and Study of Prognostic Gene Expression Signatures in Cancer

context("galgo function")

test_that("a valid galgo.Obj is generated", {
   set.seed(29042010)
  # load example dataset
   library(breastCancerTRANSBIG)
   data(transbig)
   Train <- transbig
  
   expression <- Biobase::exprs(Train)
   clinical <- Biobase::pData(Train)
   OS <- survival::Surv(time = clinical$t.rfs, event = clinical$e.rfs)
  
   # We will use a reduced dataset for the example
   expression <- expression[sample(seq_len(nrow(expression)), 100), ]
  
   # Now we scale the expression matrix
   expression <- t(scale(t(expression)))
  
   # Run galgo

   output <- GSgalgoR::galgo(generations = 2, population = 5, 
      prob_matrix = expression, OS = OS)
   expect_is(output,"galgo.Obj")
  
})


test_that("a valid galgo.Obj is generated using verbose 1", {
   set.seed(29042010)
   # load example dataset
   library(breastCancerTRANSBIG)
   data(transbig)
   Train <- transbig
   
   expression <- Biobase::exprs(Train)
   clinical <- Biobase::pData(Train)
   OS <- survival::Surv(time = clinical$t.rfs, event = clinical$e.rfs)
   
   # We will use a reduced dataset for the example
   expression <- expression[sample(seq_len(nrow(expression)), 100), ]
   
   # Now we scale the expression matrix
   expression <- t(scale(t(expression)))
   
   # Run galgo
 
   output <- GSgalgoR::galgo(generations = 2, 
                           population = 4, 
                           prob_matrix = expression,
                           OS = OS,
                           verbose = 1)
   expect_is(output,"galgo.Obj")
   
})

test_that("a valid galgo.Obj is generated using verbose 0", {
   set.seed(29042010)
   # load example dataset
   library(breastCancerTRANSBIG)
   data(transbig)
   Train <- transbig
   
   expression <- Biobase::exprs(Train)
   clinical <- Biobase::pData(Train)
   OS <- survival::Surv(time = clinical$t.rfs, event = clinical$e.rfs)
   
   # We will use a reduced dataset for the example
   expression <- expression[sample(seq_len(nrow(expression)), 100), ]
   
   # Now we scale the expression matrix
   expression <- t(scale(t(expression)))
   
   # Run galgo
   
   output <- GSgalgoR::galgo(generations = 2, 
                           population = 4, 
                           prob_matrix = expression,
                           OS = OS,
                           verbose = 0)
   expect_is(output,"galgo.Obj")
   
})

test_that("not suitable solution is correctly handled", {
   set.seed(29042010)
   # load example dataset
   library(breastCancerTRANSBIG)
   data(transbig)
   Train <- transbig
   
   expression <- Biobase::exprs(Train)
   clinical <- Biobase::pData(Train)
   OS <- survival::Surv(time = clinical$t.rfs, event = clinical$e.rfs)
   
   # We will use a reduced dataset for the example
   expression <- expression[sample(seq_len(nrow(expression)), 100), ]
   
   # Now we scale the expression matrix
   expression <- t(scale(t(expression)))
   
   # Run galgo
   
   output <- GSgalgoR::galgo(generations = 2, 
                             population = 1, 
                             prob_matrix = expression,
                             OS = OS,
                             verbose = 0)
  expect_null(output)
   
})