tests/testthat/test_gaussianSamples.R

In pneuvial/sanssouci: Post Hoc Multiple Testing Inference

context("Simulation of Gaussian samples")

m <- 123
rho <- 0.2
n <- 100
pi0 <- 0.5

test_that("output dimensions are as expected for two-sample tests", {
    sim <- gaussianSamples(m, rho, n, pi0, SNR = 1, prob = 0.5)
    
    expect_equal(ncol(sim$X), n)
    expect_equal(nrow(sim$X), m)
    expect_length(sim$H, m)
    
    tbl <- table(sim$categ)
    expect_length(tbl, 2)
    
    tbl <- table(sim$H)
    expect_length(tbl, 2)
})

test_that("output dimensions are as expected for one-sample tests", {
    sim <- gaussianSamples(m, rho, n, pi0, SNR = 1, prob = 1)
    
    expect_equal(ncol(sim$X), n)
    expect_equal(nrow(sim$X), m)
    expect_length(sim$H, m)
    
    tbl <- table(sim$categ)
    expect_length(tbl, 1)
    
    tbl <- table(sim$H)
    expect_length(tbl, 2)
})


test_that("one- and two-sample tests are generated as desired", {
    # one-sample    
    sim <- gaussianSamples(m, rho, n, pi0, SNR = 1, prob = 0)
    tbl <- table(sim$categ)
    expect_length(tbl, 1)
    expect_equal(names(tbl), "0")
    
    # one-sample    
    sim <- gaussianSamples(m, rho, n, pi0, SNR = 1, prob = 1)
    tbl <- table(sim$categ)
    expect_length(tbl, 1)
    expect_equal(names(tbl), "1")
    
    # two-sample    
    sim <- gaussianSamples(m, rho, n, pi0, SNR = 1, prob = 0.5)
    tbl <- table(sim$categ)
    expect_length(tbl, 2)
})