I am trying to use ChIPseeker to annotate and do functional enrichment analysis for H3K4me3 peaks. I am a bit confused using seq2gene function. How can I only select those genes which has H3K4me3 peak in +/- 2000bp? From the tutorial, in the arguments, there's option for "tssRegion = c(-1000, 1000), flankDistance = 3000". How can I configure this to get only those genes which has H3K4me3 peaks in TSS +/- 2000bp?
Thank you in advance for your help!