HypergeoMat

library(HypergeoMat)
library(microbenchmark)

jhpq <- hypergeomPFQ_julia()
## Starting Julia ...

microbenchmark(
  Rcpp = hypergeomPFQ(m = 30, a = c(1, 2, 3), b = c(4, 5), x = c(0.1, 0.2, 0.3)),
  Julia = jhpq(m = 30, a = c(1, 2, 3), b = c(4, 5), x = c(0.1, 0.2, 0.3)),
  times = 10
)
## Unit: milliseconds
##   expr       min        lq      mean    median       uq      max neval
##   Rcpp 3171.9094 3266.6886 3386.5241 3317.9526 3531.426 3659.559    10
##  Julia  339.2297  355.8744  755.2628  377.9874  438.013 3866.369    10