test larger data sets for correct jaccard result

- closes #204
rnabioco · Mar 19, 2017 · 43e4d85 · 43e4d85
1 parent 82266c8
commit 43e4d85
Show file tree

Hide file tree

Showing 2 changed files with 14 additions and 0 deletions.
diff --git a/NEWS.md b/NEWS.md
@@ -1,5 +1,9 @@
 # valr 0.1.2.9000
 
+## Bug fixes
+
+* Merge intervals in `bed_jaccard()` and use numeric values for calculation (fixes #204).
+
 # valr 0.1.2
 
 ## Major changes

diff --git a/tests/testthat/test_jaccard.r b/tests/testthat/test_jaccard.r
@@ -15,3 +15,13 @@ test_that("jaccard coeff is calculated correctly", {
   res <- bed_jaccard(x, y)
   expect_equal(res$jaccard, 0.25)
 })
+
+test_that("jaccard coeff is calc'd for large data sets", {
+  genome <- read_genome(valr_example('hg19.chrom.sizes.gz'))
+
+  x <- bed_random(genome, n = 1e5, seed = 10000)
+  y <- bed_random(genome, n = 1e5, seed = 20000)
+
+  res <- bed_jaccard(x, y)
+  expect_equal(res$jaccard, 0.01613906)
+})