tests/testthat/test_moments.R

In CharlotteJana/pdmppoly: Moment Approximation for Polynomial PDMPs

#======== todo =================================================================
#s1 colmoment + Abhänigkeit von LaF komplett rausnehmen

context("moments")

data(simplePdmp)
tmp <- tempfile()

test_that("'moments' works for objects of class 'multSimCsv'", {
  skip_if_not(exists('colmoment', where = asNamespace('LaF'), mode = 'function'))
    
  # simulation
  ms <- multSim(simplePdmp, seeds = 1:5)
  msCsv <- multSimCsv(simplePdmp, seeds = 1:5, prefix = tmp)
  
  # test that method 'mean' and 'moments(..., order = 1)' 
  # give equal results
  means <- mean(ms)
  meansCsv <- mean(msCsv)
  moments <- moments(ms, order = 1)
  momentsCsv <- moments(msCsv, order = 1)
  
  expect_equal(means, dplyr::select(moments, -order))
  expect_equal(means, meansCsv, tolerance = 1e-07)
  expect_equal(moments, momentsCsv, tolerance = 1e-07)
  # results are not exactly equal because 'multSimCsv' only stores 7 digits.

  # test that method 'moments' gives equal results 
  # for 'multSim' and 'multSimCsv'
  m2 <- moments(ms, order = 2)
  m2Csv <- moments(msCsv, order = 2)
  expect_equal(m2, m2Csv, tolerance = 1e-07)
  
  m10 <- moments(ms, order = 5)
  m10Csv <- moments(msCsv, order = 5)
  expect_equal(m10, m10Csv, tolerance = 1e-06)
  
})

unlink(paste0(tmp, "*")) # remove csv files