I am performing peak annotation using peakAnno() in ChIPseeker. It works well. I was wondering if there is a way in ChIPseeker to do any statistical comparison for the overlap of peaks with the different genomic regions (i.e. intron, exon, utr, distal intergenic). This would provide an accurate analysis of relative over-/under-representation since things like distal Intergenic regions make up the majority of the genome. A proper statistical test would rule out the bias due to predominance of specific elements.