tests/testthat/test_integration.R

In DistributionalDDA: Distributional Data Analysis in Bayes Spaces

test_that("test BS expectation", {

    BS1 = getBayesSpace(StdReferenceMeasureN01)
    BS3 = getBayesSpace(ReferenceMeasureDiscrete(c("A","B","C"),c(0.2,0.3,0.5)))
    BS2 = getBayesSpace(ReferenceMeasureCartesianProduct(StdReferenceMeasureN01,
                                            StdReferenceMeasureUnif01
                                            ))

    BB1 = getBayesSpaceBasis(BS1)
    BB2 = getBayesSpaceBasis(BS2)
    BB3 = getBayesSpaceBasis(BS3)
    D1 = getBayesSpaceBasisDimension(BB1)
    D2 = getBayesSpaceBasisDimension(BB2)
    D3 = getBayesSpaceBasisDimension(BB3)
    
    E1 = expectation(getReferenceMeasure(BS1),as.function(BB1),tolerance=5E-3)
    E2 = expectation(getReferenceMeasure(BS2),as.function(BB2),tolerance=5E-3)
    E3 = expectation(getReferenceMeasure(BS3),as.function(BB3),tolerance=5E-3)

    expect_equal(c(E1),rep(0,D1),tolerance=5E-2)
    expect_equal(c(E2),rep(0,D2),tolerance=5E-2)
    expect_equal(c(E3),rep(0,D3),tolerance=5E-2)

})