tests/testthat/test-basic.R

In xia-lab/MetaboAnalystR: An R Package for Comprehensive Analysis of Metabolomics Data

context("Testing MetaboAnalystR - Basic Functionality")

library(MetaboAnalystR)

test_that("Uploading CSV Files Works", {
  
  mSet <- InitDataObjects("conc", "stat", FALSE)
  mSet <- Read.TextData(mSet, "http://www.metaboanalyst.ca/resources/data/human_cachexia.csv", "rowu", "disc")
  
  expect_equal(length(mSet), 5)
  expect_match(mSet$dataSet$type, "conc")
  expect_match(mSet$dataSet$cls.type, "disc")
  expect_match(mSet$dataSet$format, "rowu")
  expect_equal(length(mSet$dataSet$cmpd), 63)
  expect_match(mSet$analSet$type, "stat")
  expect_false(mSet$dataSet$paired)
  
})

test_that("Sanity Check Works", {
  
  mSet <- InitDataObjects("conc", "stat", FALSE)
  mSet <- Read.TextData(mSet, "http://www.metaboanalyst.ca/resources/data/human_cachexia.csv", "rowu", "disc")
  mSet <- SanityCheckData(mSet)
  
  expect_equal(length(mSet$dataSet), 19)
  expect_equal(mSet$dataSet$small.smpl.size, 0)
  expect_equal(mSet$dataSet$cls.num, 2)
  
})

test_that("Replace Min Works", {
  
  mSet <- InitDataObjects("conc", "stat", FALSE)
  mSet <- Read.TextData(mSet, "http://www.metaboanalyst.ca/resources/data/human_cachexia.csv", "rowu", "disc")
  mSet <- SanityCheckData(mSet)
  mSet <- ReplaceMin(mSet)
  
  expect_equal(length(mSet$dataSet), 20)
  expect_equal(nrow(qs::qread("data_proc.qs")), 77)
  expect_equal(ncol(qs::qread("data_proc.qs")), 63)
  expect_match(mSet$msgSet$replace.msg, 
               "Zero or missing values were replaced by 1/5 of the min positive value for each variable.")
})

test_that("Normalization Works", {
  
  mSet <- InitDataObjects("conc", "stat", FALSE)
  mSet <- Read.TextData(mSet, "http://www.metaboanalyst.ca/resources/data/human_cachexia.csv", "rowu", "disc")
  mSet <- SanityCheckData(mSet)
  mSet <- ReplaceMin(mSet)
  mSet <- PreparePrenormData(mSet)
  mSet <- Normalization(mSet, "QuantileNorm", "LogNorm", "MeanCenter", ref=NULL, ratio=FALSE, ratioNum=20)  
  
  expect_equal(length(mSet$dataSet), 29)
  expect_equal(nrow(mSet$dataSet$norm), 77)
  expect_equal(ncol(mSet$dataSet$norm), 63)
  expect_match(mSet$dataSet$rownorm.method, "Quantile Normalization")
  expect_match(mSet$dataSet$trans.method, "Log10 Normalization")
  expect_match(mSet$dataSet$scale.method, "Mean Centering")
  expect_false(mSet$dataSet$combined.method)
  
})