test_that(
  "Shannon entropy. Expected = 0.86775",
  {
    expected <- 0.86775
    data(Qesem_s)
    actual <- shannon(Qesem_s$HU)
    expect_equal(round(actual,5), round(expected,5))
  }
)