tests/testthat/test-anova.R

In computational-metabolomics/structtoolbox: Data processing & analysis tools for Metabolomics and other omics

# testanova class
test_that('anova 1way TYPE 1',{
  # DatasetExperiment
  D=iris_DatasetExperiment()
  # method
  ME=ANOVA(formula=y~Species,ss_type = 'I')
  ME=model_apply(ME,D)
  # expect all true
  expect_true(all(ME$significant[,1]))

})

# testanova class
test_that('anova 1way type 2',{
    # DatasetExperiment
    D=iris_DatasetExperiment()
    # method
    ME=ANOVA(formula=y~Species,ss_type='II')
    ME=model_apply(ME,D)
    # expect all true
    expect_true(all(ME$significant[,1]))
    
})

# testanova class
test_that('anova 1way type 3',{
    # DatasetExperiment
    D=iris_DatasetExperiment()
    # method
    ME=ANOVA(formula=y~Species,ss_type='III')
    ME=model_apply(ME,D)
    # expect all true
    expect_true(all(ME$significant[,1]))
    
})


test_that('anova 2way',{
  set.seed('57475')
  # DatasetExperiment
  D=iris_DatasetExperiment()
  D$sample_meta$fake_news=sample(D$sample_meta$Species,150,replace=FALSE)
  # method
  ME=ANOVA(formula=y~Species*fake_news)
  ME=model_apply(ME,D)
  # expect all true
  expect_true(all(ME$significant[,1]))

})

# test HSDEM class
test_that('hsd 1 factor',{
  set.seed('57475')
  # DatasetExperiment
  D=iris_DatasetExperiment()
  # method
  ME=HSD(formula=y~Species)
  ME=model_apply(ME,D)
  # expect all true
  expect_true(all(ME$significant[,1]))

})

# test HSDEM class
test_that('hsd 2 factors',{
  set.seed('57475')
  # DatasetExperiment
  D=iris_DatasetExperiment()
  D$sample_meta$fake_news=sample(D$sample_meta$Species,150,replace=FALSE)
  # method
  ME=HSD(formula=y~Species*fake_news)
  ME=model_apply(ME,D)
  # expect all true
  expect_true(all(ME$significant[,1]))
})