context("Compare results with those in Baggelaar et al (2010)")

test_that("EQR-results on page 22 can be reproduced.", {
    data(EQR)
    result <- mya(EQR)
    expect_that(result$MYA, equals(0.167985609))
    expect_that(result$q05, equals(0.075350782))
    expect_that(result$q95, equals(0.333436559))
    expect_that(result$PROB_LTT, equals(0.988412826))
    expect_that(result$PROB_GTT, equals(0.011587173501))
})

test_that("DCA-results on page 24 can be reproduced.", {
    data(DCA)
    result <- mya(DCA)
    expect_that(result$MYA, equals(9.5))
    expect_that(result$q05, equals(7.5852374065))
    expect_that(result$q95, equals(11.4147625935))
    expect_that(result$PROB_LTT, equals(0.7372895))
    expect_that(result$PROB_GTT, equals(0.2627105))
})