Add `coverage()` accuracy metric

[mitchelloharawild]

coverage <- function(.dist, .actual, percentile = 80, na.rm = TRUE, ...){
  intvl <- hilo(.dist, percentile)
  within <- (.actual > intvl$lower) & (.actual < intvl$upper)
  mean(within)
}