tests/testthat/test-GLV.R

In microbiome/miaSim: Microbiome Data Simulation

test_that("simulateGLV", {

    # check simulateGLV
    set.seed(533)
    A <- miaSim::powerlawA(4, alpha = 1.01)
    tse <- simulateGLV(
        n_species = 4, A = A, t_start = 0, t_end=100,
        t_store = 100
    )

    expect_s4_class(tse, "TreeSummarizedExperiment")
    # expect_equal(dim(tse@assays@data@listData[["counts"]]), c(4, 100))
    expect_equal(dim(tse), c(4, 100))    

    expect_error(Error1 <- simulateGLV(n_species = 0.5))
    expect_error(Error2 <- simulateGLV(n_species = 5, A = 2, x = 2, b = 3))

    # check norm = TRUE
    set.seed(53)
    tse2 <- miaSim:::simulateGLV(
        n_species = 3,
        A = powerlawA(n_species = 3, alpha = 2), t_start = 0, t_end=100,
        t_store = 100, norm = TRUE
    )
    expect_type(assay(tse2, "counts"), "double")    
    expect_equal(dim(tse2), c(3, 100))    
})