This amounts to a rectangular integration over a multivariate normal, for which an efficient algorithm is available[41]. The current implementation of this integration is suitable to estimate p-values larger than 10−15. To approximate gene-wise p-values below this limit we multiply the minimum p-value of SNPs in the region with the effective number of tests within the gene (see Section Algorithmic details for gene-score calculation).
