custom annotation in annotatePeak
Entering edit mode
Last seen 6.1 years ago

Hi Guangchuang,

Is there a way to control annotations in the the annotatePeak function?  For instance I would like to ignore "downstream" and consider those peaks "intergenic".  Or for instance combine "1st exon" with "other exon" so there is one category of "exon".  

Any help would be great, thanks for doing this!


chipseeker Guangchuang Yu • 2.5k views
Entering edit mode
Guangchuang Yu ★ 1.2k
Last seen 5 months ago
China/Guangzhou/Southern Medical Univer…
> require(ChIPseeker)
Loading required package: ChIPseeker

ChIPseeker v1.15.3  For help:

If you use ChIPseeker in published research, please cite:
Guangchuang Yu, Li-Gen Wang, Qing-Yu He. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 2015, 31(14):2382-2383
> ff = getSampleFiles()[[4]]
> options(ChIPseeker.ignore_1st_exon = TRUE)
> options(ChIPseeker.ignore_1st_intron = TRUE)
> options(ChIPseeker.ignore_downstream = TRUE)
> options(ChIPseeker.ignore_promoter_subcategory = TRUE)
> x = annotatePeak(ff)
>> loading peak file...                 2018-01-09 11:36:33
>> preparing features information...         2018-01-09 11:36:33
>> identifying nearest features...         2018-01-09 11:36:34
>> calculating distance from peak to TSS...     2018-01-09 11:36:34
>> assigning genomic annotation...         2018-01-09 11:36:34
>> assigning chromosome lengths             2018-01-09 11:36:52
>> done...                     2018-01-09 11:36:52
Warning message:
In loadTxDb(TxDb) :
  >> TxDb is not specified, use 'TxDb.Hsapiens.UCSC.hg19.knownGene' by default...
> x
Annotated peaks generated by ChIPseeker
1331/1331  peaks were annotated
Genomic Annotation Summary:
     Feature  Frequency
6   Promoter 58.7528174
2     5' UTR  0.3005259
1     3' UTR  2.1036814
3       Exon  3.0052592
5     Intron  7.9639369
4 Intergenic 27.8737791

In github version, you can do it via options. Beware that annotatePeak always output detail annotations and such options only work for summarizing the result and visualization.

Entering edit mode
Last seen 6.1 years ago

Thank you for the quick response! Everything working great I appreciate it!



Login before adding your answer.

Traffic: 419 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6