Entering edit mode
Hi, I am wondering how I can find the exact definition used to do peak annotation when doing:
peakAnno <- annotatePeak(peaks, tssRegion=c(-5000, 5000),
TxDb=txdb, annoDb="org.Hs.eg.db")
For example: I know downstream is within 3kb from the gene but is it from end of gene or start of gene? etc.
If it is from org.Hs.eg.db can you let me know how I can find that out?
Thanks
Lin
Thanks for your reply. Sorry. which vignette were you referring to? this one http://bioconductor.org/packages/release/bioc/vignettes/ChIPseeker/inst/doc/ChIPseeker.html ?
What about distal intergenic?