tests/testthat/test.contourmap.R

In finnlindgren/excursions: Excursion Sets and Contour Credibility Regions for Random Fields

context("Contourmap")

test_that("Contourmap, test ind", {
  data <- integration.testdata1()
  ind1 = c(1,2,3,4)
  ind2 = c(4,3,2,1)
  ind3 = rep(FALSE,data$n)
  ind3[1:4] = TRUE

  res1 = contourmap(data$mu, data$Q, n.levels=2, ind=ind1,
                    seed=data$seed, alpha=0.1,max.threads=1)
  res2 = contourmap(data$mu, data$Q, n.levels=2, ind=ind2,
                    seed=data$seed, alpha=0.1,max.threads=1)
  res3 = contourmap(data$mu, data$Q, n.levels=2, ind=ind3,
                    seed=data$seed, alpha=0.1,max.threads=1)

  expect_equal(res1$F,res2$F,tolerance=1e-7)
  expect_equal(res2$F,res3$F,tolerance=1e-7)
})


test_that("Contourmap, P measures", {
  data <- integration.testdata1()

  res1 = contourmap(data$mu, data$Q, n.levels=4,
                    seed=data$seed, alpha=0.1,max.threads=1,
                    compute = list(F = FALSE, measures = c("P2","P1")))

  expect_equal(res1$P1,0.9217417,tolerance=1e-3)
  expect_equal(res1$P2,0.405841,tolerance=1e-3)
})