jaccard reporting invalid results

[jayhesselberth]

jaccard values should be [0, 1].

library(valr)
library(tidyverse, warn.conflicts = FALSE)
#> Loading tidyverse: ggplot2
#> Loading tidyverse: tibble
#> Loading tidyverse: tidyr
#> Loading tidyverse: readr
#> Loading tidyverse: purrr
#> Loading tidyverse: dplyr
#> Conflicts with tidy packages ----------------------------------------------
#> filter(): dplyr, stats
#> lag():    dplyr, stats

genome <- tribble(
  ~chrom, ~size,
  "chr1", 1e6,
  "chr2", 1e6
)

x <- bed_random(genome, n = 1e4, seed = 10000)

bed_jaccard(x, x)
#> # A tibble: 1 × 4
#>      len_i    len_u   jaccard      n
#>      <int>    <int>     <dbl>  <int>
#> 1 60095982 20000000 -1.498803 110200

[jayhesselberth]

also this.

library(valr)
library(tidyverse, warn.conflicts = FALSE)
#> Loading tidyverse: ggplot2
#> Loading tidyverse: tibble
#> Loading tidyverse: tidyr
#> Loading tidyverse: readr
#> Loading tidyverse: purrr
#> Loading tidyverse: dplyr
#> Conflicts with tidy packages ----------------------------------------------
#> filter(): dplyr, stats
#> lag():    dplyr, stats

genome <- tribble(
  ~chrom, ~size,
  "chr1", 1e6,
  "chr2", 1e6
)

x <- bed_random(genome, n = 1e5, seed = 10000)

bed_jaccard(x, x)
#> Warning in summarise_impl(.data, dots): integer overflow - use
#> sum(as.numeric(.))
#> # A tibble: 1 × 4
#>   len_i     len_u jaccard        n
#>   <int>     <int>   <dbl>    <int>
#> 1    NA 200000000      NA 10114884