I made an R package called ExOutBench, which is used to benchmark the results of a gene expression outlier calling method. It calculates the relative risk of having a nearby rare variant gene given outlier status. Check out the documentation here: https://jeinson.github.io/ExOutBench/index.html